【信号处理】基于窗函数实现数字信号频谱分析附matlab代码

1 简介

窗函数法是设计数字滤波器的一种方法，本文采用各种矩形窗、汉宁窗、哈明窗、布莱克曼窗等不同窗形式对滤波器性能的影响.最后运用Matlab进行仿真比较.​

2 完整代码

function y = plot_freq(Fs,L,x,s,i) 

% 绘制频域幅度波形图

% Fs：采样频率

% L：序列长度

% x：幅值

% s：标题名称

% i：标识序号

figure(i);

y_fft = fft(x);

P2 = abs(y_fft/L);

P1 = P2(1:floor(L/2)+1);

P1(2:end-1) = 2*P1(2:end-1);

f = Fs*(0:(L/2))/L;

plot(f,P1);

title(s);

xlabel('f (Hz)');

ylabel('|X(k)|');

axis([0 3000 0 0.005]); 

%{

y_fft = (fft(x)/L*2); % 幅度图 /(N/2) 恢复真实幅值

y2 = abs(fftshift(y_fft)); % 平移

f = (Fs * (0:L-1)