2048点FFT在TMS320C240x定点DSP上的实现.doc

2048点FFT在TMS320C240x定点DSP上的实现 ? ?? ?? ?? ? Accomplishment of 2048 points FFT in TMS320C240x fixed-point DSP? ?? ?? ?? ?? ?? ?? ?? ? 文其林 摘要：针对线阵CCD输出数据的处理需要，介绍了2048点FFT在TMS320C240x定点DSP上的实现。FFT算法程序采用汇编语言编写，可通过C语言函数调用，因此具有实现简单、高效的特点，且扩展性和兼容性强，可运行于嵌入式线阵CCD光谱探测系统中。关键词：FFT；DSP；劈分算法；位反转算法Abstract：Because of the demand of data processing in linear CCD，this paper introduces the realization of 2048 points FFT in TMS320C240x fixed-point DSP. Arithmetic program of FFT is programmed in assembly code，and it can be called by C function. So it has the features of easy realization，high efficiency，good expansibility and compatibility. It had run in the embedded system for linear CCD spectral detection successfully.Key words：FFT；DSP；splitting arithmetic；bit-reversed arithmetic引言傅立叶变换是一种将信号从时域转变为频域表示的变换手段，它在信号的频谱分析以及系统的分析、设计中得到了广泛的应用。在计算机系统中，实际上是以离散傅立叶变换（DFT）的方式处理数据。由于DFT的运算量比较大，实际应用中常使用DFT的快速算法－快速傅立叶变换（FFT）。快速傅立叶变换算法有基2算法、基4算法等，还有按时间抽取（DIT）和按频率抽取（DIF）的算法。目前，由于线阵CCD在光谱探测、光学传感等方面的广泛应用，常需要嵌入式芯片（如DSP）对CCD输出的数据进行实时分析处理。用于工业探测的线阵CCD与我们通常认识的面阵CCD有所不同。我们日常生活中出现的面阵CCD少说也有几十万像素，高的可达上千万像素，常用在数码相机等设备上。相比较而言，用于工业探测的线阵CCD像素不高，但在光谱响应范围、几何精度、动态范围有自己独特的优势。比如我们所用的TCD1208 CCD就是一款面向光谱探测用的线阵CCD，有2160像元，光谱响应范围400纳米-1100纳米，并且动态范围大，灵敏度高，稳定可靠。整个嵌入式线阵CCD光谱探测系统的结构框图如图1所示，在DSP上要完成CCD驱动、A/D转换、FFT运算、控制液晶显示等功能，其中FFT运算是重点。 考虑到CCD是2160像元，为了方便运用基2算法计算FFT，又不损失太多的像元，我们在DSP上计算2048点FFT。由于TI公司提供的FFT程序只能计算512点FFT，因此，就需要我们编写能够实现2048点FFT的程序。以FFT在嵌入式线阵CCD光谱探测系统中的应用为例，重点介绍2048点实数基2DIT－FFT在TMS320LF2407定点DSP中的实现。此程序也可实现1024点复数FFT。由于TMS320C240x系列DSP的硬件内核与汇编指令完全相同，此程序可以直接移植到所有TMS320C240x DSP上。根据此程序，可以方便地扩展至4096点FFT，用于5000像元CCD的数据处理。1 程序介绍与实现1.1程序流程FFT的计算量比傅立叶变换的计算量减少了，但是FFT要做到大点数，实时运算，对于普通的单片机来说还是一件比较困难的事。一方面，FFT需要对原始自然序列进行码位倒序排列；另一方面，蝶形运算是复数运算，需要多次查表相乘运算才能实现。DSP控制器就是针对这些需求而设计的专用芯片：具有专为FFT算法而设计的反序间接寻址；可实现增/减1或增/减一个变址量的间接寻址方式，为各种查表运算提供方便；能在一个指令周期完成乘和累加操作，提高了乘法运算的速度[1]。TMS320LF2407定点DSP是一款专为工业控制、电机控制和数字信号处理等用途而设计的DSP，具备单周期乘加指令，具有FFT反序间接寻址功能，最高运行速度为40MIPS。为了充分利用DSP芯片特有的反序间接寻址等功能，FFT算法程序采用汇编语言编写，主程序采用C语言，