DSP0.pptVIP

第10章 IIR滤波器 * 本章作业 2 5 13 16 18 23 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * * 2011年4月1日星期五 * 2010年5月17日星期一; 11周：2009年4月27日星期一; 12周：2009年5月4日星期一 12周：2008年5月12日星期一; 13周：2008年5月19日星期一 * 13周：2008年5月19日星期一 * Passatt：通带衰减dB Stopatt：阻带衰减dB * 第10章 IIR滤波器 * 第10章 IIR滤波器 * 第10章 IIR滤波器 * 第10章 IIR滤波器 * 10.6 脉冲响应不变法 IIR 滤波器设计 每个模拟滤波器都有冲激响应h(t)，就像每个数字滤波器都有脉冲响应h[n]一样。脉冲响应不变法选择的数字脉冲响应h[n]是满足设计要求的模拟滤波器冲激响应h(t)的采样值，如图10.20所示。即： h[n]=h(nT) (10.13) 其中T为所用的采样间隔。用这种方式得到的被采样的冲激响应将给出与原模拟滤波器非常相近的滤波器形状。 第10章 IIR滤波器 * 具有传输函数H(s)模拟滤波器的冲激响应h(t)可以求取，求取的方法超出了本书的范围。 只要用于实现脉冲响应不变法数字滤波器的采样间隔T很小，尽量减小混叠，数字和模拟滤波器的形状就会非常接近，如例10.12所示。 第10章 IIR滤波器 * 脉冲响应不变法设计步骤 设计步骤 由H(s)用拉氏反变换法计算出脉冲响应h(t)（本课程不介绍） 计算h[n]=h(nT) 见图10.20（P312 Matlab) 求出数字滤波器的传递函数H(z) Matlab方法 [bz, az] = impinvar(num, den, fs) H(z) = filt(bz, az)*fs num、den为模拟滤波器的系数 第10章 IIR滤波器 * 例10.12 用脉冲响应不变法设计截止频率为750 Hz的一阶巴特沃斯滤波器。 解： 由式(10.8)，截止频率为750 Hz的一阶巴特沃斯滤波器的传输函数为： 注意，因为不使用双线性变换，所以这里不需要频率预扭曲。由式(10.9)，得模拟滤波器的滤波器形状： 该特性由脉冲响应不变法设计的数字滤波器来逼近。由式(10.14)，数字滤波器的传输函数为： 第10章 IIR滤波器 * 由式(10.15)，数字滤波器的频率响应为： 该表达式给出了采样间隔为T时的幅度响应(滤波器形状)。图10.21比较了三种不同采样间隔下的模拟和数字滤波器形状，所有滤波器的直流增益都定标为1。 第10章 IIR滤波器 * 当采样间隔减小[如图(b)和(c)所示]，数字滤波器的形状越来越逼近模拟滤波器形状。 当采样率1/T比待求的截止频率大许多倍时，可得到最佳逼近。 第10章 IIR滤波器 * 10.7 “最佳拟合滤波器设计 正如例题中所示，递归滤波器的设计要处理非常复杂的方程。因此，完成这种设计几乎总要借助于软件的支持。 完成自动设计过程的一种方法是编制程序以实现前面节次的方程。对于待求滤波器，当方程变得过分复杂或没有方便的模拟滤波器原型相对应时，最佳的选择是通过迭代来逼近正确的滤波器。 采用这种方法的程序首先猜测一个传输函数，然后计算该滤波器形状，并与待求滤波器形状进行数学上的比较：当出现失配时，调整滤波器系数值并重新计算滤波器形状。这一过程持续进行，直到找出满足所需滤波器形状的传输函数为止。 第10章 IIR滤波器 * 最佳“拟合”设计法概述 在工具软件的帮助下，找出与期望的频率特性很相似的数字滤波器的传递函数。 猜测H(z)?比较频率特性?满足要求的H(z)。 第10章 IIR滤波器 * 10.8 带通、高通和带阻IIR滤波器 方法：先设计模拟低通滤波器，再转换为模拟带通、高通和带阻滤波器，最后离散化为数字滤波器。 若采用双线性变换法进行离散化，设计时要考虑频率的预扭曲问题。 第10章 IIR滤波器 * 截止频率为ωp的低通模拟滤波器HL(s)转换为带通、高通和带阻模拟IIR 第10章 IIR滤波器 * 带通、高通和带阻滤波器的Matlab设计 第10章 IIR滤波器 * 设计举例：例10.13 第10章 IIR滤波器 * 第10章 IIR滤波器 * 第10章 IIR滤波器 * 该滤波器的幅度响应见图10.24 第10章 IIR滤波