数字信号处理 第6章

《数字信号处理 第6章》由会员分享，可在线阅读，更多相关《数字信号处理 第6章（97页珍藏版）》请在金锄头文库上搜索。

1、第5章 无限脉冲响应(IIR)数字 滤波器的设计,5.1 学习要点 5.2 例题 5.3 教材第6章习题与上机题解答,第6章 无限脉冲响应数字滤波器的设计 习题解答1设计一个巴特沃斯低通滤波器， 要求通带截止频率fp=6 kHz，通带最大衰减ap=3 dB， 阻带截止频率fs=12kHz， 阻带最小衰减as=25 dB。 求出滤波器归一化系统函数G(p)以及实际滤波器的Ha(s)。 解: （1） 求阶数N。,将ksp和sp值代入N的计算公式, 得,所以取N=5（实际应用中， 根据具体要求， 也可能取N=4， 指标稍微差一点， 但阶数低一阶， 使系统实现电路得到简化）。,（2） 求归一化系统函数

2、G(p)。 由阶数N=5直接查教材第157页表6.2.1, 得到五阶巴特沃斯归一化低通滤波器系统函数G(p)为,或,当然， 也可以先按教材（6.2.13）式计算出极点：,再由教材（6.2.12）式写出G(p)表达式为,最后代入pk值并进行分母展开， 便可得到与查表相同的结果。（3） 去归一化（即LP-LP频率变换）， 由归一化系统函数G(p)得到实际滤波器系统函数Ha(s)。,由于本题中ap=3 dB， 即c=p=26103 rad/s， 因此,对分母因式形式， 则有,如上结果中，c的值未代入相乘， 这样使读者能清楚地看到去归一化后，3 dB截止频率对归一化系统函数的改变作用。 2 设计一个切

3、比雪夫低通滤波器， 要求通带截止频率fp=3 kHz，通带最大衰减p=0.2 dB，阻带截止频率fs=12 kHz, 阻带最小衰减s=50 dB。 求出滤波器归一化系统函数G(p)和实际的Ha(s)。 解: （1） 确定滤波器技术指标。 p=0.2 dB, p=2fp=6103 rad/s s=50 dB, s=2fs=24103 rad/s p=1,(4) 求阶数N和。,为了满足指标要求， 取N=4。,(3) 求归一化系统函数G(p),其中， 极点pk由教材（6.2.46）式求出如下：,(4) 将G(p)去归一化， 求得实际滤波器系统函数Ha(s)：,其中， sk=ppk=6103pk, k

4、=1, 2, 3, 4。 因为p4=p1*, p3=p2*, 所以， s4=s1*, s3=s2*。 将两对共轭极点对应的因子相乘， 得到分母为二阶因子的形式， 其系数全为实数。,也可得到分母多项式形式， 请读者自己计算。 3 设计一个巴特沃斯高通滤波器， 要求其通带截止频率fp=20 kHz， 阻带截止频率fs=10 kHz， fp处最大衰减为3 dB， 阻带最小衰减as=15 dB。 求出该高通滤波器的系统函数Ha(s)。,解: (1) 确定高通滤波器技术指标要求：p=20 kHz, ap=3 dBfs=10 kHz, as=15 dB(2) 求相应的归一化低通滤波器技术指标要求： 套用图

5、5.1.5中高通到低通频率转换公式， p=1, s=p/s， 得到p=1, ap=3 dB,as=15 dB,(3) 设计相应的归一化低通G(p)。 题目要求采用巴特沃斯类型， 故,所以， 取N=3， 查教材中表6.2.1， 得到三阶巴特沃斯归一化低通G(p)为,(4) 频率变换。 将G(p)变换成实际高通滤波器系统函数H(s)：,式中c=2fc=220103=4104 rad/s,4. 已知模拟滤波器的系统函数Ha(s)如下：,(1),(2),式中a、 b为常数， 设Ha(s)因果稳定， 试采用脉冲响应不变法将其转换成数字滤波器H(z)。,解: 该题所给Ha(s)正是模拟滤波器二阶基本节的两

