高通FIR数字滤波器.docVIP

课程设计报告 (2010-- 2011年度第一学期) 名 称: 数字信号处理 题 目：高通FIR数字滤波器设计 院 系： 动力系 班 级： 测控08K1 学 号： 学生姓名： 指导教师： 金秀章 设计周数： 一周 成 绩： 日期：2011年1月7日 《数字信号处理》课程设计 任 务 书 一、 目的与要求 本次课程设计是使学生们加深对数字信号处理课程内容的理解而开设的。滤波器作为数字信号处理的一个重要分支，是测控专业学生必须掌握的内容。通过FIR数字滤波器的设计，使学生掌握滤波器和一般数字信号处理系统的设计步骤及软/硬件实现手段，为以后的工作和科学研究打下坚实的基础。 二、 主要内容 掌握设计数字滤波器的基本步骤； 重点掌握利用窗口函数法设计FIR数字滤波器的设计方法，--、不同窗函数对滤波器滤波性能的影响以及滤波器单位冲激响应长度对滤波器延时特性的影响； 能够利用傅立叶变换在理论上分析滤波器的频率响应； 掌握计算机C语言的编程技巧； 三、 进度计划 序号 设计(实验)内容 完成时间 备注 1 熟悉利用窗口函数法设计FIR滤波器的方法和步骤，针对具体题目进行设计 2天 2 编制程序 2天 3 调试， 答辩 1天 四、 设计（实验）成果要求 设计报告 五、 考核方式 答辩 学生姓名：杨欣 指导教师：金秀章 2011 年01月7 日 数字信号处理课程设计 一、 目的与要求 本次课程设计是使学生加深对数字信号处理课程内容的理解而开设的．滤波器作为数字信号处理的一个重要分支，是测控专业学生必须掌握的内容。通过FIR数字滤波器的设计，使学生掌握滤波器和一般数字信号处理系统的设计步骤及软/硬件实现手段，为以后的工作和科学研究打下坚实的基础。 二、正文 1．设计题目： 用窗口函数设计一个高通FIR数字滤波器，要求设计的频率响应为： 2．设计步骤 （1）根据设计要求，确定FIR滤波器的频率响应，包括过渡带宽、通带和阻带阻尼限制： 由设计题目可得知滤波器的通带通带宽度0.6，过渡带宽度0.1，阻带宽度0.3，阻带最小衰减为 。 （2）根据前一步的过渡带宽、通带和阻带阻尼限制，确定选择的窗口函数类型和的列长度： 由于阻带最小衰减要求为-34dB，根据窗函数基本参数可查的汉宁窗的阻带最小衰减为-44dB,符合要求且在所有符合要求中最简,因此选择汉宁窗.同样可以查的汉宁窗的过渡带宽为8/N.所以0.1=8/N,N=8/(0.1)=80,取N=81. 同样可以查得汉宁窗. （3）根据第一步确定理想的频率响应，通过求的逆傅里叶变换得到，对进行右移、截短和窗口函数处理，求得因果FIR数字滤波器的单位冲激响应 : 由已知得条件可以得出理想得频率响应 ； 。 右移、截短 . 因此因果FIR数字滤波器的单位冲激响应＝* （4）求得傅里叶变换，检验是否符合要求； . 各部分MATLAB的仿真曲线如下: a). b).汉宁窗: c).实际 . d)对数幅频相频响应: 由此可以得知设计的滤波器阻带最小衰减为-44db,验证了汉宁窗的阻带最小衰减为-44db,同时说明该设计的滤波器在阻带最小衰减的性能指标满足设计要求的-34db. e).频率特性: 由频率特性可以清晰的看出过渡宽为0.1,而且当时,最大,当时,最大,符合设计要求的. （5）设计硬件电路或软件程序实现所设计的FIR数字滤波器。 (6）验证滤波效果:输入一混合信号频率分别为0.6pi和0.2pi; x=sin(0.6*pi*n)+sin(0.2*pi*n) 可以看出低频部分被虑去,高频部分完好的保存,同时可以得到群延迟了(N-1)/2=40. 三、分析总结： 通过本次高通FIR数字滤波器的设计，加深了对数字滤波器的了解，同时对其的设计方法及各种性能指标都有一定的了解。由阻带最小衰减确窗函数的类型，由过渡带宽确定窗函数的点数N，由理想频率响应通过傅立叶逆变换得到和hd(n)截断右移，由窗函数特性确定的w(n),加窗后得到实际h(n)=hd(n)*w(n) 。最后验证各项指标（包括阻带最小衰减，过渡带宽，通带阻带以及截至频率）。 四、附件: MATLAB程序代码: 理想hd(n) n=0:1:80 y=-sin(0.35*