Discussion :

[lanfou]

Bonjour,

Je génère, d'après des données que je récupère via un autre logiciel, la PSD.

Code :

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

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 :

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

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 :

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

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