(精品word)基于MATLAB的IIR数字滤波器的设计.doc

(精品word)基于MATLAB的IIR数字滤波器的设计
(精品word)基于MATLAB的IIR数字滤波器的设计
(精品word)基于MATLAB的IIR数字滤波器的设计
基于MATLAB的IIR数字滤波器的设计
字体大小：小中大作者：武卫华来源：日期：2007-05—09点击：2210
摘要：介绍IIR数字滤波器的传统设计思想与步骤.及其计算机辅助设计方法。以一数字带通滤波器为例，着重说明了基于MATLAB的三种实现手段：模拟低通原型、合适模拟带通及直接原型，为数字滤波器设计带来全新的实现手段,设计快捷方便,仿真波形直观.数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号，实质上是一个由有限精度算法实现的线性时不变离散系统.它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。数字滤波器和模拟滤波器有着相同的滤波概念，根据其频率响应特性可分为低通、高通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外，还有滤波精度高(与系统字长有关）、稳定性好（仅运行在0与l两个电平状态）、灵活性强等优点.数字滤波器按单位脉冲响应的性质可分为无限长单位脉冲响应滤波器IIR和有限长单位脉冲响应滤波器(FIR)两种。本文介绍（IIR)数字滤波器的设计与分析。
1IIR数字滤波器设计思路与步骤IIR数字滤波器可用一个n阶差分方程y（n）=Σbrx（n—r)+Σaky（n-k）,或用它的Z域系统函数：对照模拟滤波器的传递函数：不难看出，数字滤波器与模拟滤波器的设计思路相仿,其设计实质也是寻找一组系数｛b，a｝,去逼近所要求的频率响应,使其在性能上满足预定的技术要求;不同的是模拟滤波器的设计是在S平面上用数学逼近法去寻找近似的所需特性H(S）,而数字滤波器则是在Z平面寻找合适的H(z).IIR数字滤波器的单位响应是无限长的，而模拟滤波器一般都具有无限长的单位脉冲响应，因此与模拟滤波器相匹配。由于模拟滤波器的设计在理论上已十分成熟,因此数字滤波器设计的关键是将H（S）→H(Z)，即,利用复值映射将模拟滤波器离散化。已经证明，冲击响应不变法和双线性变换法能较好地担当此任，则在此基础上，数字滤波器的设计就可首先归结为模拟滤波器的设计了。数字滤波器的设计步骤如图1所示。图1数字滤波器设计步骤
2IIR数字滤波器设计方法IIR数字滤波器的设计方法有多种，可归纳为下述两种。2。1传统设计方法根据前述设计思路,首先设计一个模拟原型滤波器(截止频率为1rad/s的低通滤波器），然后在模拟域(S平面)进行频率变换,将模拟原形滤波器转换成所需类型(指定截止频率的低通、高通、带通、带阻）的模拟滤波器，再将其数字离散化，从S平面映射至Z平面,得到所需技术指标的数字滤波器。上述过程中,也可先将模拟原型离散化，得到数字原型滤波器,继而在数字域（Z平面)进行频率变换,得到所需类型的数字滤波器.模拟滤波器到数字滤波器的转换可在时域进行也可在频域实现,时域转换的关键是要使数字滤波器与模拟滤波器时域响应的采样值相等,以保持其瞬态特性不变,常用的是冲击响应不变法.频域变换法必须使得数字滤波器在-π≤ω≤π范围内的幅频特性与模拟滤波器在—π/T≤Ω≤π/T范围内的幅频特性一致,即保证S平面与z平面上幅频特性的一一单值对应关系，常用的是双线性变换法。2。2计算机辅助设计方法传统设计方法思路清晰，步骤详尽，可参阅公式、手册循章而行。但由于计算繁琐，手工计算大多只能用来进行简单低阶选频滤波器（如LP,HP,BP及BS等）的设计.计算机辅助设计方法是集电路理论、网络图论、数值分析、矩阵运算、元件建模、优化技术、高级计算机语言等多交叉学科于一身的新领域，它把计算机的快速、高精度、大存储容量、严格的逻辑判断和优良的数据处理能力与人的思维创造能力充分结合起来，极大地简化了数字滤波器的设计过程。在优秀科技应用软件MATLAB的信号处理工具箱中，提供了一整套模拟、数字滤波器的设计命令和运算函数，方便准确，简单易行，使得设计人员除了可按上述传统设计步骤快速地进行较复杂高阶选频滤波器的计算、分析外,还可通过原型变换法直接进行各种典型数字滤波器设计，即应用MATLAB设计工具从模拟原型直接变换成满足原定频域指标要求的数字滤波器.此外，MATLAB软件还提供了强大的数字滤波器优化设计功能。即预先确定一种最佳设计准则,然后直接求得在该准则下滤波器系统的传递函数{b,a}.这种最优设计方法可方便地用于任意幅频特性要求的多带通复杂滤波器系统的设计。
3各种设计方法的MATLAB实现下面以一带通数字滤波器的设计为例,说明基于MATLAB的三种实现方法。3。1基于模拟低通原型的M