数字信号正式模板.doc

南昌航空大学电子信息工程学院

数字信号处理 课程实验报告

实验名称： 用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);
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。

五、 实验总结
……