数字图像处理实验5

班级学号姓名实验组别实验日期室温报告日期成绩报告内容：(目的和要求，原理，步骤，数据，计算，小结等)实验五IIR数字滤波器的设计设计内容：1、设计低通模拟滤波器，具体要求如下：3db带宽10KHZ，阻带起始频率12KHZ，阻带应到达的最小衰减-30db代码：clear;closeall;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('幅度indb');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；阻带范围：截止频率60HZ；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('幅度indb');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('幅度indb');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/Fs];ws=[2*fsl/Fs,2*fsu/Fs];Rp=3;As=30;[n,wn]=cheb1ord(wp,ws,Rp,As);[b,a]=butter(n,wn,'stop');[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('幅度indb'