搜索
1、五种窗函数FIR数字滤波器比较N=31;wc=1/2;n=0:30;w1n=boxcar(N);k=0:1024-1;w=(pi/1024)*k;%N=2048X=w1n'*(exp(-j*pi/1024)).^(n'*k);magX=20*log10(abs(X)/max(abs(X)));figure(1);subplot(2,2,1);stem(n,w1n,'.');xlabel('n');ylabel('w1(n)');title('矩形窗时域波形（N=31）');subplot(2,2,2);plot(w/pi,magX);title('矩形窗的幅度特性');xlabel('\omega/\pi');ylabel('20lg|W1(e^j^\omega)|');gridon;%axistight;h1n=fir1(N-1,wc,w1n);[H,w]=freqz(h1n,1,1024);magH=20*log10(abs(H)/max(abs(H)));subplot(2,2,3);stem(n,h1n,'.');xlabel('n');ylabel('h1(n)');title('FIR滤波器的单位脉冲响应（N=31）');subplot(2,2,4);plot(w/pi,magH);title('FIR滤波器的幅度特性');xlabel('\omega/\pi');ylabel('20lg|Hg(e^j^\omega)|');gridon;矩形窗的四种波形描绘程序如下：三角窗的四种波形描绘程序如下：汉宁窗的四种波形描绘程序如下：哈明窗的四种波形描绘程序如下：布莱克曼窗的四种波形描绘程序如下：※结果分析：矩形窗主瓣比较集中，旁瓣较高；三角窗主瓣宽度约等于矩形窗的两倍，但旁瓣小；汉宁窗主瓣加宽并降低，旁瓣则显著减小；哈明窗旁瓣更小，旁瓣衰减速度比汉宁窗衰减速度慢；布莱克曼窗主瓣宽，旁瓣较低。2、改进升余弦窗（哈明窗）的FIR数字滤波器设计figure(1)Window=hamming(33);b=fir1(32,0.25,Window);freqz(b,1)figure(2)Window=hamming(21);b=fir1(20,0.25,Window);freqz(b,1)figure(3)Window=hamming(15);b=fir1(14,0.25,Window);freqz(b,1)※结果分析：N越大，阻带衰减越大。