快速傅里叶变换(FFT)的DSP实现.docxVIP

快速傅里叶变换(FFT)的DSP实现 快速傅里叶变换（FFT）的DSP实现 （天津大学电子信息工程学院） 摘要：本文介绍了快速傅里叶变换（FFT）的快速高效的原理及实现方法，对快速傅立叶变换（FFT）的特点进行了研究和总结。对于快速傅立叶变换（FFT） 在TMS320C54X系列数字信号处理器（DSP）实现中出现的计算溢出等问题进行了分析并提出了解决方法，同时据此使用DSP实现了快速傅立叶变换（FFT）。 关键词：数字信号处理；快速傅立叶变换；反序；计算溢出 1 引言： 傅里叶变换是一种将信号从时域变换到频域的变换方式，在语音处理、图像处理、信号处理领域中都发挥了极大的作用，是一种重要的分析工具。离散傅里叶变换（DFT）是连续傅里叶变换在离散系统中的表现形式，具有非常广泛的应用。但是由于DFT的计算量很大，因此在很长一段时间里其应用受到限制。快速傅里叶变换（FFT）是实现普通离散傅里叶变换的一种高效方法，快速傅里叶变换（FFT）的出现使得傅里叶变换在实际中得到了广泛的应用。 快速傅里叶变换并不是一种新的变换，它是离散傅里叶变换的一种快速算法。它是DSP领域中的一项重大突破。由于考虑了计算机和数字硬件实现的约束条件，研究了有利于机器操作的运算结构，使DSP的计算时间缩短了一到两个数量级，还有效的减少了计算所需的存储容量，FFT技术的应用极大的推动了DSP的理论的技术的发展。 本文中使用的是由TI公司生产的TMS320C54系列的DSP。C54x系列DSP具有很高的操作灵活性和速度。它具有一个先进的修正哈佛结构、专门硬件逻辑的CPU、片内存储器、片内外设和专用的指令集、将C54xCPU和片内存储器与外设配置组合在一起的螺旋结构。这使得该系列可以满足电子市场众多领域的应用要求。 2 DSP在数字信号处理中的优势： 数字信号处理是一门广泛应用于许多领域的新兴学科。20世纪60年代以来，随着计算机和信息技术的飞速发展，数字信号处理技术应用而生并得到迅速广泛的应用。数字信号处理算法运算量大,需要执行大量的乘累加运算。并且常具有某些特定模式，大部分处理时间花在执行相对小循环的操作上。与此同时，还要求处理芯片有专门的接口，具有很高的数据吞吐能力。 DSP芯片，也称数字信号处理器，其特殊的结构可以用来快速的实现各种数字信号处理算法。DSP芯片一般具有如下的主要特点： （1）在一个指令周期内可完成一次乘法和一次加法。 （2）程序和数据空间分开，可以同时访问指令和数据。 （3）片内具有快速RAM，通常可通过独立的数据总线同时访问两块芯片。 （4）具有低开销或无开销循环及跳转的硬件支持。 （5）快速的中断处理和硬件I/O接口支持。 （6）具有在单周期内操作的多个硬件地址产生器。 （7）可以并行执行多个操作。 （8）支持流水线结构，使取指、译码、取操作数和执行等操作可以重叠执行。 1 3 离散傅立叶变换（DFT）及快速傅里叶变换（FFT）的原理： 离散傅立叶变换可由如下的公式表示出来： X(k)= N?1n?0?x?n?WknNknN k=0,1,2,…..,N-1 式中，W?e?2???j??kn?N??cos(2?kn2?kn)?jsin() 0?k,n?N?1 NN由于计算一个X(k)值需要N次复数乘法和（N-1）次复数加法，因而计算N个X(k)的值，需要N次复数乘法和Ｎ(N-1)次复数加法。每次复数乘法包括４次实数乘法和２次实数加法，每次复数加法包括两次实数加法，因此计算Ｎ点的DFT共需要４＊N次实数乘法和2N(2N-1)次实数加法。当Ｎ很大时，运算量是很可观的，因此需要使用改进的快速DFT 算法。 快速傅立叶变换（FFT）的基本原理：FFT算法是基于可以将一个长度为Ｎ的序列的离散傅立叶变换逐次分解为较短的离散傅立叶变换来计算这一基本原理的。本文中将使用的是按时间抽取（Decimation-in-Time）的基２FFT算法。使用N/2点DFT的方式获得FFT计算公式如下式： N/2?122X(k)??r?0x1(r)WrkN/2?WkNN/2?1r?0?x(r)W2rkN/2rk?X1(k)?WN/2X2(k) k?1,1,2,......,N?1从上式可知，X1(k)和X2(k)的计算分别需要2?N/2?次复数乘法和N(N-2)/2次复数 k加法，而WNX2(k)的计算需要N/2次复数乘法，所以共需要N2?N/2次复数乘法。每 2??次复数乘法包括４次实数乘法和２次实数加