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

(数字信号处理)实验报告

实验名称 实验三 离散傅立叶变换及其特性验证 实验时间 年 月 日 专业班级 学 号 姓 名

成 绩 教师评语：

一、实验目的

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

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

二、实验原理与计算方法

1、离散时间傅立叶变换

如果序列x(n)满足绝对可和的条件，即

n

则其离散时间傅立叶变换定义为： |x(n)| ，

X(ej ) F[x(n)]

n

x(n)e

j n

(1)

假设序列x(n)在n1 n nN(即不一定在[0, N-1])有N个样本，要估计下列各点上的X(ej )：

k

M

M

k，　k 0,1,2...,M

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

j

l 1N

j

knl

x(nl)，　k 0,1,2...,M (2)

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

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

kn jM

W e;n1 n nN，　k 0,1,2...,M

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

exp jMkn

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

T

XT xT exp jnk

M

其中nTk是一个N×(M+1)维矩阵。用MATLAB实现如下： k=[0:M]; n=[n1:n2];

X=x*(exp(-j*pi/M)).^(n’*k); 2、离散傅立叶变换

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