《《数字信号处理》第六章 iir滤波器的设计方法》由会员分享，可在线阅读，更多相关《《数字信号处理》第六章 iir滤波器的设计方法（122页珍藏版）》请在金锄头文库上搜索。
1、第六章 IIR滤波器的设计,主要内容,理解数字滤波器的基本概念 了解最小相位延时系统 理解全通系统的特点及应用 掌握冲激响应不变法 掌握双线性变换法 掌握Butterworth、Chebyshev低通滤波器的特点 了解利用模拟滤波器设计IIR数字滤波器的设计过程 了解利用频带变换法设计各种类型数字滤波器的方法,6.1 引言,数字滤波器：,是指输入输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。,高精度、稳定、体积小、重量轻、灵活，不要求阻抗匹配，可实现特殊滤波功能,优点：,1、滤波器的基本概念,（1） 滤波器的功能滤波器的功能是对输入信号进行滤波
2、以增强所需信号部分，抑制不要的部分。,a） 时域说明 b） 频域说明,（2） 四种基本的滤波器,四种基本滤波器为低通（LP）、高通（HP）、带通（BP）和带阻滤波器（BRF）：,（3） 四种基本滤波器的数字表示,低通高通带通带阻,2、LP到其他滤波器的变换,由LP实现的HP,LP实现的BP,LP实现的BRF,3、 滤波器的性能指标,带宽：当幅度降低到0.707时的宽度称为滤波器的带宽（3dB带宽）,通带、阻带与过渡带：信号允许通过的频带为通带，完全不允许通过的频带为阻带，通带与阻带之间为过渡带。,滚降与滚降率：滤波器幅频特性在过渡带的衰减和衰减速度称为滚降与滚降率。,阻带衰减：输入信号在阻带的
3、衰减量,带内平坦度：通带和阻带内的平坦程度,4、数字滤波器的设计步骤,数字滤波器的设计三个步骤: (1) 按要求确定滤波器的性能参数；(2) 用一个因果稳定的离散线性移不变系统的系统函数去逼近去逼近这一性能要求；(3) 用有限精度的运算实现；实现可以采用通用计算机，也可以采用DSP。,5、数字滤波器的技术要求,选频滤波器的频率响应：,为幅频特性：表示信号通过该滤波器后各频率成分的衰减情况,为相频特性：反映各频率成分通过滤波器后在时间上的延时情况,理想滤波器不可实现，只能以实际滤波器逼近,通带最大衰减：,阻带最小衰减：,6、表征滤波器频率响应的特征参量,幅度平方响应,的极点既是共轭的，又是以单位
4、圆成镜像对称的,H(z)的极点：单位圆内的极点,相位响应,相位响应：,群延迟响应,相位对角频率的导数的负值,若滤波器通带内 = 常数，则为线性相位滤波器,7、IIR数字滤波器的设计方法,先设计模拟滤波器，再转换为数字滤波器,用一因果稳定的离散LSI系统逼近给定的性能要求：,即为求滤波器的各系数,计算机辅助设计法,s平面逼近：模拟滤波器,z平面逼近：数字滤波器,8 、将DF的技术指标转换为ALF的技术指标一、意义AF的设计有一套相当成熟的方法：设计公式； 设计图表；有典型的滤波器，如巴特沃斯，切比雪 夫等。二、一般转换方法1、2、3、4、,三、转换举例例如，一低通DF的指标：在 的通带 范围，幅
5、度特性下降小于1dB；在 的 阻带范围，衰减大于15dB；抽样频率 ； 试将这一指标转换成ALF的技术指标。解：按照衰减的定义和给定指标，则有,假定 处幅度频响的归一化值为1， 即,这样，上面两式变为,由于 ，所以当没有混叠时，根据关系式模拟filter的指标为,6.2 最小与最大相位延时系统、最小与最大相位超前系统,LSI系统的系统函数：,频率响应：,模：,相角：,当,位于单位圆内的零/极矢量角度变化为2p,位于单位圆外的零/极矢量角度变化为 0,单位圆外的零点数为mo,单位圆内的极点数为pi,单位圆外的极点数为po,则：,全部极点在单位圆内：po = 0，pi = N,因果稳定系统,1）全
6、部零点在单位圆内：,2）全部零点在单位圆外：,为最小相位延时系统,为最大相位延时系统,n 0时，h(n) = 0,最小相位延时系统的性质,1）在 相同的系统中，具有最小的相位滞后,2）最小相位延时系统的能量集中在n=0附近，而总能量相同,5）级联一个全通系统，可以将一最小相位系统转变成一相同幅度响应的非最小相位延时系统,4）在 相同的系统中， 唯一,3）最小相位序列的 最大：,6.3 全通系统,一阶全通系统：,极点：,零点：,零极点以单位圆为镜像对称,极点：,零点：,实系数二阶全通系统,两个零点（极点）共轭对称,极点：,零点：,零点与极点以单位圆为镜像对称,可以证明，一阶全通节在任何频率上，其
7、频率响应的模都为1：,N 阶数字全通滤波器,极点： 的根,零点： 的根,因为,所以,全通系统的应用,1）任一因果稳定（非最小相位延时）系统H(z)都可以表示成全通系统Hap(z)和最小相位系统Hmin(z)的级联,其中：H1(z)为最小相位延时系统，为单位圆外的一对共轭零点,而幅度响应不变：,P231 图66,2）级联一个全通系统可以使非稳定滤波器变成一个稳定滤波器,把非稳定系统的单位圆外的极点映射到单位圆内,单位圆外极点：,3）作为相位均衡器，校正系统的非线性相位，而不改变系统的幅度特性,利用均方误差最小准则求均衡器Hap(z)的有关参数,6.4 用模拟滤波器设计IIR数字滤波器,设计思想：
8、,s 平面 z 平面,模拟系统 数字系统,H(z) 的频率响应要能模仿 Ha(s) 的频率响应，即 s 平面的虚轴映射到 z 平面的单位圆,因果稳定的 Ha(s) 映射到因果稳定的 H(z) ，即 s 平面的左半平面 Res 0 映射到 z 平面的单位圆内 |z| Ws/2处衰减越快，失真越小,当滤波器的设计指标以数字域频率wc给定时，不能通过提高抽样频率来改善混迭现象,二、模拟滤波器的数字化,系数相同：,极点：s 平面 z 平面,稳定性不变：S域稳定，即极点位于s平面左半平面，故在 s 平面 ： z 平面,当T 很小时，数字滤波器增益很大，易溢出，需修正,令：,则：,试用冲激响应不变法，设计
9、IIR数字滤波器,例：设模拟滤波器的系统函数为,解：据题意，得数字滤波器的系统函数：,设T = 1s，则,模拟滤波器的频率响应：,数字滤波器的频率响应：,优点：,缺点：,保持线性关系：w=WT线性相位模拟滤波器转变为线性相位数字滤波器,频率响应混迭只适用于限带的低通、带通滤波器,h(n)完全模仿模拟滤波器的单位抽样响应ha(t) 时域逼近良好,冲激响应不变法的优缺点,6.6 阶跃响应不变法,变换原理,数字滤波器的阶跃响应g(n) 模仿模拟滤波器的阶跃响应ga(t),T 抽样周期,阶跃响应不变法同样有频率响应的混叠失真现象但比冲激响应不变法要小。,例：二阶Butterworth 归一化模拟滤波器
10、（LPF）为：,设计对应3dB截止模拟频率为50Hz的二阶Butterworth数字滤波器。设数字系统采样频率为500Hz，并采用阶跃响应不变法来设计。,解：求模拟系统函数：,最后得 （用在z-1表示）,代入T=1/500，计算ZT得,6.7 双线性变换法,冲激响应不变法、阶跃响应不变法：时域模仿逼近缺点是产生频率响应的混叠失真 为了克服这一缺点，采用双线性变换法。 使数字滤波器的频率响应与模拟滤波器的频率响应相似,一、变换原理及特点,冲激或者阶跃响应不变法的映射是多值映射,导致频率响应交叠。 改进思路：先将s域平面压缩到一个中介平面s1，然后再将s1映射到Z平面。,为使模拟滤波器某一频率与数
11、字滤波器的任一频率有对应关系，引入系数 c,2）某一特定频率严格相对应：,1）低频处有较确切的对应关系：,特定频率处频率响应严格相等，可以较准确地控制截止频率位置,二、变换常数c的选择,三、逼近情况,1）,2）,四、优缺点,优点：,避免了频率响应的混迭现象,s 平面与 z 平面为单值变换,缺点：除了零频率附近，W与w之间严重非线性,2）要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变,分段常数型模拟滤波器经变换后仍为分段常数型数字滤波器，但临界频率点产生畸变,预畸变,给定数字滤波器的截止频率 w1，则,按W1设计模拟滤波器，经双线性变换后，即可得到w1为截止频率的数字滤波器,五、模拟滤波器
12、数字化方法,可分解成级联的低阶子系统,可分解成并联的低阶子系统,6.8 常用模拟低通滤波器特性,将数字滤波器技术指标转变成模拟滤波器技术指标，设计模拟滤波器，再转换成数字滤波器,模拟滤波器,巴特沃斯 Butterworth 滤波器,切比雪夫 Chebyshev 滤波器,椭圆 Ellipse 滤波器,贝塞尔 Bessel 滤波器,1、由幅度平方函数 确定模拟滤波器的系统函数,h(t)是实函数,将左半平面的的极点归Ha(s),将以虚轴为对称轴的对称零点的任一半作为Ha(s)的零点，虚轴上的零点一半归Ha(s),Ha(s) Ha(-s)的零极点分布,由幅度平方函数得象限对称的s平面函数,对比 和 ，
13、确定增益常数,由零极点及增益常数，得,例：,解：,极点：,零点： （二阶）,零点：,的极点：,设增益常数为K0,（1） 巴特沃尔斯滤波器（Butterworth）,2、常见模拟滤波器设计,1）幅度平方函数：,当,称Wc为Butterworth低通滤波器的3分贝带宽,N为滤波器的阶数,Wc为通带截止频率,1）幅度函数特点：,3dB不变性,通带内有最大平坦的幅度特性，单调减小,过渡带及阻带内快速单调减小,当WWst（阻带截止频率）时，衰减的d1为阻带最小衰减,Butterworth滤波器是一个全极点滤波器，其极点：,2）幅度平方特性的极点分布：,极点在s平面呈象限对称，分布在Buttterwort
14、h圆上，共2N点,极点间的角度间隔为,极点不落在虚轴上,N为奇数，实轴上有极点，N为偶数，实轴上无极点,Ha(s) Ha(-s)的零极点分布 (a) N=3 (三阶） (b)N=4 （四阶）,为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(-s)。 Ha(s)的表示式为,设N=3，极点有6个（见前页图），它们分别为：（下页）,3）滤波器的系统函数：,其中分子系数 由Ha(0)=1解得。,取s平面左半平面的极点s0,s1,s2组成Ha(s)：,由前讨论知，系统函数为：,极点为：,当N为偶数，Ha(s)的极点（左半平面）皆成共轭对，记为：,这一对共轭点构成一个二阶子系统：,当N为奇数，则系统由一个一阶系统（极点s=-1）和(N-1)/2个二阶系统组成：,整个系统由N/2（N为偶数）个此二阶系统级联而成：,N为偶数,N为奇数,为归一化系统的系统函数,去归一化，得,例：试设计一个模拟低通巴特沃思滤波器，要求通带截止频率c=24000rad/s，通带最大衰减1=3dB，阻带下限截止频率st=28000rad/s，阻带最小衰减2=20dB.,
