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

快速傅立叶变换（FFT）的实现 一、 实验目的 了解FFT的原理及算法； 了解DSP中FFT的设计及编程方法； 熟悉FFT的调试方法； 二、 实验原理 FFT是一种高效实现离散付立叶变换的算法，把信号从时域变换到频域，在频域分析处 理信息。 对于长度为/V的有限长序列％（?）,它的离散傅里叶变换为： .V-I n-0 = e-j{^/N}nk ,称为旋转因子，或蝶形因子。 在x（zO为复数序列的情况下，计算汾幻： 对某个々位，需要7V次复数乘法、（7V-1）次复数加法; 对所有/V个々值，需要；V2次S数乘法和/V（AH）次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)瑚 义3⑼ 4(0) 卯） 工(2) 工⑶ -V(4) ?H5) 瑚 X(7) 图5.1时间抽取的FFT—DIT FFT信号流图 FFT算法主要分为以下四步。 第一步输入数据的组合和位倒序 把输入序列作位倒序是为了在整个运算最后的输出屮得到的序列是向然顺序。 第二步实现N点S数FFT 第一级蝶形运算；第二级蝶形运算；第三级至log2；V级蝶形运算； FFT运算中的旋转因子是一个复数，可表示： W^=e~j2^,N =cos^-jsin^- ? N N 为了实现旋转因子W;v的运算，在存储空间分別建立正弦表和余弦表，毎个表对应从0 度到180度，采用循环寻址來对正弦表和余弦表进行寻址。 第三步功率谱的计算 X⑻是由实部;幻和虚部X,(幻组成的复数：X(/c) = X/?(jt) + 7X/(/:)；汁算功率 谱时只需将FFT变换好的数据，按照实部；幻和虚部；^(幻求它们的平方和，然后对平 方和进行开平方运算。 第四步输出FFT结果。 三、实验内容与步骤 本实验要求使川FFT变换求一个时域信号的频域特性，并从这个频域特性求出该信号 的频率值，然后使用DSP汇编语言对FFT的DSP编程。 将计算机与DES-320PP—A教学实验系统连接好，打开实验箱电源，然后运行CCS软 件。 2.新建fft.pjt工程文件，并输入FFT的DSP程序及命令链接文件。如下： *1、*1、*1、I、I、*1、zj、zj、zj、zj、zj、,1、zj、zj、zj、zj、 *Radix-2,D1T,Real-input EET Program * * fft asm 氺 kIx xlx kIx xlx kIx xlx kIx xlx kIx kIx kIx kIx %2x Klx .nunregs ? global reset, start, sav_sin, sav idx, sav grp .def start, c intOO DATA data space 1024 .copy ?，fft. inc# N .set 128 LOGN .set 7 savgrp ?usect 〃tempv〃,3 sav sin .set sav grp+1 savidx ? set sav_grp+2 OUTPUT .usect OUTPUT, 256 BOS .usect stack, OEh TOS .usect 〃stack,1 .copy /ztwiddlel. inc ? copy /，twiddle2. inc c intOO ? text b start nop nop start: STM #TOS, SP LD #0, DP SSBX FRCT STM #2*N, BK STM #INPUT, AR3 STM #DATA, AR7 鵬1M AR7, AR2 STM #N-1,BRC RPTBD plend-1 STM 眠ARO LDM AR3, A READA *AR2+ ADD #1, A READA *AR2+ MAR *AR3+0B pi end: STM #0, BK LD #-l, ASM 鵬1M AR7, AR2 STM #DATA+2, AR3 STM #N/2-l,BRC LD *AR2, 16, A RPTBD slend-1 STM #3, ARO SLB 峨3, 16, A, B ADD *AR3, 16, A STH A, ASM, *AR2+ ST B, *AR3+ IlLD *AR2, A SLB 峨3, 16, A, B ADD *AR3, 16, A STH A, ASM, *AR2+0 ST B, *