快速傅里叶变换程序设计解析.docVIP

1.设计主要内容及要求； 1）掌握DSP A/D转换器使用方法。 2）研究FFT原理以及利用DSP实现的方法。 3）编写A/D采样和FFT程序，调试，观察结果。 2.对设计论文撰写内容、格式、字数的要求； （1）.课程设计论文是体现和总结课程设计成果的载体，一般不应少于3000字。 （2）.学生应撰写的内容为：中文摘要、目录正文、参考文献等文理通顺，内容正确完整，书写工整，装订整齐。 封面、任务书、、中文摘要目录正文、参考文献1.2 设计要求 1.2.1设计目的 1）理解FFT的算法以及利用DSP实现的方法。 2）能熟练的调试程序并能观察其结果。 3）熟悉TMS320F2812系列DSP芯片的软件设计方法。 1.2.2基本要求 1）掌握DSP A/D转换器使用方法。 2）研究FFT原理以及利用DSP实现的方法。 3）编写A/D采样和FFT程序，调试，观察结果。 二 设计思路 2.1 FFT算法简介 快速傅里叶变换（FFT）是一种高效实现离散傅里叶变换（DFT）的快速算法，是数字信号处理中最重要的工具之一，它在声学，语音，信号处理等领域有着广泛的应用。是将信号从时域变换到频域的一种变换形式，是信号处理领域中一种重要的分析工具。离散傅里叶变换（DFT）是连续傅里叶变换在离散系统中的表现形式。 傅里叶变换分为连续傅里叶变换和离散傅里叶变换。离散傅里叶变换简称DFT，是对离散信号进行傅里叶变换的方法，其运算量大、复杂度与变换点数的二次方成正比，因而不适用于进行实时信号处理。20世纪60年代由Cooley和Tukey提出了快速傅里叶变换（FFT）算法，它是快速计算DFT的一种高效方法，可以明显地降低运算量，大大地提高DFT的运算速度，从而使DFT在实际中得到了广泛的应用，已成为数字信号处理最为重要的工具之一。 DSP芯片的出现使FFT的实现变得更加方便。由于多数的DSP芯片都能在单指令周期内完成乘法—累加运算，而且还提供了专门的FFT指令（如实现FFT算法所必需的比特反转等），使得FFT算法在DSP芯片上实现的速度更快。本节首先简要介绍FFT算法的基本原理，然后介绍FFT算法的DSP实现。 2.2 FFT算法原理 快速傅氏变换（FFT）是离散傅氏变换的快速算法，根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅里叶变换算法进行改进获得的。对傅氏变换的理论并没有新的发现，但是对在计算机系统或者数字系统中应用离散傅里叶变换，可以说进了一大步。设为N项的复数序列，由DFT变换，任一X(M)的计算需要N次复数乘法和N-1次复数加法，而一次复数乘法等于四次实数乘法和两次实数加法，一次复数加法等于两次实数加法，即使把一次复数乘法和一次复数加法定义成一次“运算”（四次实数乘法和四次实数加法），则求出N项复数序列的X(M)，即N点DFT变换大约需要次运算。当N=1024点甚至更多的时候，需要次运算，在FFT中，利用的周期性和对称性，把一个N项序列（设N=2k，k为正整数），分为两个项的子序列，每个点DFT变换需要次运算，再N次运算把两个点的DFT变换组合成一个N点DFT变换。这样变换以后，总运算次数变成。继续上面的例子，N=1024时，总的运算次数就变成了525312次，节省了大约50%的运算量。如果将这种“一分为二”思想不断进行下去，直到分成两两一组DFT运算单元，则N点DFT变换只需要次运算，N在1024时，运算量仅有10240次，是先前直接算法的1%，点数越多，运算量节约就越大，这就是FFT的优越性。 2.3 快速傅立叶变换算法 FFT的基本思想：将大点数的DFT分解为若干个小点数DFT组合，从而减少运算量。 根据对序列分解与选取方法的不同而产生了FFT的多种算法。算法分类如图2-1所示： 图2-1 FFT的几种算法 在这次设计中，采用最基础的按时间抽取的基-2FFT算法来实现快速傅里叶变换。为了将大点数的DFT分解为小点数的DFT运算，要求序列的长度N为复合数，最常用的是的情况（M为正整数）。该情况下的变换称为基-2FFT算法。下面讨论基-2FFT情况的算法。 先将序列x(n)按奇偶项分解为两组： （2-1） 将DFT运算也相应分为两组: （2-3） 至此，一个点DFT被分解为两个点的DFT。由式（2-3）能否将全部的点的解求出来？ 分析： （2-4） 和只有个点，则由式（2-4）只能求出X(k)的前个点的DFT。要求出全部点的，需要找出、和的关系，其中k=0,1,2......(N/2)-1。由式子（2-4）可得化简得 （2-5） 这样N点DFT可全部由下式确定出来： （2-6） 上式可用一