Hallo, ich möchte ein Software-Filter (sehr langsames Filter 2 Hz Grenzfrequenz, 20000 Hz Sample Frequnez) bauen. Folgendes Filter funktioniert sehr gut mit double Zahlen: double input [SIZE]; double output [SIZE]; double gain = 3184.098757; // double weight = 0.9993718788; // xv[1] = xv[2]; xv[2] = (input[i]/*<<19*/)/gain; yv[1] = yv[2]; yv[2] = (xv[1]+xv[2])/*>>9*/ + (weight*yv[1])/*>>19*/; // output[i] = yv[2]; Das Filter hat die Kurve, die im angehängten xls file gezeigt ist Nun muß ich das Filter in integer Arithemtik umschreiben, da das Rechnen auf meinem rechenschwachen Mikrocontroller viel zu lange dauert. Mein Ergebnis: signed long input [SIZE]; signed long output [SIZE]; signed long gain = 12736; //3184.098757<<2 signed long weight = 523959; //0.9993718788<<19 xv[1] = xv[2]; xv[2] = (input[i]<<19)/gain; yv[1] = yv[2]; yv[2] = (xv[1]+xv[2])>>17 + (weight*yv[1])>>19; output[i] = yv[2]; Das Filter hat nun als Ausgangswert 0 obwohl gleich oft nach links bzw rechts geschoben wird. Woran liegt das? wo liegt der Fehler Gruß Samuel
Du wolltest wahrscheinlich die Klammern in
1 | yv[2] = (xv[1]+xv[2])>>17 + (weight*yv[1])>>19; |
etwas anders setzen, nämlich so:
1 | yv[2] = (xv[1]+xv[2]>>17) + (weight*yv[1]>>19); |
+ bindet stärker als >>. Trotzdem stimmt die Berechnung noch nicht, da ja nach Eingabewert entweder ein Überlauf passiert oder das Ergebnis 0 wird.
Bitte melde dich an um einen Beitrag zu schreiben. Anmeldung ist kostenlos und dauert nur eine Minute.
Bestehender Account
Schon ein Account bei Google/GoogleMail? Keine Anmeldung erforderlich!
Mit Google-Account einloggen
Mit Google-Account einloggen
Noch kein Account? Hier anmelden.