实验四离散傅立叶变换及其快速算法.docVIP

实验四 离散傅立叶变换及其快速算法 一．实验目的 理解信号变换的基本概念 理解离散傅立叶变换的基本概念 掌握快速傅立叶变换的应用方法 实验原理 信号是数字信号处理领域中最基本、最重要的概念。而数字信号变换技术，又是对信号进行处理操作的最基本的有效途径之一。因此，数字信号变换技术，便成为数字信号处理领域中专业人员所必须要张我的一项最基本的技能。 简单地说，数字信号变换技术就是为了处理操作上的方便和可能，通过数学变换，将一个域内的信号变换映射倒另一个域内的信号的方法。常用的数字信号变换主要有：傅立叶变换、离散余弦变换（DCT）、Z变换、Chirp z变换、Hilbert变换等。这些变换，都有着各自的理论和其应用背景。 MATLAB中的工具箱对这几种典型的变换，都提供了相对应的、具体的应用函数。这可以使得工程人员大大节省无谓的工作量，从而将主要精力放到新技术的创新和研发上面。下面将对这几种变换的含义和应用进行具体的介绍。 傅立叶变换是信号分析和处理的重要工具。有限长序列作为离散信号的一种，在数字信号处理种占有着极其重要的位置。对于有限长序列，离散傅立叶变换不仅在理论上有着重要的意义，而且有快速计算的方法－的傅立叶变换可以表示为 ＝ 逆变换为 在这里，是模拟角频率。可以看到，时域的连续函数造成频域的非周期谱，时域的非周期性造成频域的连续谱。 结论：非周期连续时间函数对应于一非周期连续频域变换函数。 2、周期连续时间信号的傅立叶级数 周期为的周期性连续时间信号傅立叶变换是离散频域函数，可表示为 逆变换为 这就是经常称之为傅立叶级数的变换形式。在这里，也是模拟角频率。可以看到，时域的连续函数造成频率域的非周期谱，频域函数的离散造成时域函数的周期性。 结论：周期连续时间函数对应于一非周期离散频域变换函数。 3、非周期离散时间信号的傅立叶变换可以表示为 逆变换为 在这里，是数字频率，它和模拟角频率的关系为。可以看到，时域的取样对应于频域的周期延拓，而时域函数的非周期性造成频域的离散谱。 结论：非周期离散时间函数对应于一周期连续频域变换函数。 4、周期离散时间信号的傅立叶变换 周期离散时间信号的傅立叶变换－离散傅立叶变换，可以表示为 逆变换为 可以看到，时域的取样对应于频域的周期延拓，而时域函数的周期性造成频域的离散谱。 结论：周期离散时间函数对应于一周期离散频域变换函数。 （二）离散傅立叶变换 离散傅立叶级数变换是周期序列，仍不便于计算机计算。但离散傅立叶级数虽是周期序列，却只有个独立的数值，所以它的许多特性可以通过有限长序列延拓来得到。对于一个长度为的有限长序列，也即只在个点上有非零值，其余皆为零，即： 把序列以为周期进行周期延拓得到周期序列，则有： 所以，有限长序列的离散傅立叶变换（DFT）为： 逆变换为 若将DFT变换的定义写成矩阵形式，则得到 X=A﹒ Dftmtx 函数：用来计算DFT变换矩阵A的函数 调用方式 A＝dftmta（n）：返回n×n的DFT变换矩阵A。若x为给定长度的行向量，则y＝x*A，返回x的DFT变换y。 Ai＝conj（dftmtx（n））/n；返回n×n的IDFT变换矩阵Ai。 应用说明 【实例5-1】 A=dftmtx(4) Ai=conj(dftmtx(4))/4 运行结果 A = 1.0000 1.0000 1.0000 1.0000 1.0000 0 - 1.0000i -1.0000 0 + 1.0000i 1.0000 -1.0000 1.0000 -1.0000 1.0000 0 + 1.0000i -1.0000 0 - 1.0000i Ai = 0.2500 0.2500 0.2500 0.2500 0.2500 0 + 0.2500i -0.2500 0 - 0.2500i 0.2500 -0.2500 0.2500 -0.2500