上海交通大学电子信息与电气工程学院dsp离散时间信号实验报告二.docVIP

设计数字滤波器实验报告 刘运鹏 5080309226 实验目的： １、掌握设计IIR和FIR滤波器的原理和具体设计方法。 ２、学会在计算机上用IIR和FIR滤波器进行数字滤波。 ３、IIR和FIR滤波器特性分析。 实验内容： １、掌握IIR和FIR滤波器原理。 ２、通过实验设计IIR和FIR滤波器。 ３、分析IIR和FIR滤波器。 ４、对所得结果加以讨论。 题目一：50.a 50.b 设计一个用在结构 中的低通数字滤波器，满足下列要求： 　采样速率为8000个样本/秒； 　通带边缘频率为1500Hz，波动为3dB； 　阻带边缘频率为2000Hz，衰减为40dB； 　通带等波动但是阻带是单调的； 　采用脉冲响应不变法。 a.在脉冲响应不变法中，选择T=1，求出并联形式的系统函数，画出对数幅度响应和脉冲响应图。 b.在脉冲响应不变法中，选择T=1/8000，求出并联形式的系统函数，画出对数幅度响应和脉冲响应图，比较上述两种设计，讨论T对设计的影响。 参数修改：1500修改为1000 程序代码： fp=1000; fs=2000; wp=2*pi*fp; ws=2*pi*fs; rp=3; rs=40; [n,wn]=cheb1ord(wp,ws,rp,rs,s); [z,p,k]=cheb1ap(n,rp); [b1,a1]=zp2tf(z,p,k); [b2,a2]=lp2lp(b1,a1,wn); [bz,az]=impinvar(b2,a2,8000); [hz,wz]=freqz(bz,az); mag=abs(hz); lmag=20*log(mag); f=(0:511)*pi/512; plot(f,lmag); grid %幅度响应 ylabel(幅值); xlabel(频率); h=impz(bz,az); plot(h); grid %脉冲响应 ylabel(幅值); xlabel(时间); fp=1000; fs=2000; wp=2*pi*fp/8000; ws=2*pi*fs/8000; rp=3; rs=40; [n,wn]=cheb1ord(wp,ws,rp,rs,s); [z,p,k]=cheb1ap(n,rp); [b1,a1]=zp2tf(z,p,k); [b2,a2]=lp2lp(b1,a1,wn); [bz,az]=impinvar(b2,a2,8000/8000); [hz,wz]=freqz(bz,az); mag=abs(hz); lmag=20*log(mag); f=(0:511)*pi/512; plot(f,lmag); grid %幅度响应 ylabel(幅值); xlabel(频率); h=impz(bz,az); plot(h); grid %脉冲响应 ylabel(幅值); xlabel(时间); 实验结果和分析： 前后两次图像相同证明T不影响滤波器设计 对数幅度图a 脉冲响应曲线a 对数幅度图b 脉冲响应曲线b 题目二：56 用cheby2函数设计一个带通数字滤波器，指标为： 低阻带边缘频率：0.3pi 高阻带边缘频率：0.6pi As=50dB 低通带边缘频率：0.4pi Rp=0.5dB 高通带边缘频率：0.5pi 画出所设计滤波器的脉冲响应和对数幅度响应（单位：dB）。 参数修改：Rs修改为40dB。 程序代码： wp=[0.4*pi,0.5*pi]; ws=[0.3*pi,0.6*pi]; rp=0.5; rs=40; [n,wn]=cheb2ord(wp/pi,ws/pi,rp,rs); [b,a]=cheby2(n,rs,wn); [h,w]= freqz(b,a); mag=abs(h); lmag=20*log(mag); f=(0:511)*pi/512; plot(f,lmag); grid %幅度响应 ylabel(幅值); xlabel(频率); h=impz(b,a); plot(h); grid %脉冲响应 ylabel(幅值); xlabel(时间); 实验结果和分析： 对数幅度图 脉冲响应曲线 第三题：59 用凯泽窗设计技术设计一个高通滤波器，技术指标为： 阻带边缘：0.4pi As=60dB 通带边缘：0.6pi Rp=0.5dB 画出设计的滤波器的脉冲响应和幅度响应（dB值）。 参数修改：As改为70dB。 程序代码： wc=0.4*pi; %截止频率