第6章IIR数字滤波器的设计方法ppt课件.ppt

《第6章IIR数字滤波器的设计方法ppt课件.ppt》由会员分享，可在线阅读，更多相关《第6章IIR数字滤波器的设计方法ppt课件.ppt（117页珍藏版）》请在三一办公上搜索。

1、第6章 IIR数字滤波器的设计方法,6.1 数字滤波器的基本概念6.2 数字滤波器的技术指标6.3 全通系统6.4 最小与最大相位延时系统、最小与最大相位超前系统6.5 模拟原型低通滤波器设计6.6 模拟频带变换法设计各种响应的模拟滤波器6.7 间接法的IIR数字滤波器设计方案6.8 模拟滤波器数字化为数字滤波器的映射方法,6.9 将样本模拟归一化低通滤波器先经模拟频带变换，再数字化方案的设计步骤6.10 将样本模拟低通滤波器直接数字化为各种频率响应数字滤波器的设计方案6.11 数字频域频带变换。将样本模拟归一化低通滤波器先数字化，再做数字频域频带变换的设计方法,理解数字滤波器的基本概念理解最

2、小相位延时系统理解全通系统的特点及应用掌握冲激响应不变法掌握双线性变换法掌握Butterworth、Chebyshev低通滤波器的特点了解利用模拟滤波器设计IIR数字滤波器的设计过程,第6章学习目标,本章作业练习,P336: 6.10(1)6.13,2要求从二阶巴特沃思模拟滤波器用双线性变换导出一低通数字滤波器，已知3dB截止频率为100Hz，系统抽样频率为1kHz。,数字滤波器：,输入输出均为数字信号，其功能是通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分。,高精度、稳定、体积小、重量轻、灵活，不要求阻抗匹配，可实现特殊滤波功能。,优点：,6.1 数字滤波器的基本概念

3、,1、数字滤波器的分类,经典滤波器：,现代滤波器：,选频滤波器,维纳滤波器卡尔曼滤波器自适应滤波器等,按功能分：低通LPF 高通HPF 带通BPF 带阻BSF 全通滤波器APF,按实现的网络结构或单位抽样响应分：,FIR滤波器,IIR滤波器（N阶）,2、数字滤波器的设计过程,用一个因果稳定的离散LSI系统的系统函数H(z)逼近此性能指标,按设计任务，确定滤波器性能要求，制定技术指标,利用有限精度算法实现此系统函数：如运算结构、字长的选择等,实际技术实现：软件法、硬件法或DSP芯片法,IIR数字滤波器的设计方法,先设计模拟滤波器，再转换为数字滤波器,用一因果稳定的离散LSI系统逼近给定的性能要求

4、：,即求滤波器的各系数：,计算机辅助设计法,s平面逼近：模拟滤波器,z平面逼近：数字滤波器,6.2 数字滤波器的技术,选频滤波器的频率响应：,为幅频特性： 表示信号通过该滤波器后各频率成分的衰减情况,为相频特性： 反映各频率成分通过滤波器后在时间上的延时情况,：通带截止频率,：阻带截止频率,：通带容限,：阻带容限,阻带：,过渡带：,通带：,理想滤波器不可实现，只能以实际滤波器逼近,通带最大衰减：,阻带最小衰减：,其中：,表征滤波器频率响应的特征参量,幅度平方响应,的极点既是共轭的，又是以单位圆成镜像对称的（h(n)为实序列）,为了系统可实现，H(z)的极点只取单位圆内的极点,相位响应,相位响应

5、：,群延迟响应,相位对角频率的导数的负值,若滤波器通带内 = 常数，则为线性相位滤波器,6.4 最小与最大相位延时系统、最小与最大相位超前系统,LSI系统的系统函数：,频率响应：,模：,相角：,当,位于单位圆内的零/极矢量角度变化为,位于单位圆外的零/极矢量角度变化为 0,单位圆外的零点数为mo,单位圆内的极点数为pi,单位圆外的极点数为po,因果稳定系统全部极点在单位圆内：po = 0，pi = N,1）全部零点在单位圆内：,2）全部零点在单位圆外：,为最小相位延时系统,为最大相位延时系统,n 0时，h(n) = 0,相位延时系统,逆因果稳定系统,1）全部零点在单位圆内：,2）全部零点在单位

6、圆外：,全部极点在单位圆外：po = N，pi = 0,为最大相位超前系统,为最小相位超前系统,相位超前系统,n 0时，h(n) = 0,最小相位延时系统的性质,1）在 相同的系统中，具有最小的相位滞后,2）最小相位延时系统的能量集中在 n = 0 附近， 而总能量相同,5）级联一个全通系统，可以将一最小相位延时系统转变成一相同幅度响应的非最小相位延时系统,4）在 相同的系统中， 唯一,3）最小相位序列的 最大：,6.3 全通系统,一阶全通系统：,极点：,零点：,零极点以单位圆为镜像对称,极点：,零点：,实系数二阶全通系统,两个零点（极点）共轭对称,极点：,零点：,零点与极点以单位圆为镜像对称

