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 !