MatlabMatlab在信号处理中的应用分析.ppt

【例10.43】设计一个带阻数字滤波器，要求 利用椭圆滤波器原型。 解： %例10.43的程序：eg10_43.m %设计IIR数字带阻滤波器 %数字带阻滤波器的技术指标 wl=0.19*pi;wsl=0.198*pi;ws2=0.202*pi;wu=0.21*pi; Rp=3;Rs=13; %设置周期为1 T=1; * %用双线性变换法转换为模拟带阻滤波器的指标 Omegau=2/T*tan(wu/2);Omegal=2/T*tan(wl/2); Omegas1=2/T*tan(ws1/2);Omegas2=2/T*tan(ws2/2); Omega0=sqrt(Omegau.*Omegal); B=Omegau-Omegal; %转换为归一化模拟低通的指标 yitas2=Omegas2/B;yita0=Omega0/B; lanmdas=yitas2/(yitas2.^2-yita0.^2);lanmdap=1; %设计归一化模拟低通滤波器 [N,Wn]=ellipord(lanmdap,lanmdas,Rp,Rs,s) [z0,p0,k0]=ellipap(N,Rp,Rs); p=p0*Wn;z=z0*Wn; k=k0*Wn^N; ba=k*real(poly(z)); aa=real(poly(p)); * %转换为模拟带阻 w0=(Omega