电子科技大学 数字信号处理 上机实验四

1、电子科技大学生命科学与技术学院标 准 实 验 报 告（实验）课程名称 数字信号处理 2012-2013-第2学期电子科技大学教务处制表电 子 科 技 大 学实 验 报 告学生姓名：张华博 学 号：2011091010004 指导教师：李永杰实验地点：清水河校区 实验时间：2013年 月 日一、实验室名称：科二 504 机房 二、实验名称：FIR滤波器的特性及应用三、实验学时：2学时四、实验原理：五、实验目的：（详细填写） 1. 熟悉相位信息的特性2. 测试FIR滤波器的滤波特性六、实验内容：（详细填写） 1. 相位信息的作用2. 测试线性相位FIR滤波器的滤波特性七、实验器材（设备、元器件）：

2、 八、实验步骤： 九、实验数据及结果分析：（详细填写）（包括程序、图、结果等） 1、n=0:99;m=2/100:2/100:2;y=sin(0.1*pi*n)+2*cos(0.2*pi*n);XXX = fft(y);plot(m,abs(XXX);title('信号频谱图');figure(2); plot(m,unwrap(angle(XXX); title('信号相位谱');x1=cos(0.85*pi*n);x2=0.2*n;Pha1=unwrap(angle(XXX)+x1;Pha2=unwrap(angle(XXX)+x2;figure(3)plo

3、t(m,Pha1);title('加随机相位的相位谱');figure(4)plot(m,Pha2);title('加线性相位的相位谱'); cc1 = abs(XXX).*exp(i*Pha1);dd1 = ifft(cc1, 'symmetric'); cc2 = abs(XXX).*exp(i*Pha2); dd2 = ifft(cc2, 'symmetric'); figure(5)plot(n,y,'b',n,dd1,'r',n,dd2,': b');title('

4、fft反变换与原信号对比图');legend('原信号','加随机相位的信号','加线性相位的信号');2、h=-1 2 -3 6 -3 2 -1; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1;zplane(h,a)

5、;M = 100;H,w = freqz(h,a,M);figure(2); plot(w,abs(H);xlabel('omega');ylabel('Amplitude');figure(3); plot(w,unwrap(angle(H); xlabel('omega');ylabel('Phase'); h1=-1 2 -3 -3 2 -1; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa1=1; %µ&

6、#177;ΪFIRÂ˲¨Æ÷ʱ£¬È¡a=1;figure(4); zplane(h1,a1);M = 100;H1,w1 = freqz(h1,a1,M);figure(5); plot(w1,abs(H1);xlabel('omega');ylabel('Amplitude');figure(6); plot(w1,unwrap(angle(H1); xlabel('omega');ylabel(&

7、#39;Phase'); h2=-1 2 -3 0 3 -2 1; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa2=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1;figure(7); zplane(h2,a2);M = 100;H2,w2 = freqz(h2,a2,M);figure(8); plot(w2,

8、abs(H2);xlabel('omega');ylabel('Amplitude');figure(9); plot(w2,unwrap(angle(H2); xlabel('omega');ylabel('Phase'); h3=-1 2 -3 3 -2 1; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa3=1; %µ±ÎªFIRÂ˲¨

9、Æ÷ʱ£¬È¡a=1;figure(10); zplane(h3,a3);M = 100;H3,w3 = freqz(h3,a3,M);figure(11); plot(w3,abs(H3);xlabel('omega');ylabel('Amplitude');figure(12); plot(w3,unwrap(angle(H3); xlabel('omega');ylabel('Phase');h=-1 2 -3 6 -3 2 -1; a=1;

10、coef = sum(abs(h)/sum(a);h = h/coef; zplane(h,a);M = 50;H,w = freqz(h,a,M);figure(2); plot(w/pi,abs(H);figure(3); plot(w/pi,unwrap(angle(H); n=0:49;xs=sin(0.95*pi*n);xn=sin(0.5*pi*n);x=xs+xn;y = filter(h, a, x);h2,w2=freqz(y)figure(4)plot(n,xs,n,y,'r');legend('Ô­ÐÅ

11、86;Å','Â˲¨ºó');figure(5)plot(w1/pi,abs(h1),w2/pi,abs(h2),'r');title('ƵÆ×ͼ¶Ô±È');legend('Ô­ÐźÅ','Â˲¨ºó'); 十、实验结论：（详细填写） 1、给一信号相位谱加一线性相位后，做fft反变换所得图像为原图平移一定相位；若加一随机相位，做反变换后，不仅相位发生平移，而且幅度也发生变化。2、1、2、3、4型FIR滤波器均为线性相位滤波器，滤波之后图形有相位平移。