用MATLAB设计IIR数字滤波器-实验七

;set(gca,'YTickMode','manual','YTick',[0,Attn,ripple,1]);gridsubplot(2,2,2);plot(w/pi,angle(H)/pi);ylabel('\phi');title('相位响应');axis([0,1,-1,1]);set(gca,'XTickMode','manual','XTick',[0,0.2,0.35,1]);set(gca,'YTickMode','manual','YTick',[-1,0,1]);gridsubplot(2,2,3);plot(w/pi,dbH);title('幅度响应(dB)');ylabel('dB');xlabel('频率(\pi)');axis([0,1,-40,5]);set(gca,'XTickMode','manual','XTick',[0,0.2,0.35,1]);set(gca,'YTickMode','manual','YTick',[-50,-15,-1,0]);gridsubplot(2,2,4);zplane(bd,ad);axis([-1.1,1.1,-1.1,1.1]);title('零极点图');程序运行结果如下:n=3Omgc=7.7867bd=0.5221-1.35311.3531-0.5221ad=1.0000-1.45101.0283-0.2711由频率特性曲线可知,该设计结果再通阻带截止频率处能满足Rp≤1dB、As≥15dB的设计指标要求,系统的极点全部在单位圆内,是一个稳定系统。由n=3可知,由3阶的模拟低通滤波器原型用双线性变换法设计出来的切比雪夫Ⅱ型数字高通滤波器是一个3阶系统。其传递函数为:设计一个切比雪夫Ⅱ型数字带通滤波器,要求:;阻带,滤波器采样周期Ts=0.001s。列出传递函数并作频率响应曲线和零极点分布图。