第六章 IIR滤波器的设计方法

DSP黄海哈尔滨理工大学2014年5月*6.1引言数字滤波器： 是指输入输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。高精度、稳定、体积小、重量轻、灵活，不要求阻抗匹配，可实现特殊滤波功能优点：*1、滤波器的基本概念（1）滤波器的功能滤波器的功能是对输入信号进行滤波以增强所需信号部分，抑制不要的部分。a）时域说明b）频域说明*（2）四种基本的滤波器四种基本滤波器为低通（LP）、高通（HP）、带通（BP）和带阻滤波器（BRF）：*（3）四种基本滤波器的数字 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示低通高通带通带阻*2、滤波器的性能指标 带宽：当幅度降低到0.707时的宽度称为滤波器的带宽（3dB带宽）* 通带、阻带与过渡带：信号允许通过的频带为通带，完全不允许通过的频带为阻带，通带与阻带之间为过渡带。* 滚降与滚降率：滤波器幅频特性在过渡带的衰减和衰减速度称为滚降与滚降率。* 阻带衰减：输入信号在阻带的衰减量* 带内平坦度：通带和阻带内的平坦程度*3、数字滤波器的设计步骤数字滤波器的设计四个步骤:(1)按要求确定滤波器的性能参数；(2)用一个因果稳定的离散线性移不变系统的系统函数去逼近去逼近这一性能要求；(3)用有限精度的算法来实现这个系统函数，包括选择运算结构（第五章的基本结构）、选择合适的字长（输入、中间变量、输出）、有效数字的处理方法（舍入、截尾）；(4)实现可以采用通用计算机，也可以采用DSP。*4、数字滤波器的技术要求 选频滤波器的频率响应：*理想滤波器不可实现，只能以实际滤波器逼近**5、表征滤波器频率响应的特征参量 幅度平方响应H(z)的极点：单位圆内的极点* 相位响应* 群延迟响应相位对角频率的导数的负值*6、IIR数字滤波器的设计方法 先设计模拟滤波器，再转换为数字滤波器用一因果稳定的离散LSI系统逼近给定的性能要求： 计算机辅助设计法s平面逼近：模拟滤波器z平面逼近：数字滤波器*用模拟滤波器设计IIR数字滤波器 设计思想： H(z)的频率响应要能模仿Ha(s)的频率响应， 即s平面的虚轴映射到z平面的单位圆 因果稳定的Ha(s)映射到因果稳定的H(z)， 即s平面的左半平面Re[s]<0 映射到z平面的单位圆内|z|<1* 设计方法：-冲激响应不变法-阶跃响应不变法-双线性变换法*冲激响应不变法一、变换原理数字滤波器的单位冲激响应h(n)模仿模拟滤波器的单位冲激响应ha(t)T—抽样周期*设则：*从频率响应来看： 数字滤波器的频率响应是模拟滤波器频率响应的周期延拓，周期为2p/T 只有当模拟滤波器的频率响应是带限的，且带限于折叠频率以内时，即 才能使数字滤波器的频响在折叠频率以内重现模拟滤波器的频响而不产生混迭失真* 实际系统不可能严格限带，都会混迭失真，在|W|>Ws/2处衰减越快，失真越小 当滤波器的设计指标以数字域频率wc给定时，不能通过提高抽样频率来改善混迭现象*二、模拟滤波器的数字化**当T很小时，数字滤波器增益很大，易溢出，需修正令：则：*试用冲激响应不变法，设计IIR数字滤波器例：设模拟滤波器的系统函数为解：据题意，得数字滤波器的系统函数：设T=1s，则*模拟滤波器的频率响应：数字滤波器的频率响应：* 优点： 缺点： 保持线性关系：w=WT 线性相位模拟滤波器转变为线性相位数字滤波器 频率响应混迭 只适用于限带的低通、带通滤波器 h(n)完全模仿模拟滤波器的单位抽样响应ha(t)时域逼近良好冲激响应不变法的优缺点*阶跃响应不变法 变换原理数字滤波器的阶跃响应g(n) 模仿模拟滤波器的阶跃响应ga(t) T—抽样周期** 阶跃响应不变法同样有频率响应的混叠失真现象但比冲激响应不变法要小。*双线性变换法冲激响应不变法、阶跃响应不变法：时域模仿逼近缺点是产生频率响应的混叠失真为了克服这一缺点，采用双线性变换法。使数字滤波器的频率响应与模拟滤波器的频率响应相似*一、变换原理及特点脉冲响应不变法的映射是多值映射,导致频率响应交叠。改进思路：先将s域平面压缩到一个中介平面s1，然后再将s1映射到Z平面。** 为使模拟滤波器某一频率与数字滤波器的任一频率有对应关系，引入系数c*2）某一特定频率严格相对应：1）低频处有较确切的对应关系： 特定频率处频率响应严格相等，可以较准确地控制截止频率位置二、变换常数c的选择*三、逼近情况1）2）*四、优缺点 优点：避免了频率响应的混迭现象s平面与z平面为单值变换* 缺点：除了零频率附近，W与w之间严重非线性2）要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变分段常数型模拟滤波器经变换后仍为分段常数型数字滤波器，但临界频率点产生畸变*预畸变 给定数字滤波器的截止频率w1，则按W1设计模拟滤波器，经双线性变换后，即可得到w1为截止频率的数字滤波器*五、模拟滤波器数字化方法 可分解成级联的低阶子系统* 可分解成并联的低阶子系统*常用模拟低通滤波器特性 将数字滤波器技术指标转变成模拟滤波器技术指标，设计模拟滤波器，再转换成数字滤波器 模拟滤波器 巴特沃斯Butterworth滤波器 切比雪夫Chebyshev滤波器 椭圆Ellipse滤波器 贝塞尔Bessel滤波器*1、由幅度平方函数确定模拟滤波器的系统函数h(t)是实函数将左半平面的的极点归Ha(s)将以虚轴为对称轴的对称零点的任一半作为Ha(s)的零点，虚轴上的零点一半归Ha(s)Ha(s)Ha(-s)的零极点分布* 由幅度平方函数得象限对称的s平面函数*例：解：极点：零点：（二阶）零点：的极点：设增益常数为K0*1）巴特沃尔斯滤波器（Butterworth）2、常见模拟滤波器设计幅度平方函数：当称Wc为Butterworth低通滤波器的3分贝带宽N为滤波器的阶数Wc为通带截止频率*1）幅度函数特点：当W＝Wst（阻带截止频率）时，衰减的d1为阻带最小衰减*Butterworth滤波器是一个全极点滤波器，其极点：2）幅度平方特性的极点分布：* 极点在s平面呈象限对称，分布在Buttterworth圆上，共2N点 极点间的角度间隔为 极点不落在虚轴上 N为奇数，实轴上有极点，N为偶数，实轴上无极点Ha(s)Ha(-s)的零极点分布(a)N=3(三阶）(b)N=4（四阶）*3）滤波器的系统函数：为归一化系统的系统函数去归一化，得*4）滤波器的设计步骤： 根据技术指标求出滤波器阶数N： 确定技术指标：由得：同理：令则：* 求出归一化系统函数： 或者由N，直接查表得其中技术指标Wc给出或由下式求出： 其中极点： 去归一化阻带指标有富裕或通带指标有富裕* 例：设计Butterworth数字低通滤波器，要求在频率低于0.2prad的通带内幅度特性下降小于1dB。在频率0.3p到p之间的阻带内，衰减大于15dB。分别用冲激响应不变法和双线性变换法。1、用冲激响应不变法设计1）由数字滤波器的技术指标：2）得模拟滤波器的技术指标：选T=1s*a）确定参数 用通带技术指标，使阻带特性较好，改善混迭失真3）设计Butterworth模拟低通滤波器*b)求出极点（左半平面）c)构造系统函数c’)去归一化*4）将Ha(s)展成部分分式形式:变换成Butterworth数字滤波器：*用冲激响应不变法设计出的Butterworth滤波器*2、用双线性变换法设计1）由数字滤波器的技术指标：2）考虑预畸变，得模拟滤波器的技术指标：*a）确定参数 用阻带技术指标，使通带特性较好，因无混迭问题3）设计Butterworth模拟低通滤波器*b)求出极点（左半平面）c)构造系统函数*c’)去归一化*4）将Ha(s)变换成Butterworth数字滤波器：*2）切贝雪夫滤波器（Chebyshev）N：滤波器的阶数Wc：截止频率，不一定为3dB带宽0<e<1，表示通带波纹大小，e越大，波纹越大CN(x)：N阶Chebyshev多项式TypeIChebyshev* N为偶数 N为奇数*Chebyshev滤波器的三个参量： Wc：通带截止频率，给定 e：表征通带内波纹大小 N：滤波器阶数，等于通带内最大最小值的总数由通带衰减决定阻带衰减越大所需阶数越高Ws为阻带截止频率*TypeIIChebyshevfilter通带内：单调特性阻带内：等波纹起伏* 例：用双线性变换法设计Chebyshev数字低通滤波器，要求在频率低于0.2prad的通带内幅度特性下降小于1dB。在频率0.3p到p之间的阻带内，衰减大于15dB。1）由数字滤波器的技术指标：2）考虑预畸变，得模拟滤波器的技术指标：*a）确定参数3）设计Chebyshev模拟低通滤波器*b)求左半平面极点*c)构造系统函数*c’)去归一化*设计的四阶Chebyshev滤波器*椭圆滤波器(Ellipticfilter) 带内均匀波动 最快的滚降*贝塞尔滤波器（Bessel）*最大相位平坦特性*小结：利用模拟滤波器设计IIR数字滤波器的步骤 将数字滤波器的技术指标转变成模拟滤波器的技术指标通带截止频率阻带截止频率通带截止频率阻带截止频率 确定数字滤波器的技术指标： 冲激响应不变法 双线性变换法* 按模拟滤波器的技术指标设计模拟低通滤波器 Butterworth低通滤波器 Chebyshev低通滤波器 将模拟低通滤波器转换成数字低通滤波器 冲激响应不变法 双线性变换法**********************************************************************************