FFT dans MATLAB: mauvaise fréquence 0Hz

Question

Je veux utiliser fft dans MATLAB pour analyser certaines données expérimentales sauvegardées dans un fichier excell. mon code:

A=xlsread('Book.xls'); G=A'; x=G(2, :); 
N=length(x); 
F=[-N/2:N/2-1]/N; 
X = abs(fft(x-mean(x),N)) 
X = fftshift(X); 
plot(F,X) 

Mais il trace un graphique avec un grand mauvais composant de 0 Hz, ma vraie fréquence est d'environ 395Hz et il ne figure pas dans le graphique tracé. S'il vous plaît dites-moi ce qui ne va pas.

Toute aide serait appréciée.

Answer 1

Supposons que l'on lit le signal du fichier:

G = xlsread('Book.xls'); 
t = G(:,1); 
x = G(:,2); 
N = length(x); 

D'abord, nous estimons que la fréquence d'échantillonnage de l'axe du temps, et de construire le vecteur de fréquence:

Fs = 1 ./ abs(t(2)-t(1)); 
F = (-N/2:N/2-1)*Fs/N; 

Calculons alors la FFT et l'intrigue:

X = abs(fft(x-mean(x),N)); 
X = fftshift(X); 
stem(F,X) 

enfin trouver le pic et la fréquence correspondante:

>> [~,ind] = max(X); 
>> F(ind) 
ans = 
     -400 

vous voudrez peut-être effectuer un zoom avant près de l'origine pour voir les choses plus clairement:

xlim([-1000 1000])