Discussion :

[Shagocode]

Bonjour,

je suis en train d'essayer d'extraire la phase d'un signal pour avoir la fréquence en fonction du temps, pour simplifier mon problème je vais mettre un cas équivalent en plus simple qui donne le même résultat que je ne comprends pas :

- J'ai un signal sinusoidal de fréquence 5MHz (que j'appelle testo),
- Je calcule la fft de ce signal
- j'en prends la phase avec "angle" (phaseff)
- je fais la ifft de ça pour repasser en domaine temporel (phaseft)
- comme la phase est l'intégrale de la fréquence, je dérive ce résultat pour trouver la fréquence (GFfilt1 )

-> le problème : si je plot GFfilt1 je m'attends à avoir les 5MHz de mon signal et donc une droite horizontale d'ordonnée 5e6 en fonction de t, or ce n'est pas du tout le cas, j'ai juste un courbe qui varie aux alentours de 0, pourtant en regardant abs(fft) de mon signal j'ai bien un pic à 5MHz

-> parallèlement à ça j'ai fait la transfo de Hilbert pour en extraire l'angle aussi, ça marche pas mieux.

voici le code, quelqu'un voit le soucis ? je pense que je m'y prend mal quelque part mais je sais pas où..

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
fs = 1e9;
N = 100000;
t = 1/fs*[1:N];
t = t';
testo = sin(2*pi*5e6*t);
 
Lf = length(testo);
NFFTf = 2^nextpow2(Lf);
Ytest = fft(testo,NFFTf);
f = fs/2*linspace(0,1,NFFTf/2+1);
ampff=abs(Ytest);
plot(ampff);
plot(f,2*abs(Ytest(1:NFFTf/2+1)));
phaseff = unwrap(angle(Ytest)); %phase domaine freq
 
Xtest=ifft(phaseff,NFFTf);
phaseft = unwrap(angle(Xtest)); %phase domaine temp
 
Hfilt1 = hilbert(testo);
Hphase = angle(Hfilt1);
 
GHfilt1 = gradient(Hphase);
GFfilt1 = gradient(phaseft);
 
 
 figure
plot(GHfilt1);hold on;plot(GFfilt1,'r');

[tachmou]

Salut,

Je ne m'y connais en rien... Cependant, tu peux regarder par là.

Ensuite, je ne sais pas d'où tu sors que la fréquence est la dérivée de la phase, mais pour moi, c'est n'importe quoi... sauf d'un point de vue unité, ok.

A+