matlabで高速フーリエ変換し,振幅スペクトルを表示する.
function showfft(sig,fs) %sig:信号.1次元配列 %fs:信号のサンプリング周波数 N=length(sig); % 信号のサンプル数 plot(fs*(0:N-1)/N,abs(fft(sig)/N)) %振幅スペクトルの描画:縦軸が振幅,横軸が周波数 title('amplitude spectrum') xlabel('frequency [Hz]') ylabel('amplitude [a.u.]') end