西南交通大学信息科学与技术学院数字信号处理课件 第九章.pptVIP

第9章IIR 数字滤波器设计 9.1.1 数字滤波器指标 通常，需要设计的数字滤波器的幅度或相位（延时）响应是确定的。 在有些情况下，可能会指定滤波器的单位抽样响应或是阶跃响应。 在大多数实际应用中, 关键问题是用一个可实现的传输函数去逼近给定的滤波器幅度响应指标。 9.1.1 数字滤波器指标 在本章中，我们将仅讨论幅度逼近问题。 在这里我们指出四种基本的理想滤波器类型，它们的幅度响应分别如下图所示： 9.1.1 数字滤波器指标 这些理想滤波器都是不可实现的。 实际上, 数字滤波器在通带和阻带的幅度响应指标以一定容限的允许误差来表征。 另外，过渡带就是通带和阻带间的部分。 9.1.1 数字滤波器指标 例如, 给定理想低通滤波器的幅度响应 |G(ej?)| 如下图所示： 归一化的数字低通滤波器幅度响应指标如下图所示： 映射的非线性导致了频率轴的失真，称为频率畸变。 畸变的效果如右图所示： 例9.3 高通IIR数字滤波器的设计 切比雪夫1 型IIR 数字高通滤波器的设计 指标: FP=700 Hz, FS=500 Hz, αp=1 dB, αs=32 dB, FT= 2 kHz 确定归一化截止角频率为： 利用MATLAB 的代码来设计： [N, Wn] = cheb1ord(1, 1.9626105, 1, 32, ’s’) [B, A] = cheby1(N, 1, Wn, ’s’); [BT, AT] = lp2hp(B, A, 1.9626105); [num, den] = bilinear(BT, AT, 0.5); 原型模拟巴特沃兹低通滤波器的指标为: Ωp=1, Ωs=2.3617627 , αp=1dB, αs=40dB 利用MATLAB 编程代码设计如下： [N, Wn] = buttord(1, 2.3617627, 1, 40, ’s’) [B, A] = butter(N, Wn, ’s’); [BT, AT] = lp2bp(B, A, 1.1805647, 0.777771); [num, den] = bilinear(BT, AT, 0.5); 阻带带宽为 由 带阻IIR数字滤波器的设计 利用MATLAB 编程代码设计如下 [N, Wn] = ellipord(0.4234126, 1, 1, 40, ’s’); [B, A] = ellip(N, 1, 40, Wn, ’s’); [BT, AT] = lp2bs(B, A, 1.1805647, 0.777771); [num, den] = bilinear(BT, AT, 0.5); IIR滤波器的谱变换 IIR滤波器的谱变换 因此 IIR滤波器的谱变换 回想对于一个稳定的全通函数 A(z) 满足以下条件 IIR滤波器的谱变换 低通滤波器到低通滤波器的谱变换 上式中 的α 是关于给定的两个截止频率的函数。 低通滤波器到低通滤波器的谱变换 在单位圆上，上式变换可简化为： 低通滤波器到低通滤波器的谱变换 求解上式可得： 低通滤波器到低通滤波器的谱变换 这里 低通滤波器到低通滤波器的谱变换 低通滤波器到低通滤波器的变换如下所示： 低通滤波器到高通滤波器的谱变换 所求的变换为： 低通滤波器到高通滤波器的谱变换 例如 – 低通数字传输函数如下： 低通滤波器到高通滤波器的谱变换 可得到期望的高通滤波器的传输函数如下： 低通滤波器到高通滤波器的谱变换 低通滤波器到带通滤波器的谱变换 所求的变换如下所示： 低通滤波器到带通滤波器的谱变换 参数 α 和 β 由下式给出 低通滤波器到带通滤波器的谱变换 修正后为 低通滤波器到带阻滤波器的谱变换 所求变换如下式所示： 低通滤波器到带阻滤波器的谱变换 参数 α 和 β 由下式给出 作业 Read textbook from 427 to 454 Problems 9.2, 9.3, 9.6, 9.9, 9.15, 9.22, 9.25 M9.1, M9.5 作业 阅读教材 427 to 454 习题 9.2, 9.3, 9.6, 9.9, 9.15, 9.22, 9.25 M9.1, M9.5 为了避免混淆，我们用z-1 来表示原型低通数字滤波器GL(z) 的单位延时，并用 来表示变换后的滤波器 目的：将给定低通数字数字IIR的传输函数GL(z)转换成另一个低通、高通、带通或带阻滤波器传输函数 分别定义z平面和 平面如下所示： 我们将从z域到 域的变换表示如下 由 , 得 , 因此可得 因此函数