[信息与通信]第6章 无限脉冲响应数字滤波器的设计.ppt

图6.6.1 例6.6.1图 (a)零极点分布； (b)幅度特性 * 解 (1)数字带通滤波器技术指标为 通带上截止频率 ωu=0.4πrad 通带下截止频率 ωl=0.3πrad 通带内最大衰减αp=3dB， 阻带上截止频率ωs2=0.5πrad 阻带下截止频率ωs1=0.2πrad 阻带内最小衰减αs=18dB。 * (2) 模拟带通滤波器技术指标如下： 设T=1，则有 (通带中心频率) (带宽) * 将以上边界频率对带宽B归一化，得到 ηu=3.348,ηl=2.348; ηs2=4.608,ηs1=1.498; η0=2.804 (3) 模拟归一化低通滤波器技术指标： 归一化阻带截止频率 归一化通带截止频率 λp=1 αp=3dB,αs=18dB * (4) 设计模拟低通滤波器： 查表6.2.1，得到归一化低通传输函数G(p), * (5) 将归一化模拟低通转换成模拟带通： (6)通过双线性变换法将Ha(s)转换成数字带通滤波器H(z)。 下面将(5)、(6)两步合成一步计算： 将上式代入(5)中的转换公式，得 将上面的p等式代入G(p)中，得 * 带通数字滤波器的MATLAB实现例6.5.2wl=0.3*pi; wu=0.4*pi; ap=3; %通带ws1=0.2*pi; ws2=0.5*pi; as=18; %阻带T=1;wl = (2/T)*tan(wl/2);%数字滤波器指标转换为模拟滤波器指标wu = (2/T)*tan(wu/2);ws1= (2/T)*tan(ws1/2);ws2= (2/T)*tan(ws2/2);wo = sqrt(wl*wu);B=wu-wl;[N,wc] = buttord([wl wu],[ws1 ws2],ap,as,s)[z,p,k] = buttap(N);[b,a] = zp2tf(z,p,k);[bt,at] = lp2bp(b,a,wo,B)figure;freqs(bt,at);title(模拟带通滤波器);[bd,ad] = bilinear(bt,at,1/T);figure;freqz(bd,ad);title(双线性不变法设计数字带通滤波器); * 例 设计一个数字带阻滤波器，通带下限频率ωl=0.19π,阻带下截止频率ωs1=0.198π，阻带上截止频率ωs2=0.202π，通带上限频率ωu=0.21π，阻带最小衰减αs=13dB，ωl和ωu处衰减αp=3dB。采用巴特沃斯型。 * (2) 模拟带阻滤波器的技术指标： 设T=1，则有 阻带中心频率平方为 Ω20=ΩlΩu=0.421 阻带带宽为 ?B=Ωu-Ωl=0.07rad/s 解 (1) 数字带阻滤波器技术指标： ωl=0.19πrad,ωu=0.21πrad,αp=3dB; ωs1=0.198πrad,ωs2=0.202πrad,αs=13dB * 将以上边界频率对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)。 * 带阻数字滤波器的MATLAB实现例6.5.3 wl=0.19*pi; wu=0.21*pi; ap=3; %通带ws1=0.198*pi; ws2=0.202*pi; as=13; %阻带T=1;wl = (2/T)*tan(wl/2);%数字滤波器指标转换为模拟滤波器指标wu = (2/T)*tan(wu/2);ws1= (2/T)*tan(ws1/2);ws2= (2/T)*tan(