数字信号实验报告4.doc

南京信息工程大学 实验（实习）报告 实验（实习）名称用窗函数法设计FIR数字滤波器日期2012.6.1 指导老师 行鸿彦 院（系）电信 专业电信 年级09级 班次2班 姓名徐宏岭学号20091305090 一、实验目的 1、掌握用窗函数法设计FIR数字滤波器的原理和方法。 2、熟悉线性相位FIR数字滤波器的特性。 3、了解各种窗函数对滤波特性的影响。 二、 实验原理 如果所希望的滤波器的理想频率响应函数为（），则其对应的单位脉 冲响应为 （4－1） 窗函数设计法的基本原理时用有限长单位脉冲响应序列h(n)逼近。用窗函数将截断，并进行加权处理，得到： （4－2） 就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为 （4－3） 式中，N为所选窗函数的长度。 由书本第七章可知，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。 这样选定窗函数类型和长度N后，求出单位脉冲响应，并按式（4－3）求出。是否满足要求，要进行验算。一般在h(n)尾部加零使长度满足2的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果补满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。 如果要求线性相位特性，则h(n)还必须满足： 根据上式中的正负号和长度N的奇偶性又将线性相位FIR滤波器分为四类。要根据所设计的滤波器特性正确选择其中一类。例如，要设计线性相位低通特性，可选择一类，而不能选择一类。 三、 实验内容 1、用升余弦窗设计一线性相位低通FIR数字滤波器，截止频率。窗口长度N＝15，33。要求在两种窗口长度情况下，分别求出h(n),打印出相应的幅频特性和相频特性曲线，观察3dB带宽和20dB带宽。总结窗口长度N对滤波特性的影响。 设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函数，即 其中 n＝33，，用四种窗函数设计线性相位低通滤波器。绘制相应的 幅频特性曲线，观察3dB和20dB带宽以及阻带最小衰减，比较四种窗函数对滤波器特性的影响。 实验结果 1.N=15时的矩形窗： N=33时的矩形窗： 2.N=15时的汉明窗： N=33时的汉明窗： N=15时的汉宁窗： N=33时的汉宁窗： N=15时的布莱克曼窗： N=33时的布莱克曼窗： 1、矩形窗的频率响应主瓣宽度为4π/N，第一副瓣比主瓣低13dB。 2、哈明窗是改进的升余弦窗，能量更加集中在主瓣中，主瓣的能量约占99.96%，第一旁瓣的峰值比主瓣小40dB,但主瓣宽度和汉宁窗相同，为8π/N。 3、汉宁窗的幅度函数由三部分相加，使能量更集中在主瓣中，但代价是主瓣宽度加宽到8π/N。 4、布莱克曼窗的幅度函数由五部分组成，它们都是移位不同，且幅度也不同的函数，使旁瓣再进一步抵消，阻带衰减进一步增加，过渡带是矩形窗过渡带的3倍。 5、调整窗口长度N可以有效低控制过渡带的宽度，减少带内波动以及加大阻带的衰减只能从窗函数的形状上找解决方法，如果能找到的窗函数形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就减小了，旁瓣的减小可使通带、阻带波动减小，从而加大阻带衰减，但这样总是以加宽过渡带为代价的。 五、实验代码 主程序： s = -1; while(s0) clc; N=input(请输入窗函数长度 N=); s=1; end close all; i=0;wc=pi/4; while(s) n=0:N-1; hd=ideal(wc,N); k=input(请选择窗口类型:\n[1](boxcar)\n[2](hamming)\n[3](hanning)\n[4](blackman)\n请选择：,s); k=str2num(k); if(k==1) B=boxcar(N); string=[Boxcar,N=,num2str(N)]; els