MATLAB中的FFT：錯誤的0Hz頻率

Question

我想在MATLAB中使用fft來分析保存爲excell文件的一些exprimental數據。 我的代碼：

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) 

但它繪製了大錯爲0Hz分量的圖，我的真實頻率是大約的395Hz，它沒有在繪製的圖所示。 請告訴我什麼是錯的。

任何幫助，將不勝感激。

Answer 1

假設我們從文件中讀取信號：

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

首先我們估計從時間軸上的取樣頻率，並生成頻率向量：

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

然後計算FFT和情節：

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

終於找到峯值和相應的頻率：

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

您可能希望縮小在原點附近，看東西更清楚：

xlim([-1000 1000])