7、,N 阶数字全通滤波器,极点：,零点：,全通系统的应用,1）任一因果稳定系统H(z)都可以表示成全通系统 Hap(z)和最小相位系统Hmin(z)的级联,其中：H1(z)为最小相位延时系统， 为单位圆外的一对共轭零点,而幅度响应不变：,2）级联一个全通系统可以使非稳定滤波器变成一个稳定滤波器,把非稳定系统的单位圆外的极点映射到单位圆内,单位圆外极点：,3）作为相位均衡器，校正系统的非线性相位，而不改变系统的幅度特性,利用均方误差最小准则求均衡器Hap(z)的有关参数,例：如下图是一个线性时不变因果系统，（1）求系统函数H(z)；（2）若a1=-2，b0=-2，b1=1，画出系统的零极点图；（3

8、）用几何法确定系统的幅频响应；（4）讨论系统的稳定性，若系统不稳定，则在不改 变系统的滤波特性的情况下，重新设计系统。,w(n),w(n+1),6.8 模拟滤波器数字化为数字滤波器的映射方法,设计思想：,s 平面 z 平面,模拟系统 数字系统,H(z) 的频率响应要能模仿 Ha(s) 的频率响应，即 s 平面的虚轴映射到 z 平面的单位圆,因果稳定的 Ha(s) 映射到因果稳定的 H(z),即 s 平面的左半平面 Res 0 映射到 z 平面的单位圆内 |z| 1,设计方法：,- 冲激响应不变法,- 阶跃响应不变法,- 双线性变换法,6.8.1 冲激响应不变法,数字滤波器的单位冲激响应 模仿模

9、拟滤波器的单位冲激响应,T抽样周期,1、变换原理,2、混迭失真,仅当,数字滤波器的频响在折叠频率内重现模拟滤波器的频响而不产生混迭失真：,数字滤波器的频率响应是模拟滤波器频率响应的周期延拓，周期为,实际系统不可能严格限带，都会混迭失真，在 处衰减越快，失真越小,当滤波器的设计指标以数字域频率 给定时，不能通过提高抽样频率来改善混迭现象,3、模拟滤波器的数字化方法,系数相同：,极点：s 平面 z 平面,稳定性不变：s 平面 z 平面,当T 很小时，数字滤波器增益很大，易溢出，需修正,令：,则：,解：据题意，得数字滤波器的系统函数：,设T = 1s，则,模拟滤波器的频率响应：,数字滤波器的频率响应

10、：,4、优缺点,优点：,缺点：,保持线性关系：线性相位模拟滤波器转变为线性相位数字滤波器,频率响应混迭只适用于限带的低通、带通滤波器,h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好,6.8.2 双线性变换法,1、变换原理,使数字滤波器的频率响应与模拟滤波器的频率响应相似。,冲激响应不变法：时域模仿逼近缺点是产生频率响应的混叠失真,为使模拟滤波器某一频率与数字滤波器的任一频率有对应关系，引入系数 c,2、变换常数c的选择,2）某一特定频率严格相对应：,1）低频处有较确切的对应关系：,特定频率处频率响应严格相等，可以较准确地控制截止频率位置,3、逼近情况,1）,2）,4、优缺点,优点：,避免

11、了频率响应的混迭现象,s 平面与 z 平面为单值变换,缺点：,除了零频率附近， 与 之间严重非线性,预畸变,给定数字滤波器的截止频率 ，则,按1设计模拟滤波器，经双线性变换后，即可得到 1为截止频率的数字滤波器。,5、模拟滤波器的数字化方法,可分解成级联的低阶子系统,可分解成并联的低阶子系统,6.5 模拟原型低通滤波器设计,将数字滤波器技术指标转变成模拟滤波器技术指标，设计模拟滤波器，再转换成数字滤波器,模拟滤波器,巴特沃斯 Butterworth 滤波器,切比雪夫 Chebyshev 滤波器,椭圆 Ellipse 滤波器,贝塞尔 Bessel 滤波器,1、由幅度平方函数 确定模拟滤波器的系统

12、函数,h(t)是实函数,将左半平面的极点归,虚轴上的零点一半归,由幅度平方函数得象限对称的s平面函数,将 因式分解，得到各零极点,对比 和 ，确定增益常数,由零极点及增益常数，得,例：,解：,极点：,零点： （二阶）,零点：,的极点：,设增益常数为K0,2、模拟Butterworth 低通滤波器,幅度平方函数：,当,称 为Butterworth低通滤波器的3分贝带宽,N 为滤波器的阶数,为通带截止频率,1）幅度函数特点：,3dB不变性,通带内有最大平坦的幅度特性，单调减小,过渡带及阻带内快速单调减小,当 （阻带截止频率）时，衰减 ，为阻带最小衰减,Butterworth滤波器是一个全极点滤波器

