数字巴特沃斯滤波器 19.pptVIP

数字信号处理 陈鹏 3299340 数字巴特沃斯滤波器 Digital Butterworth Filter Lesson 19 复习提问 数字滤波器设计的三大步骤？ 冲激响应不变法的设计思想是？ 冲激响应不变法哪个频段失真较大，为什么？它适用于设计哪种滤波器？ 冲激响应不变法适合哪种类型的系统函数？ 冲激响应不变法的稳定性如何？ 巴特沃斯滤波器 Butterworth Filter 巴特沃斯滤波器设计公式、步骤 冲激响应不变法设计数字巴特沃斯滤波器 双线性变换法设计数字巴特沃斯滤波器 设计公式 巴特沃斯(Butterworth)滤波器都幅度响应在通带内具有 最平坦的特性，且在通带和阻带内幅度特性是单调变化 的。模拟巴特沃斯滤波器的幅度平方函数为 幅度特性 幅度特性 幅度特性 现在来分析Butterworth滤波器极点的分布特点。如果用 s 代替 ，即经解析延拓，得到 由此得到极点 极点分布 用极点构建系统函数 用极点构建系统函数 用极点构建系统函数 Ａ值的确定 用极点构建系统函数 Ｎ为偶数时，模拟Butterworth滤波器的系统函数为 式中， 为左半平面的极点， 为 的共轭极点 N为奇数时，模拟Butterworth滤波器的系统函数为 　　　　为负实轴上的极点。 设计步骤 总结设计数字Butterworth滤波器的步骤如下 根据实际需要规定滤波器的数字截止频率 处的衰减，单位为 dB 由数字截止频率 　处的衰减计算模拟巴特沃斯滤波器的阶数N和频率 设计步骤 求模拟巴特沃斯滤波器的极点，并由 s 平面左半平面 　的极点构成系统函数 左半平面的极点： 系统函数： 使用冲激不变法或双线性变换法将 　转换成数字滤波器的系统函数 举例—冲激响应不变法 举例—冲激响应不变法 举例—冲激响应不变法 举例—冲激响应不变法 举例—冲激响应不变法 举例—双线性变换法 举例—双线性变换法 举例—双线性变换法 现在可以认为模拟频率是归一化的频率，即令 T=1，得 举例—双线性变换法 举例—双线性变换法 举例—双线性变换法 举例—双线性变换法 *五邑大学 信息学院 * 如下图所示为Butterworth滤波器的幅度特性，其中 为角频率，在 处幅度响应的平方为 0.5，N为滤波器的阶数，当 时，幅度响应为1。 从式子和图都可以看出，随着N的增大，幅度响应曲线在截止频率附近变得越来越陡峭，即在通带内有更大部分的幅度接近于1，在阻带内以更快的速度下降至零。Butterworth滤波器存在极点，而零点在 由此看出，巴特沃斯滤波器的极点分布特点：在 s 平面上共有2N个极点等角距地分布在半经为 的圆周上． 这些极点对称于虚轴，而虚轴上无极点；N为奇数时，实轴上有两个极点；N为偶数时，实轴上无极点；各个极点间的角度为 。图示为 N=3时各极点的分布情况。 知道巴特沃斯滤波器的极点分布后，便可以由 s 平面左半平面的极点构成系统函数 ，根据极点分布，可以得到 上式中， 　　是 s 平面左半平面的极点， 　　是右半平面的极点， 　A 和 B 都为常数。 Butterworth滤波器有 2N 个极点，且对称于虚轴，所以可将左半平面的极点分配给 ，以便得到一个稳定的系统，把右半平面的极点分配给 ， 不是所需要的，可以不管它，于是有巴特沃斯滤波器系统函数： N为偶数，A由滤波器在 处的单位冲激响应来确定，即 于是得到 N为奇数可得到一样的结果。 例 4.2 设计一个数字巴特沃斯滤波器，在通带截止频率 　　 　　　　　处衰减不大于 1dB，在阻带截止频率 处衰减不小于 15dB。 解 (1) 根据滤波器的指标得 (2) 设 T=1，将数字域指标转换成模拟域指标得 代入巴特沃斯滤波器的幅度平方函数得 解这两个方程得 按此值设计的滤波器满足通带指标要求，阻带指标将超 过给定值。 (3) 把 　代入式子 得到 s 平面左半平面的3对极点分别为： 由这3对极点构成的滤波器的系统函数为