输入信号序列和采样率,该子函数可以画出该信号的频谱图。
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
请发表评论