Goertzel Algorithmus (Frequenz + Phase)

[stoxx]

Hallo liebe Delphianer


wisst ihr zufällig, wo man im Goertzel Algorithmus die Phase (und Amplitude) herbekommt?

hab den Algorithmus mal in Delphi umgesetzt, die Frequenzbestimmung funktioniert soweit wie gewünscht. (Projekt im Anhang)
In diesem PDF steht leider nur, dass man die Phase bei der reinen Frequenzbestimmung ja nicht bräuchte .. hmm

http://www.ti.com/lit/an/spra096a/spra096a.pdf

ich bin jetzt in Mathe so tief leider nicht gleich so fit, das zu verstehen.
Vielleicht kann ja jemand helfen ?

Delphi-Quellcode:

function Goertzel(const x : TDoubleArray; frequency : Double; SampleRate : Integer) : Extended;
var Skn, Skn1, Skn2 : Extended;
    i : Integer;
    WNk : Extended;
    coeff : Extended;
    power : Extended;
begin

  Skn := 0; Skn1 := 0; skn2 := 0;
  coeff := 2 * Cos( 2 * Pi * frequency / SampleRate);
  for I := 0 to Length(x)- 1  do begin

    Skn :=  x[i] + (coeff * Skn1) - Skn2 ;
    Skn2 := Skn1;
    Skn1 := Skn;

  end; // for I 
 result := Skn2 * Skn2 + Skn1 * Skn1 - coeff * Skn1*skn2;
end;

Danke !