数字信号处理21-6.2IIR DF设计(一).pptVIP

* * * * Butterworth模拟滤波器幅频的最高值为1，故在画幅频图时没有做归一化操作。 * * 6.2 IIR数字滤波器设计 基本概念 模拟低通滤波器的设计 冲激响应不变法 双线性变换法 设计IIR滤波器的频率变换法 * 一、 基本概念 数字滤波器的理想幅频特性 连续 周期 理想 * 1 滤波器的技术指标 理想特性 实际特性 通带容限 阻带容限 通带纹波 通带截止频率?p 阻带截止频率?s 通带 阻带 过渡带 ? 低通滤波器频率响应幅度特性的容限图 * 在通带内，1 - ? 1 ≤ | H (e j? )| ≤ 1+ ? 1 ，| ? | ≤ ?p 在阻带内，| H (e j? )| ≤ ? 2，?s ≤ | ? | ≤ ? 在具体技术指标中往往使用通带允许最大衰减(波纹) ? p和阻带应达最小衰减 (波纹) ? s描述。 1 滤波器的技术指标 * 2 滤波器的主要参数 幅度平方响应 当单位冲激响应 h (n) 为实序列时 显然 H (z) H (z -1) 的极点 (零点) 关于单位圆 | z | = 1镜像对称。 相位响应 群延迟响应 * 3 IIR滤波器设计方法 DF设计步骤 (1) 确定滤波器设计指标 (2) 用因果稳定系统逼近设计指标 (3) 用有限精度运算实现这个系统 (4) 实际的技术实现。 * 3 IIR滤波器设计方法 间接法：利用模拟滤波器的理论来设计数字滤波器. (1) 将数字滤波器的设计转换为模拟滤波器的设计。 (2) 设计满足技术指标的模拟滤波器。 (3) 将模拟滤波器转换为数字滤波器。 ?p, ?s ?p,?s H(s) H(z) 频率 变换 设计模拟 滤波器 冲激响应 不变法 双线性变换法 IIR DF实现方法 * 3 IIR滤波器设计方法 直接法：最优化设计法 [步骤] (1) 选择一种最优准则； (2) 求在此最优准则下滤波器系统函数的系数 ak，bm 。 最优化方法又称为计算机辅助设计法。 * 二、模拟低通滤波器的设计 实际低通滤波器的技术指标 幅度平方函数的确定 由幅度平方函数确定系统函数 巴特沃思低通滤波器的实际设计步骤 MATLAB实现 重点：巴特沃思低通滤波器设计 * 1 实际低通滤波器的技术指标 1 0 通带 过渡带 阻带 ) ( W j e H p W s W s ? p ? - 1 W Ws: 阻带截止频率 ? p: 通带波动 ? s: 阻带波动 通带衰减(dB) 阻带衰减(dB) Wp: 通带截止频率 * 1 实际低通滤波器的技术指标 模拟低通滤波器的设计步骤： 1.由给定的技术指标确定幅度平方函数 2.由幅度平方函数确定系统函数 * 2 幅度平方函数的确定 幅度平方函数的几种常用类型 巴特沃思(Butterworth)滤波器：滤波器具有单调下降的幅频特性。 切比雪夫(Chebyshev)滤波器：幅频特性在通带或者在阻带有波动，可以提高选择性。 贝塞尔(Bessel)滤波器：通带内有较好的线性相位特性。 椭圆(Ellipse)滤波器：选择性相对前三种是最好的，但在通带和阻带内均为等波纹幅频特性。 * 2 幅度平方函数的确定 (1)巴特沃思低通滤波器（BWF）的幅度平方函数定义 式中，N = 1, 2, 3, … 为滤波器阶数； ?c ：3dB截止频率。 当 ? = ?c 时 * 2 幅度平方函数的确定 (2)频率响应的特点 | Ha(j?)|单调下降； 阶次 N 越高，特性越接近矩形； BWF幅度特性 N=8 N=4 N=2 | Ha(j?)|2 在 ? = 0 点的1至2N-1阶导数都为零，称为最大平坦性； 任意阶次 N 下都有 | Ha(j?c)| = 0.707。 * 2 幅度平方函数的确定 (3)参数的确定： ①已知 对于巴特沃思低通滤波器： * 2 幅度平方函数的确定 (3)参数的确定： ①已知 对N向上取整 * 2 幅度平方函数的确定 (3)参数的确定： ①已知 将N带入原方程： 整理化简得： 阻带或通带衰减公式 * 2 幅度平方函数的确定 ②已知 得： 对N向上取整 (3)参数的确定： 3dB截止频率 * 3 由幅度平方函数确定系统函数 令 得： 零极点关于虚轴对称 (1) 系统函数确定 * 3 由幅度平方函数确定系统函数 由 幅度平方函数确定系统函数的步骤： ? 2 = -s2 ， 1.对 作代换 得 2.对 分子、分母分别进行因式分解， 得其零、极点； 3.选取s左