数字信号处理实验七FIR滤波器的设计.docVIP

学生实验报告 开课学院及实验室 电子楼317 2013年05月27日 学院 机械与电气 工程学院 年级、专 业、班 姓名 学号 实验课程名称 数字信号处理实验 成绩 实验项目名称 实验七FIR滤波器的设计 指导 老师 一、实验目的 复习巩固窗函数法设计FIR滤波器的概念和方法，观察用儿种常用窗函数设计的FIR数宁滤波 器技术指标。 实验原理 FIR滤波器的设计问题在于寻求一系统函数H(z),使其频率响应H(en逼近滤波器要求的理 想频率响应Hd(ei(t))，其对应的单位脉冲响应為(町o 用窗函数设计FIR滤波器的基本方法 设计思想从时域从发，设计瓜砒逼近理想hd (n) o设理想滤波器叽(严)的单位脉冲响应为 hd (n)。以低通线性相位FIR数字滤波器为例。 乙(严)=￡鵝(心 心(力=丄『￡ (严严de 2兀— 打⑺)一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。耍想得到一 个因果的有限长的滤波器h(n),最直接的方法是截断〃(/?) = hd(n)wM ,即截取为有限长因果序列, 并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的耍求，h(n)必须 是偶对称的。对称中心必须筹于滤波器的延时常数，即 (h(n) = hd(n)w(n) [a = (N-l)/2 用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象， 且最大波纹大约为幅度的9%,这个现象称为吉布斯(Gibbs)效应。为了消除吉布斯效应，一般采用 其他类型的窗函数。 典型的窗函数 (1)矩形窗(Rectangle Window) w(〃)= Rn(7?) 其频率响应和幅度响应分别为: / = sin(NQ/2)F 咛 sin ⑹ / 2) (2)三角形窗 (2)三角形窗(Bartlett Window) vv(n)= vv(n)= 2/1 N-1 2-^- N_\ N-l 0n——- 2 N-l nN-l 2 其频率响应为叭宀諾[將笋宀 哈宁(Hanning)窗，又称升余弦窗 ⑷is(笔叽(“) 其频率响应和幅度响应分别为: 叫?｛0.5WZ)+ 0.25叫3涪)+ W +語川「沙 %)5初+ 0.25叽定严如語)] 汉明(Hamming)窗，又称改进的升余弦窗 w(n) = [0?54-0.46cos(2^)]Rn (n) N-l 其幅度响应为W(a)) = 0.54% (cd) + 0.23[W? {co一 - ) + 叫(q+仝 )] TOC \o 1-5 \h \z N-] N 布莱克曼(Blankman)窗，又称二阶升余弦窗 2/7 4-77 7T w(n) = [042 - 0.5 cos( ) + 0.08 cos( )]/?v (n) N-l N-l 1 Itt 2tt W (ty) = 0.42% (CD) + 0.25 卩竹⑹-亍〒)+ 肌(ty+—-)] 其幅度响应为 NT N-\ + 0.04[叫(e - 丁〒)+ 肌⑷+-—)] N N-I (6)凯塞(Kaiser)窗 丽,加1-[1-2]/(叱理2,0“纫-1 其屮B是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的交换关系，一般说来，B越大，过渡 带越宽，阻带越小衰减也越大。Io(?)是第一类修正零阶贝塞尔函数。 若阻带最小衰减表示为人=-201og10①，P的确定可采用下述经验公式: ° 人5 21 0 = J 0.5842（人 一 21严 +0.07886（4—21） 21 A$ 50 0.1102（人一8.7） 50 若滤波器通带和阻带波纹相等即6 p= 6 s时，滤波器阶数可通过下式确定: 人-7.95 14.36AF 式中:空旦 式中: 空旦』厂叭 171 2 兀 三、使用仪器、材料 1、 硬件计算机 2、 软件:Matlab 四、实验步骤 知识准备 在实验编程之前，认真复习有关FIR滤波器设计的有关知识，尤其是窗函数的有关内容，阅读本 次实验指导，熟悉窗函数及四种线性相位FIR滤波器的特性，掌握窗函数设计滤波器的具体步骤。 编制窗函数设计FIR滤波器的程序。绘制滤波器的幅频和相位曲线，注意长度N对曲线的影 响。至少选两种不同长度，每种长度选两种窗函数。 设计一线性相位FIR数字低通滤波器，截止频率?=0.2龙，过渡带宽度AOV0.4%,阻带 衰减 40dBo 设计一线性相位FIR数宁高通滤波器，要求通带截止频率叫=0.5兀,阻带截止频率 ? = 0.25〃，通带最大衰减舛=\dB ,阻带最小衰减4 = 40dB。 对一段白噪声信号分别用上述设计好的滤波器进行滤波，分析滤波以后信号的频谱。 对一段频率在滤波器通带范围内的正弦波进行滤波，对比滤波前后信号的时域特性。 五、实验过