实验报告1FFT及DSP实现.docVIP

实验报告 实验名称 FFT的DSP实现 课程名称 高速DSP原理与应用 任课老师 刘建国 姓名 郭乾 学号 2009300817 班级 3491 日期 2012年6月日 一、实验目的 1、进一步熟悉C28X的汇编指令 2、加深对FFT算法的理解，体会利用的周期性和对称性来减少DFT的运算次数的原理。 二、实验要求 用汇编语言实现FFT算法，并画流程框图。 三、实验原理 设序列x(n)的长度为N，且满足,M为自然数。 利用旋转因子周期性，可以得到快速算法（FFT）。点的FFT共进行M级运算，每级由N/2个典型运算组成。同一级中，每一个蝶形运算的两个输入数据只对计算本蝶形运算有用，而且每个蝶形运算的输入、输出数据结点又同在一条水平线上，这就意味着计算完一个蝶形运算后，所得输出数据可立即存入原输入数据所占用的存储单元。这样，经过M级运算后，原来存放输入序列数据的N个存储单元（数组A）中便依次存放X（k）的N个值。N点DIT-FFT运算流图中，每级有N/2个蝶形。每个蝶形都要乘以因子，称其为旋转因子，p为旋转因子的指数。但各级的旋转因子和循环方式都有所不同。 N点序列的DFT和IDFT变换定义式如下： 编程思想: 第L级中，每个蝶形的两个输入数据相距B=个点；每级有B个不同的旋转因子；同一旋转因子对应着间隔为点的个蝶形。总结上述运算规律，可采用下述运算方法。先从输入端开始，逐级进行，共进行M级运算。在进行第L级运算时，依次求出B个不同的旋转因子，每求出一个旋转因子，就计算完它对应的所有个蝶形。这样我们可以用如下的三重循环程序实现DIT-FFT运算。 DIT-FFT运算流程框图 另外，DIT-FFT算法运算流图的输出X（k）为自然顺序，但为了适应原位运算，其输入序列不是安x（n）的自然序列排列，这种经过M次偶奇抽选后的排列成为序列x（n）的倒序。因此，在运行M级蝶形之前应先对序列x（n）进行倒序。 四、实验环境 软件环境：CCS3.3 硬件环境：N 五、实验过程、数据记录、处理及结论 用matlab产生一个单频信号的1024点采样数据 程序段如下： clear all; clc t=1:1:1024; y=500*sin(t); y=fix(y); sprintf( .word %d\n,y) 然后作为FFT的输入 信号的时域波形 输入信号的频谱: FFT的结果 六、讨论 FFT（快速傅里叶变换）是《数字信号处理》课程的一个核心内容。有一次，数字信号处理课老师在课堂上布置了预习时域抽取算法FFT原理的作业。我通过画8点DFT的三次时域抽取分解运算流图，加深了自己对蝶形运算的理解。在这个基础上通过推导，我比较了用FFT和DFT计算长度 的序列的乘法和加法次数，直观的明白了DIF一FFT算法相对于DFT算法在计算速度上的优越性。 七．附录（原程序代码） 程序代码如下： .sect .sd var_x: .copy VAR.asm .data .global mystart,_c_int00 M: .byte 00h L: .byte 9h B: .word 00h J: .word 00h P1: .word 00h K: .word 00h SAVE_Y: .copy FFT_Input.asm YINZI: .copy YINZI.txt .text _c_int00: mystart: SETC OBJMODE C28ADDR .C28_AMODE N .set 1024 YZ .set 512 MOVL XAR2,#FirstIn MOVL XAR3,#VAR_FIRST MOV @AR0,#N MOV AR1,#N-1 NOP*,ARP2 SETC AMODE .LP_AMODE LOOP: MOV AL,*++,ARP3 MOV *BR0++,AL,ARP1 XBANZ LOOP,*--,ARP2 LOOP1: MOV @@L,#1 ;add ar1,#1 MOV AL,#2 ;;b=2l-1 MOV T,@@L-1 LSL AL,T MOV @@B,AL LOOP2: MOV @@J,#0 MOV AL,@@M SUB AL,@@L MOV T,AL MOV AL,#2 LSL AL,T MOV T,AL MPY ACC,T,@@J MO