实验八 傅立叶变换.pdf

实验八 傅立叶变换--第1页 实验八 傅立叶变换 一、实验目的 利用MATLAB 进行离散时间傅立叶变换；分析时域采样频率与频谱混叠；由离散序列 恢复模拟信号。 二、实验原理及实验内容 无限长离散序列的傅立叶变换称为离散时间傅立叶变换，它的正逆变换的形式如下：    F()  f (n)ej n n 1     j n  f (n)  F( )e d 2  可见它在时域上是离散序列，而在频域上是连续函数，即具有连续的频谱。如果时域 序列是有限长的，并把它做周期延拓，它的频谱就向一些等间隔的点靠拢。随着延拓周期数 的无限增加，其连续频谱就收敛于周期性的离散点，此时就要用离散傅立叶变换来处理。 一般规则如表所示 表 傅立叶变换一般规则 时域信号（傅立叶反变换） 频谱曲线（傅立叶变换） 变 换 名 称 连续信号 连续频谱 傅立叶变换 离散信号（有限样本点） 周期性连续频谱 离散时间傅立叶变换 多周期离散信号 连续频谱向离散点集中 离散时间傅立叶变换 周期性离散信号 周期性离散频谱 离散傅立叶变换 周期性连续信号 离散频谱 傅立叶级数 1．离散时间傅立叶变换 取一个周期的正弦信号，作8 点采样，求它的连续频谱。然后对该信号进行N 个周期 延拓，再求它的连续频谱，比较分析其结果。 建模： 为了求离散信号的连续频谱，不能直接用fft 函数，可从离散傅立叶变换的定义出发。    X()  x(n)ej n n 设置一系列较密频率点wi，就可求出一系列X(wi)=x*exp(-j*wi*n’)，因为x 是行向量， n’以及exp(-j*wi*n’)为同长的列向量，两者的向量点乘就包含了逐项相乘后的连加演算，得 出一个标量X(wi) ，即该频点上的频率响应。不同的wi 可用for 循环来解决。但MATLAB 中往往可以用元素群运算来代替一个for 循环，为此把w 也设成行向量，放在n’之后，n’*w 及exp(-j*n’*w)就成为一个矩阵，其行数与x 相同，列数与w 的长度相同。用x 左乘它以后， 各列分别得出相应wi 处的响应，最后得到X=x*ex