9FIR设计.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * DSP原理及应用 * 第6章 应用程序设计 2．快速傅氏变换FFT 由于WN2=[e-j(2?/N)]2=[e-j2?/(N/2)]=WN/2，则有： k = 0,1,…,N/2-1 (6.3.3) X1(k)和X2(k)分别为x1(n)和x2(n)的N/2点的DFT。 由于对称性，WNk+N/2 = -WNk， 则X(k+N/2)=X1(k) - WNkX2(k)。 N点X(k)可分为两部分： 前半部分：X(k)=X1(k) + WNkX2(k) k=0, 1, … N/2-1 (6.3.4) 后半部分：X(k+N/2)=X1(k)-WNkX2(k) k=0, 1, …N/2-1 (6.3.5) * DSP原理及应用 * 第6章 应用程序设计 2．快速傅氏变换FFT 基2 DIT FFT的蝶形运算： WNk -1 xm-1(p) xm-1(q) xm(p) xm(q) 蝶形算法： xm(p)= xm-1(p)+ xm-1(q)WNk (6.3.6) xm(q)= xm-1(p) - xm-1(q)WNk (6.3.7) * DSP原理及应用 * 第6章 应用程序设计 2．快速傅氏变换FFT 在基数为2的FFT中，设N=2M，共有M级运算，每级有N/2个2点FFT蝶形运算，因此，N点FFT总共有(N/2)log2N个蝶形运算。 例如：在N点FFT中，当N=8时，共需要3级，12个基2 DIT FFT的蝶形运算。 * DSP原理及应用 * 第6章 应用程序设计 2．快速傅氏变换FFT 信号流程如图： x(0) x(4) x(2) x(6) x(1) x(5) x(3) x(7) WN0 x(7) x(6) x(5) x(4) x(3) x(2) x(1) x(0) -1 -1 -1 WN0 WN0 WN0 -1 -1 WN0 WN2 WN0 WN2 -1 -1 -1 WN0 WN1 WN2 WN3 -1 -1 -1 -1 * DSP原理及应用 * 第6章 应用程序设计 2．快速傅氏变换FFT 输入是经过比特反转的倒位序列——位码倒置。 排列顺序：x(0)、x(4)、x(2)、x(6)、x(1)、x(5)、x(3)、x(7) 输出是按自然顺序排列。 排列顺序：x(0)、x(1)、x(2)、x(3)、x(4)、x(5)、x(6)、x(7) * DSP原理及应用 * 第6章 应用程序设计 6.3 快速傅里叶变换(FFT)的DSP实现 下面以8~1024复数点FFT算法为例，介绍用DSP芯片实现的方法。 6.3.2 FFT算法的DSP实现 1．FFT运算序列的存储分配 FFT运算时间是衡量DSP芯片性能的一个重要指标,在用DSP芯片实现FFT算法时，合理地安排好所使用的存储空间是十分重要的。 * DSP原理及应用 * 第6章 应用程序设计 1．FFT运算序列的存储分配 (1) 程序空间的分配 程序存储器 ┇ sine1 E000 ┇ E1FF 正弦系数表 cosine1 E200 ┇ E3FF 余弦系数表 fft_prg E400┇ E4A2 程序代码 ┇ .vectors FF80┇ FFFF 复位向量和 中断向量表 E000h~E1FFh：正弦系数表 E200h~E3FFh：余弦系数表 E400h~E4A2h：程序代码 FF80h~FFFFh：复位向量和 中断向量表 * DSP原理及应用 * 第6章 应用程序设计 1．FFT运算序列的存储分配 (2) 数据空间的分配 数据存储器 0000 ┇ 005F 存储映像 寄存器 .bss 0060 0061 0062 暂存单元 stack 0063 ┇ 006C 堆 栈 ┇ sine 0400 ┇ 05FF 正弦系数表 数据存储器 ┇