第 1 页共 10 页实验 7 用 MATLAB设计 IIR 数字滤波器一、实验目的：1、加深对IIR 数字滤波器的基本设计方法的理解。2、掌握用模拟滤波器原型设计IIR 数字滤波器的方法。3、了解 MATLAB有关 IIR 数字滤波器设计的子函数的调用方法。二、实验内容及步骤1、阅读并输入实验原理中介绍的例题程序，观察输出的数据和图形，结合基本原理理解每一条语句的含义。2、用双线性变换法设计切比雪夫Ⅱ型数字滤波器，列出传递函数并描绘模拟和数字滤波器的幅频响应曲线。① 设计一个数字低通滤波器，要求：ω p=0.2П，Rp=1dB ；阻带： ω s=0.35П，As=15dB ，滤波器采样频率Fs=10Hz。wp=0.2*pi; %滤波器的通带截止频率ws=0.35*pi; %滤波器的阻带截止频率Rp=1;As=15; %滤波器的通阻带衰减指标ripple=10^(-Rp/20); % 滤波器的通带衰减对应的幅度值Attn=10^(-As/20); %滤波器的阻带衰减对应的幅度值%转换为模拟滤波器的技术指标Fs=10;T=1/Fs; Omgp=(2/T)*tan(wp/2);%原型通带频率的预修正Omgs=(2/T)*tan(ws/2);% 原型阻带频率的预修正%模拟原型滤波器计算[n,Omgc]=buttord(Omgp,Omgs,Rp,As,'s') %计算阶数 n 和截止频率[z0,p0,k0]=buttap(n); %设计归一化的巴特沃思模拟滤波器原型ba1=k0*real(poly(z0)); %求原型滤波器的系数b aa1=real(poly(p0)); %求原型滤波器的系数a 课程名称：数字信号处理实验成绩：指导教师：王丽霞实 验 报 告院系：信息工程学院班级：11 通信专升本学号：030100068 姓名：周海霞日期：2011 第 2 页共 10 页[ba,aa]=lp2lp(ba1,aa1,Omgc); %变换为模拟低通滤波器%也可将以上4 行替换为 [bb,aa]=butter(n,Omgc,'s'); 直接求模拟滤波器系数%用双线性变换法计算数字滤波器系数[bd,ad]=bilinear(ba,aa,Fs) [sos,g]=tf2sos(bd,ad) %转换成级联型%求数字系统的频率特性[H,w]=freqz(bd,ad); dbH=20*log10((abs(H)+eps)/max(abs(H))); subplot(2,2,1);plot(w/pi,abs(H)); ylabel('|H|');title(' 幅度响应 ');axis([0,1,0,1.1]); set(gca,'XTickMode','manual','XTick',[0,0.25,0.4,1]); set(gca,'YTickMode','manual','YTick',[0,Attn,ripple,1]);grid subplot(2,2,2);plot(w/pi,angle(H)/pi); ylabel('\phi');title(' 相位响应 ');axis([0,1,-1,1]); set(gca,'XTickMode','manual','XTick',[0,0.25,0.4,1]); set(gca,'YTickMode','manual','YTick',[-1,0,1]);grid subplot(2,2...
