FFT反变换的DSP实现.docVIP

数字信号处理算法的DSP实现 课程报告 题目：FFT反变换的DSP实现 班 级：信息工程1班 姓 名：张庆 学 号：200720101127 日 期：2010.12.25 一、FFT反变换相关理论 对于离散傅里叶变换(DFT)的数字计算，FFT是一种有效的方法。一般假定输入序列是复数。当实际输入是实数时，利用对称性质可以使计算DFT非常有效。 一个优化的实数FFT算法是一个组合以后的算法。原始的2N个点的实输入序列组合成一个N点的复序列，之后对复序列进行N点的FFT运算，最后再由N点的复数输出拆散成2N点的复数序列，这2 N点的复数序列与原始的2N点的实数输入序列的DFT输出一致。 使用这种方法，在组合输入和拆散输出的操作中，FFT运算量减半。这样利用实数FFT算法来计算实输入序列的DFT的速度几乎是一般FFT算法的两倍。 二、FFT反变换DSP实现原理 实数FFT运算序列的存储分配 如何利用有限的DSP系统资源，合理的安排好算法使用的存储器是一个比较重要的问题。本文中，程序代码安排在0x3000开始的存储器中，其中0x3000~0x3080存放中断向量表，FFT程序使用的正弦表和余弦表数据(.data段)安排在0xc00开始的地方，变量(.bss段定义) 存放在0x80开始的地址中。另外，本文中256点实数FFT程序的数据缓冲位0x2300~0x23ff , FFT变换后功率谱的计算结果存放在0x2200~0x22ff中。 三、FFT反变换DSP实现程序说明 连续定位.cmd文件程序如下： MEMORY { PAGE 0: PRGROM : origin=00080h length=00200h VECT: origin = 0xff80, len = 0x80 PAGE 1: S_PRAM : origin=00060h length=00020h INTRAM1: origin=00400h length=00200h INTRAM2: origin=00800h length=00200h INTRAM3: origin=00c00h length=01000h EXTRAM1: origin=01c00h length=00400h EXTRAM2: origin=02000h length=010h EXTRAM3: o= 2100h l=500h } SECTIONS { vectors : {} VECT PAGE 0 .text : {} PRGROM PAGE 0 .data : {} INTRAM3 PAGE 1 power : {} EXTRAM1 PAGE 1 sin_tbl : {} INTRAM1 PAGE 1 cos_tbl : {} INTRAM2 PAGE 1 fft_vars: {} S_PRAM PAGE 1 stack : {} EXTRAM2 PAGE 1 .stack : {} EXTRAM3 PAGE 1 .bss : {} EXTRAM3 PAGE 1 } 2. 基2实数FFT运算的算法，该算法主要分为以下四步进行： 输入数据的组合和位排序 实现数据位倒序存储的具体程序段如下： .mmregs .include fft_size.inc .def bit_rev .ref _real_fft_input, fft_data .asg AR2,REORDERED_DATA .asg AR3,ORIGINAL_INPUT .asg AR7,DATA_PROC_BUF .text bit_rev: SSBX FRCT STM #_real_