《IIR滤波器设计》PPT课件.ppt

第 4 章 IIR滤波器设计,葛青,4.1滤波器的概念及分类,数字滤波器：通过对数字信号的运算处理，改变信号频谱，完成滤波作用的算法或装置。数字滤波器具有稳定、适应性强、线性相位等优点,避免了模拟滤波器固有的电压漂移、温度漂移和噪声等问题。,2.滤波原理,线性时不变系统 输入输出关系可用如下数学函数来描述 时域： 频域：,幅频特性：表示信号通过该滤波器后各频率成分的衰减情况 相频特性：反映各频率成分通过滤波器后在时间上的延时情况（ 联系无失真传输条件）,2.滤波原理,滤波原理,2.滤波原理,滤波器的设计就是用一个因果稳定的线性移不变系统的系统函数去逼近理想滤波器的性能。 数字滤波器设计:确定传输函数H(z)的过程。,滤波器的设计原理,表征滤波器频率响应的特征参量（3个参量） （1）幅度平方响应：幅度的平方与频率的关系 （2）相位响应 （3）群延迟响应（物理意义）:相位对角频率的导数的负值。若滤波器通带内群延迟响应=常数，则为线性相位滤波器，联系无失真传输问题。,无失真传输系统,如果系统响应仅是输入信号在时间上的延时和幅 度上的放大或缩小，就认为信号在传输过程中没 有失真，这种系统称为无失真传输系统。,无失真传输系统,3.滤波器的分类,经典滤波器：用于提取有用的信号，抑制不需要的干扰；或将信号分成不同的频率成分； 现代滤波：信号检测和信号估计,经典滤波器和现代滤波器,模拟和数字滤波器的理想幅频特性曲线,模拟滤波器幅频特性曲线,数字滤波器幅频特性曲线,IIR(infinite impulse response): 有限冲激响应滤波器,单位冲激响应h(n)是无限长序列 FIR(finite impulse response): 无限冲激响应滤波器,单位冲激响应h(n)是有限长序列,数字滤波器的分类,IIR系统与FIR系统,从系统函数的构造来区分 IIR系统：至少有一个极点。包括全极点系统（分子只有常数项）和零极点系统（分子不止常数项）；有反馈环路，采用递归型结构。 FIR系统：收敛域内无极点，是全零点系统。无反馈环路，多采用非递归结构。,IIR滤波器及FIR滤波器的系统函数,有限冲激响应滤波器的传输函数为 无限冲激响应滤波器的传输函数为,不全为零,4.数字滤波器的性能要求,一个理想滤波器，要求所在通频带内幅频响应是一常数；相位频率相应为零或是频率的线性函数。但一个实际的滤波器要是不可能得到上述幅频和相频响应。以低通滤波器为例，频率响应有通带、过渡带及阻带三个范围。,4.数字滤波器的性能要求,p：通带衰减,s：阻带衰减,p：通带截止频率 s：阻带截止频率,p-s ：过渡带,通带允许的最大衰减 阻带最小衰减：,4.数字滤波器的性能要求,式中假定 被归一化为1,按设计任务，确定滤波器性能要求，制定技术指标 一个因果稳定的离散LSI系统的系统函数H(z)逼近此性能指标 利用有限精度算法实现此系统函数：如运算结构、字长的选择等 实际技术实现：软件法、硬件法或DSP芯片法,数字滤波器的设计过程,数字滤波器的设计,IIR滤波器设计主要内容包括： 巴特沃思、切比雪夫模拟低通滤波器设计； 脉冲响应不变法和双线性变换法的数字化变换方法； 数字高通、带通和带阻滤波器的设计。 而FIR滤波器是直接采用的数字式设计方法。针对FIR滤波器特征，首先介绍了其线性相位的实现条件，然后介绍了窗函数法和频率抽样法的设计方法。,设计IIR数字滤波器一般有以下两种方法： 1. 模拟滤波器:首先设计一个合适的模拟滤波器，然后将它转换成满足给定指标的数字滤波器，这种方法适合于设计幅频特性比较规则的滤波器，例如低通、高通、带通、带阻等。 2.直接在频域或者时域中进行数字滤波器设计，由于要联立方程，设计时需要计算机作辅助设计。,4.2 IIR数字滤波器的设计,IIR滤波器以模拟低通滤波器为基础的设计方法，为了设计其他的选频滤波器（高通，带通，带阻等），需要对低通滤波器进行频率转换，在设计过程中有两种不同的变换,频带变换和模拟/数字变换。根据这两种变换的先后次序，引出两种设计方法。,4.2 IIR数字滤波器的设计,其它非低通数字滤波器的设计过程,IIR的设计方法借用模拟滤波器,先设计一个合适的模拟滤波器，然后变换成满足给定指标的数字滤波器。,设计步骤 将给定的数字滤波器的技术指标转换为模拟滤 波器的技术指标； 根据转换后的技术指标设计模拟原型滤波器； 按照一定规则将模拟滤波器转换为数字滤波器。,4.2.1 模拟低通滤波器的设计,常用模拟低通滤波器及其特点 1.巴特沃思(Butterworth)滤波器: 特点是从通带中心向两边幅频特性单调下降； 2.切比雪夫型I(Chebyshev)滤波器:通带内有等波纹幅频特性函数，阻带内有单调下降的幅频特性函数； 3.切比雪夫型II(Chebyshev)滤波器:通带内有单调下降的幅频特性函数，阻带内有等波纹幅频特性函数； 4.椭圆滤波器(又称考尔Cauer滤波器或双切比雪夫滤波器):通带和阻带内均有等波纹幅频特性； 5.贝塞尔(Bessel)滤波器（又称最大平坦群时延滤波器）:在中心频率附近逼近线性相位特性；,模拟滤波器的设计过程,（1）根据要求确定设计指标； （2）选择滤波器类型； （3）计算滤波器阶数； （4）通过查表或计算确定滤波器的系统 函数Ha(s); （5）综合实现并装配调试；,4.2.1 模拟低通滤波器的设计,1.巴特沃思低通滤波器 （1）基本性质 巴特沃思滤波器以巴特沃思函数来近似滤波器的系统函数。巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器。,N为滤波器阶数 c为3dB截止频率,巴特沃思滤波器的低通模平方函数表示,巴特沃思低通滤波平方幅频特性函数,巴特沃斯滤波器的主要特征 a.对所有的N:,即,b. 是 的单调下降函数。,c. 随着阶次N的增大而更接近于理想低通滤波器。,归一化原型滤波器是指截止频率 已经归一化成 的低通滤波器。对于截止频率为 某个低通滤波器，则令 替归一化原型滤波器系统函数中的 ， 对于其他高通、带通、带阻滤波器，可应用后面讨论到的频带变换法，由其变换得出。,巴特沃斯低通滤波器的幅度平方函数为：,N为滤波器的阶次， 为3dB截频。,(2)巴特沃思低通滤波器设计过程,通带最大衰减和阻带最小衰减可由定义式求出: (设=0处幅度归一化到1), 则:,(2)巴特沃思低通滤波器设计过程,1)求出N,Nmin应取向上取整。,2)如技术指标未给出 ，则可由下式计算： 或 采用上式求 ，阻带指标有富裕量； 下式则通带指标有富裕量。,查表求得归一化传输函数 , 令 代替归一化原型滤波器系统函数中 的 , 即得到实际滤波器传输函数Ha(s)。,3)求系统函数Ha(s),低通巴特沃斯滤波器设计步骤：,（1）由 求滤波器阶次N； （2）由N查表，求出归一化传输函数 ； （3）令 代入 ，得实际 滤波器传输函数Ha(s) 。,Nmin应取向上取整。,N及 计算,或,通带边频: 通带的最大衰减为: 阻带边频为: 阻带的最小衰减为:,例:设计一巴特沃思低通滤波器，使其满足以下 指标:,解:根据滤波器技术指标,取N=4, 查表得四阶巴特沃思多项式，得归一化系统函数表达式,确定阶次N,巴特沃斯归一化低通滤波器参数,巴特沃斯归一化低通滤波器参数,巴特沃斯归一化低通滤波器参数,代入,1.巴特沃思滤波器设计步骤：,(1)根据技术指标求 求出滤波器的阶数N和3dB截止频率,Nmin应取向上取整。,或,1.巴特沃思滤波器设计步骤：,(2)求系统函数Ha(s) 1)查表求 Ha(s) 查表求得归一化传输函数 ,令 代替归一化原型滤波器系统函数中的 , 得到实际滤波器传输函数。,求滤波器的极点，并由s平面左半平面的极点构成 。 幅度平方函数极点为： 将s平面左半平面的N个极点分配给,2)极点法求系统函数Ha(s),例：已知fp=5kHz，p=2dB，fs=12kHz， s=30dB，设计巴特沃斯低通滤波器。 解：（1）确定阶次N：,取N=5，其极点为（归一化极点）,（2）查表求极点： P1,2=-0.3090j0.9511; p3,4=-0.8090 j0.5857; p0=-1.000 所以,（3）求c： 所以 s比题目中给出的指标小，过渡带符合要求。,（4）去归一化，求Ha(S) 将p=s/c代入Ha(P)，得：,切比雪夫滤波器引入原因： Butterworth滤波器频率特性，无论在通带与阻带都随频率单调变化，因此如果在通带边缘满足指标，则在阻带内肯定会有富裕量，也就是会超过指标的要求，因此并不经济。 更有效的方法是将指标的精度要求均匀的分布在通带内，或均匀分布在阻带内，或同时均匀分布在通带与阻带内，这时就可设计出阶数较低的滤波器。这种精度均匀分布的办法可通过选择具有等波纹特性的逼近函数来完成。,4.2.2切比雪夫滤波器设计,通带内等波纹，阻带内单调-切比雪夫1型 通带内单调，阻带内等波纹-切比雪夫2型 切比雪夫1型的幅度平方函数：,是小于1的正数，称为纹波参数，是表示通带内纹波大小的一个参数， 愈大，纹波也愈大。 p称为有效通带截止频率。 频率通常对p归一化，,在通带内，在1和 之间起伏变化；在阻带内是单调下降的； 当N为奇数时，滤波器在 处的幅度响应为1；当N为偶数时，滤波器 在 处的幅度响应为 。当 时，,4.2.2切比雪夫滤波器设计,基本性质,的确定,设允许的通带纹波为，那么,阶数N 的确定,(1)查表得 (2)查表得幅度平方函数的极点,其中,的确定,极点分布,2N个等角度间隔（间隔为/N）分布在 为长半轴， 为短半轴的椭圆 上的点。,切比雪夫滤波器的传输函数,系数A由s=0时滤波器幅度响应的值确定: 当N为奇数时， 当N为偶数时，,去归一化得到滤波器的传输函数Ha(s),切比雪夫低通滤波器的设计步骤,(1)由给定的设计指标确定切比雪夫滤波器的参数、N和p,(2)确定Ha(s) 1)查表求得归一化传输函数 ,令 代替归 一化原型滤波器系统函数中的 ，即得到实际滤波 器传输函数。,切比雪夫低通滤波器的设计步骤,归一化Chebyshev低通滤波器的传输函数为,切比雪夫低通滤波器的设计步骤,归一化Chebyshev低通滤波器的传输函数为,为了得到单位峰值增益，一般选,切比雪夫低通滤波器的设计步骤,(2)确定Ha(s) 2）滤波器的极点，并由s平面左半平面的极点构成Ha(s)。,切比雪夫低通滤波器的设计步骤,系数A由s=0时滤波器幅度响应的值确定: 当N为奇数时， 当N为偶数时，,切比雪夫低通滤波器的设计步骤,去归一化得到滤波器的传输函数Ha(s),切比雪夫滤波器是一种全极型滤波器，它是由切比雪夫多项式的正交函数 推导出来的，采用了在通带内等波动，在通带外单调衰减的准则去逼近理想特性。 它比巴特沃斯滤波器在通带内有更均匀的特性，是在所有全极型滤波器中过渡带最窄的最优滤波器。,4.2.2切比雪夫滤波器设计,N=4,1,c,1/(1+2),(2)设计过程 1)根据要求的滤波器指标确定波纹参数和N。,滤波器的阶数N的确定：,的确定：,待定参数可通过通带内的最大衰减来求得,【例】 设计一个满足下列指标的Chebyshev低通滤波器。当 时， ；当 时， 。 【解】 取N3,归一化Chebyshev低通滤波器的传输函数为,归一化Chebyshev低通滤波器的传输函数为,为了得到单位峰值增益，一般选,2)查表求得归一化传输函数 ,令 代替归 一化原型滤波器系统函数中的 ，即得到实际滤波 器传输函数。,(2)设计过程,基本思路： 转换要求： （1）确保稳定性； （2）模仿的频响特性基本一致,4.3 IIR数字滤波器的实现,由模拟滤波器设计IIR数字滤波器,常用方法,工程上常用: 脉冲响应不变法 双线性变换法,1.变换原理,数字滤波器的单位脉冲响应序列h(n)正好等于模拟滤波器的冲激响应ha(t)的采样值，即 h(n)=ha(nT), T为采样周期.如以 Ha(s) 及 H（z）分别表示 ha(t) 的拉氏变换及 h(n) 的Z变换，即,数字滤波器的系统函数h(n)的Z变换便是H(z)。,4.3.1脉冲响应不变法( Impulse Invariance ),下面我们分析从模拟滤波器到数字滤波器S平面和Z平面之间的映射关系。,抽样信号的傅里叶变换,时域离散非周期信号,频域连续周期的频谱,抽样信号的傅里叶变换,数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。可能有混叠失真,-3S/2,冲激不变法s平面与z平面的映射关系,s平面到z平面的多值映射,但是，实际的模拟滤波器频率响应都不是严格限带的， 变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。 当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。,脉冲响应不变法中的频响混叠现象,1)滤波器指标以模拟域形式给出，此时,解决混叠的方法：,已确定，采样频率增加，令采样时间间隔（T=1/fs）减小混叠减小。,2)滤波器指标以数字域形式给出，此时,未确定，但,已定，采样频率增加，为保证omegas,不变，必有,增加，增加滤波器的阶数N，混叠减小。,脉冲响应不变法的基本思路,设模拟滤波器的系统函数若只有单阶极点，且分母的阶数高于分子阶数 NM，则可表达为部分分式形式 其拉氏反变换为 对ha(t)采样得到数字滤波器的单位脉冲响应序列,脉冲响应不变法模拟到数字的转换,再对h(n)取Z 变换，得到数字滤波器的传递函数,脉冲响应不变法模拟到数字的转换,（1）S平面的每一个单极点s=sk变换到Z平面上z=eskT处的单极点。 （2） Ha(s)与H(z)的部分分式的系数相同，都是Ak。 （3）如果模拟滤波器是因果稳定的，则所有极点sk位于S平面的左半平面，即Resk0， 则变换后的数字滤波器 的全部极点在单位圆内，即|eskT|=eReskT1，因此数字滤波器也是因果稳定的。,Ha(s)和H(z) 比较,练习：设模拟滤波器的系统函数为,试用脉冲响应不变法将Ha(s)转换成IIR数字滤波器的系统函数H(z)。,解 直接利用公式可得到数字滤波器的系统函数为,设T=1，则有,冲激响应不变法的实现步骤 (1)确定数字滤波器的性能指标（ ) 和( ) (2)利用 ，将数字滤波器的性能指标变换 成模拟滤波器的性能指标( )和( ）(3)利用成熟的模拟滤波器设计方法设计模拟滤波器，求出 (4)将模拟滤波器系统函数 转变成数字滤波器的系统函数,4.冲激响应不变法的优缺点,优点：一个线性相位的模拟滤波器可以映射成为一个线性相位的数字滤波器。 缺点：根据奈奎斯特抽样定理，冲激响应不变法只适用于限带的模拟滤波器，即只适用低通、带通滤波器的设计，不适合高通、带阻滤波器的设计，否则会出现混叠，而且 这种频率混叠现象的出现不会因抽样周期的减少而消失。,例: IIR低通滤波器的设计指标如下： 1）通带截止频率P0.1rad， 阻带起始频率s0.25rad 2）通带最大衰减p=3dB， 阻带最小衰减s=15dB,设计举例,解：,1) 根据DF指标, 将DF指标转换为归一化LPF指标,1.巴特沃思滤波器设计步骤：,(1)根据技术指标求 求出滤波器的阶数N和3dB截止频率,Nmin应取向上取整。,或,4) 求H(z),（1）频率坐标变换是线性的，数字滤波器在 无频率混叠时能较好地重现原模拟滤波器的频率特性； （2）由于数字域的冲激响应完全模仿模拟域的冲激响应，故时域特性逼近好； （3）存在频率效应，对于非带限的模拟滤波器不适用。,脉冲响应不变法的优缺点,-3S/2,4.3.2 双线性变换法,（ Bilinear Transformation method ） 为了克服脉冲响应不变法的频谱混叠 现象，采用非线性频率压缩方法 -双线性变换法。,从s到z平面的映射： 映射时应满足的条件： 如原模拟滤波器是因果稳定的，则 数字滤波器也应是因果稳定的， 并且性能指标符合设计要求。,4.3.2双线性变换法（Bilinear Transform),1.变换原理,s平面到z平面的映射关系 二次映射法S S1 Z 为了将S平面的j轴压缩到S1平面j1轴上的 到 一段上，可通过以下的正切变换实现：,这里C是待定常数，下面会讲到用不同的方法确定C 。,当 由 时, 由 经过变化到 ，即S平面的整个 轴被压缩到S1平面的 一段. 将这一关系解析扩展至整个S平面， 则得到S平面到S1平面的映射关系： 再将 S1 平面通过标准变换关系映射到Z平面，即令,令,(1)整个s平面到s1平面上的转换为：,(2)s1平面到z平面的变换则为：,(3)从s平面到z平面的变换为：,或,（）,（）,4.3.2 双线性变换法,2.变换常数的选择,常数的选择可以使模拟滤波器的频响特性和数字滤波器的频响特性在不同的频率范围有对应的关系，起到调节二者频带间关系的作用。选择的方法有两种。,1)使模拟滤波器和数字滤波器的频响特性在低频部分有较确切的对应关系，即当 较小的时候，有,由此得,2)使数字滤波器的某一特定频率(例如截止频率 与模拟原型滤波器的特定频率 严格对应，即,2.变换常数的选择,由于在待定的模拟频率和待定的数字频率处频率响应应严格相等，因而可以较准确的控制截止频率位置。,3.双线性变换可以消除频率混叠的原因,S/2,-S/2,如令,则得,3. 模拟角频率和数字角频率的映射关系,3. 模拟角频率和数字角频率的映射关系,双线性变换的频率对应关系,双线性变换法虽然避免了“频率混叠效应”，但出现了模拟频率与数字频率为一种非线性的关系情形即：,可见：模拟滤波器与数字滤波器的响应在对应的频率关系上发生了“畸变”，也造成了相位的非线性变化，这是双线性变换法的主要缺点。具体而言，在上刻度为均匀的频率点映射到上时变成了非均匀的点，而且随频率增加越来越密。,双线性变换法除了不能用于线性相位滤波器设计外，仍然是应用最为广泛的设计IIR数字滤波器的方法。,（1）确定数字低通技术指标： 通带截止频率 、通带衰减 阻带截止频率 、阻带衰减 （2）预畸变：将数字低通指标转换成模拟低通指标（ 和 不变 ） 边界频率的变换关系：,4.双线性变换法设计数字低通滤波器的步骤,4.双线性变换法设计数字低通滤波器的步骤,（3）设计模拟低通滤波器传输函数Ha(s)； （4）将Ha(s)转换成数字低通滤波器： 这里的采样间隔T可任意选取， 通常取 T=1。,例：用双线性变换法将图示RC低通滤波器转换成数字低通滤波器。,解： 传输函数,数字滤波器的传递函数为：,例2：,设计数字低通滤波器，要求在通带内频率 低于0.2时，容许幅度误差在1dB以内； 在频率0.3到之间的阻带衰减大于15dB。 解：（1）数字低通指标,（2）模拟低通指标：,（3）设计巴特沃斯模拟低通滤波器：,查表后，得,得模拟滤波器传输函数Ha(s),取N=6,通带指标有富裕量,去归一化,令,代入,（4）用双线性变换法转换成数字低通滤波器,例3:,用双线性变换法和模拟低通巴