资源描述：
《实验报告1FFT的DSP实现》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验报告实验名称FFT的DSP实现课程名称任课老师高速DSP原理与座用刘建国姓学班日名号级期200930081734912012年6月19日一、实验目的1、进一步熟悉C28X的汇编指令2、加深对FFT算法的理解，体会利用吟的周期性和对称性来减少DFT的运算次数的原理。二、实验要求用汇编语言实现FFT算法，并画流程框图。设序列x(n)的长度为N,且满足N=2M,M为自然数。"的周期性和对称性来减小DFT运算次数，可以得到三、实验原理利用旋转因子快速算法(FFT)。N=疋点的FFT共进行M级运算，每级由
2、N/2个典型运算组成。同一级中,每一个蝶形运算的商个输入数据只对计算木蝶形运算有用，而但每个蝶形运算的输入、输出数据结点乂同在一条水平线上，这就意味着计算完一个蝶形运算后，所得输出数据可立即存入原输入数据所占用的存储单元。这样，经过M级运算后，原來存放输入序列数据的N个存储单元(数组A)中使依次存放X(k)的N个值。N点DIT-FFT运算流图中，每级有N/2个蝶形。每个蝶形都要乘以因子W#,称其为旋转因子，p为旋转因子的指数。但各级的旋转因子和循环方式都有所不同。N点序列的DFT和IDFT变换定义式
3、如下:x［幻=ZX［诃弃斤=01八_人编程思想：第L级中，每个蝶形的两个输入数据和距B=2L~l个点；每级有B个不同的旋转因子；同一旋转因子对应着间隔为2厶点的2旳"个蝶形。总结上述运算规律，可采用下述运算方法。先从输入端开始，逐级进行，共进行M级运算。在进行第L级运算时，依次求出B个不同的旋转因子，每求出一个旋转因子，就计算完它对应的所有2‘w"个蝶形。这样我们可以用如下的三重循环程序实现DIT-FFT运算。DIT-FFT运算流程框图读入X5),MN=2^M倒丿芋L=l,MB=2"(L-1)^―
4、・K=J,N-l,2^L']/T<=A(K)十A*(k+B)*W“pA*(k+B)<=A(K)-A*(k+B)*W"pA(K)<=T输山结束另外，DIT-FFT算法运算流图的输出X(k)为口然顺序，但为了适应原位运算，其输入序列不是安x(n)的自然序列排列，这种经过M次偶奇抽选后的排列成为序列x(n)的倒序。因此，在运行M级蝶形Z前应先对序列x(n)进行倒序。四.实验环境软件环境:CCS3.3硬件环境：N五、实验过程、数据记录、处理及结论用matlab产生一
5、个单频信号的1024点采样数据程序段如下：clearall;clct=l:1:1024;y=500*sin(t);y=fix(y);sprintf.word%d',y)然后作为FFT的输入信号的时域波形EFT的结果23eP18e*8-1.4e*8-92e*7-46e*7-0-46e*792e*7•28e・8.011.122233344455.66677738891001111221331441561671781®919909,46276)LinAuto118e*6-23e七六、讨论FFT（快速傅
6、里叶变换）是《数字信号处理》课程的一个核心内容。有一次，数字信号处理课老师在课堂上布置了预习时域抽取算法FFT原理的作业。我通过画8点DFT的三次时域抽取分解运算流图，加深了口己对蝶形运算的理解。在这个基础上通过推导，我比较了用FFT和DFT计算长度N=2"的序列的乘法和加法次数，直观的明白了DIF一FFT算法相对于DFT算法在计算速度上的优越性。七.附录（原程序代码）程序代码如下：•sect〃•sd〃var_x:.copy〃VAR.asm"•data.globalmystart,_c_int00M
7、:.byteOOhL:.byte9hB:.wordOOhJ:.wordOOhPl:-wordOOhK:•wordOOhSAVE_Y:・copy"FFT_Input.asmzzYINZI:.copy〃YINZI・lxt〃•textcintOO:mystart:SETCOBJMODEC28ADDR.C28JVMODEN.set1024YZ.set512MOVLXAR2,ttFirstlnMOVLXAR3,#VARFIRSTMOV@ARO,#NMOVAR1,#N-1NOP*,ARP2SETCAMODEMO
8、V*BRO++,AL,ARP1XBANZLOOP,,ARP2MOV@@L,#1;addarl,#1MOVAL,#2;;b=21-lMOVT,@@L-1LSLAL,TMOV@@B,ALMOV@@J,#0MOVAL,@@MSUBAL,@@LMOVT,ALMOVAL,#2LSLAL,TMOVT,ALMPYACC,T,@@JMOVL@@P1,ACCMOVAL,@@JMOV@@K,AL•LP_AMODELOOP:MOVAL,*++,ARP3L00P1:L00P2:L00P3:M
