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

《基于dsp的快速傅立叶变换(fft)的实现(汇编语言)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、快速傅立叶变换（FFT）的实现一、实验目的1.了解FFT的原理及算法；2.了解DSP中FFT的设计及编程方法；3.熟悉FFT的调试方法；二、实验原理FFT是一种高效实现离散付立叶变换的算法，把信号从时域变换到频域，在频域分析处理信息。对于长度为/V的有限长序列％（«）,它的离散傅里叶变换为：.V-In-0=e-j{^/N}nk,称为旋转因子，或蝶形因子。在x（zO为复数序列的情况下，计算汾幻：对某个々位，需要7V次复数乘法、（7V-1）次复数加法;对所有/V个々值，需要；V2次S数乘法和/V（A

2、H）次S数加法。对于N相当大时（如1024）来说，直接计算它的DFT所作的计算量是很大的，FFT的基本思想在于：利用W;v=ej（+、nk的周期性即：W^+N=对称性：=-W^,+N/2将原有的N点序列分成两个较短的序列，这些序列的DFT可以很简单的组合起来得到原序列的DFT.按时间抽取的FFT——DITFFT信号流图如图5.1所示：义3⑼4(0)卯）工(2)工⑶-V(4)•H5)瑚X(7)图5.1时间抽取的FFT—DITFFT信号流图FFT算法主要分为以下四步。第一步输入数据的组合和位倒序把输

3、入序列作位倒序是为了在整个运算最后的输出屮得到的序列是向然顺序。第二步实现N点S数FFT第一级蝶形运算；第二级蝶形运算；第三级至log2；V级蝶形运算；FFT运算中的旋转因子是一个复数，可表示：W^=e~j2^,N=cos^-jsin^-•'NN为了实现旋转因子W;v的运算，在存储空间分別建立正弦表和余弦表，毎个表对应从0度到180度，采用循环寻址來对正弦表和余弦表进行寻址。第三步功率谱的计算X⑻是由实部;幻和虚部X,(幻组成的复数：X(/c)=X/?(jt)+7X/(/:)；汁算功率谱时只需将

4、FFT变换好的数据，按照实部；幻和虚部；^(幻求它们的平方和，然后对平方和进行开平方运算。第四步输出FFT结果。三、实验内容与步骤1.本实验要求使川FFT变换求一个时域信号的频域特性，并从这个频域特性求出该信号的频率值，然后使用DSP汇编语言对FFT的DSP编程。2.将计算机与DES-320PP—A教学实验系统连接好，打开实验箱电源，然后运行CCS软件。2.新建fft.pjt工程文件，并输入FFT的DSP程序及命令链接文件。如下：*1、*1、*1、'I、'I、*1、zj、zj、zj、zj、zj、

5、,1、zj、zj、zj、zj、*Radix-2,D1T,Real-inputEETProgram**fft

6、rp+1savidx•setsav_grp+2OUTPUT.usect"OUTPUT",256BOS.usect"stack",OEhTOS.usect〃stack",1.copy/ztwiddlel.inc"•copy/，twiddle2.inc"cintOO•textbstartnopnopstart:STM#TOS,SPLD#0,DPSSBXFRCTSTM#2*N,BKSTM#INPUT,AR3STM#DATA,AR7鵬1MAR7,AR2STM#N-1,BRCRPTBDplend-1STM

7、眠AROLDMAR3,AREADA*AR2+ADD#1,AREADA*AR2+MAR*AR3+0Bpiend:STM#0,BKLD#-l,ASM鵬1MAR7,AR2STM#DATA+2,AR3STM#N/2-l,BRCLD*AR2,16,ARPTBDslend-1STM#3,AROSL'B峨3,16,A,BADD*AR3,16,ASTHA,ASM,*AR2+STB,*AR3+IlLD*AR2,ASL'B峨3,16,A,BADD*AR3,16,ASTHA,ASM,*AR2+0STB,*AR3+0%

8、IlLD*AR2,Aslcnd:鵬1MAR7,AR2STM#DATA+4,AR3STM#N/4-l,BRCLD*AR2,16,ARPTBDs2end-lSTM#5,AROSL'B峨3,16,A,BADD*AR3,16,ASTHA,ASM,*AR2+STB,*AR3+IlLD*AR2,ASL'B峨3,16,A,BADD*AR3,16,ASTHA,ASM,*AR2+STHB,ASM,*AR3+MAR峨3+ADD*AR2,*AR3,ASL'B*AR2,*AR3-,BSTHA,ASM,*AR2+SIB*