2

これは入力信号です: オリジナル

    plot(abs(fft(ecg)))

fft

私も試してみました

    fvtool(x_vals)

私に与えた:

別

ただし、x 軸を Hz で表示したい。したがって、基本的には、この信号の周波数スペクトルを Hz で表示したいと考えています。

ありがとう!

4

1 に答える 1

2
function [f amp] = getspectrum( Mdata, Mf )

%  Mdata    data 
%  Mf       sampling rate / frequency (Hz)

NFFT = 2 ^ nextpow2(length(Mdata)); 
Y = fft(double(Mdata), NFFT) / length(Mdata);
f = (double(Mf) / 2 * linspace(0, 1, NFFT / 2))'; % Vector containing frequencies in Hz
amp = 2 * abs(Y(1:(NFFT / 2))); % Vector containing corresponding amplitudes

これが役立つことを願っています。

于 2012-12-16T15:28:25.200 に答える