无限脉冲响应数字滤波器的设计.ppt

例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第157页,共183页，星期六，2024年，5月(2)模拟带阻滤波器的技术指标：设T=1，则有阻带中心频率平方为Ω20=ΩlΩu=0.421阻带带宽为?B=Ωu-Ωl=0.07rad/s第158页,共183页，星期六，2024年，5月将以上边界频率对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第159页,共183页，星期六，2024年，5月(4)设计模拟低通滤波器：(5)将G(p)转换成模拟阻带滤波器Ha(s)：第160页,共183页，星期六，2024年，5月(6)将Ha(s)通过双线性变换，得到数字阻带滤波器H(z)。第161页,共183页，星期六，2024年，5月wp=[0.19*pi,0.21*pi];ws=[0.198*pi,0.202*pi];rp=3;rs=13;[N,wpo]=buttord(wp,ws,rp,rs,s);[b,a]=butter(N,wpo,stop);freqz(b,a);第162页,共183页，星期六，2024年，5月wp=0.7*pi;ws=0.5*pi;Rp=1;Rs=32;Ts=2;Wp=(2/Ts)*tan(wp/2);Ws=(2/Ts)*tan(ws/2);[N,Wc]=cheb1ord(Wp,Ws,Rp,Rs,s);[z,p,k]=cheb1ap(N,1)[Bp,Ap]=zp2tf(z,p,k);[b,a]=lp2hp(Bp,Ap,1.9261);[bz,az]=bilinear(b,a,1/Ts);freqz(bz,az);第163页,共183页，星期六，2024年，5月第164页,共183页，星期六，2024年，5月fpl=2025,fpu=2225;fsl=1500;fsu=2700;fs=8000;rp=1;rs=40wp=[2*fpl/fs,2*fpu/fs];ws=[2*fsl/fs,2*fsu/fs];[N,wpo]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N,rp,rs,wpo);freqz(B,A)第165页,共183页，星期六，2024年，5月第166页,共183页，星期六，2024年，5月fsl=2025,fsu=2225;fpl=1500;fpu=2700;fs=8000;rp=1;rs=40wp=[2*fpl/fs,2*fpu/fs];ws=[2*fsl/fs,2*fsu/fs];[N,wpo]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N,rp,rs,wpo,stop);freqz(B,A)第167页,共183页，星期六，2024年，5月第168页,共183页，星期六，2024年，5月解(1)用脉冲响应不变法设计数字低通滤波器。①数字低通的技术指标为ωp=0.2πrad,αp=1dB;ωs=0.3πrad,αs=15dB②模拟低通的技术指标为T=1s,Ωp=0.2πrad/s,αp=1dB;Ωs=0.3πrad/s,αs=15dB第125页,共183页，星期六，2024年，5月③设计巴特沃斯