数字滤波器研究课程论文..docVIP

研 究 生 课 程 论 文 (2010-2011学年) 数字滤波器的研究 研究生： 提交日期： 研究生签名： 学 号 学 院 电子与信息学院 课程编号 S0809074 课程名称 现代数字信号处理 学位类别 硕士 任课教师 教师评语： 成绩评定： 分 任课教师签名： 年 月 日 中文题目：数字滤波器的研究 研究生姓名 摘要：无限冲激响应（IIR）数字滤波器和有限冲激响应（FIR）数字滤波器广泛应用于数字信号处理系统中。FIR 滤波器具有很好的线性相位特性，因此线性相位特性在工程实际中如数据通信、图像处理等应用领域具有比较重要的意义。IIR 数字滤波器具结构简单和实现容易的特点, 这对于相位要求不是特别严格的场合有很重要的意义。利用MATLAB 实现IIR 低通数字滤波器的设计,并对被干扰的信号进行滤波,达到了预期结果。 关键词: IIR；FIR；线性相位；数字滤波器；MATLAB 1　引言 数字滤波器是现代数字信号处理的不可或缺的一部分，数字滤波器按其响应形式可分为无限脉冲响应（IIR）滤波器和有限脉冲响应（FIR）滤波器和两种。当系统无严格相位要求时，IIR 数字滤波器与FIR 数字滤波器相比，可用较低的阶数获得较高的选择性[1]。对于相同的设计指标，FIR 滤波器所要求的阶数比IIR滤波器高5～10 倍，而且信号的延迟也较大[2]。MATLAB 作为一种矩阵运算为基础的交互式程序语言，是进行科学研究常用且必不可少的工具。MATLAB着重针对科学计算、工程计算和绘图的需求。它用解释方式工作，键入程序立即得出结果，人机交互性能好，为科学人员所乐于接受。MATLAB 提供了丰富的函数，其中fir1 函数实现了加窗线性相位FIR 滤波器设计的经典方法，fir1 主要用于常用的标准通带滤波器设计，包括低通、带通、高通和带阻数字滤波器[3]。一般来讲，数字滤波器的设计要经过三步：确定设计指标、模拟逼近和数字转换。通常在设计滤波器之前，先根据具体的应用确定一些技术指标，然后就可以根据数学知识和滤波器的基本原理提出一个滤波器的模型来逼近给定的指标，逼近的结果通常是得到以差分方程或脉冲响应描述的滤波器，最后可以根据这个描述用硬件或软件实现，至此完成一个滤波器设计的全过程[3]。 2.数字滤波器的设计 2.1模拟椭圆滤波器对加噪语音信号的滤波： 设计模拟滤波器的一般方法是：首先根据技术指标确定滤波器的传输函数H(s） （2） 其中，是模拟的频率，是数字的归一化频率，T是积分步长。 用公式（2）预畸所求数字滤波器的的数字频率指标（数字滤波器通带截止频率0.2π，阻带截止频率0.6π），得到一个等价的模拟低通滤波器的频率指标。 3： 得到原型模拟低通滤波器的指标后作频率响应，从而得模拟的传输函数H(s)，然后对传输函数进行双线性变换，从而得到数字滤波器的传输函数H(z). 结果分析： 图3.模拟原型与数字滤波器幅度相应比较 图4. 模拟原型与数字滤波器相位相应比较 由图1可以看出模拟原型低通滤波器的通带截止频率为10KHz对应着数字低通滤波器的截止频率0.2π rad与双线性变换结果一致。由图2可以看出，相位响应就有很大差别了，这是由双线性变换的频率畸变所致。 2.3 FIR滤波器的结构与设计 对于IIR滤波器必须明确要求推导出来的传输函数是稳定的，而另一方面，对于FIR滤波器的设计，由于滤波器的传输函数是以1 /z的多项式表示的，所以FIR数字滤波器都是稳定的[11]。 利用傅里叶加窗级数法设计线性相位低通FIR滤波器（没有使用到fir1函数） 设计指标：通带截止频率4rad/s，阻带截止频率6rad/s，最大通带衰减为0.2dB,最小阻带衰减为42dB，抽样率为18rad/s. 图5.由截短得到的理想FIR低通滤波器的冲击响应 图6.各种窗函数设计的低通滤波器 结果分析： Hamming窗函数的长度是30 通带波纹是Rp = 0.0394dB 最小阻带衰减时As = 52dB Hanning窗函数的长度是29 通带波纹是Rp= 0.0711dB 最小阻带衰减时As = 44dB Blackman窗函数的长度是52 通带波纹是Rp= 0.0027dB 最小阻带衰减时As = 75dB 均满足滤波器的性能指标。 结论 在电子系统中，具有严格的线性相位特性的FIR 数字滤波器被广泛使用，是信号处理的基本构成之一。利用matlab的强大运算功能，基于matlab信号处理工具箱（signal processing toolbox）的