数字信号处理实验报告_王炜_20150518湖南工业大学..docx

电气与信息工程学院数字信号处理实验报告老师： 舒小华姓名： 王炜学号： 12401720207班级： 电子信息1202时间： 2015年5月实验一 快速傅里叶变换及其应用一、实验目的在理论学习的基础上，通过本实验，加强对FFT的理解，熟悉MATLAB中有关函数。应用FFT对典型信号进行频谱分析。了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT.应用相关FFT实现序列的线性卷积和相关。二、实验内容实验用到的信号序列：高斯序列衰减正弦序列三角波序列反三角波序列上机实验内容：观察高斯序列的时域和幅频特性，固定信号中参数p=8，改变q的值，使q分别等于2、4、8，观察他们的时域和幅频特性，了解当q去不同值时，对信号序列的时域和频域特性的影响；固定q=8,改变p,使p分别为等于8、13、14，观察参数p变化对信号序列的时域以及频域特性的影响，注意p等于多少时，会发生明显的泄露现象，混叠是否也随之出现？记录实验中观察到的现象，绘制相应的时域序列和幅频特性曲线。代码如下：clear all;%p=8,q=2%[Xa1,Fa1]=gauss(8,2);k=0:15;subplot(5,2,1);plot(k,Xa1);Xlabel(n);Ylabel(时域特性);text(10,0.5,p=8,q=2);subplot(5,2,2);plot(k,Fa1);Xlabel(n);Ylabel(幅频特性);text(8,3,p=8,q=2);%p=8,q=4%[Xa1,Fa1]=gauss(8,4);k=0:15;subplot(5,2,3);plot(k,Xa1);Xlabel(n);Ylabel(时域特性);text(10,0.5,p=8,q=4);subplot(5,2,4);plot(k,Fa1);Xlabel(n);Ylabel(幅频特性);text(8,3,p=8,q=4);%p=8,q=8%[Xa1,Fa1]=gauss(8,8);k=0:15;subplot(5,2,5);plot(k,Xa1);Xlabel(n);Ylabel(时域特性);text(10,0.5,p=8,q=8);subplot(5,2,6);plot(k,Fa1);Xlabel(n);Ylabel(幅频特性);text(8,3,p=8,q=8);%p=8,q=8%[Xa1,Fa1]=gauss(13,8);k=0:15;subplot(5,2,7);plot(k,Xa1);Xlabel(n);Ylabel(时域特性);text(10,0.5,p=13,q=8);subplot(5,2,8);plot(k,Fa1);Xlabel(n);Ylabel(幅频特性);text(8,3,p=13,q=8);%p=8,q=8%[Xa1,Fa1]=gauss(14,8);k=0:15;subplot(5,2,9);plot(k,Xa1);Xlabel(n);Ylabel(时域特性);text(10,0.5,p=14,q=8);subplot(5,2,10);plot(k,Fa1);Xlabel(n);Ylabel(幅频特性);text(8,3,p=14,q=8);观察衰减正弦序列的时域和幅频特性，a=0.1，f=0.0625，检查峰值出现的位置是否正确，注意频谱的形状，绘出频谱特性曲线，改变发，使分别等于0.4375和0.5625，观察这两种情况下，频谱的形状和频谱出现的位置，有无混叠、泄露现象？说明产生现象的原因。代码如下：clear all;a=0.1;f=0.0625;i=0:160x=exp(-a*(i)).*sin(2*pi*f*(i));subplot(3,2,1);plot(i,x);xlabel(n);title(a=0.1,f=0.0625 时域特性)subplot(3,2,2);G=fft(x);plot(i,abs(G));xlabel(k);title(a=0.1,f=0.0625 幅频特性)f=0.4375;x=exp(-a*(i)).*sin(2*pi*f*(i));subplot(3,2,3);plot(i,x);xlabel(n);title(a=0.1,f=0.4375 时域特性)subplot(3,2,4);G=fft(x);plot(i,abs(G));xlabel(k);title(a=0.1,f=0.4375 幅频特性)%f=0.5625;f=0.5625;x=exp(-a*(i)).*sin(2*pi*f*(i));subplot(3,2,5);plot(i,x);xlabel(n);title(a=0.1,f=0.5625 时域特性)subplot(3,2,6);G