信号与系统仿真实验报告(5)

Columns 11 through 15

1.1009 - 0.0673i 1.0959 - 0.1034i 1.0820 - 0.1556i 1.0458 - 0.2476i 0.8900 - 0.4403i

Columns 16 through 20

0.0052 - 0.0000i 0.8900 + 0.4403i 1.0458 + 0.2476i 1.0820 + 0.1556i 1.0959 + 0.1034i

Columns 21 through 25

1.1009 + 0.0673i 1.1032 + 0.0375i 1.1052 + 0.0118i 0.0375i

Columns 26 through 30

1.1009 - 0.0673i 1.0959 - 0.1034i 1.0820 - 0.1556i 0.4403i

Column 31

0.0052 - 0.0000i

源程序：

stem(abs(X),'.') %绘制幅频特性曲线 stem(angle(X),'.') %绘制相频特性曲线

运行结果： 幅频特性：

- 21 -

1.1052 - 0.0118i 1.0458 - 0.2476i 1.1032 - 0.8900 -

相频特性：

- 22 -

五、结果分析及实验体会 结果分析：

问题（1）中，从y的图形可以看出，x和h卷积的结果序列值随n的增大呈现先增大后减小的趋势。问题（2）中，从响应y的图形可以看出当自变量n取奇数时，y的值为零，当n取偶数时，y的值从某一负数值逐渐趋于零。通过对单位抽样响应的幅频特性进行分析可以得出，当频率在图形的一个周期以内的相当宽的一段范围内，幅度值较高，只有在图形的一个周期的起始和终止部分，幅度值较低。由相频特性曲线可以看出，在图形的一个周期开始的一段频率内，相位会出现一个正的峰值，在图形的一个周期快结束的一段频率内，相位会出现一个负的峰值。在图形的一个周期内中间的一段频率内，相位变化较为平缓，并且其绝对值普遍较低。

- 23 -

实验四：信号的DFT分析

一、实验题目：信号的DFT分析

计算余弦序列x(n)?cos(n)RN(n)的DFT。分别对N=10、16、22时计算

8?DFT，绘出X(k)幅频特性曲线，分析是否有差别及产生差别的原因。 二、实验目的：

（1）通过实验复习巩固课本中有关DFT的知识 （2）能够应用matlab绘制不同点的DFT图形 （3）能够根据所得图形化结果对信号进行分析 三、解题分析：

在matlab信号处理工具箱中，matlab提供了4个内部函数用于计算DFT和IDFT，它们分别是：fft(x)，fft(x,N)，ifft(X)，ifft(X,N)。 fft(x) 计算M点的DFT。M是序列x的长度，即M=length(x)。 fft(x,N) 计算N点的DFT。若M>N，则将原序列截短为N点序列，再计算其N点DFT；若M

ifft(X) 计算M点的IDFT。M是序列X的长度。

ifft(X,N) 计算N点IDFT。若M>N，则将原序列截短为N点序列，再计算其N点IDFT；若M

- 24 -

四、实验程序

计算余弦序列x(n)?cos(?8n)RN(n)的DFT。分别对N=10、16、22时计算DFT，

绘出X(k)幅频特性曲线，分析是否有差别及产生差别的原因。

程序： 当取N=10时

n=0:9;

x=cos(n*(pi/8)); y=fft(x);

stem(abs(y)); grid on;

当取N=16时

n=0:15;

x=cos(n*(pi/8)); y=fft(x);

stem(abs(y)); grid on;

- 25 -

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库信号与系统仿真实验报告(5)在线全文阅读。