用窗函数法设计数字滤波器.docVIP

用窗函数法设计FIR数字滤波器 一、试验目标 1.掌握用窗函数法设计FIR数字滤波器原理和方法。 2.熟悉线性相位FIR数字滤波器特征。 3.了解多种窗函数对滤波特征影响。 二、试验仪器 微型计算机 matlab软件 三、试验原理和方法 假如所期望滤波器理想频率响应函数为 ，则其对应单位脉冲响应为 = （2-1） 窗函数设计法基础原理是用有限长单位脉冲响应序列迫近。因为往往是无限长序列，且是非因果，所以用窗函数将截断，并进行加权处理，得到： ＝ (2-2) 就作为实际设计FIR数字滤波器单位脉冲响应序列，其频率响应函数为： ＝ (2-3) 式中，N为所选窗函数长度。 由第七章可知，用窗函数法设计滤波器性能取决于窗函数类型及窗口长度N取值。设计过程中，要依据对阻带最小衰减和过渡带宽度要求选择适宜窗函数类型和窗口长度N。多种类型窗函数可达成阻带最小衰减和过渡带宽度见第七章。 这么选定窗函数类型和长度N后，求出单位脉冲响应＝·，并按式（2-3）求出。是否满足要求，要进行验算。通常在尾部加零使长度满足于2整数次幂，方便用FFT计算。假如要观察细节，补零点数增多即可。假如不满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。 假如要求线性相位特征，则还必需满足 (2-4)依据上式中正负号和长度N奇偶性又将线性相位FIR滤波器分成四类。要依据设计滤波特征正确选择其中一类。比如，要设计线性低通特征，可选择一类，而不能选一类。 四、试验内容 1.复习用窗函数法设计FIR数字滤波器一节内容，阅读本试验原理掌握设计步骤。 2.编写程序 ① 编写能产生矩型窗、哈明窗、汉宁窗、莱克曼窗窗函数子程序。 ② 编写主程序。主程序框图图 开始 开始 读入窗口长度N 计算hd（n） 调窗函数子程序 求ω（n） 计算h（n）=hd（n）*ω（n） 调FFT子程序 对h（n）进行DFT 调绘图子程序绘制 幅频特征曲线 相位特征曲线 结束 计算幅度特征 和相位特征 设 ： (2-4) (2-5) (2-6) 画图时，20lg打印幅度特征。第k点对应频率。为使曲线包络更 靠近幅度特征曲线，DFT变换区间要选大些。比如窗口长度N=33时，可经过在末尾补零方法，使长度变为64，再进行4点DFT，则可得到更正确幅度衰减特征曲线。 3.上机试验内容 用四种窗函数设计线性相位低通FIR数字滤波器，截止频率rad，N=256。绘 制对应幅频特征曲线，观察3dB带宽和20dB带宽和阻带最小衰减，比较四种窗函数对滤波器特征影响。 设计低通FIR数字滤波器时，通常以理想低通滤波特征为迫近函数，即 其中 (2-8) (2-9) 五、思索题 1.假如给定通带截止频率和阻带截止频率和阻带最小衰减，怎样用窗函数法设计线性相位低通滤波器？写出设计步骤。 2.定性说明本试验中3dB截止频率理论值在什么位置？是否等于理想低通截止频率ω？ 3.假如要求用窗函数法设计带通滤波器，且给定上、下边带截止频率为和，试求理想带通单位脉冲响应。 六、试验汇报要求 1.简述试验目标及原理。 2.根据试验步骤及要求，比较多种情况下滤波性能，说明窗口长度N和窗函数类型对滤波特征影响。 3.总结用窗函数法设计FIR滤波器关键特点。 4. 简明回复思索题。 参考程序1： N=input(窗宽度N=);k=input(窗型：1.矩形窗,2.hanning(升余弦窗),3.hamming（改善升余弦窗）,4.Blackman 请选择：);subplot(2,2,1);w=pi/5;a=(N-1)/2;n=0:(N-1);m=n-a+eps;h=sin(w*m)./(pi*m);if k==1??? B=bartlett(N);else if k==2??????? B=hanning(N);??? else if k==3??????????? B=hamming(N);??????? else if k==4??????????????? B=blackman(N);