4.2基2FFT算法-数字信号处理.ppt

图4.2.8 倒序规律 4.2 基2FFT算法 图4.2.9 倒序程序框图 4.2.5 频域抽取法FFT(DIF―FFT) 4.2 基2FFT算法 设序列x(n)长度为N=2M，首先将x(n)前后对半分开，得到两个子序列，其DFT可表示为如下形式： 偶数 奇数 将X(k)分解成偶数组与奇数组，当k取偶数 (k=2r,r=0,1,…,N/2-1)时 (4.2.14) 4.2 基2FFT算法 当k取奇数(k=2r+1,r=0,1,…,N/2-1)时 (4.2.15) 将x1(n)和x2(n)分别代入(4.2.14)和(4.2.15)式，可得 (4.2.16) 4.2 基2FFT算法 图4.2.10 DIF―FFT蝶形运算流图符号 4.2 基2FFT算法 图4.2.11 DIF―FFT一次分解运算流图(N=8) 4.2 基2FFT算法 图4.2.12 DIF―FFT二次分解运算流图(N=8) 4.2 基2FFT算法 图4.2.13 DIF―FFT运算流图(N=8) 4.2 基2FFT算法 图4.2.14 DIF―FFT的一种变形运算流图 4.2 基2FFT算法 图4.2.15 DIF―FFT的一种变形运算流图 4.2 基2FFT算法 上述FFT算法流图也可以用于离散傅里叶逆变换(Inverse Discrete Fourier Transform，简称IDFT)。比较DFT和IDFT的运算公式： IDFT的高效算法 4.2 基2FFT算法 如果希望直接调用FFT子程序计算IFFT，则可用下面的方法： 对上式两边同时取共轭，得 由于 4.2 基2FFT算法 4.3 进一步减少运算量的措施 充分利用无关紧要的旋转因子 ，减少复数乘法运算。如： 多类蝶形单元运算 可以使复数乘法次数降至 4.3 进一步减少运算量的措施 充分利用一些特殊的复数运算，减少复数乘法运算。如： 从实数运算考虑，计算N=2M点DIT―FFT所需实数乘法次数为 (4.3.4) 4.3 进一步减少运算量的措施 在FFT运算中，旋转因子WmN=cos(2πm/N)-jsin(2πm/N)，求正弦和余弦函数值的计算量是很大的。直接计算效率低。 可以采用预先建立旋转因子表，提高运算速度，但不足之处是占用内存较多。 4.3 进一步减少运算量的措施 旋转因子的生成 设x(n)为N点实序列，取x(n)的偶数点和奇数点分别作为新构造序列y(n)的实部和虚部，即 对y(n)进行N/2点FFT，输出Y(k)，则 根据DIT―FFT的思想及式(4.2.7)和(4.2.8)，可得到 实序列的FFT算法 4.3 进一步减少运算量的措施 由于x(n)为实序列，所以X(k)具有共轭对称性，X(k)的另外N/2点的值为 4.3 进一步减少运算量的措施 4.4 其他快速算法简介 分裂基FFT算法 离散哈莱特变换(DHT) 基4FFT 基8FFT 第四章快速付里叶变换(FFT)Fast Fourier Transforming 快速付里叶变换FFT 有限长序列通过离散傅里叶变换 (DFT)将其频域离散化成有限长序列。但其计算量太大(与N的平方成正比），很难实时地处理问题，因此引出了快速傅里叶变 换(FFT)。 FFT并不是一种新的变换形式，它只是DFT的一种快速算法，并且根据对序列分解与选取方法的不同而产生了FFT的多种算法。 FFT在离散傅里叶反变换、线性卷积和线性相关等方面也有重要应用。 4.1 引言 1965年库利(Cooley J.W)和图基(J.W.Turkey)在《Mathematic of Computation》 杂志上发表了著名的“机器计算付里级数的一种算法”文章，提出一种快速计算DFT的方法和计算机程序--揭开了FFT发展史上的第一页。 1984年法国的杜哈梅尔(P. Dohamel) 和霍尔曼(H. Hollmann)提出分裂基算法，使运算效率进一步提高。 FFT产生故事 4.1 引言 本章主要内容 直接计算DFT算法存在的问题及改进途径。 多种DFT算法(时间抽取算法DIT算法，频率抽取算法DIF算法，线性调频Z变换即CZT法) FFT的应用 4.1 引言 问题提出： 设有限长序列x(n),非零值长度为N,计算对x(n)进行一次DFT运算，共需多大的运算工作量? DFT计算存在的问题及改进途径 4.2 基2FFT算法 4.2 基2FFT算法 有限长序列x(n)的N点DFT为 可得，计算X(k)的所有N个值，共需