Kleines Problem bei Koefizientenbestimmung bei FIR

Jan R. Gast #3682993
Lesenswert?
• ▲
▼ 
Hallo,

möchte in ein DSP ein FIR Tiefpassfilter "einbauen". Dafür Wollte ich 
die Koeffizienten bestimmen, indem ich mit der Inversen ZDFT die 
Dirac-Impulsantwort eines Recteckfensters Berechene.



Natürlich habe ich jetzt die Sincfunktion herausbekommen, Allerdings, 
mit viel zu großen Amplituden, was zu einer Irren Verstärkung führt.

Was ist falsch?

Habe mich dann im Internet mal auf den Weg gemacht um etwas zu finden, 
wo diese Koeffizienten berechnet werden.

Habe das hier gefunden,
http://www.fh-kl.de/~hans.neuschwander/download_skripte/digfil/DigFilter.pdf

Hier wurden die Rohkoefizienten ebenfalls mit der Inversen FT berechnet, 
und die Lösungsfunktion entspricht bis auf die Multiplikation mit 2*fa 
genau der Meinen.

Was habe ich falsch gemacht? Danke schonmal.
Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Beitrag melden
Jan R. Gast #3683107
Lesenswert?
• ▲
▼ 
Detlef Kunz schrieb:
> Wieviel Taps (Stützstellen bzw. Koeffiziente) verwendest Du und was
> kommt dabei heraus, wenn Du alle Koeffizienten deines FIR addierst?

Mit Samplerate 1Khz und Fg= 100Hz

ist die Summe aller Koefizienten 1146,56

Das habe ich mit der Formel, die ich hergeleitet habe berechnet, nicht 
mit der aus dem Script.
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Michael A. #3683258
Lesenswert?
• ▲
▼ 
Jan R. schrieb:
> Habe das hier gefunden,
> http://www.fh-kl.de/~hans.neuschwander/download_skripte/digfil/DigFilter.pdf

Gl. 31 auf Seite 34 ist korrekt.

Bei dir steht aufm Bild
was wahrscheinlich
 heißen soll, das muss
 sein.

Desweiteren ist dein Sinus-Argument falsch (si-Funktion ist stets 
sin(x)/x, das seh ich bei dir nicht). Vielleicht auch nochmal eine 
anderen Quelle anschauen, da die Fenstermethode meiner Meinung nach in 
deiner Quelle schlecht erklärt wird.
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Jan R. Gast #3683276
Lesenswert?
• ▲
▼ 
bei der formel im skrip kürzt sich fg weg. da Im Zähler und  im nenner 
.keine ahnung warum das da noch steht. Bei mir bleibt aber Fa stehen, 
was seltsam ist.

Ist mein Vorgehen der integration Richtig? Ist FA gleich Fa da einmal 
großes A einmal kleines.  Ist es auch korekt, die grenzen auf +\- Fg zu 
legen, da die grenzfrequenz kleiner der Samplerate sein muss. und 
integral von 0, 0 ist.
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Detlef K. #3683344
Lesenswert?
• ▲
▼ 
Jan R. schrieb:
> Detlef Kunz schrieb:
>> Teile alle Koeffizienzen durch 1146,56
>> Die Summe aller Koeffizienten muss 1 ergeben.
>
> Ja das ergibt 1.

Hmm, lol
Natürlich ergibt das 1.
Du hast mich falsch verstanden.
Denn die Summe der Koeffizienten representiert die Verstärkung deines 
Tiefpasses, also im Moment G=1146,56 (61.2dB)
Wenn Du eine Verstärkung von 1 erwartest, dann musst Du jeden einzelnen 
Koeffizienten durch die Summer aller Koeffizienten Teilen.
Die Summe der neuen Koeffizienten sollte dann 1 ergeben.
:)
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Michael A. #3686323
Lesenswert?
• ▲
▼ 
Jan R. schrieb:
> Aber meine Ursprüngliche Frage, warum bei mir da fa sthenbleibt, ist
> noch nicht geklärt. Wo mein Fehler ist.

Das ist natürlich ein bisschen Kaffeesatzleserei, wenn uns die ganzen 
Zwischenschritte der Rechnung vorenthalten werden.

Ich kann mich nur wiederholen, such dir ein anderes Skript, die 
Erklärung der Fenstermethode deiner Quelle ist meiner Meinung nach 
schlecht gemacht.

Dein Ansatz mit der iDTFT von -fg bis fg ist aber korrekt. Wobei dein 
Frequenzgang F(f) abhängig von f sein muss, wenn du über f integrierst. 
Das mit dem e^N/2-Term erschwert die Sache nur, ich würde das ganze 
erstmal für ein nicht-kausales Filter berechnen und dann um N/2 
verschieben.

Was laut Glaskugel dein Fehler sein könnte, ist die Normierung auf fA im 
e-Term, die gehört da nicht hin, sonst müsstest du Integrationsvariable 
und -grenzen anpassen.
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Thread beobachten

Antwort schreiben

Bitte melde dich an, um einen Beitrag zu schreiben.

Bestehender Account

oder

Mit Google-Account einloggen

Noch kein Account?

Die Registrierung ist kostenlos und dauert nur eine Minute.

Jetzt registrieren