matlabで振幅スペクトルの表示(fft)

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