-第四章-快速傅立叶变换-Fast-Fourier-Transfo.ppt

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) 回顾：N点DIT―FFT运算流图(N=8) * bjhk 如上所述，N点DIT―FFT运算流图中，每级都有N/2个蝶形。每个蝶形都要乘以因子WNP，称其为旋转因子，p称为旋转因子的指数。 2)旋转因子的变化规律 观察FFT运算流图发现，第L级共有2L-1个不同的旋转因子。N=23=8时的各级旋转因子表示如下： L=1时，WNp=WN/4J， N/4 =21 =2L， J=0 L=2时， WNp =WN/2J， N/2 =22 =2L， J=0，1 L=3时， WNp =WNJ， N =23 =2L， J=0，1，2，3 * bjhk 对N=2M的一般情况，第L级的旋转因子为： * bjhk 设序列x(n)经时域抽选(倒序)后，存入数组X中。如果蝶形运算的两个输入数据相距B个点(B=2L-1)，应用原位计算，则蝶形运算可表示成如下形式： 下标L表示第L级运算，XL(J)则表示第L级运算后数组元素X(J)的值。 * bjhk 3) 编程思想及流程图 开始 送入x(n)和N=2M 调整输入x(n)的顺序 for(L=1; L=M; L++) B = 2L-1 for(J=0; J=B-1; J++) p = J·2M-L for(k = J; k= N-1; k=k+2L) 输出结果 结束 * bjhk 4）码位倒序 由N=8蝶形图看出：原位计算时，FFT输出的X(k)的次序正好是顺序排列的，即X(0)…X(7),但输入x(n)都不能按自然顺序存入到存储单元中，而是按x(0),x(4),x(2),x(6) ,x(1),x(5),x(3),x(7)的顺序存入存储单元，即为乱序输入，顺序输出。 这种顺序看起来相当杂乱，然而它是有规律的。即码位倒读规则。 * bjhk 自然顺序n 二进制码表示 码位倒读 码位倒置顺序n’ 以N=8为例： 0 1 2 3 4 5 6 7 000 001 010 011 100 101 110 111 000 100 010 110 001 101 011 111 0 4 2 6 1 5 3 7 看出：码位倒读后的顺序刚好是数据送入计算机内的顺序。 * bjhk 倒序规律 * bjhk 对于数N，在其二进制最高位加1，等于加N/2。 若已知某个反序号为J，为求下一个反序号，可先判J的最高位： 1) 若为0，则把该位变成1（即加N/2）就得到下 一个反序号， 2) 若为1，则需判断次高位： ① 若次高位为0，则把最高位变0（相当减去 N/2）后，再把次高位变1（即加N/4）。 ② 若次高位为1，则需判断次次高位…… 分析： * bjhk 倒 序 排 列 算 法 的 流 程 图 正序序列已在数组A[ ] 中，输 入 N LH= N/2 , j=LH , N1 = N-2 j=j-k k=k/2 k=LH jk j=j+k T = A(j) A(i) = A(j) A(j) = T for(i=1; i=N1; i++) i ≥ j N Y Y N * bjhk 第四节 按频率抽选的基2-FFT算法 在基2快速算法中，频域抽取法FFT也是一种常用的快速算法，简称DIF―FFT。 设序列x(n)长度为N=2M，首先将x(n)前后对半分开，得到两个子序列，其DFT可表示为如下形式 * bjhk * bjhk 第四章 快速傅立叶变换 Fast Fourier Transform * bjhk 第一节 直接计算DFT的问题及改进途径 1、问题的提出 设有限长序列x(n)，非零值长度为N，若对x(n)进行一次DFT运算，共需多大的运算工作量？ 计算成本? 计算速度? * bjhk 2. DFT的运算量 回忆DFT和IDFT的变换式： 1）x(n)为复数， 也为复数。 2）DFT与IDFT的计算量相当。 注意： * bjhk 计算机运算时（编程实现）： N次复乘，N-1次复加 N个点 以DFT为例： * bjhk 复数乘法 复数加法 一个X(k) N N – 1 N个X(k) (N点DFT) N 2 N (N – 1) 实数乘法 实数加法 一次复乘 4 2 一次复加 2