数字信号处理（西电版） 第五章 IIR数字滤波器的设计方法 复习ppt课件

1、第5章 IIR数字滤波器的设计方法 5.1 根本概念 5.2 IIR滤波器设计的特点 5.3 常用模拟低通滤波器的设计方法 5.4 用脉冲呼应不变法设计IIR数字滤波器 5.5 用双线性变换法设计IIR数字滤波器 5.6 设计IIR滤波器的频率变换法 5.7 Z平面变换法 5.1 根本概念 数字滤波器是对数字信号实现滤波的LTI系统, 将输入的数字序列经过特定运算转变为输出的数字序列。 设计数字滤波器就是根据义务要求，确定滤波器的性能参数，用一个因果稳定的LTI离散时间系统的系统函数去逼近性能要求5.1 根本概念 1 数字滤波器的分类 经典滤波器与现代滤波器 经典：输入信号中有用频率成分与希望

2、滤除的频率各占有不同的频带 现代： 维纳滤波器、卡尔曼滤波器、自顺应滤波器等 IIR滤波器与FIR滤波器 低通、高通、带通与带阻等滤波器 数字滤波器的理想幅频特性 p s理想滤波器是非因果的，需求设计一个因果可实现的滤波器实现2 滤波器的技术目的 滤波器的性能要求通常以频率呼应幅度特性的允许误差来表征 以右图低通为例, 频率呼应包括 通带、过渡带与阻带 1(2)为通(阻)带的容限, p(s)为通(阻)带截止频率p s 通带允许的最大衰减波纹Ap和阻带应到达的最小衰减As式中 |H(ej0)|=1 (归一化) IIR滤波器普通采用递归型的实现构造3 FIR型滤波器和IIR型滤波器 IIR滤波器的

3、系统函数H(z)在Z平面上不仅有零点，而且有极点 FIR滤波器普通采用非递归型，系统函数极点只能在Z平面原点 按照实践义务要求， 确定滤波器的性能目的 用因果稳定的离散LTI系统的(FIR、IIR)系统函数逼近性能要求 利用有限精度算法实现系统函数。包括选择运算构造，适宜的字长以及有效数字的处置方法等 4 滤波器的设计步骤5.2 IIR滤波器设计的特点 将IIR滤波器的系统函数用极、零点表示： 对系统函数的设计就是确定各系数ak, bk或零极点ck，dk和A，使滤波器满足给定的性能要求MN 设计方法 最优化设计法(计算机辅助设计法)从统计的概念出发，对要提取的有用信号估计，在统计目的最优的意义

4、下，使估计值最优逼近有用信号，减弱或消除噪声 最优化设计法的设计步骤1选择一种最优准那么，如最小均方误差准那么 2求在此最正确准那么下滤波器系统函数的系数ak, bk。经过迭代运算，找到使为最小时的一组系数ak, bk 设计方法 利用模拟滤波器的实际来设计数字滤波器根据性能参数，先设计一个模拟滤波器，再按一定算法(脉冲呼应不变法和双线性变换法)转换为满足预定目的的数字滤波器 本质：一个由S平面到Z平面的复变函数的映射变换该映射变换须满足的条件下面先讨论模拟滤波器的设计方法，再引见由模拟变换为IIR数字的两种算法：脉冲呼应不变法和双线性变换法 S平面虚轴 j必需映射到Z平面的单位圆ej上 S平面

5、的左半平面必需映射到Z平面单位圆的内部|z|15.3 模拟低通滤波器的设计方法 常用的模拟原型滤波器有巴特沃思Butterworth滤波器、切比雪夫Chebyshev滤波器、椭圆Ellipse滤波器、贝塞尔Bessel滤波器等 模拟原型滤波器的特点Butterworth具有单调下降的幅频特性Chebyshev的幅频特性在通(阻)带有动摇, 可以提高选择性Bessel通带内有较好的线性相位特性Ellipse的选择性相对前三种是最好的, 但在通带和阻带内均为等波纹幅频特性5.3.1 幅度平方函数确定系统函数(模拟滤波器的设计方法)滤波器的技术目的给定后，需求设计一个传输函数Ha(s) ，希望其幅度

