Übersetzung Programm in integer Arithmetik

Gast #928057
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
Gast #928147
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.

Antwort schreiben

Eine eigene Frage? Dieser Thread ist über 18 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.