MATLAB中心电信号的FFT

Question

这是输入信号：

    plot(abs(fft(ecg)))

我也试过了

    fvtool(x_vals)

这给了我：

但是我想要以赫兹为单位的 x 轴。所以本质上我想看到这个信号的频谱，以赫兹为单位。

谢谢！

Answer 1

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

我希望这可能会有所帮助。