基于dsp的快速傅立叶变换fft的实现汇编语言.docVIP

快速傅立叶变换（FFT）的实现 一、实验目的 1.了解FFT的原理及算法； 2.了解DSP中FFT的设计及编程方法； 3.熟悉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——DIT FFT信号流图如图5.1所示： 图5.1 时间抽取的FFT—DIT FFT信号流图 FFT算法主要分为以下四步。 第一步 输入数据的组合和位倒序 把输入序列作位倒序是为了在整个运算最后的输出中得到的序列是自然顺序。 第二步 实现N点复数FFT 第一级蝶形运算；第二级蝶形运算；第三级至log2N级蝶形运算； FFT运算中的旋转因子是一个复数，可表示： 为了实现旋转因子的运算，在存储空间分别建立正弦表和余弦表，每个表对应从0度到180度，采用循环寻址来对正弦表和余弦表进行寻址。 第三步 功率谱的计算 X(k)是由实部和虚部组成的复数：；计算功率谱时只需将FFT变换好的数据，按照实部和虚部求它们的平方和，然后对平方和进行开平方运算。 第四步 输出FFT结果。 三、实验内容与步骤 １. 本实验要求使用FFT变换求一个时域信号的频域特性，并从这个频域特性求出该信号的频率值，然后使用DSP汇编语言对FFT的DSP编程。 ２. 将计算机与DES-320PP—A教学实验系统连接好，打开实验箱电源，然后运行CCS软件。 2.新建fft.pjt工程文件，并输入FFT的DSP程序及命令链接文件。如下： ************************************** *Radix-2,DIT,Real-input FFT Program * * fft.asm * ************************************* .mmregs .global reset,start,sav_sin,sav_idx,sav_grp .def start,_c_int00 .data DATA .space 1024 .copy fft.inc N .set 128 LOGN .set 7 sav_grp .usect tempv,3 sav_sin .set sav_grp+1 sav_idx .set sav_grp+2 OUTPUT .usect OUTPUT,256 BOS .usect stack,0Fh TOS .usect stack,1 .copy twiddle1.inc .copy twiddle2.inc .text _c_int00 b start nop nop start: 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+0B plend: 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