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

Question

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

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