第7章有限脉冲响应数字滤波器设计汇编.docVIP

成 绩： 《数字信号处理》 作业与上机实验 （第七章） 班 级： 学 号： 姓 名： 任课老师： 完成时间： 2014.10.28 信息与通信工程学院 2013—2014学年第 2 学期 第7章 有限脉冲响应数字滤波器设计 1、教材p238： 19. (1)MATLAB代码如下： Fs=80000; T=1/Fs; M=1024; fp=15000;fs=20000;rs=40; wp=2*pi*fp*T;ws=2*pi*fs*T; Rp=-20*log10(1-0.02);As=40; [N1,wp1]=ellipord(wp/pi,ws/pi,Rp,As); [B,A]=ellip(N1,Rp,As,wp1); [Hk,wk1]=freqz(B,A,1000); mag=abs(Hk);pah=angle(Hk); Bt=ws-wp; alph=0.5842*(rs-21)^0.4+0.07886*(rs-21); N=ceil((rs-8)/2.285/Bt); wc=(wp+ws)/2/pi; hn=fir1(N,wc,kaiser(N+1,alph)); Hk=fft(hn,M); k=0:M/2-1; wk=(2*pi/M)*k; figure(1); plot(wk/pi,20*log10(abs(Hk(k+1))),:,linewidth,2.5); hold on plot(wk1/pi,20*log10(mag),linewidth,2); hold off legend(FIR滤波器,IIR滤波器); title(损耗函数); figure(2) plot(wk/pi,angle(Hk(k+1))/pi,:,linewidth,2.5); hold on plot(wk1/pi,pah/pi,linewidth,2); hold off legend(FIR滤波器,IIR滤波器); title(相频特性曲线); 两种数字滤波器的损耗函数和相频特性的比较分别如图1.1、1.2所示： 图1.1损耗函数比较图 图1.2 相频特性比较图 20. (1)MATLAB代码如下： wc=pi/4; N=21; M=1024; n=0:N-1; h_boxcar=fir1(N-1,wc/pi,boxcar(N)); h_hanning=fir1(N-1,wc/pi,hanning(N)); h_hamming=fir1(N-1,wc/pi,hamming(N)); h_blackman=fir1(N-1,wc/pi,blackman(N)); plot(n,h_boxcar); hold on plot(n,h_hanning,:,linewidth,2); plot(n,h_hamming,+,linewidth,2); plot(n,h_blackman,o); hold off legend(矩形窗,汉宁窗,哈明窗,布莱克曼窗); title(单位冲激响应); Hk=fft(h_boxcar,M);k=0:M/2-1; wk=(2*pi/M)*k; figure(); plot(wk/pi,20*log10(abs(Hk(k+1))),linewidth,2); Hk=fft(h_hanning,M) hold on plot(wk/pi,20*log10(abs(Hk(k+1))),:,linewidth,3); Hk=fft(h_hamming,M) plot(wk/pi,20*log10(abs(Hk(k+1))),o); Hk=fft(h_blackman,M) plot(wk/pi,20*log10(abs(Hk(k+1))),*); hold off legend(矩形窗,汉宁窗,哈明窗,布莱克曼窗); title(损耗函数); 四种窗函数设计的单位脉冲响应和损耗函数的比较如图2.1 2.2所示： 图2.1 单位脉冲响应比较图 图2.2 损耗函数比较图 21. MATLAB代码如下： wc=pi/4; N=21; M=1024; n=0:N-1; h_boxcar=fir1(N-1,