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 -