数字信号处理实验报告书(实验四).doc

实验报告 课程名称： 数字信号处理 实验题目： 实验四 模拟原型滤波器的设计 院 系： 物理与电子信息系 班 级： 姓 名： 学 号： 指导教师： 实验时间： 2012年 12月12日 实验四 模拟原型滤波器的设计 一、实验目的 学会运用MATLAB设计模拟低通滤波器原型的设计方法。 二、实验原理 由于IIR数字滤波器是在已知的归一化的低通模拟滤波器的基础上设计的，主要包括巴特沃斯低通滤波器、切比雪夫低通滤波器、椭圆低通滤波器，因此把这些低通滤波器称为滤波器原型。 三、实验内容及结果 设计一个模拟原型低通滤波器，通带截止频率为6khz，通带最大衰减为1db，阻带截止频率为15khz，阻带最小衰减为30db。要求：分别实现符合以上指标的巴特沃斯滤波器、切比雪夫I型滤波器、椭圆滤波器，绘制幅频特性、幅频特性分布图、相频特性、零极点分布图，并列写系统函数表示式，注意：是所求系统的系统函数，而不是归一化的系统函数。 巴特沃斯模拟原型低通滤波器 MATLAB源程序与运行结果图形： fp=6000;wp=2*pi*fp %输入滤波器的通带截止频率 fs=15000;ws=2*pi*fs %输入滤波器的阻带截止频率 Rp=1;As=30; %输入滤波器的通阻带衰减指标 %计算滤波器的阶数和3dB截止频率 [n,wc]=buttord(wp,ws,Rp,As,s) %计算n阶模拟低通原型 [z0,p0,k0]=buttap(n) %字母后加0表明这是原型滤波器的各指标，而不是所求的滤波器的 b0=k0*poly(z0) %求归一化滤波器分子系数b0 a0=poly(p0) %求归一化滤波器分母系数a0 [H,w]=freqs(b0,a0); %求归一化系统的频率特性 dbH=20*log10(abs(H)/max(abs(H))); %将归一化系统的幅频特性化为分贝值，注意此时的分贝值为负的，为了使分贝图和幅频响应图一致 subplot(2,2,1),plot(w*wc/(2*pi),abs(H)),grid %画所求滤波器的幅频响应图 axis([0,20000,0,1.1]);ylabel(幅度);xlabel(f(Hz)); subplot(2,2,2),plot(w*wc/(2*pi),angle(H)),grid %画所求滤波器的相频响应图 axis([0,20000,-4,4]);ylabel(相位);xlabel(f(Hz)); subplot(2,2,3),plot(w*wc/(2*pi),dbH),grid %画所求滤波器的幅频响应分贝图 axis([0,20000,-30,2]); ylabel(幅度(dB));xlabel(f(Hz)); subplot(2,2,4),plot(p0*wc,x); %画所求滤波器的极点图，用x号标出 axis square,axis equal,grid on %使x、y轴等比例显示 %以下程序可不写 wx0=[wp,ws]/wc; %设置频率向量 Hx=freqs(b0,a0,wx0); %计算通带阻带该两点的频率特性，即H(jwp)、H(jws) dbHx=-20*log10(abs(Hx)/max(abs(H))) %将H(jwp)、H(jws)化为分贝值，即通带最大衰减 、阻带最小衰减 wp = 3.7699e+004 ws = 9.4248e+004 n = 5 wc = 4.7241e+004 z0 = [] p0 = -0.3090 + 0.9511i -0.3090 - 0.9511i -0.8090 + 0.5878i -0.8090 - 0.5878i -1.0000 k0 = 1 b0 = 1 a0 = 1.0000 3.2361 5.2361 5.2361 3.2361 1.0000 dbHx = 0.4327 30.0000 系统函数： 切比雪夫I型模拟原型低通滤波器 MATLAB源程