数字信号DSP第四章2.ppt

* 输入序列的重排 * 倒序与重排程序框图 倒序程序 模块 数据交换 J (N/2,N/4,N/8, …，2，1) * 3) 旋转因子的变化规律 N=2M 点DIT-FFT运算流图中，共有M级蝶形运算。 每级都有N/2个蝶形，每个蝶形都要乘以旋转因子： P：旋转因子指数 找到旋转因子与运算级数L之间的关系！！ 由于各级旋转因子和循环方式不同，为编写程序，则 * 旋转因子与运算级数 L (L=1,2, …,M)之间的关系： N=23=8时，各级旋转因子表示： L=1时， L=2时， L=3时， 第L级运算共有2 L-1个不同的旋转因子。 注意：编程时，将运算级数L作为循环变量。 * 第L级旋转因子 的指数 对N=2M的一般情况，第L级的旋转因子为 第L级有2L-1个不同的旋转因子，每个旋转因子分别对应着2M-L个蝶形。 * 4) 蝶形运算规律 其中， p=J·2 M-L； J=0,1,…，2 L-1-1； L=1,2,…，M L：第L级运算； XL (J) :第L级运算后的值 第L级蝶形运算的两个输入数据相距B=2 L-1 个点， 蝶形运算可表示为： XL-1(J)+X L-1(J+B) XL (J) XL-1(J) - X L-1(J+B) XL(J+B) * 实数运算的蝶形算法： 设 * 5) 编程思想及程序框图 先从输入端（第1级）开始，逐级进行，共进行M级运算。 在进行第L级运算时，依次求出B=2L-1个不同的旋转因 子，每求出一个旋转因子，就计算出它对应的所有2M-L 个蝶形。这样，可以利用三重循环程序实现DIT-FFT运 算，程序框图如下 * DIT―FFT运算和程序框图 共M 级运算 每级的旋转因子数目,以及 蝶形运算两输入数据的距离 同一旋转因子p对应间隔2L的2M-L个蝶形 * 对于任何流图，只要保持各节点所连的支路及传输系数不变，则不论节点位置怎么排列所得流图总是等效的，所得最后结果都是x(n)的DFT的正确结果，只是数据的提取和存放的次序不同而已。 将和x(4)水平相连的所有节点和x(1)水平相连的所有节点位置对调，再将和x(6)水平相连的所有节点与和x(3)水平相连的所有节点对调，其余诸节点保持不变。 4.2.4 DIT算法的其他形式流图 * 输入顺序、 输出倒序的FFT流图 特点：旋转因子规律性强，适用于大型数据处理系统 * ① 数据存放的方式不同； ②取用系数的顺序不同。 经过简单变换，也可得输入与输出都是按自然顺序排列的流图以及其他各种形式的流图 。 * 输入、输出均为顺序的FFT流图 特点：输入和输出均为顺序排列，适用于硬件处理系统。 谢谢大家！ 感谢您的观看！ * 利用同一组存储单元存储输入、输出数据的方法，称为原位（址）计算 Company LOGO * 4.2 时域抽取基-2FFT(DIT-FFT) 将序列x(n)按n的奇偶分成两组： 序列x(n)的N点DFT 4.2.1 算法原理 * 可约性 * 这样就把 N点DFT化成两个N/2点DFT的组合，其中 * 后半部分 前半部分 蝶形运算 旋转因子 * 蝶形运算流图 每个蝶形运算量：1次复乘，2次复加 * 每个碟形： 2输入2输出 * 时域抽取将一个N点DFT分解为2个N/2点DFT(N=8) N/2个碟形 * 经一次分解后的运算量： 复数乘法 复数加法 一个N / 2点DFT (N / 2)2 N / 2 (N / 2 –1) 两个N / 2点DFT N 2 / 2 N (N / 2 –1) 一个蝶形 1 2 N / 2个蝶形 N / 2 N 总计 经一次分解，运算量减少了近一半 * N / 2仍为偶数，进一步分解：N / 2 → N / 4 * * 由两个N/4点DFT组合成N/2点DFT * 这样逐级分解，直到2点DFT，当N = 8时，即分解到 X3(k)，X4(k)，X5(k)，X6(k)，k = 0, 1。 * * 时域抽取将一个N点DFT分解为4个N/4点DFT(N=8) * 这种方法的每一步分解，都是按输入序列在时间上 的次序是属于偶数还是属于奇数来分解为两个更短的 子序列，所以称为“按时间抽取法”。 经过逐次分解，分解到计算2点DFT（2点DFT计算 不需要乘法运算）。最后将N点DFT运算全部转化为 蝶形运算，对于N=2M点的DFT，共需要M级蝶形运算。 * N点DIT-FFT运算流图(N=8) 第一级蝶形 第二