6、种典型形式。 所以， 求解该题具有代表性， 解该题的过程， 就是导出这两种典型形式的Ha(s)的脉冲响应不变法转换公式。 设 采样周期为T。,(1),Ha(s)的极点为s1=a+jb, s2=ajb 将Ha(s)部分分式展开（用待定系数法）：,比较分子各项系数可知， A1、 A2应满足方程：,解之得， A1=1/2， A2=1/2， 所以,套用教材（6.3.4）式， 得到,按照题目要求， 上面的H(z)表达式就可作为该题的答案。 但在工程实际中， 一般用无复数乘法器的二阶基本节结构来实现。 由于两个极点共轭对称， 所以将H(z)的两项通分并化简整理， 可得,这样， 如果遇到将,用脉冲响应不变法

7、转换成数字滤波器时， 直接套用上面的公式即可， 且对应结构图中无复数乘法器， 便于工程实际中实现。,(2),Ha(s)的极点为s1=a+jb, s2=ajb 将Ha(s)部分分式展开：,套用教材（6.3.4）式， 得到,通分并化简整理， 得到,5 已知模拟滤波器的系统函数如下：,(1),(2),试采用脉冲响应不变法和双线性变换法将其转换为数字滤波器。 设T=2 s。 解: . 用脉冲响应不变法,(1),方法一 直接按脉冲响应不变法设计公式， Ha(s)的极点为,将T=2代入上式， 得,方法二 直接套用4题（2）所得公式。 为了套用公式， 先对Ha(s)的分母配方， 将Ha(s)化成4题中的标准

8、形式：,c为一常数,由于,所以,对比可知， ， 套用公式， 得,(2),或通分合并两项得, 用双线性变换法 （1）,(2),6 设ha(t)表示一模拟滤波器的单位冲激响应， 即,用脉冲响应不变法， 将此模拟滤波器转换成数字滤波器（用h(n)表示单位脉冲响应， 即 h(n)=ha(nT)）。 确定系统函数H(z)， 并把T作为参数， 证明： T为任何值时， 数字滤波器是稳定的， 并说明数字滤波器近似为低通滤波器还是高通滤波器。,解: 模拟滤波器系统函数为,Ha(s)的极点s1=0.9， 故数字滤波器的系统函数应为,H(z)的极点为z1=e0.9T, |z1|=e0.9T,题6解图,所以， T0时

9、， |z1|1， H(z)满足因果稳定条件。 对T=1和T=0.5， 画出H(ej)曲线如题6解图实线和虚线所示。,由图可见， 该数字滤波器近似为低通滤波器。 且T越小， 滤波器频率混叠越小， 滤波特性越好（即选择性越好）。 反之，T越大， 极点 离单位圆越远， 选择性越差， 而且频率混叠越严重， =附近衰减越小， 使数字滤波器频响特性不能模拟原模拟滤波器的频响特性。 7 假设某模拟滤波器Ha(s)是一个低通滤波器， 又知， 数字滤波器H(z)的通带中心位于下面哪种情况？并说明原因。 （1） =0（低通）。 （2） =（高通）。 （3） 除0或以外的某一频率（带通）。,解： 方法一 按题意可写

10、出,故,即,原模拟低通滤波器以=0为通带中心， 由上式可知， =0时， 对应于=， 故答案为(2)。,方法二 找出对应于=0的数字频率的对应值即可。 令z=1， 对应于ej=1， 应有=0， 则对应的不是模拟低通中心频率， 所以， 答案（1）即=0（低通）不对。 令z=1， 对应于ej=1， 应有=， 则即将=0映射到=处， 所以答案为(2)。 ,方法三 直接根据双线性变换法设计公式及模拟滤波器由低通到高通频率变换公式求解。 双线性变换设计公式为,当T=2时，,则H(z)亦为低通。 如果将Ha(s)变换为高通滤波器：,， 这时，如果Ha(s)为低通，,则可将Hah(s)用双线性变换法变成数字高

