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

电子科技大学生命科学与技术学院 标 准 实 验 报 告 （实验）课程名称 数字信号处理 2012-2013-第2学期 电子科技大学教务处制表 电 子 科 技 大 学 实 验 报 告 学生姓名：张华博 学 号：2011091010004 指导教师：李永杰 实验地点：清水河校区 实验时间：2013年 月 日 一、实验室名称：科二 504 机房 二、实验名称：FIR滤波器的特性及应用 三、实验学时：2学时 四、实验原理： 五、实验目的：（详细填写） 1. 熟悉相位信息的特性 2. 测试FIR滤波器的滤波特性 六、实验内容：（详细填写） 1. 相位信息的作用 2. 测试线性相位FIR滤波器的滤波特性 七、实验器材（设备、元器件）： 八、实验步骤： 九、实验数据及结果分析：（详细填写）（包括程序、图、结果等） 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) plot(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(fft反变换与原信号对比图); legend(原信号,加随机相位的信号,加线性相位的信号); 2、 h=[-1 2 -3 6 -3 2 -1]; %×??o11?ìh(n)DòáD a=1; %μ±?aFIR??2¨?÷ê±￡?è?a=1; zplane(h,a); 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]; %×??o11?ìh(n)DòáD a1=1; %μ±?aFIR??2¨?÷ê±￡?è?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(Phase); h2=[-1 2 -3 0 3 -2 1]; %×??o11?ìh(n)DòáD a2=1; %μ±?aFIR??2¨?÷ê±￡?è?a=1; figure(7); zplane(h2,a2); M = 100; [H2,w2] = freqz(h2,a2,M); figure(8); plot(w2,abs(H2)); xlabel(\omega); ylabel(Amplitude); figure(9); plot(w2,unwrap(angle(H2))); xlabel(\omega);ylabel(Phase); h3=[-1 2 -3 3 -2 1]; %×??o11?ìh(n)DòáD a3=1; %μ±?aFIR??2¨?÷ê±￡?è?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))); x