FFT的计算机实现.docxVIP

快速傅立叶变换(FFT )的计算机实现 邓凯电气0706 1.FFT运算的原理 DFT运算过程中如果有限长序列的长度很长时，即 N很大时，在运算过程 中所做的乘法和加法运算将很多。即便是采用高速的计算机进行运算，也很难达 到信息实时处理的要求。由库利、图基等科学家提出的快速傅里叶变换 FFT算 法大大减少了运算过程中的乘法和加法次数， 适合信息处理对实时性的要求，从 而得到广泛应用。 按时间抽取的FFT算法的基本思想是将输入的有限长序列首先分成奇数序 列和偶数序列，分别计算出奇、偶序列的 DFT，然后根据DFT的周期性和对称 性质，将其化简，接着将已分成的奇、偶序列再次分别划分成奇、偶序列，即前 面的奇序列按其长度再划分成奇数序列和偶数序列；前面的偶序列按其长度再划 分成奇数序列和偶数序列。分别计算其 DFT后，再按上述方法进行化简。如此 反复，直至被划分成的奇、偶序列长度为 1为止。 1.1基二FFT算法原理 若设X(k)=DFT[x(n)]，且0W n, k N-1，N为偶数(如果N为奇数，则添上 一个零值点使长度N为偶数)。把它分为奇序列和偶序列： x1(i) x(2i) x2(i) x(2i 1) (0 i — 1) 2 又 X1 (k) DFT[x 1 (i)] X 2 (k) DFT[x 2 (i)] (0 i N 1) 2 X(k) X1(k) X2(k)WNk 则有X(kN) X1(k) X 则有 X(k N) X1(k) X2(k)WN (0 i 1) 其中，X1[k] X2[k]分别是 xi[i]和 X2[i]的-点 DFT 2 即要求X[k]的值仅仅只需要求X1[k] X2[k]在(0 i - 1)部分的值即可 2 这样，我们就可以将Xn[k] —直分为奇序列和偶序列来求其值，直到奇偶序列长 度为1为止。 1.2蝶形图 对于FFT运算，我们通常用蝶形图来表示 X(k) Xi(k) X2(k)WNk 比如 n X(k —) Xi( k) X2(k)WN 我们表示为 X(k) Xi(k) X2 X(k) Xi(k) X2(k)WN 图 1.21 X(k —) Xi(k) X2(k)血 2 下面以N=8为例，运用FFT算法原理计算DFT，如图1.22、图1.23、图 1.24所示。 图1.22按时间抽取，将一个N=8点DFT分解为两个4点DFT