11、通；,这正是题中所给变换关系，所以数字滤波器,通带,中心位于=， 故答案（2）正确。8 题8图是由RC组成的模拟滤波器， 写出其系统函数Ha(s)， 并选用一种合适的转换方法， 将Ha(s)转换成数字滤波器H(z)， 最后画出网络结构图。,题 8 图,解: 模拟RC滤波网络的频率响应函数为,显然， Ha(j)具有高通特性， 用脉冲响应不变法必然会产生严重的频率混叠失真。 所以应选用双线性变换法。 将Ha(j)中的j用s代替， 可得到RC滤波网络的系统函数：,用双线性变换法设计公式, 可得,H(z)的结构图如题8解图所示。,题8解图,由图可见， 在模拟域由一个R和一个C组成的RC滤波网络， 用双

12、线性变换法转换成数字滤波器后， 用两个乘法器、 两个加法器和一个单位延迟器实现其数字滤波功能。 也可用软件实现该数字滤波功能。 由滤波器差分方程编写程序较容易。 为此， 由H(z)求出差分方程。,编程序实现差分方程中的计算， 即可实现对输入信号序列x(n)的高通滤波。 9 设计低通数字滤波器， 要求通带内频率低于0.2 rad时， 容许幅度误差在1 dB之内； 频率在0.3到之间的阻带衰减大于10 dB。 试采用巴特沃斯型模拟滤波器进行设计， 用脉冲响应不变法进行转换， 采样间隔T=1 ms。,解: 本题要求用巴特沃斯型模拟滤波器设计， 所以， 由巴特沃斯滤波器的单调下降特性， 数字滤波器指标

13、描述如下： p=0.2 rad, p=1 dBs=0.3 rad, s=10 dB采用脉冲响应不变法转换， 所以， 相应的模拟低通巴特沃斯滤波器指标为,(1) 求滤波器阶数N及归一化系统函数G(p)：,取N=5。 查教材6.1节的表6.2.1(第157页)， 可知模拟滤波器系统函数的归一化低通原型为,将G(p)部分分式展开：,其中， 系数为 A0=0.1382+j0.4253, A1=0.8091j1.1135, A2=1.8947 A3=0.8091+j1.1135, A4=0.1382j0.4253(2) 去归一化求得相应的模拟滤波器系统函数Ha(s)。 我们希望阻带指标刚好， 让通带指标

14、留有富裕量， 所以按教材（6.2.20）式求3 dB截止频率c为,其中, Bk=cAk, sk=cpk。 (3) 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z):,我们知道， 脉冲响应不变法的主要缺点是存在的频率混叠失真， 使设计的滤波器阻带指标变差。 另外， 由该题的设计过程可见。 当N较大时， 部分分式展开求解系数Ak或Bk相当困难， 所以实际工作中用得很少， 主要采用双线性变换法设计， 见第10题。,10 要求同题9， 试采用双线性变换法设计数字低通滤波器。 解: 已知条件如下： 数字滤波器指标： p=0.2 rad, p=1 dBs=0.3 rad, s=10 dB采用

15、双线性变换法， 所以要进行预畸变校正， 确定相应的模拟滤波器指标（为了计算方便， 取T=1 s）：,(1) 求相应模拟滤波器阶数N：,其中， ksp与题9相同（因为p、 s相同）， 即,(2) 查教材表6.2.1， 得,(3) 去归一化， 求出Ha(s)：,(4) 用双线性变换法将Ha(s)转换成H(z)：,请读者按T=1 ms进行设计， 比较设计结果。,11 设计一个数字高通滤波器， 要求通带截止频率p=0.8 rad， 通带衰减不大于3 dB，阻带截止频率s=0.5 rad， 阻带衰减不小于18 dB。 希望采用巴特沃斯型滤波器。 解: （1） 确定数字高通滤波器技术指标： p=0.8 rad, p=3 dB s=0.5 rad, s=18 dB （2） 确定相应模拟高通滤波器技术指标。 由于设计的是高通数字滤波器， 所以应选用双线性变换法， 因此进行预畸变校正求模拟高通边界频率（假定采样间隔T=2 s）：,