数字信号第六章.ppt

第六章 IIR数字滤波器的理论与设计 6.1 数字滤波器设计概述 IIR数字滤波器的设计方法 用模拟滤波器来设计数字滤波器的设计原理 6.2 模拟滤波器的数字仿真 时域的数字仿真 频域的数字仿真 6.3 脉冲响应不变法 模拟滤波器的系统函数: 模拟滤波器的复频率s与数字滤波器的复频率z之间的关系 因此，若给定模拟滤波器的系统函数，可求得数字滤波器的系统函数： S平面与Z平面的映射关系 分析数字滤波器的角频率： 6.5 模拟滤波器的逼近 问题的提出 1. 用模拟滤波器来设计数字滤波器的3步骤 2. 低通滤波器的平方幅度波形及其逼近 3. 根据平方幅度特性给出的性能指标设计模拟滤波器的系统函数Ha(s) 2. 最平坦函数 4. 由平方幅度函数求得模拟滤波器系统函数 由平方幅度函数得极点： 5. 一般情况下的低通滤波器 根据指标求模拟滤波器的系统函数（即确定N， ） 两个未知数，建立两个方程来求解 幅度平方函数： 通带边界频率处（归一化， ）： 阻带边界频率处： 低通滤波器平方幅度特性的一般形式： C型滤波器的思想：将零点分散开来，在通带内多个位置出现最大值，从而改善性能。 1.Chebyshev多项式 定义： （1）关于奇偶性 （2） 关于分界点 CN(1)=1 分界点连续 N：0～8时的Chebyshev多项式 （3）关于零点分布 2. C 型低通滤波器的幅频响应 3. C型滤波器的极点和系统函数 可以证明，这些根共2N个，它们成复共轭对出现； 关于虚轴对称，却没有一个在虚轴上； 将左半平面的极点分给Ha(S)，这N个极点为： 由此求得系统函数的分母，Q（S） 思想： 若在阻带内有有限s平面的传输零点，使其靠近通 带，就会使过渡带的衰减特性变陡。 N为偶数： 6.5.5 三种滤波器的比较 6.5.6 滤波器图表法设计 低通滤波器的指标 2.用列线图求滤波器的阶次N 查表求得系统函数 目的：由低通原型滤波器导出其它类型的滤波器。 频率响应函数：周期2*pi 幅频响应： 四类数字滤波器的理想幅频特性： 定义从一个z平面到另一个Z平面的映射： 可以证明，满足要求的映射函数的一般形式： 通过两个典型点的映射确定 和alpha 根据A点、B点的映射，和6.90式建立方程组，可求得两个参数: 椭圆上的极点还可以由作图法确定 再根据最高次幂作个修正，得C型滤波器的系统函数 平方幅度特性： Cauer滤波器简介 系统函数： N为奇数： 1. 关于振幅频率特性 2. 关于过渡带的陡度 3. 设计复杂性和滤波器频响特性对参数 变化的灵敏度。 通带允许的最大衰减 阻带要求的最小衰减 衰减： 幅度响应特性经常用衰减表示 分母多项式： 分子：下表的系数以 为基准频率，故 。因此分子为1 B型 C型 分母多项式：查表确定 分子： C型的系数以 为基准频率，故 。 因此分子为 6.6 数字滤波器的变换 任意的数字滤波器的设计过程 对数字滤波器： 当单位取样响应为实序列时，为偶函数 因此，对一数字滤波器的频率特性，可只考虑0～pi区间的情况 使得： 从原型数字低通滤波器到其它滤波器的变换： 映射的要求： 稳定因果的低通数字滤波器的原型，映射得到的数字滤波器必须也是稳定、因果的。 具体来说： 1.这个映射必须是一个有理函数 2.Z平面的单位圆内部必须映射到Z平面的单位圆内部 S平面与z平面的映射关系： 令： 非线性映射，整个虚轴映射到z平面单位圆的一周之上，一一对应，不存在频域混叠失真的问题 例：已知模拟滤波器的传输函数如下，试采用双线性变换法将其转换成数字滤波器，设T=2s 解： 5.3.2 巴特沃思低通逼近 巴特沃思低通滤波器幅度平方函数定义为 : 1. 3dB带宽：当Ω=0时，|Ha(j0)|=1； 当Ω=Ωc时，|Ha(jΩc)|2=1/2，10lg|Ha(jΩc)|2= -3dB， Ωc为3dB截止频率。当Ω=Ωc时，不管N为多少，所有的特性曲线都通过－3 dB点。 3. N的影响 即：角频率为0及很大时，幅频特性接近理想情况 这两处曲线趋于平坦，因此B型特性又叫最平坦特性 通带内，N越大，平方幅度特性随角频率增大下降越慢 阻带内，N越大，平方幅度特性随角频率增大下降越快 因此，N越大，B型滤波器幅频特性越接近理想的矩形 离截止频率越近，幅频特性与理想特性相差越大 平方幅度函数 极点为： 这2N个极点就是Ha(s)Ha(-s)的极点 而若sk是Ha(s)的极点，则-sk就是Ha(-s)的极点 将左半平面的N个极点sk分给Ha(s)，有 求出N， 后，按照公式6.33求解 标称化频率 （