用巴特莱特窗函数法设计数字FIR带通滤波器设计.doc

课 程 设 计 课程设计名称：数字信号处理课程设计 专 业 班 级 ： 电信 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 乔丽红 课程设计时间： 6.16-6.20 电子信息工程 专业课程设计任务书 学生姓名 专业班级 学号 题 目 用巴特莱特窗函数法设计数字FIR带通滤波器 课题性质 其他 课题来源 自拟 指导教师 乔丽红 同组姓名 主要内容 用巴特莱特窗函数法设计一个数字FIR带通滤波器，要求通带边界频率为400Hz，500Hz，阻带边界频率为350Hz，550Hz，通带最大衰减1dB，阻带最小衰减40dB，抽样频率为2000Hz，用MATLAB画出幅频特性，画出并分析滤波器传输函数的零极点； 信号经过该滤波器，其中450Hz，600Hz，滤波器的输出是什么？用Matlab验证你的结论并给出的图形。 任务要求 1．掌握用巴特莱特窗函数法设计数字FIR带通滤波器的原理和设计方法。 2．掌握用Kaiser方程估计FIR数字滤波器长度的方法。 3．求出所设计滤波器的Z变换。 4．用MATLAB画出幅频特性图并验证所设计的滤波器。 参考文献 1、程佩青著，《数字信号处理教程》，清华大学出版社，2001 2Sanjit K. Mitra著，孙洪，余翔宇译，《数字信号处理实验指导书（MATLAB版）》，电子工业出版社，2005年1月 3郭仕剑等，《MATLAB 7.x数字信号处理》数字信号处理 理论算法与实现数字信号处理是把许多经典的理论体系作为自己的理论基础，同时又使自己成为一系列新兴学科的理论基础。FIR数字滤波器的主要优点有：一、具有严格的线性相位特性；二、不存在稳定性问题；三、可利用DFT来实现。这些优点使FIR数字滤波器得到了广泛应用。窗函数法是一种设计FIR数字滤波器的基本方法，但它不是最佳设计方法，在满足同样设计指标的情况下，用这种方法设计出的滤波器的阶数通常偏大。在窗函数法的基础上，以所定义的逼近误差最小为准则来进行优化设计的算法，由于其中的逼近误差可根据不同的设计要求进行定义，故此算法适应性强，它可用于设计选频型滤波器，又适用于非选频型滤波器的设计常用的窗函数有矩形窗函数、三角窗函数、汉宁（Hann）窗函数、海明（Hamming）窗函数、布莱克曼（Blackman）窗函数、凯塞（Kaiser）窗函数等。FIR系统不像IIR系统那样易取得较好的通带和阻带衰减特性,要取得较好的衰减特性,一般要求H（z）阶次要高,也即M要大。FIR系统有自己突出的优点：系统总是稳定的;易实现线性相位;允许设计多通带（或多阻带）滤波器,后两项都是IIR系统不易实现的。FIR数字滤波器的设计方法有多种,如窗函数设计法、频率采样法和Chebyshev逼近法等。随着Matlab软件尤其是Matlab的信号处理工作箱的不断完善,不仅数字滤波器的计算机辅助设计有了可能,而且还可以使设计达到最优化。经过该滤波器，其中450Hz，600Hz，滤波器的输出是什么？用Matlab验证你的结论并给出的图形。 二　设计原理及设计思路 1.设计FIR数字滤波器的基本方法，系统频率响应可写成：，令=，H(w)为幅度函数，称为相位函数。这与模和辐角的表示方法不同，H(w)为可为正可为负的实数，这是为了表达上的方便。如某系统频率响应 =sin4w,如果采用模和幅角的表示方法，sin4w的变号相当于在相位上加上（因-1= ），从而造成相位曲线的不连贯和表达不方便，用则连贯而方便。窗函数法又称傅里叶级数法，其设计是在时域进行的。?一般是无限长且非因果的，设计时需用一个合适的窗函数把截成有限长的因果序列，使对应的频率响应（的傅里叶变换）尽可能好地逼近理想频率响应。窗函数法的主要缺点是：一、不容易设计预先给定截止频率的滤波器；二、满足同样设计指标的情况下所设计出的滤波器的阶数通常偏大。 4/(2M+1) 0.92/M 20.9dB 巴特利特 4/(2M+1) 2.1/M 25dB 汉宁 8/(2M+1) 3.11/M 43.9dB 海明 8/(2M+1) 3.32/M 54.5dB 布莱克曼 12/(2M+1) 5.56/M 75.3dB 　FIR数字滤波器设计的基本步骤如下（1）确定技术指标在设计一个滤波器之前,必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中,数字滤波器常被用来实现选频操作。因此,指标的形式一般在频域中给出幅度和相位响应。幅度指标主要