数字图像处理实验5.docVIP

班级 学号 姓名 实验组别 实验日期 室温 报告日期 成绩 报告内容：(目的和要求，原理，步骤，数据，计算，小结等) 实验五 IIR数字滤波器的设计 设计内容： 1、设计低通模拟滤波器，具体要求如下： 3db带宽10KHZ，阻带起始频率12KHZ，阻带应达到的最小衰减-30db 代码： clear;close all;clc; % 开始准备 fp=10000;fs=12000;Fs=20000;Rp=1;As=30;T=1/Fs; % 设计指标 wp=fp/Fs*2; ws=fs/Fs*2; % 求归一化频率 [n,wn]=buttord(wp,ws,Rp,As,s) %模拟原型滤波器的阶数 [z,p,k]=buttap(n); %求出零极点型系统函数 [bap,aap]=zp2tf(z,p,k) %将其转化为b/a型 [b,a]=lp2lp(bap,aap,wn); %模拟滤波器（低通) [db,mag,pha,w]=freqs_m(b,a,0.5*pi); %求出滤波器的幅频 [ha,x,t]=impulse(b,a); %求出滤波器的相频 subplot(2,2,1),plot(w/pi,mag);title(幅度); subplot(2,2,2),plot(w/pi,db);title(幅度 in db); subplot(2,2,3),plot(w/pi,pha/pi);title(相位); subplot(2,2,4),plot(t,ha),title(冲激响应); figure,freqs(b,a); （b）结果： 2、设计数字滤波器 （1）带通滤波器：通带范围：100HZ---200HZ ； 阻带范围：截止频率60 HZ；240HZ ，通带允许最大衰减3dB, 阻带应达到的最小衰减-30dB (a)代码： fpl=100;fpu=200;fsl=60;fsu=240;Fs=1000; wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];Rp=3;As=30; [n,wn]=cheb1ord(wp,ws,Rp,As); [b,a]=butter(n,wn,bandpass); [db,mag,pha,grd,w]=freqz_m(b,a); subplot(2,2,1),plot(w/pi,mag);title(幅度); subplot(2,2,2),plot(w/pi,db);title(幅度 in db); subplot(2,2,3),plot(w/pi,pha/pi);title(相位); subplot(2,2,4),plot(w/pi,grd),title(群延时); figure,freqz(b,a); （b）结果： （2）低通滤波器wp=pi/3, ws=pi/2通带允许最大衰减3dB, 阻带应达到的最小衰减-30dB。 （a）代码： fpl=100;fpu=200;fsl=60;fsu=240;Fs=1000; wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];Rp=3;As=30; [n,wn]=cheb1ord(wp,ws,Rp,As); [b,a]=butter(n,wn,bandpass); [db,mag,pha,grd,w]=freqz_m(b,a); subplot(2,2,1),plot(w/pi,mag);title(幅度); subplot(2,2,2),plot(w/pi,db);title(幅度 in db); subplot(2,2,3),plot(w/pi,pha/pi);title(相位); subplot(2,2,4),plot(w/pi,grd),title(群延时); figure,freqz(b,a); （b）结果： （3）带阻滤波器：阻带范围：100HZ---200HZ 通带截止频率60HZ和240HZ。通带允许最大衰减3dB, 阻带应达到的最小 （a）代码： fpl=100;fpu=200;fsl=60;fsu=240;Fs=1000; wp=[2*fpl/Fs,2*fpu/F