iir数字滤波器设计matlab代码

《iir数字滤波器设计matlab代码》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、大连理工大学实验报告学院(系)：电子信息与电气T程学部专业:电子信息T程班级:姓名：学号：组：实验时间：实验室：大黑楼C221实验台:指导教师签字：成绩:实验名称实验三IIR数字滤波器设计一、实验目的和要求(1)掌握双线性变换法设计IIR数字滤波器。(2)掌握脉冲响应不变法设计IIR数字滤波器。二、实验原理和内容1.巴特沃思低通滤波器巴特沃思低通滤波器幅度平方函数定义为IH,側22NQ式中N为正整数，代表滤波器的阶次，.称为截止频率。当Q=.时，有?1即=3dB=201g所以又称.为巴特沃思低通滤波器的3分贝带宽。

2、巴特沃思低通滤波器的幅度特性如图所示。图1巴特沃思低通滤波器的幅度特性)1211/21.切贝雪夫低通滤波器切贝雪夫I型滤波器的幅度平方函数为其中，s是表示通带波纹大小的一个参数，s越大，波纹也越大。2■为0对0(.的归一化频率，Q为截止频率，C/V(x)是N阶切贝雪夫多项式。切贝雪夫I型低通滤波器的幅度特性如图所示图2切贝雪夫I型低通滤波器的幅度特性三、主要仪器设备Matlab四、实验内容及结果1.用双线性变换法设计一个巴特沃斯IIR低通数字滤波器。设计指标参数为••在通带内频率低于0.4*pi时，最大衰减小于2dB

3、;在阻带内［0.6*pi,pi］频率区间上，最小衰减大于15dB。要求：•(1)用buttord，buttap，zp2tf，lp21p，bilinear等指令实现•(2)用buttord，butter指令实现•并分别绘制两种设计方法的频响曲线(1)双线形变换法设计巴特沃斯IIR低通数字滤波器，代码如下:clcclearwp=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]=butt

4、ord(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，w1]=freqz(Bz，Az，w);H=20*logl0(abs(H));plot(w/pi,H);gridon;titlef双线形变换法设计巴特沃斯IIR低通数字滤波器’)；xlabelf频率/pi’);ylabclC幅度/dB1);（2）直接法设计巴特沃斯IIR低通数字滤波器，

5、代码如下：:clcclearwp=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，w1]=freqz（Bz，Az，w）;H=20*logl0（abs（H））;plot（w/pi，H）;gridon;title（’直接法设计巴特沃斯IIR低通数字滤波器xlabelf频率/pi'）;ylabelC幅度/dB）;S接法设讣巴特沃斯IIR低逍数卞滤波器5000.10.20.30.40.50.60.7

6、0.80.91频率/pi(1)设计一个巴特沃兹低通滤波器，给定技术指标是A=100//Z九=3⑽Hzap=3dBax=20dB，抽样频率力6=1000Hz。代码如下：clcclearfp=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低通数字滤波器’);gpep3c6ez(S

7、)