基于DSP的快速傅立叶变换(FFT)的实现(汇编语言).doc

快速傅立叶变换(FFT)的实现一、;;;二、实验原理FFT是一种高效实现离散付立叶变换的算法,把信号从时域变换到频域,在频域分析处理信息。对于长度为N的有限长序列x(n),它的离散傅里叶变换为:,称为旋转因子,或蝶形因子。在x(n)为复数序列的情况下,计算X(k):对某个k值,需要N次复数乘法、(N-1)次复数加法;对所有N个k值,需要次复数乘法和N(N-1)次复数加法。对于N相当大时(如1024)来说,直接计算它的DFT所作的计算量是很大的,FFT的基本思想在于:利用的周期性即:对称性:将原有的N点序列分成两个较短的序列,这些序列的DFT可以很简单的组合起来得到原序列的DFT。按时间抽取的FFT——:—DITFFT信号流图FFT算法主要分为以下四步。第一步输入数据的组合和位倒序把输入序列作位倒序是为了在整个运算最后的输出中得到的序列是自然顺序。第二步实现N点复数FFT第一级蝶形运算;第二级蝶形运算;第三级至log2N级蝶形运算;FFT运算中的旋转因子是一个复数,可表示:为了实现旋转因子的运算,在存储空间分别建立正弦表和余弦表,每个表对应从0度到180度,采用循环寻址来对正弦表和余弦表进行寻址。第三步功率谱的计算X(k)是由实部和虚部组成的复数:;计算功率谱时只需将FFT变换好的数据,按照实部和虚部求它们的平方和,然后对平方和进行开平方运算。第四步输出FFT结果。三、,并从这个频域特性求出该信号的频率值,然后使用DSP汇编语言对FFT的DSP编程。-320PP—A教学实验系统连接好,打开实验箱电源,S软件。,并输入FFT的DSP程序及命令链接文件。如下:***************************************Radix-2,DIT,Real-inputFFTProgram**************************************** .mmregs .globalreset,start,sav_sin,sav_idx, start, .space ""N .set 128LOGN .set 7sav_grp .usect "tempv",3sav_sin .set sav_grp+1sav_idx .set sav_grp+2OUTPUT .usect "OUTPUT",256BOS .usect "stack",0FhTOS .usect "stack",1 .copy "" .copy "" .text_c_int00 bstart nop nopstart: STM #TOS,SP LD #0,DP SSBX FRCT STM #2*N,BK STM #INPUT,AR3 STM #DATA,AR7 MVMM AR7,AR2 STM #N-1,BRC RPTBD plend-1 STM #N,AR0 LDM AR3,A READA*AR2+ ADD#1,A READA*AR2+ MAR *AR3+0Bplend: STM #0,BK LD #-1,ASM MVMM AR7,AR2 STM #DATA+2,AR3 STM #N/2-1,BRC LD *AR2,16,A RPTBD s1end-1 STM #3,AR0 SUB *AR3,16,A,B ADD *AR3,16,A STH A,ASM,*AR2+ ST B,*AR3+ ||LD *AR2,A SUB *AR3,16,A,B ADD *AR3,16,A STH A,ASM,*AR2+0 ST B,*AR3+0% ||LD *AR2,As1end: MVMM AR7,AR2 STM #DATA+4,AR3 STM #N/4-1,BRC LD *AR2,16,A RPTBD s2end-1 STM #5,AR0 SUB *AR3,16,A,B ADD *AR3,16,A STH A,ASM,*AR2+ ST B,*AR3+ ||LD *AR2,A SUB *AR3,16,A,B ADD *AR3,16,A STH A,ASM,*AR2+ STH B,ASM,*AR3+ MAR *AR3+ ADD *AR2,*AR3,A SUB *AR2,*AR3-,B STH A,ASM,*AR2+ SUB *AR2,*AR3,A