DFT算法原理FFT的算法原理.docVIP

《数字信号处理》课程论文 离散傅里叶变换(DFT) 算法简介及其应用举例 姓 名：安昱 学 号：12011243986 专 业：通信工程 班 级：2011级（1）班 指导老师：武永峰 学 院：物理电气信息学院 完成日期：2013年11月 11 日 离散傅里叶变换(DFT)算法简介及其应用举例 （安昱 12011243986 2011级1班） [摘要]：离散傅立叶变换（DFT）实现了信号首次在频域表示的离散化，使得频域也能够用计算机进行处理。并且这种DFT变换可以有多种实用的快速算法。使信号处理在时、频域的处理和转换均可离散化和快速化。最后就该项目做了总结。 [关键词]DFT算法 Matlab语言 频域采样 DFT应用 DFT算法原理 快速傅氏变换（FFT）是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。它对傅氏变换的理论并没有新的发现，但是对于在计算机系统或者说数字系统中应用离散傅立叶变换，可以说是进了一大步。 设x(n)为N项的复数序列，由DFT变换，任一X（m）的计算都需要N次复数乘法和N-1次复数加法，而一次复数乘法等于四次实数乘法和两次实数加法，一次复数加法等于两次实数加法，即使把一次复数乘法和一次复数加法定义成一次“运算”（四次实数乘法和四次实数加法），那么求出N项复数序列的X（m）,即N点DFT变换大约就需要N2次运算。当N=1024点甚至更多的时候，需要N2=1048576次运算，在FFT中，利用WN的周期性和对称性，把一个N项序列（设N=2k,k为正整数），分为两个N/2项的子序列，每个N/2点DFT变换需要（N/2）2次运算，再用N次运算把两个N/2点的DFT变换组合成一个N点的DFT变换。这样变换以后，总的运算次数就变成N+2（N/2）2=N+N2/2。继续上面的例子，N=1024时，总的运算次数就变成了525312次，节省了大约50%的运算量。而如果我们将这种“一分为二”的思想不断进行下去，直到分成两两一组的DFT运算单元，那么N点的DFT变换就只需要Nlog2N次的运算，N在1024点时，运算量仅有10240次，是先前的直接算法的1%，点数越多，运算量的节约就越大，这就是FFT的优越性。 二、FFT的算法原理 FFT算法的输出X(K)为自然顺序，但为了适应原位计算，其输入序列不是按x(n)的自然顺序排序，这种经过M-1次奇偶抽选后的排序为序列的倒序。因此，在运算之前应先对序列x(n)进行倒序。倒序的规律就是把顺序数的二进制位倒置，即可得到倒序值。倒序数是在M位二进制数最高位加一，逢2向右进位。M位二进制数最高位的权值为N/2，且从左到右二进制位的权值依次为你N/4，N/8，…，2，1。因此，最高位加一相当于十进制运算J+N/2。（J表示当前倒序数的十进制数值） 三、DFT的定义 设序列x(n)长度为M，定义x(n)的N点DFT为 式中，N称为离散傅里叶变换区间长度，要求N ≥ M。为 书写简单，令，因此通常将N点DFT表示为 定义X(k)的N点离散傅里叶逆变换(IDFT)为 四、DFT算法举例 1、用DFT(FFT)对时域离散信号进行频谱分析 给定参考信号如下： 用以8为周期进行周期性延托形成的周期序列 分别以变换区间N=8,16,32 对进行DFT(FFT),画出相应的幅频特性曲线 分别以变换区间N=8,16 对分别进行DFT(FFT),画出相应的幅频特性曲线 分别以变换区间N=4,8,16对分别进行DFT(FFT),画出相应的幅频特性曲线 对进行频谱分析, 请自己选择变换区间,要求画出幅频特性曲线 2、MATLAB函数简介 MATLAB是一门计算机编程语言，取名来源于Matrix Laboratory，本意是专门以矩阵的方式来处理计算机数据，它把数值计算和可视化环境集成到一起，非常直观，而且提供了大量的函数，使其越来越受到人们的喜爱，工具箱越来越多，应用范围也越来越广泛。本题要用到一维快速傅里叶变换函数fft、取模（实数取绝对值）函数abs和求相位函角函数angle、绘图函数plot和stem。这些函数都属于MATLAB基本函数 ● fft和ifft：一维快速傅里叶变换和一维快速傅里叶逆变换函数。调用各式如下： Xk=fft(xn