数字信号处理Matlab实验三-IIR数字滤波器的设计.docVIP

数字信号处理Matlab实验三-IIR数字滤波器的设计 XXXX大学XXXX学院 课程名称： 数字信号处理 班级： 姓名： 学号： 实验地点： 日期： 实验名称 IIR数字滤波器的设计 实验目的： 加深理解IIR数字滤波器的时域特性和频域特性，掌握IIR数字滤波器的设计原理与设计方法，以及IIR数字滤波器的应用。 实验内容： IIR数字滤波器一般为线性移不变的因果离散系统，N阶IIR数字滤波器的系统函数可以表达为的有理多项式，即 式中：系数至少有一个非零。对于因果IIR数据滤波器，应满足。 IIR数字滤波器的设计主要通过成熟的模拟滤波器设计方法来实现。 首先在频域将数字滤波器设计指标转换为模拟滤波器设计指标，然后将任意的模拟滤波器为原型模拟低通滤波器指标，根据模拟滤波器的设计指标来设计出模拟低通滤波器，然后又经过相应的复频域转换得到H(s),最后又H(s)经过脉冲响应不变法或双线性变换法得到所需要的IIIR数字滤波器H(z)。 由此可见，IIR数字滤波器设计的重要环节是模拟滤波器的设计。设计模拟低通滤波器的主要方法有Butterwort、Chebyshev、和椭圆等滤波器设计方法。 实验步骤 Butterwort数字滤波器设计 (1) Butterwort滤波器是通带阻带都单调衰减的滤波器。调用buttord函数可以确定巴特沃斯滤波器的阶数，其格式为：[N,Omegac]=buttord(Omegap,Omegas,Rp,As,’s’)。 其中，输入参数Rp,As分别为通带最大衰减和阻带最小衰减，以dB为单位； Attn=1/(10^(As/20)); %阻带衰减 %模拟巴特沃斯原型滤波器计算 [N,OmegaC]=buttord(OmegaP,OmegaS,Rp,As,s); %原型的阶数和截止频率的计算 [z0,p0,k0]=buttap(N); %归一化巴特沃斯原型设计函数 p=p0*OmegaC;z=z0*OmegaC; %将零点极点乘以OmegaC， 得到非归一化零极点 k=k0*OmegaC^N; %将k（）乘以Omegac^N，得到非归一化k ba=k*real(poly(z)); %有零点计算分子系数向量 aa=real(poly(p)); %有极点计算分母系数向量 [bd,ad]=impinvar(ba,aa,Fs); %调用脉冲响应不变法函数 %检验频率响应 [H,w]=freqz(bd,ad,1000,whole); %计算数字系统频率响应 H=(H(1:1:501));w=(w(1:1:501)); %取其前一半值，并化为列向量 mag=abs(H); %求其幅特性 db=20*log10((mag+eps)/max(mag)); %化为分贝值 pha=angle(H); %求其相特性 grd=grpdelay(bd,ad,w); %求其群延迟 subplot(2,2,1);plot(w/pi,mag); title(幅度响应); xlabel();ylabel(|H|);axis([0,1,0,1.1]); set(gca,XTickMode,manual,XTick,[0,0.4,0.6,1]);grid on; %画刻度线 set(gca,XTickMode,manual,YTick,[0,Attn,Ripple,1]); subplot(2,2,3);plot(w/pi,db); title(幅度(单位：dB)); xlabel(频率(单位：pi));ylabel(分贝);axis([0,1,-100,50]); set(gca,XTickMode,manual,XTick,[0,0.4,0.6,1]); %画刻度线 set(gca,YTickMode,manual,YTick,[-100,-50,0,50]);grid on; set(gca,YTickLabelMode,manual,YTickLabels,[-100,-50,0,50]); subplot(2,2,2);plot(w/pi,pha/pi); title(相位响应); xlabel();ylabel(单位：pi);axis([0,1,-1,1]); set(gca,XTickMode,manual,XTick,[0,0.4,0.6,1]); %画刻度线 set(gca,YTickMode,manual,YTick,[-1,0,1]);grid on; subplot(2,2,4);plot(w/pi,grd); title(群延迟); xlabel