matlab程序加窗

wlp=0.35*pi;whp=0.65*pi;wc=[wlp/pi,whp/pi];N=ceil(8/0.15);n=0:N-1;window=hanning(N);[h1,w]=freqz(window,1);subplot(411);stem(window);%画窗函数axis([06001.2]);grid;xlabel('n');title('Hanning窗函数');subplot(412);plot(w/pi,20*log(abs(h1)/abs(h1(1))));%画Hanning窗频谱图axis([01-3500]);grid;xlabel('w/pi');ylabel('幅度(dB)');title('Hanning窗函数的频谱');hn=fir1(N-1,wc,hanning(N));%基于窗函数的FIR滤波器[h2,w]=freqz(hn,1,512);%freqz是用离散傅里叶变换的标准公示计算的fft使用快速傅里subplot(413);stem(n,hn);axis([060-0.250.25]);grid;xlabel('n');ylabel('h(n)');title('Hanning窗函数的单位脉冲响应');subplot(414);plot(w/pi,20*log(abs(h2)/abs(h2(1))));grid;xlabel('w/pi');ylabel('幅度(dB)');title('滤波器相对幅度响应')结果：