【MATLAB】画信号频谱的子函数

function [f,spec,NFFT]=spec_fft_plot(sample,L,Fs)
% 输入数据说明：
% sample：信号序列；
% L：信号序列的长度；
% Fs：该信号的采样频率。
% 输出数据说明：
% f：频率；
% spec：频谱图纵坐标；
NFFT = 2^nextpow2(L); 
% NFFT = L;
spec=abs(fft(sample,NFFT)./L);
flag = 1;
if flag ==1
    f = Fs/2*linspace(0,1,NFFT/2+1);
    spec=spec(1:length(f));
else
    f = Fs/2*linspace(-1,1,NFFT);
    spec = fftshift(spec);
end