数字信号处理实验 实验八.doc

实 验 报 告 课程:数字信号处理实 班级:姓 名:学号: 实验源程序代码及截图: function hd=ideal_lp(wc,M; alpha=(M-1/2; n=[0:1:(M-1]; m=n-alpha+eps; hd=sin(wc*m./(pi*m; wp=0.2*pi;ws=0.3*pi; tr_width=ws-wp; M=ceil(6.6*pi/tr_width+1; n=[0:1:M-1]; wc=(ws+wp/2; hd=ideal_lp(wc,M; w_ham=(hamming(M; h=hd.*w_ham; [H,w]=freqz(h,1,1000,whole; H=(H(1:1:501;w=(w(1:1:501; mag=abs(H; db=20*log10((mag+eps/max(mag; pha=angle(H; grd=grpdelay(h,1,w; delta_w=2*pi/1000; Rp=-(min(db(1:1:wp/delta_w+1; As=-round(max(db(ws/delta_w+1:1:501; subplot(1,1,1; subplot(2,2,1;stem(n,hd;title(理想脉冲响应 ; axis([0 M-1 -0.1 0.3];xlabel(n;ylabel(hd(n; subplot(2,2,2;stem(n,w_ham;title(Hamming窗 ; axis([0 M-1 0 1.1];xlabel(n;ylabel(w(n; subplot(2,2,3;stem(n,h;title(实际脉冲响应 ; axis([0 M-1 -0.1 0.3];xlabel(n;ylabel(h(n; subplot(2,2,4;plot(w/pi,db;title(相对标尺的幅度响应 ;grid; axis([0 1 -100 10];xlabel(frequency in pi units;ylabel(Decibels; set(gca,XTickMode,manual,XTick,[0,0.2,0.3,1]; set(gca,YTickMode,manual,YTick,[-50,0]; set(gca,YTickLabelMode,manual,YTickLabels,[50;0]; 20 40 60 理 想 脉 冲 响 应 n h d (n 020 4060 Hamming 窗 n w (n 20 4060 实 际 脉 冲 响 应 n h (n 0.20.3 1 -50 相 对 标 尺 的 幅 度 响 应 frequency in pi units D e c i b e l s ws1=0.2*pi;wp1=0.35*pi; wp2=0.65*pi;ws2=0.8*pi; As=60; tr_width=min((wp1-ws1,(ws2-wp2 M=ceil(11*pi/tr_width+1 n=[0:1:M-1]; wc1=(ws1+wp1/2;wc2=(wp2+ws2/2; hd=ideal_lp(wc2,M-ideal_lp(wc1,M; w_bla=(blackman(M; h=hd.*w_bla; [H,w]=freqz(h,1,1000,whole; H=(H(1:1:501;w=(w(1:1:501; mag=abs(H; db=20*log10((mag+eps/max(mag; pha=angle(H; grd=grpdelay(h,1,w; delta_w=2*pi/1000; Rp=-(min(db(1:1:wp/delta_w+1; As=-round(max(db(ws/delta_w+1:1:501; subplot(1,1,1; subplot(2,2,1;stem(n,hd;title(理想脉冲响应 ; axis([0 M-1 -0.1 0.3];xlabel(n;ylabel(hd(n; subplot(2,2,2;stem(n,w_bla;title(blackman窗 ; axis([0 M-1 0 1.1];xlabel(n;ylabel(w(n; subplot(2,2,3;stem(n,h;title(实际脉冲响应 ; axis([0 M-1 -0.1 0.3];xlabel(n;ylabel(h(n; subplot(2,2,4;plot(w/pi,db;title(相对标尺的幅度响应 ;grid; axis([0 1 -100 10];xlabel(frequency in pi units;ylabel(Decibels; set(gca,XTickMode,manual