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])