《数字信号处理》第六章 IIR滤波器的设计方法【荐】.ppt

为使模拟滤波器某一频率与数字滤波器的任一频率有对应关系，引入系数 c 2）某一特定频率严格相对应： 1）低频处有较确切的对应关系： 特定频率处频率响应严格相等，可以较准确地控制截止频率位置 二、变换常数c的选择 三、逼近情况 1） s平面虚轴 z平面单位圆 2） 左半平面 单位圆内 s平面 z平面 右半平面 单位圆外 虚轴 单位圆上 四、优缺点 优点： 避免了频率响应的混迭现象 s 平面与 z 平面为单值变换 缺点：除了零频率附近，W与w之间严重非线性 2）要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变 1）线性相位模拟滤波器 非线性相位数字滤波器 分段常数型模拟滤波器经变换后仍为分段常数型数字滤波器，但临界频率点产生畸变 预畸变 给定数字滤波器的截止频率 w1，则 按W1设计模拟滤波器，经双线性变换后，即可得到w1为截止频率的数字滤波器 五、模拟滤波器数字化方法 可分解成级联的低阶子系统 可分解成并联的低阶子系统 6.8 常用模拟低通滤波器特性 将数字滤波器技术指标转变成模拟滤波器技术指标，设计模拟滤波器，再转换成数字滤波器 模拟滤波器 巴特沃斯 Butterworth 滤波器 切比雪夫 Chebyshev 滤波器 椭圆 Ellipse 滤波器 贝塞尔 Bessel 滤波器 1、由幅度平方函数 确定模拟滤波器的系统函数 h(t)是实函数 将左半平面的的极点归Ha(s) 将以虚轴为对称轴的对称零点的任一半作为Ha(s)的零点，虚轴上的零点一半归Ha(s) Ha(s) Ha(-s)的零极点分布 由幅度平方函数得象限对称的s平面函数 对比 和 ，确定增益常数 由零极点及增益常数，得 例： 解： 极点： 零点： （二阶） 零点： 的极点： 设增益常数为K0 （1） 巴特沃尔斯滤波器（Butterworth） 2、常见模拟滤波器设计 1）幅度平方函数： 当 称Wc为Butterworth低通滤波器的3分贝带宽 N为滤波器的阶数 Wc为通带截止频率 1）幅度函数特点： 3dB不变性 通带内有最大平坦的幅度特性，单调减小 过渡带及阻带内快速单调减小 当W＝Wst（阻带截止频率）时，衰减的d1为阻带最小衰减 Butterworth滤波器是一个全极点滤波器，其极点： 2）幅度平方特性的极点分布： 极点在s平面呈象限对称，分布在Buttterworth圆上，共2N点 极点间的角度间隔为 极点不落在虚轴上 N为奇数，实轴上有极点，N为偶数，实轴上无极点 Ha(s) Ha(-s)的零极点分布 (a) N=3 (三阶） (b)N=4 （四阶） 为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(-s)。 Ha(s)的表示式为 设N=3，极点有6个（见前页图），它们分别为：（下页） 3）滤波器的系统函数： 其中分子系数 由Ha(0)=1解得。 取s平面左半平面的极点s0,s1,s2组成Ha(s)： 由前讨论知，系统函数为： 极点为： 当N为偶数，Ha(s)的极点（左半平面）皆成共轭对，记为： 这一对共轭点构成一个二阶子系统： 当N为奇数，则系统由一个一阶系统（极点s=-1）和(N-1)/2个二阶系统组成： 整个系统由N/2（N为偶数）个此二阶系统级联而成： N为偶数 N为奇数 为归一化系统的系统函数 去归一化，得 例：试设计一个模拟低通巴特沃思滤波器，要求通带截止频率Ωc=2π×4000rad/s，通带最大衰减δ1=3dB，阻带下限截止频率Ωst=2π×8000rad/s，阻带最小衰减δ2=20dB. 解： （1）求阶数N： 联立求解，（过程见教材）得： 因此可取N=4； （2）求极点： 由式： （3）求系统函数Ha(s)： 实际上，求出N=4时，可以直接查表（表6-4），得到归一化的（Ωc=1）4阶巴特沃思低通滤波器的系统函数，然后再用s/Ωc代替其中的s即可得到同样的结果。 4）滤波器的设计步骤： 根据技术指标求出滤波器阶数N： 确定技术指标： 由 得： 同理： 令 则： 求出归一化系统函数： 或者由N，直接查表得 其中技术指标Wc给出或由下式求出： 其中极点： 去归一化 阻带指标有富裕 或 通带指标有富裕 逆因果稳定系统 1）全部零点在单位圆内：