3 pièce(s) jointe(s)
Revenir en temporel via IFFT depuis une dsp
Bonjour,
Je génère, d'après des données que je récupère via un autre logiciel, la PSD.
Code:
1 2 3
| [sigmaphi] = kourou; %appel fonction kourou
Sigmaphi_dsp = psd(spectrum.welch,sigmaphi); |
La variable sigmaphi contient mes données
Je sais grosso modo à quoi m'attendre, donc je plot voir si c'est correct.
Code:
1 2 3
| figure
loglog(Sigmaphi_dsp.Frequencies,Sigmaphi_dsp.Data)
grid on |
Ok, c'est ce à quoi je m'attends, pas de problème.
Maintenant, je voudrai revenir en temporel. Donc, comme mon signal est réel, je prends la racine carrée de ma DSP et j'en fait l'IFFT.
Code:
1 2 3 4 5
| %Repassage en temporel
TF_Sigmaphi = sqrt(Sigmaphi_dsp.Data);
Sigmaphi_temporel = ifft(TF_Sigmaphi); |
Je me retrouve donc en sortie avec un signal complexe et lorsque je plot soit la partié réelle, soit le module, je me retrouve avec :
1) un truc tout moche qui n'est pas ce à quoi je m'attends
2) une échelle en abscisse qui correspond à mon nombre de points de la FFT ou de mon IFFT (longueur du vecteur de départ) et pas une échelle temporelle.
Concernant mes données, je n'ai aucune donnée pouvant m'informer de quoi que ce soit au niveau temporel.
Ci joint les screenshots de mes figures.
Est-ce que mon raisonnement est erroné ou autre ?
Merci