第六章作业1).docVIP

《数字信号处理》 作业与上机实验 （第六章） 班 级： N电信12-1F 学 号： 24122201242 姓 名： 曹福利 任课老师： 李宏民 完成时间： 11.16 信息与通信工程学院 2013—2014学年第 2 学期 第6章 无限脉冲响应数字滤波器设计 1、教材p195： 13,14,15,16,17,18,19 2、某信号为：，其中设计最低阶数IIR数字滤波器，按下图所示对进行数字滤波处理，实现： 1）将频率分量以高于50dB的衰减抑制，同时以低于2dB的衰减通过和频率分量； 2）将和频率分量以高于50dB的衰减抑制，同时以低于2dB的衰减通过频率分量； 要求：按数字滤波器直接型与级联型结构图编写滤波程序，求得；IIR滤波器采用双线性（高通）与冲击响应不变法（低通）设计；画出所设计的滤波器频率特性图、信号时域图；给出滤波器设计的MATLAB代码与滤波器实现的代码；选择合适的信号采样周期T与两种滤波器设计方法中对模拟滤波器冲击响应的采样周期T’, 注意区分二者的异同。 1、 13 . 解：t=1; fs=4000; wpu=0.45*pi; wpl=0.25*pi; wsu=0.55*pi; wsl=0.15*pi; wpz=[0.25,0.45]; wsz=[0.15,0.55]; wp=2/t*tan(wpz/2); ws=2/t*tan(wsz/2); rp=3; as=40; [n,wc]=buttord(wp,ws,rp,as,s); [b,a]=butter(n,wc,s); [bz,az]=impinvar(b,a,fs); [nd,wdc]=buttord(wpz,wsz,rp,as); [bd,adz]=butter(nd,wdc); hk=freqz(bd,adz); plot(abs(hk)) 图一.Hz t=1; fs=4000; wpu=0.45*pi; wpl=0.25*pi; wsu=0.55*pi; wsl=0.15*pi; wpz=[0.25,0.45]; wsz=[0.15,0.55]; wp=2/t*tan(wpz/2); ws=2/t*tan(wsz/2); rp=3; as=40; [n,wc]=buttord(wp,ws,rp,as,s); [b,a]=butter(n,wc,s); [bz,az]=impinvar(b,a,fs); [nd,wdc]=buttord(wpz,wsz,rp,as); [bd,adz]=butter(nd,wdc); hk=freqz(bd,adz); figure(1) plot(angle(hk)) figure(2) [hp,w]=freqz(bd,adz,4000); plot(w/pi,20*log10(hp)); 图二.相频特性曲线 图三.损耗函数 14. 解：T=1/80000; wp=4000*2*pi*T; ws=20000*2*pi*T; rp=0.5; rs=45; [n,wc]=buttord(wp,ws,rp,as,s); [b,a]=butter(n,wc,s); [bz,az]=impinvar(b,a); figure(1) freqz(b,a); [bz,az]=impinvar(b,a); figure(2) freqz(bz,az); 图四 15. 解：T=1/80000; wp=4000*2*pi*T; ws=20000*2*pi*T; Rp=0.5; As=45; [N,wp1]=cheblord(wp,ws,Rp,As,s); [b,a]=chebyl(N,wp1,Rp,s); [bz,az]=impinvar(b,a,fs); [db,mag,pha,grd,w]=freqz_m(bz,az); plot(w/pi,db); 图五 16． 解：T=1/2500000; wp=325000*2*pi*T; ws=225000*2*pi*T; rp=1; rs=40; [N,wc]=ellipord(wp,ws,rp,rs); [b,a]=ellip(N,rp,rs,wc); figure(1) freqz(b,a); [bz,az]=impinvar(b,a); figure(2) freqz(bz,az); 图六 17． 解： fsl=560000; fsu=7