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