1、南昌航空大学电子信息工程学院数字信号处理 课程实验报告实验名称: 用 FFT 作谱分析实验时间: *年*月*日指导教师: 班 级 : *学 号 : *姓 名 : *成 绩 : *一、 实验目的(1)(2)(3)二、 实验内容(1)(2)(3)(4) 三、 实验原理四、 实验步骤、相应的实验结果、相应的结果分析(1) 编写一个 M 文件,以文件名 sy3_1.m 保存在当前路径,具体代码如下:N=10;n1=0:N-1;x1=ones(1,6),zeros(1,N-6);Xk1=dft(x1,N);figure(1);subplot(2,1,1);stem(n1,x1);subplot(2,1,
2、2);stem(n1,abs(Xk1);N=20;n2=0:N-1;x2=ones(1,6),zeros(1,14);Xk2=dft(x2,N);figure(2);subplot(2,1,1);stem(n2,x2);subplot(2,1,2);stem(n2,abs(Xk2);在命令窗口中键入 sy3_1 并回车运行,可以得到如图 3-1 所示的图形。图 3-1由图 3-1 可知,对一个信号进行时域补零,再进行 FFT 运算得到的离散谱,包络不变,但是间隔更小,相当于进行了频域的插值。(3) 编写一个 M 文件,以文件名 sy3_3.m 保存在当前路径,具体代码如下:。在命令窗口中键入sy3_3并回车运行,可以得到如图 3-3所示的图形。图 3-3由图3-3 可知,第一个图中N=128,两个频率成分无法区分;第二个图中N=128,然后再补零为512 ,由实验3-2 知,补零无法提高频谱分辨率,所以两个频率成分还是无法区分;第三个图,是截取了序列的512个数据点,这是得到的频谱图可以清楚看到有两个频率成分。因为F=2.05-2=0.05 ;fc=2.05, 所以Nmin=2fc/F=820。五、 实验总结