数字信号处理实验 实验八

1、实 验 报 告课程:数字信号处理实 班级:姓 名:学号: 实验源程序代码及截图: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:50

2、1'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' subp

3、lot(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('fre

4、quency 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'204060 理 想 脉 冲

5、 响 应nh d (n 020 4060Hamming 窗nw (n 204060 实 际 脉 冲 响 应nh (n 0.20.31-50相 对 标 尺 的 幅 度 响 应 frequency in pi unitsD e c i b e l sws1=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(w

6、c1,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

7、,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

8、9;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,'YTic

9、kLabelMode','manual','YTickLabels','50''0' >> sy8_2 tr_width =0.4712 M =75 204060理 想 脉 冲 响 应nh d (n blackman 窗nw (n 204060实 际 脉 冲 响 应nh (n 0.20.31-50相 对 标 尺 的 幅 度 响 应frequency in pi unitsD e c i b e l s窗函数法设计 FIR 滤波器的主要特点分析:FIR 滤波器目前常用的设计方法有窗函数法和频率采样法, 窗函数法是从 时域进行 设计,而频率采样法是从频域进行设计。窗函数法由于简单、物理意义清晰,因而得到了 较为广泛的应用。 窗函数法设计的基本思想是:首先根据技术指标要求, 选取合适的阶数 N 和窗函数的类型 w(n,使其幅频特性逼近理想滤波器幅频特性。其次,因为理想滤波器的 hd(n是无限长的,所以需要对 hd(n 进行截断,数学上称这种方法为窗函数法。 简而言 之,用窗函数法设计 FIR 滤波器是在时域进行的,先用傅里叶变换求出理想滤波器单位抽 样相应 hd(n