Phasenverschiebung zw. 2 Signalen bestimmen

>>Gibt es quasi eine DFT mit nur einer Frequenz (diese stimmt ja überein)???

Ja, das ist der Goertzel Algorithmus.

>>Hmmm...Nulldurchgänge angucken? ;)
Schlecht bei schwachem SNR und nicht subsamplegenau.

Cheers
Detlef

Ja, es scheint irgendwas quadratisches zurückzukommen. Die Phase hat mit 
dem Verhältnis der s_prev1 und s_prev2 zu tun. Code ist nicht 
übersetzbar und Zeug wie s_prev2 = s_prev1; s_prev1 = s; ist nicht 
akzeptabel, such Dir lieber ne andere Implementation. Englische 
Wikipedia erklärt Goertzel relativ erschöpfend.

Cheers
Detlef

Hallo Moppi,

durch eine Quadraturdemodulation bekommst Du abtastwertgenau die 
Amplituden- und Phaseninformation.

Eines Deiner Sinussignale scheint ja Dein Stimulus zu sein, nennen wir 
es S_out und dürfte eine sehr gute Qualität haben. Dein in Phase und 
Amplitude verändertes Detektionssignal nennen wir S_in. Für die 
Ermittlung der Amplituden- und Phasenänderung von S_in zu S_out 
benötigen wir noch ein um 90° verschobenes S_out, was ich mit 
S_out(-90°) bezeichne. Nun folgt die sogenannte Quadraturdemodulation:

x = S_out * S_in
y = S_out(-90°) * S_in

Da bei dieser Demodulation in beiden Pfaden auch ein Signalanteil mit 
doppelter Frequenz entsteht, ist ein Tiefpaßfilter zur Unterdrückung 
dieses Anteils nötig. Das kann ganz einfach als IIR-Filter 1.Ordnung 
implementiert werden (zwei Multiplikationen und eine Addition, bei 
geschickter Implementierung auch mit Bitshifts statt Multiplikationen 
realisierbar).

Da das Ergebnis nun allerdings kartesisch vorliegt (rechtwinklige 
Koordinaten x und y), mußt Du dann die Umrechnung in polare 
Werte(Amplitude und Phase) vornehmen (Hast Du einen Mikrocontroller oder 
DSP im Einsatz?):

Amplitude = wurzel(x^2+y^2)
Phase = artan(y/x).

Diese Auswertung ist Subsample-genau (nahezu beliebig, wenn bei der 
Rechnung aufgepaßt wird) und funktioniert auch mit extrem kleinen 
Signal-Rausch-Abständen.
Auf diese Weise arbeiten nämlich Funkempfänger sehr präzise mit sehr 
kleinen SNRs.

Viele Grüße!

Gerrit, DL9GFA

Hallo Moppi,

am besten kann man es in der numerischen Simulation zeigen und 
optimieren. Ich verwende gern Scilab dazu (http://www.scilab.org).

Anbei findest Du ein Scilab-Skript, das die nötige Mathematik umsetzt, 
sowie das Ergebnis als Grafik. Ich habe auf S_out einen Amplitudenfaktor 
von 0.01 und einen Phasenfehler von 30°angewendet, die auch prompt 
detektiert werden ;o). Mein Skript ist massiv dokumentiert; spiele doch 
einfach etwas damit herum. Dazu mußt Du Dir nur das freie Scilab 
installieren und im Dateimanager einen Doppelklick auf moppi01.sce 
durchführen. Mit <Strg>+<l> oder über Menü des gestarteten Editors 
kannst Du das Skript ausführen lassen.

In der Grafik ist das Einschwingen auf die korrekten Werte zu sehen. Man 
kann durch Vergrößerung des Filterkoeffizienten die Grenzfrequenz nach 
oben schieben und damit ein deutlich schnelleres Einschwingen erzielen. 
Die Unterdrückung der unerwünschten Mischprodukte wird dann aber 
geringer, was zu größeren Fehlern führt. Das hängt direkt von Deinem 
Verhältnis Abtastfrequenz zur Trägerfrequenz ab. Mit Filtern höherer 
Ordnung, die eine hohe Tiefpaß-Grenzfrequenz ermöglichen und trotzdem 
gleichzeitig eine hohe Unterdrückung des Mischproduktes bei der 
doppelten Frequenz ermöglichen, kann ein Einschwingen innerhalb weniger 
Abtastwerte erreicht werden. Zu dieser Optimierung mußt Du Dich dann 
aber intensiver mit digitalen Filtern auseinander setzen.

Ich hoffe das hilft. Viele Grüße!

Gerrit, DL9GFA

Hallo Gerrit, lang nix gelesen!


(Gast) Moppi =/= moppi. Btw.: irgendwer hat sich heut mein "vergessenes" 
Passwort zusenden lassen....
Schon seltsam.



73 de Jochen!

Hallo Jochen,

sorry, ich dachte ich hätte hier noch nen Account... .
Dass es dich mit dem gleichen Nickname gibt wusste ich nicht.... .

Grüße

Karsten

PS: schöner Nickname!

Ich bleibe mal beim Goertzel: DA komme ich einfach nicht weiter.
Hier mal mein Delphi-Code:

1

     amplitude:=StrToInt(Edit1.Text);  //Amplitude abfragen

2

     scan_freq:=StrToInt(Edit2.Text);  //Scan-Frequenz = Soll-Frequenz

3

     sample_freq:=StrToInt(Edit3.Text); //Sample-Frequenz

4

     N:=round(sample_freq/scan_freq);   //Anzahl Messwerte

5

6

     SetLength(soll,5*N+1);             //Arraylänge bestimmen

7

8

     for i:=0 to 5*N do

9

     begin

10

          soll[i]:=amplitude*sin(2*pi*i/N);  //5Sinus-Schwingungen erzeugen

11

     end;

12

13

14

    //goertzel!

15

     k:=0.5+(N*scan_freq/sample_freq);

16

     w:=(2*pi/N)*k;

17

     cosine:=cos(w);

18

     sine:=sin(w);

19

     coeff:=2*cosine;

20

     q1:=0;

21

     q2:=0;

22

     for i:=0 to 5*N do

23

     begin

24

          q0:=coeff*q1-q2+soll[i];

25

          q2:=q1;

26

          q1:=q0;

27

     end;

28

     real_value:=q1-q2*cosine;

29

     imag_value:=q2*sine;

30

     measured_amplitude:=sqrt(real_value*real_value+imag_value*imag_value);

Ist quasi die Implementierung nach 
http://www.embedded.com/story/OEG20020819S0057
Komischerweise bekomme ich nie die richtige Amplitude angezeigt... . 
Erkennt jemand einen Fehler?


Grüße


Karsten