IIR数字滤波器设计matlab代码.docVIP

大连理工大学实验报告 学院(系)：电子信息与电气T程学部专业:电子信息T程班级: 姓 名： 学号： 组： 实验时间： 实验室：大黑楼C221 实验台: 指导教师签字： 成绩: 实验名称 实验三IIR数字滤波器设计 一、实验目的和要求 掌握双线性变换法设计IIR数字滤波器。 掌握脉冲响应不变法设计IIR数字滤波器。 二、实验原理和内容 巴特沃思低通滤波器 巴特沃思低通滤波器幅度平方函数定义为 IH,側2\2N IH,側2 \2N Q 式中N为正整数，代表滤波器的阶次，.称为截止频率。当Q = .时，有 ? 1 即 = 3dB=201g = 3dB 所以又称.为巴特沃思低通滤波器的3分贝带宽。 巴特沃思低通滤波器的幅度特性如图所示。 图1巴特沃思低通滤波器的幅度特性 )12 11/2 切贝雪夫低通滤波器 切贝雪夫I型滤波器的幅度平方函数为 其中，s是表示通带波纹大小的一个参数，s越大，波纹也越大。2■为0对0(.的归一化频率，Q 为截止频率，C/V(x)是N阶切贝雪夫多项式。 切贝雪夫I型低通滤波器的幅度特性如图所示 图2切贝雪夫I型低通滤波器的幅度特性 三、 主要仪器设备 Matlab 四、 实验内容及结果 1.用双线性变换法设计一个巴特沃斯IIR低通数字滤波器。设计指标参 数为??在通带内频率低于0.4*pi时，最大衰减小于2dB ;在阻带内［0.6*pi, pi］频率区间上，最小衰减大于15dB。要求： ?(1)用 buttord，buttap，zp2tf，lp21p，bilinear 等指令实现 ?(2)用 buttord，butter 指令实现 ?并分别绘制两种设计方法的频响曲线 (1)双线形变换法设计巴特沃斯IIR低通数字滤波器，代码如下: clc clear wp = 0.4*pi; ws = 0.6*pi; rp = 2; rs = 15; Ts =0.0001; Fs=l/Ts; %采样频率 wp = 2/Ts*tan(wp/2); ws = 2/Ts*tan(ws/2); [N，wc]=buttord(wp，ws，rp，rs，’s’)； [z,p,k]=buttap(N); [b,a]=zp2tf(z,p,k); [B，A]=lp21p(b，a，wc); [Bz，Az]=bilinear(B，A，Fs); w=0:0.001 :pi; [H，w 1 ]=freqz(Bz，Az，w); H=20*logl0(abs(H)); plot(w/pi,H); grid on; titlef双线形变换法设计巴特沃斯IIR低通数字滤波器’)； xlabelf 频率/pi’); ylabclC 幅度/dB1); （2）直接法设计巴特沃斯IIR低通数字滤波器，代码如下：: clc clear wp=O ? 4; ws=O ? 6; rp=2; rs= 15; [N，wc]=buttord（wp，\vs，rp，rs）; [Bz，Az]=butter（N，wc）; w=0:0.001:pi; [H，w 1 ]=freqz（Bz，Az，w）; H=20*logl0（abs（H））; plot（w/pi，H）; grid on; title（’直接法设计巴特沃斯IIR低通数字滤波器 xlabelf 频率/pi）; ylabelC 幅度/dB）; S接法设讣巴特沃斯IIR低逍数卞滤波器 S接法设讣巴特沃斯IIR低逍数卞滤波器 50 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 频率/pi 设计一个巴特沃兹低通滤波器，给定技术指标是 A = 100//Z 九=3⑽Hzap=3dBax= 20dB，抽样频率力 6 = 1000Hz。 代码如下： clc clear fp=100; fs=300; Fs=1000; af = 3; as = 20; wp = 2氺fp/Fs; ws = 2氺fs/Fs; [n, Wn] =buttord (wp, ws, af, as); [b，a]=buttcr (n, Wn); figure(l) freqz (b, a, 512, l); titleC巴特沃斯IIR低通数字滤波器’); gp ep3c6ez gp ep3c6ez (SD① (SD① JO)① p)Ds05qd 五、实验体会 本次实验中用了双线性变换法、冲击响应不变法及直接设计法设计了 IIR滤波器。通过本次实 骑，我更直观的看到了滤波器的特性，对理论有了更深刻的认识。Matlab是一个十分强大的工具， 有很多很好用的函数，比起人工计算设计滤波器来简中.了很多很多，很直观也很实用。