用matlab验证验证DTFT的特性,实验三离散傅立叶变换及其特性验证.doc

实验三离散傅立叶变换及其特性验证

实验名称：实验三 离散傅立叶变换及其特性验证

一、实验目的

1、掌握离散时间傅立叶变换(DTFT)的计算方法和编程技术。

2、掌握离散傅立叶变换(DFT)的计算方法和编程技术。

3、理解离散傅立叶变换(DFT)的性质并用Matlab进行验证。

二、实验原理与计算方法

1、离散时间傅立叶变换

如果序列x(n)满足绝对可和的条件，即，则其离散时间傅立叶变换定义为：

(1)

如果x(n)是无限长的，则不能直接用Matlab由x(n)计算X(ej()，但可以用它来估计X(ej()表达式在[0,π]频率区间的值并绘制它的幅频和相频(或实部和虚部)曲线。

如果x(n)是有限长的，则可以用Matlab对任意频率(处的X(ej()进行数值计算。如果要在[0,π]间按等间隔频点估计X(ej()，则(1)式可以用矩阵－向量相乘的运算来实现。假设序列x(n)在(即不一定在[0, N-1])有N个样本，要估计下列各点上的X(ej()：

它们是[0,π]之间的(M+1)个等间隔频点，则(1)式可写成：

(2)

将{x(nl)}和{X(ej(k)}分别排列成向量x和X，则有：

X=Wx (3)

其中W是一个(M+1)×N维矩阵：

将{k}和{n}排成列向量，则

在Matlab中，把序列和下标排成行向量，对(3)式取转置得：

2、离散傅立叶变换

一个有限长序列的离散傅立叶变换对定义为：

(4)