6、平方函数满足给定的目的，典型滤波器的幅度平方函数都有本人的表达式，可以直接援用模拟滤波器幅度呼应常用幅度平方函数|Ha(j)|2表示 冲激呼应ha(t)为实函数模拟滤波器的设计通常总是先设计出低通，再经过频率变换转换成要求类型(高通,带通等)，而低通是经过幅度平方函数来设计的式中Ha(s)是系统函数; Ha(j)是频率呼应特性； |Ha(j)|是幅度特性 根据 |Ha(j)|2 求 Ha(s) Ha(s)Ha(-s)的极(零)点分布是成象限对称的 Ha(s)必需稳定，因此极点必定落在s左半平面 零点的分布只和滤波器的相位特征有关，如无特殊要求，可将对称零点的任一半应为共轭对取为Ha(s)的零点

7、 根据Ha(s)极(零)点可确定增益常数，最终得到系统函数Ha(s) 巴特沃思低通滤波器幅度平方函数定义 式中，N为正整数，代表滤波器的阶数。c为3dB截止频率。当=c时，衰减为 3 dB5.3.2 巴特沃思低通逼近 最平幅度逼近巴特沃思低通滤波器在通带内有最大平坦的幅度特性，因此又称为最平幅度特性滤波器|Ha(j)|2单调减小，N越大，通带内特性越平坦，过渡带越窄 巴特沃思滤波器的极(零)点分布 (公式法求解低通Ha(s)零点全在s=处，有限S平面内只需极点，属“全极点型滤波器 Ha(s)Ha(-s)的2N个极点等间隔分布在半径为c的圆巴特沃思圆上，极点间的角度间隔为/N rad 令=s/j

8、 k=1, 2, , 2N Ha(s)Ha(-s)的极点为 例N3及N4时，Ha(s)Ha(-s)的极点分布分别如以下图示N为奇(偶)数时，实轴上有(无)极点。且极点不会落在虚轴上 取S左半平面的N个极点为Ha(s)的极点，即得到Ha(s)的表示式 k=1, 2, , N 低通巴特沃思滤波器的设计步骤阶数N一样的模拟低通经归一化后，其Ha(s)表示式都是一样的，可以经过查表获得，因此设计巴特沃思滤波器的本质就是求滤波器阶次N和截止频率c (归一化) 阶数N的计算 由目的p, Ap，s和As得到 巴特沃思低通滤波器技术目的关系式阶数N的大小主要影响幅度特性下降速度知代入上式整理得到由上式可解出滤

9、波器阶数N和截止频率c或 根据阶数N，查表得到归一化系统函数HaN(s) 根据c将HaN(s)去归一化，得到实践要求的系统函数Ha(s) 低通巴特沃思滤波器设计步骤总结 知p, Ap，s和As，计算滤波器阶数N和截止频率c或 根据N查表得到归一化系统函数HaN(s) 根据c将HaN(s)去归一化，得到实践要求的系统函数Ha(s) 例 导出三阶巴特沃思模拟低通滤波器的系统函数设c2 rad/s解1 由知得到幅度平方函数为Ha(s)的表示式为P157k=1, 2, , N 当N=3可求得Ha(s)的极点skk=1, 2, 3，从而得到系统函数 例 导出三阶巴特沃思模拟低通滤波器的系统函数设c2 r

10、ad/s解2 由P158表5可知归一化的三阶巴特沃思模拟低通为 知 c2 rad/s，那么阐明：巴特沃斯归一化低通滤波器参数表有多种方式例设计一个满足下面要求的模拟低通巴特沃思滤波器： (1) 通带截止频率：p=0.2；通带最大衰减：Ap=7 dB。 (2) 阻带截止频率：s=0.3；阻带最小衰减：As=16dB。 知p, Ap，s和As，计算滤波器阶数N和截止频率c 根据N查表得到归一化系统函数HaN(s) 根据c，将HaN(s)去归一化取 N = 3或 知p, Ap，s和As，计算滤波器阶数N和截止频率c 根据c，将HaN(s)去归一化，得到实践要求的系统函数Ha(s) 根据N查表得到归一

11、化系统函数HaN(s)或代入上式将 该映射变换须满足以下要求： S平面虚轴 j必需映射到Z平面的单位圆ej上 S平面的左半平面必需映射到Z平面单位圆的内部|z|1 利用模拟滤波器来设计数字滤波器 常用的变换方法有：脉冲呼应不变法和双线性变换法使稳定的模拟系统映射为具有一样幅频特性的稳定的数字系统5.4 用脉冲呼应不变法设计IIR数字滤波器使数字滤波器的单位脉冲呼应序列h(n)与相应的模拟滤波器的冲激呼应ha(t)，在采样点处的量值相等，即h(n)=ha(nT) 5.4.1 变换原理第1章中采样序列的Z变换与模拟信号拉氏变换的关系P48脉冲呼应不变法将模拟滤波器的S平面变换成数字滤波器的Z平面

12、脉冲呼应不变法的映射关系 数字滤波器与模拟滤波器的频率呼应间的关系为 (P49) 数字滤波器的频率呼应是模拟滤波器频率呼应的周期延拓 只需当模拟滤波器的频率呼应是限带的，且满足采样定理时，数字滤波器的频率呼应才不产生混叠失真 提高采样频率为fs，可减小频率呼应的混叠效应 5.4.2 混叠失真脉冲呼应不变法中的频响混叠景象 脉冲呼应不变法设计中，由一个较为复杂的模拟系统函数求出数字滤波器系统函数是一个很费事的变换过程 由于乘积的z变换并不等于各部分变换的乘积，所以不宜采用级联分解 但各项和的z变换是线性关系，因此用部分分式表达系统函数，特别适宜于对复杂模拟系统函数的变换5.4.3 模拟滤波器的数

13、字化方法设模拟滤波器的系统函数Ha(s)只需单阶极点，且分母的阶次大于分子的阶次，将Ha(s)用部分分式表示 数字滤波器的单位脉冲呼应等于对ha(t)的采样，因此 将Ha(s) 拉氏反变换u(t)是单位阶跃函数 模拟滤波器与数字滤波器系统函数的关系对h(n)求Z变换，得数字滤波器的系统函数 S平面的单极点s=sk对应于Z平面上z=eskT处的单极点 Ha(s)与H(z)的部分分式的系数一样，都是Ak 模拟滤波器因果稳定，那么数字滤波器也因果稳定 脉冲呼应不变法只能保证极点有这种代数对应关系与模拟滤波器的系统函数Ha(s) 相比 Ha(s)与H(z)间的变换关系为 脉冲呼应不变法的设计步骤，可直

14、接将 Ha(s)写成单极点的部分分式和的方式，然后将各部分分式用上式关系替代，即得到H(z) 由于数字滤波器频率呼应幅度与采样间隔T成反比(p49) | 为使数字滤波器增益不随采样频率变化，通常令 h(n)=Tha(nT)此时: 将上式代入到例 设模拟滤波器的系统函数为 试利用脉冲呼应不变法将Ha(s)转换成IIR数字滤波器的系统函数H(z) 解 数字滤波器的系统函数为 设T=1，那么有 模拟滤波器的频响Ha(j)以及数字滤波器的频响H(ej)分别为: 由于Ha(j)不是限带的，所以H(ej)产生了频谱混叠失真 T取不同值时的幅频特性图数字滤波器的单位脉冲呼应完全模拟模拟滤波器的单位冲激呼应，

15、因此时域特性逼近好 频率坐标变换是线性的，即=T。因此该数字滤波器可以很好重现原模拟滤波器的频率特性 会产生频率混叠景象从S平面到平面为多值映射关系，适宜低通、带通滤波器的设计，而不适宜高通、带阻的设计 5.4.4 优缺陷5.5 用双线性变换法设计IIR数字滤波器 5.5.1 变换原理为抑制频率混叠景象，可采用非线性频率紧缩方法，将整个频率轴的频率范围紧缩到-/T/T之间，再用z=esT转换到Z平面上使S平面与Z平面建立起一一对应的单值关系将此关系解析延拓到整个S平面和S1平面，令j=s，j1=s1，得 再将S1平面经过规范变换z=es1T映射到Z平面 实现步骤得到两平面的单值映射关系 利用正

16、切变换实现频率紧缩，将S平面的虚轴紧缩至S1平面的-/T到/T段双线性变换符合映射变换应满足的两点要求：5.5.2 逼近的情况 S平面的虚轴映射到Z平面的单位圆 S平面的左半平面映射到Z平面的单位圆内稳定模拟滤波器经双线性变换后所得的数字滤波器也是稳定的 阐明：所谓“双线性变换，仅仅是指变换公式中s与z的关系无论是分子部分，还是分母部分都是“线性的 由于S平面与Z平面的单值对应关系，防止了频率混叠景象5.5.3 优缺陷 S平面与Z平面的严重非线性关系影响DF频响对模拟频响的模拟 该方法适于分段常数型的滤波器设计，普通设计滤波器通带与阻带均要求是分段常数(非线性的正切关系) 双线性变换法幅度和相

17、位特性的非线性映射 分段常数的滤波器，在双线性变换后各个分段边缘的临界频率点将会产生畸变，这种畸变可以预处置加以校正 数字滤波器的系统函数可经过s到z之间的简单代数关系得到 频率呼应也可用直接代换的方法得到 5.5.4 模拟滤波器的数字化方法 阶数较高时，H(z)的方式可经过查表得到设模拟系统函数的表达式为 数字滤波器的系统函数双线性变换法中a(s)的系数与H(z)的系数间关系 双线性变换时频率的预畸变 总结：双线性变换法设计IIR数字滤波器的步骤假设设计带通滤波器(右图) 频率预畸变 知带通的数字域频率与采样频率时利用下式计算相应的模拟域频率双线性变换时频率的预畸变 总结：双线性变换法设计I

18、IR数字滤波器的步骤 频率预畸变 知带通的模拟域频率与采样频率时先将模拟频率转为数字域频率=2fT 12再按先前步骤处置 总结：双线性变换法设计IIR数字滤波器的步骤 按模拟频率目的设计模拟滤波器的系统函数Ha(s) 将 代入Ha(s)，得H(z) 与频率呼应 阐明 双线性变换法后，模拟滤波器与数字滤波器仍为同种类型低通、高通、带通等的滤波器当强调模拟滤波器的瞬态呼应时，通常采用脉冲呼应不变法; 其他那么采用双线性变换法 例 设计一阶数字低通滤波器，3 dB截止频率为c=0.25，将双线性变换运用于模拟巴特沃思滤波器 频率预畸变 模拟滤波器的系统函数设计 知 数字低通c=0.25，那么巴特沃思

19、模拟低通c为 转换为数字滤波器系统函数由上题可知，T不参与设计，即双线性变换法中用 设计与用 设计得到的结果一致。 例 用双线性变换法设计三阶巴特沃思数字低通，采样频率为fs=4 kHz，其3dB截止频率为fc=1 kHz。 三阶模拟巴特沃思滤波器为 知截止频率 fc=1 kHz解：1. 频率预畸变2. 确定预畸变的模拟滤波器的截止频率 先转化成数字频率c=2fcT=0.5将c代入三阶模拟巴特沃思滤波器Ha(s)，得到 模拟滤波器的系统函数设计 转换为数字滤波器系统函数滤波器Ha(s)不是截止频率为fc=1 kHz的实践滤波器，只是一个“样本函数 用双线性变换法设计得到的三阶巴特沃思数字低通滤

20、波器的频响 取 N = 25.6 设计IIR滤波器的频率变换法 设计各种IIR数字滤被器，通常是先设计一个归一化截止频率的模拟低通原型滤波器，然后用如下两种方法转换为数字滤波器： 模拟频带变换法 数字频带变换法将归一化的模拟低通原型滤波器经模拟频带变换成所需类型的模拟滤波器，然后用冲激呼应不变法或双线性变换法转换为所需类型的数字滤波器将归一化的模拟低通原型滤波器经过冲激呼应不变法或双线性变换法转换为数字低通滤波器，然后用数字频带变换法，转换为所需类型数字滤波器5.6 设计IIR滤波器的频率变换法 模拟频带变换法 两种等效的设计方法 模拟域频率变换关系可经过查表得到截止频率为c的模拟低通滤波器到

21、其它频率选择性滤波器的转换公式 普通直接用归一化原型转换，取c=1, 可使设计过程简化 数字化可采用脉冲呼应不变法与双线性变换法 把第一种方法中的两步合成一步来实现 把数字滤波器的性能要求转换为对应的模拟滤波器的性能要求 根据性能要求设计模拟滤波器 经过脉冲呼应不变法或双线性变换法，转化为数字滤波器5.6.1 模拟低通滤波器变换成数字低通滤波器 设计步骤由于采用脉冲呼应不变法设计，因此模拟低通的 fc=1 kHz 例 用脉冲呼应不变法设计三阶巴特沃思数字低通滤波器，采样频率为fs=4 kHz，其3 dB截止频率为fc=1 kHz 知所求数字巴特沃思低通滤波器 fs=4 kHz，fc=1 kHz

22、 设计模拟低通滤波器查表得归一化三阶模拟巴特沃思原型低通滤波器的传送函数 去归一化 ， s=s/c = s/(2fc )代入上式 经过脉冲呼应不变法转化为数字滤波器将上式因式分解将部分分式系数代入下式 H(z)只与c有关，即只与截止频率fc与采样频率fs的相对值有关，而与它们的绝对大小无关整理后得到 c=cT=2fcT=0.5 例如fs=4kHz，fc=1 kHz与fs=40 kHz，fc=10kHz的数字滤波器将具有同一个系统函数。将c=cT=2fcT=0.5代入得到 该系统函数可采用并联型的运算构造来实现 以下图为三阶巴特沃思数字低通滤波器的频响幅度特性脉冲呼应不变法存在微小的混淆景象由上

23、表知，低通模拟原型到模拟高通的变换关系为 5.6.2 模拟低通滤波器变换成数字高通滤波器截止频率为c的模拟低通滤波器到其它频率选择性滤波器的转换公式 c为模拟低通截止频率c为实践高通截止频率利用双线性变换方法设计时，模拟高通与数字高通变换关系为 将两式结合，得到直接从模拟低通原型变换成数字高通表达式c为模拟低通截止频率c为实践高通截止频率式中 由此得到数字高通系统函数为 数字高通和模拟低通的极点数目或阶次一样根据双线性变换，模拟高通与数字高通频率间的关系为 于是 有又因 Ha(s)为模拟低通传送函数 系数C的求解 模拟低通与数字高通间的关系或 其变换关系曲线如右图示令s=j，z=ej， 代入

24、，得到 =0 映射到=, 即z= -1上;=映射到=0, 即z= 1上经过上述频率变换，即可将模拟低通直接变换为数字高通或 阐明：数字高通指从到经过例 设计一个巴特沃思高通数字滤波器，其通带截止频率-3dB点处为fc=3kHz，阻带上限截止频率fst=2kHz，通带衰减不大于3 dB， 阻带衰减不小于 14dB，采样频率fs=10kHz 1 求对应的各数字域频率:2求常数C 3求低通原型st采用归一化c=1原型低通作为变换的低通原型，那么低通到高通的变换中所需的C为4求阶次N取N=3 按阻带衰减求原型归一化模拟低通滤波器的阶次N5 求归一化巴特沃思低通原型的Ha(s)。查表可得Ha(s)为 6

25、 求数字高通滤波器的系统函数H(z)， 有 由表5-3知，低通模拟原型到模拟带通的变换关系为 5.6.3 模拟低通滤波器变换成数字带通滤波器截止频率为c的模拟低通滤波器到其它频率选择性滤波器的转换公式 由双线性变换，模拟带通与数字带通整理后得到 式中 将参数D、E进一步整理后得到 式中:0、2和1分为数字带通中心频率与上、下边带截止频率定义: 0为带通滤波器通带的中心频率，B为带通滤波器的通带宽度0、2和1间的关系c的求解公式当D和E经过中心频率和带宽或者是中心频率和边带频率确定后，数字带通系统函数那么也确定下来Ha(s)为模拟低通原型传送函数 数字带通的极点数或阶数是模拟低通极点数的两倍 模

26、拟低通与数字带通频率间的关系其映射关系为： 低通的通带(=0 附近) 带通的通带(=0附近)低通的阻带(=) 带通的阻带(=0, ) 令s=j，z=ej代入变换式中 得到 频率间关系曲线=0=0=-=0 模拟低通变换到数字带通 例 采样频率fs=100kHz, T=10s，设计三阶巴特沃思数字带通，其上、下边带的 3dB截止频率分为f2=37.5 kHz，f1=12.5 kHz1=2f1T=212.51031010-6=0.252=2f2T=237.51031010-6=0.75 确定模拟低通到数字带通的变换公式 确定待设计的数字目的再代入变换公式得 确定模拟低通原型 3 dB截止频率为c=4/T 三阶巴特沃思滤波器的系统函数为 由N=3，查表得到归一化原型系统函数为 转换为数字带通 阐明： 如采用归一化低通原型(c=1) 设计数字带通，结果是一样的。由于 s/c中的c和D中的c相互抵消 设计时，只需用c=1 的归一化原型HaN(s)设计 其他类型的滤波器，也可直接利用归一化原型设计 模拟低通原型到模拟带阻的变换关系为 5.6.4 模拟低通滤波器变换成数字带阻滤波器由双线性变换，模拟带通与数字带通整理后得