数字信号处理课程设计之带阻滤波器.doc

课程设计报告 (2013-- 2014年度第一学期) 名 称: 数字信号处理 题 目：带阻FIR数字滤波器设计 院 系： 动力系 班 级： 测控11K2 学 号： 学生姓名： 指导教师： 白康 设计周数： 一周 成 绩： 日期：2014年01月13日 《数字信号处理》课程设计 任 务 书 一、 目的与要求 能够运用本课程中学到的知识，设计基于窗口函数法的FIR 数字滤波器。要求掌握数字信号处理的基本方法；FIR滤波器的设计步骤和方法；能够熟练采用C 语言或MATLAB语言进行计算机辅助设计和仿真验证设计内容的合理性。 二、 主要内容 掌握设计数字滤波器的基本步骤； 重点掌握利用窗口函数法设计FIR数字滤波器的设计方法，不同窗函数对滤波器滤波性能的影响以及滤波器单位冲激响应长度对滤波器延时特性的影响； 能够利用傅立叶变换在理论上分析滤波器的频率响应； 掌握计算机C语言或MATLAB的编程技巧； 三、 进度计划 序号 设计(实验)内容 完成时间 备注 1 熟悉利用窗口函数法设计FIR带阻滤波器的方法和步骤，针对具体题目进行设计 2天 2 编制程序 2天 3 调试， 答辩 1天 四、 设计（实验）成果要求 设计报告 五、 考核方式 答辩+平时表现 学生姓名： 指导教师：白康 20014年1月13 日 数字信号处理课程设计 一、设计目的与要求 能够运用本课程中学到的知识，设计基于窗口函数法的FIR 数字滤波器。要求掌握数字信号处理的基本方法；FIR滤波器的设计步骤和方法；能够熟练采用C 语言或MATLAB语言进行计算机辅助设计和仿真验证设计内容的合理性。 二、正文 1．设计题目： 采用窗口函数法设计一个带阻FIR数字滤波器，要求设计的频率响应为 2．设计步骤 （1）根据设计要求，确定FIR滤波器的频率响应，包括截止频率、过渡带、阻带最小衰减： 由设计题目可得知滤波器的通带0.6通带宽度0.4和0，通带宽度0.2，过渡带宽度0.1和，过渡带宽0.1，阻带宽度0.2，阻带最小衰减为。 （2）根据前一步的截止频率、过渡带、阻带最小衰减，确定选择的窗口函数类型 和的列长度： 由于阻带最小衰减要求为-26dB，根据窗函数基本参数可查的汉宁窗的阻带最小衰减为-44dB,符合要求且在所有符合要求中最简,因此选择汉宁窗.同样可以查的汉宁窗的过渡带宽为8/N.所以0.1=8/N,N=8/(0.1)=80,取N=81. 同样可以查得汉宁窗 （3）根据第一步确定理想的频率响应，通过求的逆傅里叶变换得到，对进行右移、截短和窗口函数处理，求得因果FIR数字滤波器的单位冲激响应 : 由已知的条件可以得出理想得频率响应 。 右移、截短 . 因此因果FIR数字滤波器的单位冲激响应＝* （4）求得傅里叶变换，检验是否符合要求； . 3.各部分设计结果和MATLAB的仿真波形如下: 1).理想 2).汉宁窗w(n): 3).实际， 4).幅度响应(db): 由此可以得知设计的滤波器阻带最小衰减为-44db,验证了汉宁窗的阻带最小衰减为-44db,同时说明该设计的滤波器在阻带最小衰减的性能指标满足设计要求的-26db. 检验: 1).f1=sin(0.1*pi*n) 2).f2=2*sin(0.8*pi*n) 3).f3=4*sin(0.4*pi*n) 4)叠加后的图像f4=+f1+f2+f3 5).通过滤波器后的图像f4*h(n) 6).f5=f1+f3 分析：通过滤波器后，只有f1和f3的图形能通过滤波器，所以f5=f1+f3的图像应该和f4*h(n)的图像差不多。 三、收获和体会： 通过本次带阻FIR数字滤波器的设计，加深了对数字滤波器的了解，同时对其的设计方法及各种性能指标都有一定的了解。基本掌握了MATLAB的用法 ，可以根据给定的条件写出需要的程序，得到理想的波形。在课程设计的初期，由于对matlab软件的不熟悉，在算完数据编辑程序运行的时候，会出现很多的问题，例如程序的错误，多或少写一行程序，多加一个括号之类的，导致程序不能运行，不能输出波形，还有就是X轴和Y轴长度的选取不合适，导致输出波形之后无法完全显示或者出现重叠的现象，这些都是在课程设计时常遇到的问题，通过查找相关资料，请教老师和同学，我一一解决了这些困难，