实验5IIR数字滤波器的设计15生物医学工程曹美颖20151507023一、实验目的（1）加深对脉冲响应不变法、双线性变换法的了解；（2）掌握MATLAB进行滤波器设计的子函数。二、实验内容1.利用脉冲响应不变法，用巴特沃斯滤波器原型设计一个低通滤波器，满足：采样频率为10KHz，通带截止频率、通带最大衰减、阻带截止频率和阻带最小衰减分别为代码：wp=0.2*pi;ws=0.3*pi;rp=1;as=15;fs=10000;omgp=wp*fs;omgs=ws*fs;[n,omgc]=buttord(omgp,omgs,rp,as,'s');[z,p,k]=buttap(n);[bap,aap]=zp2tf(z,p,k);[ba,aa]=lp2lp(bap,aap,omgc);[bd,ad]=impinvar(ba,aa,fs);[H,wn]=freqz(bd,ad);plot(wn/pi,20*log10(abs(H)))gridon10.10.20.30.4-15-14-13-12-11-10-9-8-7-6-5-4-3-2-10结果：0.10.20.30.4-15-14-13-12-11-10-9-8-7-6-5-4-3-2-102.利用双线性变换法设计巴特沃斯高通数字滤波器，满足：通带截止频率为400Hz，阻带截止频率为200Hz，通带衰减小于3dB，阻带衰减大于15dB，采样频率为1000Hz。（选作）代码：wp=0.8*pi;ws=0.4*pi;rp=3;2as=15;fs=1000;omgp=wp*fs;omgs=ws*fs;[n,omgc]=buttord(omgp,omgs,rp,as,'s');[z,p,k]=buttap(n);[bap,aap]=zp2tf(z,p,k);[ba,aa]=lp2hp(bap,aap,omgc);[numd,dend]=bilinear(ba,aa,fs);[H,wn]=freqz(numd,dend);plot(wn/pi,20*log10(abs(H)))gridonxlabel('频率(Hz)')ylabel('幅度')title('数字滤波器幅频曲线')0.30.40.50.60.70.80.91-15-14-13-12-11-10-9-8-7-6-5-4-3-2-10(Hz)频率幅度数字滤波器幅频曲线3.设计巴特沃斯低通数字滤波器，满足：通带截止频率、通带最大衰减、阻带截止频率和阻带最小衰减分别为。3代码：wp=0.2;ws=0.3;rp=1;as=15;[n,wn]=buttord(wp,ws,rp,as);[b,a]=butter(n,wn);[H,wn]=freqz(b,a);plot(wn/pi,20*log10(abs(H)));gridonxlabel('频率(Hz)')ylabel('幅度')title('数字滤波器幅频曲线')结果：0.10.20.30.4-15-14-13-12-11-10-9-8-7-6-5-4-3-2-10(Hz)频率幅度数字滤波器幅频曲线4.设计巴特沃斯高通数字滤波器，满足：通带截止频率为400Hz，阻带截止频率为200Hz，通带衰减小于3dB，阻带衰减大于15dB，采样频率为1000Hz。代码：wp=0.8;ws=0.4;4rp=3;as=15;[n,wn]=buttord(wp,ws,rp,as);[b,a]=butter(n,wn,'high');[H,wn]=freqz(b,a);plot(wn/pi,20*log10(abs(H)))gridonxlabel('频率(Hz)')ylabel('幅度')title('数字滤波器幅频曲线')结果：0.30.40.50.60.70.80.91-15-14-13-12-11-10-9-8-7-6-5-...
