基于TMS320LF2407的FFT算法的实现

       1 快速傅里叶变换的原理

　　非周期性连续时间信号x(t)的傅里叶变换可以表示为

　　式中计算出来的是信号x(t)的连续频谱。但是，在实际的控制系统中能够得到的是连续信号x(t)的离散采样值x(nT)。因此需要利用离散信号x(nT)来计算信号x(t)的频谱。

　　有限长离散信号x(n)，n=0，1，…，N-1的DFT定义为：

　　可以看出，DFT需要计算大约N2次乘法和N2次加法。当N较大时，这个计算量是很大的。利用WN的对称性和周期性，将N点DFT分解为两个N/2点的 DFT，这样两个N/2点DFT总的计算量只是原来的一半，即(N/2)2+(N/2)2=N2/2，这样可以继续分解下去，将N/2再分解为N/4点 DFT等。对于N=2m 点的DFT都可以分解为2点的DFT，这样其计算量可以减少为(N/2)log2N次乘法和Nlog2N次加法。图1为FFT与DFT-所需运算量与计算点数的关系曲线。由图可以明显看出FFT算法的优越性。

　　将x(n)分解为偶数与奇数的两个序列之和，即

　　x1(n)和x2(n)的长度都是N/2，x1(n)是偶数序列，x2(n)是奇数序列，则

　　其中X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点DFT。由于X1(k)和X2(k)均以N/2为周期，且WN k+N/2=-WN k，所以X(k)又可表示为：

　　上式的运算可以用图2表示，根据其形状称之为蝶形运算。依此类推，经过m-1次分解，最后将N点DFT分解为N/2个两点DFT。图3为8点FFT的分解流程。

　　FFT算法的原理是通过许多小的更加容易进行的变换去实现大规模的变换，降低了运算要求，提高了与运算速度。FFT不是DFT的近似运算，它们完全是等效的。

　2 快速傅里叶算法在TMS320LF2407上的实现

　　根据FFT算法的特点，处理器要在一个指令周期内完成乘和累加的工作，因为复数运算要多次查表相乘才能实现。其二就是间接寻址，可以实现增/减1个变址量，方便各种查表方法。再次，FFT变换的输入序列x(n)是按所谓的码位倒序排列的，处理器要有反序间接寻址的能力。DSP控制器专门设计了特有的反序间接寻址，并能在一个指令周期内完成乘和累加的运算。因此，对数字信号的分析处理，DSP比其它的处理器有绝对的优势。本文采用TI公司C2000系列TMS320LF2407芯片来实现FFT算法。

　　TMS320LF2407定点DSP是一款专为工业控制、电机控制和数字信号处理等用途而设计的DSP，具备单周期乘加指令，具有FFT反序间接寻址功能，最高运行速度为40MIPS。为了充分利用DSP芯片特有的反序间接寻址等功能，FFT算法程序采用汇编语言编写，主程序采用C语言，因此程序具有良好的兼容性和可扩展性。

　　主程序流程图如图4所示。系统初始化主要完成DSP的系统控制和状态寄存器、等待状态发生器控制寄存器、中断寄存器等的必要设置。

　　本程序采样函数为：x=sin(20πt)，采样频率为640Hz。

　　输入数据波形如图5所示。一般情况下，我们只关心信号频域的幅度谱。幅度谱|X(k)|2的计算：X(k)=XR(k)+jX(k)，|X(k)2|=|Xr(k)|2+|Xi(k)|2。FFT计算结果的信号幅度谱|X(k)|2如图6所示。

　　输入信号频率是10Hz，根据公式f=kfs/N，f是原始信号的频率，k表示峰值出现的位置，fS是采样频率，N是计算的点数，从幅度谱中看出，峰值出现在k=1处，那么，f=1×640/64=10，与原始信号的实际频率一致，说明计算结果正确。