БПФ в MATLAB: неправильная частота 0 Гц - PullRequest
Новая
       16

БПФ в MATLAB: неправильная частота 0 Гц

1 голос
/ 28 мая 2010

Я хочу использовать fft в MATLAB для анализа некоторых экспериментальных данных, сохраненных в виде файла 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)

Но он строит график с большим неправильным компонентом 0 Гц, моя истинная частота составляет около 395 Гц, и он не показан на графике. Пожалуйста, скажите мне, что не так.

Любая помощь будет оценена.

1 Ответ

8 голосов
/ 28 мая 2010

Предположим, мы читаем сигнал из файла: 

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;

Затем вычислите БПФ и построите график: 

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

наконец найдите пик и соответствующую частоту: 

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

Возможно, вы захотите приблизиться к исходной точке, чтобы лучше видеть вещи: 

xlim([-1000 1000])
...