数字信号处理--第6章无限脉冲响应数字滤波器的设计精品.pptVIP

* 数字信号处理 将上式代入(5)中的转换公式，得 将上面的p等式代入G(p)中，得 * 数字信号处理 例6.5.3设计一个数字带阻滤波器，通带下限频率ωl=0.19π,阻带下截止频率ωs1=0.198π，阻带上截止频率ωs2=0.202π，通带上限频率ωu=0.21π，阻带最小衰减αs=13dB，ωl和ωu处衰减αp=3dB。采用巴特沃斯型。 解 (1) 数字带阻滤波器技术指标： ωl=0.19πrad,ωu=0.21πrad,αp=3dB; ωs1=0.198πrad,ωs2=0.202πrad,αs=13dB * 数字信号处理 (2) 模拟带阻滤波器的技术指标： 设T=1，则有 阻带中心频率平方为 Ω20=ΩlΩu=0.421 阻带带宽为 ?B=Ωu-Ωl=0.07rad/s * 数字信号处理 将以上边界频率对B归一化： ηl=8.786,ηu=9.786, ηs1=9.186,ηs2=9.386; η20=ηlηu=85.98 (3) 模拟归一化低通滤波器的技术指标： 按照(6.2.48)式，有 λp=1,αp=3dB * 数字信号处理 (4) 设计模拟低通滤波器： (5) 将G(p)转换成模拟阻带滤波器Ha(s)： * 数字信号处理 (6) 将Ha(s)通过双线性变换，得到数字阻带滤波器H(z)。 * 数字信号处理 6.6 IIR 数字滤波器的直接设计法 1. 零极点累试法 称为零极点累试法。在确定零极点位置时要注意： (1)极点必须位于z平面单位圆内，保证数字滤波器因果稳定； (2)复数零极点必须共轭成对，保证系统函数有理式的系数是实的。 * 数字信号处理 图6.6.1 例6.6.1图 (a)零极点分布； (b)幅度特性 * 数字信号处理 2.在频域利用幅度平方误差最小法直接设计IIR数字滤波器 设IIR滤波器由K个二阶网络级联而成，系统函数用H(z)表示， (6.6.1) 式中，A是常数；ai,bi,ci,di是待求的系数；Hd(e jω) 是希望设计的滤波器频响。如果在(0，π)区间取N点数 字频率ωi,i=1,2,:,N，在这N点频率上，比较|Hd(e jω)|和 |H(e jω)|，写出两者的幅度平方误差E为 (6.6.2) * 数字信号处理 而在(6.6.1)式中共有(4K+1)个待定的系数，求它们的原则是使E最小。下面我们研究采用(6.6.1)式网络结构，如何求出(4K+1)系数。 按照(6.6.2)式，E是(4K+1)个未知数的函数，用下式表示： 上式θ表示4K个系数组成的系数向量。为推导公式方便，令 (6.6.3) * 数字信号处理 为选择A使E最小，令 (6.6.4) 设θk是 θ 的第k个分量(ak或bk或ck或dk)， (6.6.5) 因为 ,式中H*i表示对Hi函数共轭。 (6.6.6) * 数字信号处理 将上式具体写成对ak,bk,ck,dk的偏导，得到： (6.6.7) 式中，k=1,2,3,:,K；i=1,2,3,:,N。 同理求得 (6.6.8) (6.6.9) (6.6.10) * 数字信号处理 由于系统函数是一个有理函数，极、零点均以共轭成对的形式存在，对于极点z1,一定有下面关系： (6.6.11) * 数字信号处理 图6.6.2 例6.6.2图 (a)要求的幅度特性；(b)k=1，2时的幅度特性 例6.6.2 设计低通数字滤波器，其幅度特性如图 6.6.2(a)所示。截止频率ωs=0.1πrad。 * 数字信号处理 解 考虑到通带和过渡带的重要，在0~0.2π区间，每隔