Goertzel Algorithmus (Frequenz + Phase)

[stoxx]

so .. diesen Fehler hab ich jetzt gefunden .. interessant.
Dein altes Delphi7 compiliert (zum Glück) anders .. sonst hätte man das Problem wohl nicht gleich gefunden .. bzw lag es genau daran, was im Buch besprochen wurde, dass es kein "guter" Algorithmus wäre, wegen Fehleranfällgikeit. (hier Fehlerfortpflanzung)

soo .. Deine Normierung werde ich mir also jetzt mal genauer ansehen !

hab vielen Dank ! .. meine Nacht ist gerettet
(Die klammer in der For Schleife)


Mit einer Klammer, die Du sicherheitshalber auch noch setzen kannst, geht es auch jetzt bei mir. erstaunlich..

Delphi-Quellcode:

  q1 := 0;
  q2 := 0;
  for i := Fenstergroesse downto 0 do  
  begin                                
    q3 := Signal[i]+ ( alpha*q1-q2 );
    q2 := q1;
    q1 := q3;
  end;


  re := q1-(alpha*q2)/2;
  imag := q2 * sin(2 * Pi * Frequenz);

  re := re / Fenstergroesse * 2;
  imag := imag / Fenstergroesse * 2;

im Prinzip kann man sich aber bei diesem Algo die Fenstergröße sparen ...