MATLAB中的FFT：错误的0Hz频率

【问题标题】：FFT in MATLAB: wrong 0Hz frequencyMATLAB中的FFT：错误的0Hz频率
【发布时间】：2012-06-15 19:19:37
【问题描述】：

我想在 MATLAB 中使用 fft 来分析一些保存为 Excel 文件的实验数据。我的代码：

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，并且未显示在绘制的图表中。请告诉我有什么问题。

任何帮助将不胜感激。

【问题讨论】：

【解决方案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])

【讨论】：