3-4离散傅立叶变换及其应用-2.ppt

4 快速傅氏变换算法FFT FFT是Fast Fourier Transform的缩写，换句话说，FFT是实现 DFT(Discret Fourier Transform)的一种快速运算手段。FFT算法可分为时间抽选法和频率抽选法。 一、直接计算DFT的问题及改进途径 运算量 FFT算法分类: 时间抽选法 DIT: Decimation-In-Time 频率抽选法 DIF: Decimation-In-Frequency 二 时间抽选法FFT 时间抽选法，就是在时域内逐次将序列分解成奇数子序列和偶数子序列，通过求子序列的DFT而实现整个序列的DFT，将计算DFT的运算量从N2复乘减少到(N2)log2N次复乘。 按时间抽选的基-2FFT算法 1、算法原理 设序列点数 N = 2L，L 为整数。 若不满足，则补零 则x(n)的DFT: 再利用周期性求X(k)的后半部分 三、蝶形结 蝶形结描述的另一种方法 分解后的运算量： 这样逐级分解，直到2点DFT 当N = 8时，即分解到X3(k)，X4(k)，X5(k)，X6(k)，k = 0, 1 四、FFT算法中一些概念 （1）“级”概念 将N 点DFT先分成两个N/2点DFT，再是四个N/4点DFT…直至N/2个两点DFT.每分一次称为“一”级运算。 因为N=2M所以N点DFT可分成M级 如上图所示依次m=0,m=1….M-1共M级 （2）“组”概念 (3) 因子的分布 2、运算量 当N = 2L时，共有L级蝶形，每级N / 2个蝶形，每个蝶形有1次复数乘法2次复数加法。 3、算法特点 1）原位计算 2）倒位序 3）蝶形运算 对N = 2L点FFT，输入倒位序，输出自然序， 第m级运算每个蝶形的两节点距离为 2m–1 第m级运算： 蝶形运算两节点的第一个节点为k值，表示成L位二进制数，左移L – m位，把右边空出的位置补零，结果为r的二进制数。 4）存储单元 输入序列x(n) : N个存储单元 4、DIT算法的其他形式流图 输入倒位序输出自然序 输入自然序输出倒位序 输入输出均自然序 相同几何形状 输入倒位序输出自然序 输入自然序输出倒位序 三 、按频率抽选的基-2FFT算法 1、算法原理 按k的奇偶将X(k)分成两部分： 令 N /2仍为偶数，进一步分解：N /2 N /4 逐级分解，直到2点DFT 设序列点数N=2L，L为整数。 将X(k)按k的奇偶分组前，先将输入x(n)按n的顺序分成前后两半： 则X(2r)和X(2r+1)分别是x1(n)和x2(n)的 N / 2点DFT，记为X1(k)和X2(k) x1(0) x1(1) -1 x1(2) x1(3) -1 x2(0) x2(1) -1 x2(2) x2(3) -1 N/2点 DFT N/2点 DFT x(0) x(7) x(1) x(2) x(3) x(4) x(5) x(6) X1(0)=X(0) X2(0)=X(1) X1(1)=X(2) X1(2)=X(4) X1(3)=X(6) X2(1)=X(3) X2(2)=X(5) X2(3)=X(7) * N (N – 1) N 2 N个X(k) (N点DFT) N – 1 N 一个X(k) 复数加法 复数乘法 2N+2 (N – 1)=2 (2N – 1) 4N 一个X (k) 2N (2N – 1) 4N 2 N个X (k) (N点DFT) 2 一次复加 2 4 一次复乘 实数加法 实数乘法 将序列x(n)按n的奇偶分成两组： N为2的整数幂的FFT算法称基-2FFT算法。 即蝶式计算结构也即为蝶式信号流图 上面频域中前/后半部分表示式可以用蝶形信号流图表示。 作图要素： (1)左边两路为输入 (2)右边两路为输出 (3)中间以一个小圆表示加、减运算（右上路为相加输出、右下路为相减输出) (4)如果在某一支路上信号需要进行相乘运算，则在该支路上标以箭头，将相乘的系数标在箭头旁。 (5)当支路上没有箭头及系数时，则该支路的传输比为1。 X1(k) X2(k) X1(k) X2(k) N (N / 2 –1) N 2 / 2 两个N / 2点DFT N N / 2 N / 2个蝶形 总计 2 1 一个蝶形 N / 2 (N / 2 –1) (N / 2)2 一个N / 2点DFT 复数加法 复数乘法 运算量减少了近一半 N / 2仍为偶数，进一步分解：N / 2 N / 4 同理： 其中： 每一级都有N/2个蝶形单元，例如：N=8,则每级都有4个蝶形单元。每一级的N/2个蝶形单元可以分成若干组，每一组具有相同的结构，相同的 因子分布，第m级的组数为： 例：N=8