13、，其极点：,2）幅度平方特性的极点分布：,极点在s平面呈象限对称，分布在Buttterworth圆上，共2N点,极点间的角度间隔为,极点不落在虚轴上,N为奇数，实轴上有极点，N为偶数，实轴上无极点,左半平面的极点：,3）滤波器的系统函数：,为归一化系统的系统函数,去归一化，得,4）滤波器的设计步骤：,确定技术指标：,由,得：,同理：,令,则：,根据技术指标求出滤波器阶数N及,求出归一化系统函数：,或者由N，直接查P278表6.2得,其中极点：,去归一化,阻带指标有富裕,或,通带指标有富裕,例6.4：P276,例：设计Butterworth数字低通滤波器，要求在频率低于 rad的通带内幅度特性下

14、降小于1dB。在频率 到 之间的阻带内，衰减大于15dB。分别用冲激响应不变法和双线性变换法。,1、用冲激响应不变法设计,1）由数字滤波器的技术指标：,2）得模拟滤波器的技术指标：选T = 1 s,P302例6.8,a）确定参数,用通带技术指标，使阻带特性较好，改善冲激响应不变法引起的混迭失真。,3）设计Butterworth模拟低通滤波器,b) 求出极点（左半平面）,c) 构造系统函数,c) 去归一化,4）将 展成部分分式形式:,变换成Butterworth数字滤波器：,设T=1s,2、用双线性变换法设计,1）由数字滤波器的技术指标：,2）考虑预畸变，得模拟滤波器的技术指标：,a）确定参数,

15、因双线性变换法无混迭问题，所以用阻带技术指标，使通带特性较好。,3）设计Butterworth模拟低通滤波器,b) 求出极点（左半平面）,c) 构造系统函数,c) 去归一化,4）将 变换成Butterworth数字滤波器：,3、Chebyshev低通逼近,N：滤波器的阶数,幅度平方函数：,：通带截止频率，不一定为3dB带宽,，表示通带波纹大小， 越大，波纹越大,：N阶Chebyshev多项式,1）幅度函数特点：,通带外：迅速单调下降趋向0,N为偶数,N为奇数,通带内：在1和 间等波纹起伏,2）Chebyshev滤波器的三个参量：,：通带截止频率，给定,：表征通带内波纹大小,N：滤波器阶数，等于

16、通带内最大最小值的总数,由通带衰减决定,阻带衰减越大所需阶数越高,为阻带截止频率,3）幅度平方特性的极点分布：,4）滤波器的系统函数：,其中：,5）滤波器的设计步骤:,确定技术指标：,根据技术指标求出滤波器阶数N及 ：,求系统函数：,或者由N和 ，直接查表6.5得,其中极点由下式求出：,去归一化,例6.5：P277,例：用双线性变换法设计Chebyshev数字低通滤波器，要求在频率低于 rad的通带内幅度特性下降小于1dB。在频率 到 之间的阻带内，衰减大于15dB。,1）由数字滤波器的技术指标：,2）考虑预畸变，得模拟滤波器的技术指标：,P305例6.9,a）确定参数,3）设计Chebysh

17、ev模拟低通滤波器,b) 求左半平面极点,c) 构造系统函数,c) 去归一化,4）将 变换成Chebyshev数字滤波器：,将数字滤波器的技术指标转变成模拟滤波器的技术指标,小结：利用模拟滤波器设计IIR数字滤波器的步骤,通带截止频率 、通带衰减,阻带截止频率 、阻带衰减,通带截止频率,阻带截止频率,通带截止频率,阻带截止频率,确定数字滤波器的技术指标：,冲激响应不变法,双线性变换法,按模拟滤波器的技术指标设计模拟低通滤波器,Butterworth低通滤波器,Chebyshev低通滤波器,将模拟低通滤波器转换成数字低通滤波器,冲激响应不变法,双线性变换法,6.7 间接法的IIR数字滤波器设计方

18、案,6.9 将样本模拟归一化低通滤波器先经模拟频带变换，再数字化方案的设计步骤,归一化：,1、模拟低通模拟低通,当通带截止频率 时，相当于去归一化,2、模拟低通模拟带通,变换关系：,B为通带带宽,归一化变换：,其中：,例：设计一个数字带通滤波器，通带范围为0.3 rad到0.4 rad，通带内最大衰减为3dB，0.2 rad以下和0.5 rad以上为阻带，阻带内最小衰减为18dB。采用Butterworth模拟低通滤波器、双线性变换法。,解：1）确定数字带通滤波器的技术指标：,2）转换为模拟带通滤波器的技术指标：,（选T=1s）,归一化：,3）转换为归一化模拟低通滤波器技术指标：,取小者：,又有,4）设计归一化模拟低通滤波器：,查表得：,5）将归一化模拟低通滤波器转换成模拟带通滤波器：,6）通过双线性变化法将模拟带通滤波器转换成数字带通滤波器：,模拟低通到数字带通的直接变换：,其中：,3、模拟低通模拟带阻,变换关系：,B为阻带带宽,归一化变换：,其中：,模拟低通到数字带阻的直接变换：,其中：,4、模拟低通模拟高通,变换关系：,归一化变换：,其中：,模拟低通到数字高通的直接变换：,其中：,P284例6-10,