Digitaler Filter steigt nicht aud Istwert

Gast #3692192
Hi zusammen!

Ich wollte mich in meinen Ferien zum ersten Mal mit digitalen Filtern 
befassen. Dabei bin ich über das Programm Winfilter 
(www.winfilter.20m.com) gestolpert. Das Programm ist ja supertoll, da 
man sich den C-Code gleich ausgeben kann.

Zum Test habe ich mal folgenden Filter entworfen:
IIR, Low Pass, Butterworth, Fsample=1000Hz, Fcut1=2Hz, Float

Das Programm spuckt mir dann folgenden C-Code aus:



/**************************************************************
WinFilter version 0.8
http://www.winfilter.20m.com
akundert@hotmail.com

Filter type: Low Pass
Filter model: Butterworth
Filter order: 2
Sampling Frequency: 1000 Hz
Cut Frequency: 2.000000 Hz
Coefficents Quantization: float

Z domain Zeros
z = -1.000000 + j 0.000000
z = -1.000000 + j 0.000000

Z domain Poles
z = 0.991114 + j -0.008807
z = 0.991114 + j 0.008807
***************************************************************/
#define NCoef 2
float iir(float NewSample) {
    float ACoef[NCoef+1] = {
        0.00003865708800424168,
        0.00007731417600848337,
        0.00003865708800424168
    };

    float BCoef[NCoef+1] = {
        1.00000000000000000000,
        -1.98222892979252840000,
        0.98238545061412497000
    };

    static float y[NCoef+1]; //output samples
    static float x[NCoef+1]; //input samples
    int n;

    //shift the old samples
    for(n=NCoef; n>0; n--) {
       x[n] = x[n-1];
       y[n] = y[n-1];
    }

    //Calculate the new output
    x[0] = NewSample;
    y[0] = ACoef[0] * x[0];
    for(n=1; n<=NCoef; n++)
        y[0] += ACoef[n] * x[n] - BCoef[n] * y[n];

    return y[0];
}



Zum Test bin ich jetzt hingegangen und habe die Funktion in meinem 
dsPIC33 1000 mal pro Sekunde folgendermaßen aufgerufen:

float filteredValue = iir(790);

Ich darf doch dann erwarten, dass der gefilterte Wert nach 2 Sekunden 
den Wert 790 annimmt. Tut er aber nicht. Der Wert bleibt bei 779,5757 
stehen. Woran kann das liegen? Hat das schon mal jemand gehabt?

Grüße euer Torti
#3692280
Hast Du es nachgerechnet?

Eventuell mit einem anderen Programm?

Wie kommst Du auf die Idee das 790 herauskommen soll?

Für mich schaut das Ergebnis vernünftig aus mit diesen Koeffizienten.

Im Anhang die Sprungantwort Deines Filters => ist übereinstimmend mit 
Deinen Ergebnissen (mit maximal kleinen Rundungsfehlern).
Angehängte Dateien:
#3692282
Torti schrieb:
> Der Wert bleibt bei 779,5757
> stehen. Woran kann das liegen?

Das liegt daran, dass das Filter im Durchlassbereich eine Dämpfung von 
A=0.1057dB aufweist.

Also kommt das gut hin mit den 780.

Bei mir siehts mit

3.913020539914440e-5
2*3.913020539914440e-5
3.913020539914440e-5

für die Zähler-Koeffizienten besser aus. Wahrscheinlich hat dein 
Programm nicht genau genug gerechnet.

Hab auch nochmal die beiden Sprungantworten geplottet. Blau deine 
Koeffizienten, rot meine.
Angehängte Dateien:

Antwort schreiben

Eine eigene Frage? Dieser Thread ist über 12 Jahre alt. Stell sie in einem eigenen Thread – dort lesen mehr Leute mit, und sie wird öfter und schneller beantwortet.

Neuen Thread starten

Eine Lösung, Ergänzung oder Rückmeldung zu diesem Thread? Dann antworte hier.

Formatierung
  • [c]C-Code[/c]
  • [code]Code in anderen Sprachen, ASCII-Zeichnungen[/code]
  • [math]Formel in LaTeX-Syntax[/math]
  • [[Titel]] – Link zu einem Artikel
  • Verweis auf einen anderen Beitrag: Rechtsklick auf dessen Beitrags-ID (#8078628), „Adresse kopieren“, in den Text einfügen

Mehr zur Formatierung

Längeren Sourcecode als Anhang, nicht in den Text. JPG nur für Fotos und Scans; Zeichnungen und Screenshots als PNG oder GIF. Siehe Bildformate.

Schon registriert? Anmelden

Wird nicht angezeigt. Wir schicken dir einen Link; sobald du damit bestätigst, erscheint dein Beitrag.

Mit dem Abschicken legst du einen kostenlosen Account an und erkennst die Nutzungsbedingungen und die Datenschutzerklärung an.