DSP技术-05电子通信10.pptVIP

一、FFT运算的DSP实现 时域 —— 频域 离散信号 周期谱 周期信号 离散谱 离散周期信号 离散周期谱 DFT的定义 DFT和FFT计算量对比 2、N点蝶2FFT x1[n] = x[2n] x2[n] = x[2n+1] 上式的最后一步是因为 我们已经将一个N点的DFT分解成为两个N/2点的DFT。 但是，X(k)有N点，但X1(k)和X2(k)都只有N/2点，因此，前面计算的只是X(k)的前一半项的结果。 对于后一半X(k)，有 这样，只要计算出（0，N/2-1）区间的X1(k)和X2(k)，也就可以很方便地计算整个（0，N-1）区间的全部X(k)，从而大大地节省了运算量。 将N点DFT分为两个N/2点DFT 将N/2点DFT分为两个N/4点DFT 2点DFT 8点蝶2FFT 按频率抽取（DIF） x1[n] = x[n] x2[n] = x[n+N/2] n=0，1，……，N/2-1 DIT与DIF对比 DIT输入是混序的，频域的输出是顺序的； DIF输入是顺序的，频域的输出是混序的。 DIT的复数乘法出现在加减之前； DIF的复数乘法出现在加减之后。 定点DSP计算1024点FFT的时间 DSP 时间（ms） TMS320C25 10.9 TMS320C6201 0.067 DSP56001 1.65 浮点DSP计算1024点FFT的时间 DSP 时间（ms） TMS320C30 3.87 TMS320C40 1.02 ADSP21060 0.46 DSP96001 0.6 ADSP21160 0.45?s 顺序、混序与位倒序 DIT与DIF总有一边是混序的； 所有DSP都提供了位倒序（bit reverse）寻址指令； 3、FFT运算的DSP实现 主要内容： 1 基二实数FFT运算的算法 2 FFT运算模拟信号的产生及输入 3 实序列FFT汇编源程序及链接命令文件 4 观察信号时域波形及其频谱 1. 基二实数FFT运算的算法 该算法主要分为以下四步： 第一步，输入数据的组合和位倒序 把输入序列作位倒序是为了在整个运算最后的输出中得到的序列是自然顺序。 第二步，N点复数FFT 在DATA数据处理缓冲器里进行N点复数FFT运算。 第三步，分离复数FFT的输出为奇部分和偶部分 分离FFT输出为RP、RM、IP和IM四个序列，即偶实数、奇实数、偶虚数和奇虚数四部分。 第四步，产生2N点的复数FFT输出序列 产生2N=256个点的复数输出，它与原始的256个点的实输入序列的DFT一致。 2. FFT运算模拟信号的产生及输入 产生模拟输入信号可以有两种方法： 第一种方法是用C语言程序产生mdata.inc文件，然后，在DSP汇编语言程序中通过.copy汇编命令将生成的数据文件mdata.inc复制到汇编程序中。 第二种方法是建立与输入信号对应的数据流文件mdata.dat，具体操作步骤见教材内容。 3. 实序列FFT汇编源程序及链接命令文件 【例】256点实序列FFT的DSP实现。 256点实序列FFT汇编源程序如下： ************************************* *Radix-2,DIT,Real-input FFT Program * * fft.asm * ************************************* .mmregs .global reset,start,sav_sin,sav_idx,sav_grp .def start,_c_int00 .data DATA .space 1024 .copy mdata1.inc ; mdata1.inc 为模拟输入信号数据，起始地址标号INPUT N .set 128 ;复数点数 LOGN .set 7 ;蝶形级数 sav_grp .usect tempv,3 ;