Discussion :

[angelmn]

bonjour à tous
cà fait longtemps que je peine à réaliser la fft d'un signal avec matlab, est ce que quelqu'un peut m'aider svp??
j'ai un signal semi périodique avec une période de 294 point et une fréquence d'echantillonnage de 360 hz. Donc mon signal a une fréquence de 1.2244 hz si j'ai bien compris.
alors j'essai de realiser la fft de ce signal sauf que je tombe sur des resultats complètement bizarres. je ne trouve pas de pic à f0 ni une période fixe séparant les autres pics.
voici mon code

Nx=length(x);
>> tf=fft(x,Nx);
>> plot(abs(tf));

Merci de m'aider

[FLB]

As tu bien pris en compte ton échelle?
En reprenant la doc:

Code :

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

1
2
3
4
5
6
7
8
9
10
doc fft;
Nx=length(x);
Fs = 360;
NFFT = 2^nextpow2(Nx); % Next power of 2 from length of y
Y = fft(x,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
plot(f,2*abs(x(1:NFFT/2+1))) 
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

[angelmn]

Je vous remercie pour votre réponse; j'ai essayé de voir dans le help mais ce que je ne comprends toujours pas c'est le choix de l’échelle et la signification de ses valeurs ; comment reconnaître à quoi correspondent les abscisses des pics est ce que ça correspond à la valeur de f0 ou pas.
voici ce que j'ai obtenu je ne sais pas si c'est juste ou pas

et voici mon signal


et ici le code que j'ai utilisé

Code :

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

1
2
3
4
5
6
7
8
9
 
Fs = 360;                    % Sampling frequency
T = 1/Fs;                     % Sample time
L = length(x);
y=x;
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
plot(f(10:length(f)),2*abs(Y(10:NFFT/2+1)))

Merci de m'aider