Discussion :

[cpalperou]

Bonjour,
je suis un absolute beginner en traitement du signal!!
Je dois extraire la phase du signal reçu par une antenne. Ce signal étant initialement émis par une autre antenne émettrice.
Le signal reçu est sous la forme d'un fichier .wav
Je fais ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
fln='myaudiofile.wav';
[Y,sample_rate]=wavread(fln);
Yfft=fft(Y);
Yfft_real=real(Yfft);
Yfft_ima=imag(Yfft);
phi=atan(Yfft_imag./Yfft_real);

Hormis la distinction entre les 4 cadrans du plan complexe (qui impose d'ajouter ou de retrancher pi selon le signe des parties réels ou imaginaires), est-ce correct?

le truc qui me paraît étrange est que phi(2)="complexe conjuguée de phi(end)" , phi(3)="complexe conjugué de phi(end-1)", phi(4)="complexe conjugué de phi(end-2)", ...........

Ne devrais je pas plutôt utiliser la fonction angle(Yfft)?

bien à vous

[cpalperou]

Bonjour,
je me permets d'insister un peu. Personne pour m'aider?
Merci

[cpalperou]

re,
J'ai compris que la fft est symétrique.
Maintenant, j'ai besoin de la phase de mon signal en fonction du temps. Mon souci est que dès que je fais fft, je passe dans le domaine spectrale et plus dans le domaine temporel
Je fais ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
fln='myaudiofile.wav';
[y,sample_rate]=wavread(fln);
Lf=length(y);
T=1/sample_rate;
t=(0:Lf-1)*T;
yfft=fft(y,Lf);
yfft=yfft(1:Lf/2);
phi=angle(yfft);
f=sample_rate*(0:Lf-1);

Alors, tracer la phase en fonction de la fréquence: OK! mais en fonction du temps, est-ce correct?
Bien à vous

[Lord Bouillon]

Salut, oui la représentation de la phase sur le temps est tout a fait légitime. Dans ton premier code atan(Q/I) représente un détecteur de phase et non la phase elle même.

Mais de manière globale, tu n'as pas besoins de passer par une fft à chaque fois pour obtenir la phase, voici un code très simple qui te le prouve:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
sigRead  = wavread("/home/matlab/test.wav");
signal = (complex(sigRead(:,1), sigRead(:,2))).';
 
Fnb = 1024;     %coef du filtre
filtBW = 150/192000;  % BW/fs
 
lpf = fir1(Fnb,filtBW);   %filtre fir1
y=filter(lpf, 1,signal) ; 
phase=angle(y);

En testant ça avec une wave j'obtiens le diagramme suivant:


Alors je ne suis pas précisément en fonction du temps mais en fonction du nombre de samples, mais cela n'est pas bien grave tu à juste à diviser par la fréquence d’échantillonnage pour retomber sur ton temps. La phase est vraiment pas belle car je travail sur des waves avec un SNR faible donc le bruit me fait bien varier la phase et c'est pour ça que tu vois des piques ainsi qu'un saut de phase qui missmatch.

Cela étant ces problèmes sont justes spécifiques à ma wave et c'est d'ailleurs pour cela que je l'ai filtré avec un filtre fir1. Si ta wave n'est pas trop bruité tu pourras t'en passer ! Bon courage !

[cpalperou]

Bonjour, merci Lord et désolé de répondre un peu tard (vacances!)
Je comprends bien votre réponse mais je n'arrive pas à l'appliquer car la sortie de mon wavread est un vecteur de nombres réels de dimension (168815*1). Ainsi, je ne peux pas faire votre 2ème commande!