Goertzel Algorithmus (Frequenz + Phase)

[stoxx]

Hallo Peter,

ich hab nochmal eine Frage zu Deinem Code.


Also wenn ich einfach Deine Demo so in Delphi starte, und einen Phasenwinkel von 0 einstelle. die Frequenz erzeuge, dann suche.

und wenn ich dann auf Goertzel Button gehe.

dann kommt ein realteil von 0 und ein Imaginäranteil von 7500 raus.
Das entspricht einem Winkel von 90.

und dann korrigierst Du diesen Winkel aber um 90 Grad, so dass ein (zufällig) richtiger Phasenwinkel von 0 rauskommt.

phas := 90 - 180 / Pi * (arctan(imag/re));
Das ist sehr eigenartig. warum stimmt der imaginäre Anteil nicht?
Darf ich fragen, wo Du die Formel her hast?
der Code funktioniert soweit, erstmal .. bei einfachen Frequenzen.
Aber bei vielen Überlagerungen, stimmt das dann auch noch?
Der Imaginärteil bei einem Phasenwinkel von 0 müsste doch Null sein, und der Realteil irgendein positiver Wert, oder bin ich da falsch gewickelt?
Welche Formel hast Du dafür jetzt verwendet?
irgendwie stimmt da was noch nicht, oder was denkst Du?