信号、系统与数字信号处理 第8章 无限冲激响应(IIR)数字滤波器的设计.ppt

所以 由此得到： 也可以查表3.9-1得到 的归一化 为 去归一化后 (3)求 从以上求解过程可见 第(2)、 (3)步的计算量相当于大， 借助MATLAB我们可以得到所需的结果。 方法二 在第(1)步求出 后，从六阶巴氏模拟原型低通 、 开始用MATLAB程序设计。 用模拟低通原型滤波器的指标，全部设计由MATLAB 方法三 程序完成。 通带截止频率 例8-10设计指标同例8-9用双线性变换法设计一个切比 ，通带最大衰减 阻带边缘频率 ，阻带最小衰减 令 解： ，得 (1)求 、 、 雪夫数字低通滤波器， ，与巴特沃思滤波器相比，指标相同，但阶数 取 减少。 (2) (3) 确定零极点 (4)确定模拟低通系统函数 可以利用表3.9-1得到归一化 ，再去归一化得到 (5)由模拟低通系统函数经双线性变换法确定数字低通系 统函数 骤一样，就是模拟低通原型的 与用冲激不变法设计切比雪夫数字滤波器相比，基本步 由图8-13与图8-11巴特沃思滤波器振幅频响相比，阻 带指标有所改善。这是因为切比雪夫数字滤波器实际 阶数N取4，比计算所需的理论阶数3.016有较大的富 余量。 。 ， §8.4 原型变换法 前面两节讨论的是由巴特沃思、切比雪夫模拟原型低通 滤波器设计数字低通的方法。而实际待求的数字滤波器 由模拟低通原型设计任意的DF，有三种方法可以选择。 有各种不同的低通、高通、带通、带阻滤波器，本节着 重讨论由模拟低通原型设计实际数字滤波器的方法。 如图8-14所示 即 （只适用双线性变换） 所需数字滤波器 三、模拟低通原型 即 所需数字滤波器 所需模拟 二、模拟低通原型 即 一、模拟低通原型 数字低通 所需数字滤波器 由图可见，由模拟低通原型滤波器出发，设计数字滤波 器的一种方法是由模拟低通原型滤波器 激不变法或双线性变换法得到数字低通滤波器的 再由数字低通得到所需的数字滤波器的 ，用冲 。这种 方法的第一步，实际就是前面两节讨论的s平面与z平面 的映射变换。关键是第二步，第二步实质是数字域z平 面之间的变换，这种变换也称z平面变换法。 设计数字滤波器的第二种方法是由模拟低通原型滤波器 ，设计所需的模拟低通、高通、带通、带阻滤 波器 ，再由 经冲激不变法或双线性变 换法得到所需的数字滤波器的 。这种方法的第 二步，是前面两节讨论的s平面与z平面的映射变换。 换，这种变换也称s平面变换法。 关键是第一步，第一步的实质是模拟域s平面之间的变 设计数字滤波器的第三种方法是由模拟低通原型滤波器 ，直接设计所需的数字滤波器的 。由 于冲激不变法频谱的混叠效应，应用受到一定限制，所 以此法适用双线性变换法。因为这种方法直接由模拟原 型滤波器变换到数字滤波器所以也称原型变换法。 分别讨论与这三种方法相关的z平面变换法、 s平面变 换法以及原型变换法。 z平面变换或s平面变换的作用都是改变原滤波器的频率 1、z平面变换法——数字域的频率变换 特性，其实质是频率变换 ，所以统称频率变换法。 的方法。因为 前面讨论了由模拟原型低通设计数字低通的方法。现在 为原型，通过适当的频率 讨论以已知数字低通 变换，设计其它所需数字滤波器 的自变量都是数字域变量， 与 这种变换的 所以称z平面变换法。 §8.4.1 频率变换法 z平面变换的作用是改变原滤波器的频率特性，其实质 是数字域的频率变换 ，所以属于频率变换法。 为了区分变换前后的两个不同的z平面，设变换前为小 平面，变换后为大 平面， 的映射关 系为 则 应满足的条件： 对映射关系的要求，或 且 由条件ⅱ） 可得 因为： Z平面的单位圆 z平面的单位圆 ⅱ）具有相应的频率关系 Z平面的单位圆内； z平面的单位圆内 ⅰ）稳定性不变 （稳定系 a.只要有一个极点 ，就有一个零点 。零、极点 成对出现，并且极点一定在单位圆内 统）。 在单位圆上恒为1，称为全通函数。 函数 任意一个全通函数可以表示为 全通函数特点： b.当 由 时，全通函数的相位 变化 量为 ， 是全通函数的阶数。 时 当 ——模恒为1 由 变化量为 可知， 应为一阶全通函数。 (1)低通到低通 均为低通，但截止频率不同。 其中： 当 为实数，且 时，可以实现所要求的映射关系 ， 所以 因为 以及 因为 即有： 验证： 频率变换关系：由 可以解出 关系曲线如图 8-15示。 若数字低通原型截止频率为 ，变换后为 ，代入 （8-31）式可得 确定 ，现在 时，变换频率扩展（原 ） 时，变换频率压缩（原 ，现在 ） 关系曲线可知 由 解出 上式分子、分母同乘 可得 只要将低通 低通变换的 ，单位圆上的 改为 频响就旋转了 ，数字低