plot(abs(fft(ecg)))
私も試してみました
fvtool(x_vals)
私に与えた:
ただし、x 軸を Hz で表示したい。したがって、基本的には、この信号の周波数スペクトルを Hz で表示したいと考えています。
ありがとう!
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
これが役立つことを願っています。