专题八之三-通用数字信号处理方法的DSP实现v2013.ppt

实数FFT运算序列的存储分配 基二实数FFT运算的算法 第一步，输入数据的组合和位倒序：（1）把输入序列作位倒序，是为了在整个运算最后的输出中得到的序列是自然顺序。（2）原始的输入的2N = 256个点的实数序列复制放到标记有“d_input_addr”的相邻单元，当成N = 128点的复数序列d[n]。其中，奇数地址是d[n]的实部，偶数地址是d[n]的虚部。复数序列经过位倒序，存储在数据处理缓冲器中 “fft-data”。 bit_rev: STM #d_input_addr,ar3 ；在AR3中放入输入地址 STM #fft_data,ar7 ；在AR7中放入处理后输出的地址 MVMM DATA_PROC_BUF,ar2 ; AR2中装入第一个位倒序 ; 数据指针 STM #K_FFT_SIZE-1,BRC STM #K_FFT_SIZE,ar0 ; AR0=输入数据数目的一半(128) RPTB bit_rev_end MVDD *ar3+,*ar7+ ; 将原始输入缓冲中的数据放入到 ; 位倒序缓冲中去之后输入缓冲 ; (AR3)指针加1, 位倒序缓冲(AR2)指 ; 针也加1 MVDD *ar3-,*ar7+ ; 将原始输入缓冲中的数据放入到 ; 位倒序缓冲中去之后输入缓冲(AR3) ; 指针减1以保证位倒序寻址正确 MAR *ar3+0B ；按位倒序寻址方式修改AR3 bit_rev_end: 基二实数FFT运算的算法 第二步，N点复数FFT：（1）在数据处理缓冲器里进行N点复数FFT运算。由于在FFT运算中要用到旋转因子WN，它是一个复数。我们把它分为正弦和余弦部分，用Q15格式将它们存储在两个分离的表中。每个表中有128项，对应从0度到180度。因为采用循环寻址来对表寻址，128 = 27 28，因此每张表排队的开始地址就必须是8个LSB位为0的地址。 基二实数FFT运算的算法 我们把128点的复数FFT分为七级来算，第一级是计算两点的FFT蝶形结，第二级是计算四点的FFT蝶形结，然后是八点、十六点、三十二点六十四点、一百二十八点的蝶形结计算。最后所得的结果表示为： 其中，R[k]、I[k]分别是D[k]的实部和虚部。 D[k] = F{d[n]} = R[k] + j I[k] 这一步中，实现FFT计算的具体程序如下： ;stage1 ：计算FFT的第一步，两点的FFT .asg AR1,GROUP_COUNTER ；定义FFT计算的组指针 .asg AR2,PX ；定义AR2为指向参加蝶形运算第一个数据的指针 .asg AR3,QX ；定义AR3为指向参加蝶形运算第二个数据的指针 .asg AR4,WR ；定义AR4为指向余弦表的指针 .asg AR5,WI ；定义AR5为指向正弦表的指针 .asg AR6,BUTTERFLY_COUNTER ；定义AR6为指向蝶形结的指针 .asg AR7,DATA_PROC_BUF ；定义在第一步中的数据处理缓冲指针 .asg AR7,STAGE_COUNTER ；定义剩下几步中的数据处理