快速傅里叶变换全面.pptVIP

如上所述，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 ., * 对N=2M的一般情况，第L级的旋转因子为： ., * 设序列x(n)经时域抽选(倒序)后，存入数组X中。如果蝶形运算的两个输入数据相距B个点(B=2L-1)，应用原位计算，则蝶形运算可表示成如下形式： 下标L表示第L级运算，XL(J)则表示第L级运算后数组元素X(J)的值。 ., * 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) 输出结果 结束 ., * 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)的顺序存入存储单元，即为乱序输入，顺序输出。 这种顺序看起来相当杂乱，然而它是有规律的。即码位倒读规则。 ., * 自然顺序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 看出：码位倒读后的顺序刚好是数据送入计算机内的顺序。 ., * 倒序规律 ., * 对于数N，在其二进制最高位加1，等于加N/2。 若已知某个反序号为J，为求下一个反序号，可先判J的最高位： 1) 若为0，则把该位变成1（即加N/2）就得到下 一个反序号， 2) 若为1，则需判断次高位： ① 若次高位为0，则把最高位变0（相当减去 N/2）后，再把次高位变1（即加N/4）。 ② 若次高位为1，则需判断次次高位…… 分析： ., * 倒 序 排 列 算 法 的 流 程 图 正序序列已在数组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 ., * 第四节 按频率抽选的基2-FFT算法 在基2快速算法中，频域抽取法FFT也是一种常用的快速算法，简称DIF―FFT。 设序列x(n)长度为N=2M，首先将x(n)前后对半分开，得到两个子序列，其DFT可表示为如下形式 ., * 第一节 直接计算DFT的问题及改进途径 1、问题的提出 设有限长序列x(n)，非零值长度为N，若对x(n)进行一次DFT运算，共需多大的运算工作量？ 计算成本? 计算速度? ., * 2. DFT的运算量 回忆DFT和IDFT的变换式： 1）x(n)为复数， 也为复数。 2）DFT与IDFT的计算量相当。 注意： ., * 计算机运算时（编程实现）： N次复乘，N-1次复加 N个点 以DFT为例： ., * 复数乘法 复数加法 一个X(k) N N – 1 N个X(k) (N点DFT) N 2 N (N – 1) 实数乘法 实数加法 一次复乘 4 2 一次复加 2 一个X (k) 4N 2N+2 (N – 1)=2 (2N – 1) N个X (k) (N点DFT) 4N 2 2N (2N – 1) 运算量 (a+jb)(c+jd)=(ac-bd)+j(bc+ad) ., * 例：计算一个 N点DFT ，共需N2次复乘。以做一次 复乘1μs计，若N =4096，所需时间为 例：石油勘探，有24个通道的记录，每通道波形记 录长度为5秒，若每