10 快速付里叶变换FFT-DIT.pptVIP

快速傅立叶变换(FFT)Fast FourierTransforming 快速傅立叶变换FFT 0 引言 1 直接计算DFT存在的问题及改进途径 2 基--2按时间抽取的FFT算法 3 基--2按频率抽取的FFT算法 4 离散反傅立叶变换的快速计算方法 5 线性卷积的FFT算法 0 引 言 0.1 快速付立叶变换FFT DFT计算量太大, 很难实时地处理问题, 因此引出了快速傅里叶变换(FFT) 。 特别说明：FFT并不是一种新的变换形式,它只是DFT的一种快速算法.并且根据对序列分解与选取方法的不同而产生了FFT的多种算法。 FFT在离散傅里叶反变换、线性卷积和线性相 关等方面也有重要应用。 0.2 主要内容 直接计算DFT算法存在的问题及改进途径 多种FFT算法(时间抽取算法DIT算法，频率抽取算法DIF算法，线性调频Z变换即CZT法) FFT的应用 1.1 直接计算DFT计算量 问题提出： 设有限长序列x(n),非零值 长度为N，计算对x(n)进行一次DFT运 算，共需多大的运算工作量? 1.2 改善DFT运算效率的基本途径 两个途径： 1.合并法：合并DFT运算中的某些项。 2. 分解法：将长序列DFT利用对称性和周期性，分解为短序列DFT。 方法： 2 基--2按时间抽取的FFT算法 2.1 算法原理 按时间抽取----将该序列按时间顺序的奇偶分解为越来越短的子序列。 基数2----N=2M，M为正整数。若不满足这个条件，可以补零加长。 2.2 算法步骤 步骤2.代入DFT中 步骤3.求出子序列的DFT 2.3 蝶形运算 (2)N/2(4点)--N/4(2点)FFT 将x1(r)和x2(r)分解成奇、偶两个N/4点(2点)点的子序列。 (d)另一个2点的DFT蝶形流图 (b)2个1点的DFT蝶形流图 (4)完整的N=8的按时间抽取FFT的运算流图 2.4 FFT算法中一些概念 将N 点DFT先分成两个N/2点DFT，再是四个N/4点DFT…直至N/2个两点DFT。每分一次称为“一”级运算。 因为N=2M所以N点DFT可分成M级 （2）“组”的概念 (3)旋转因子的分布 （4）按时间抽取法 2.5 按时间抽取的FFT算法与直接计算DFT运算量的比较 每级都由N/2个蝶形单元构成，因此每一级运算都需要N/2次复乘和N次复加（每个结加减各一次）。这样（N=2M)M级运算共需要： 复乘次数： 复加次数： 2.6 按时间抽取FFT算法的特点 原位运算（in-place) 码位倒读规则 例：求 N=23=8点FFT变换 频域：X(0)~X(3)，由X(k)给出 X(4)~X(7), 由X(k+N/2)给出 时域：x(0),x(2),x(4),x(6)为偶子序列 x(1),x(3),x(5),x(7)为奇子序列 （1）先按N=8--N/2=4，做4点的DFT： 2.3 蝶形运算 将N=8点分解成2个4点的DFT的信号流图 X(4)~X(7) 同学们自已写 (a)先将4点分解成2点的DFT： (2)N/2(4点)--N/4(2点)FFT (b)求2点的DFT： (c)一个2点的DFT蝶形流图 (2)N/2(4点)--N/4(2点)FFT (2)N/2(4点)--N/4(2点)FFT (3)将N/4（2点）DFT再分解成2个1点的DFT 取x(0)、x(4)为例: (a)求2个一点的DFT 代入上面流图可知: 这是一个蝶形结 这里用到对称性： 进一步简化： (3)将N/4（2点）DFT再分解成2个1点的DFT （1）“级”的概念 每一级都有N/2个蝶形单元，例如：N=8,则每级都有4个蝶形单元。每一级的N/2个蝶形单元可以分成若干组，每一组具有相同的结构，相同的 因子分布，第m级的组数为： 例：N=8=23，分3级。 m=0级,分成四组,每组系数为 m=1级,分成二组,每组系数为 m=2级,分成一组,每组系数为 2.4 FFT算法中一些概念 结论：每由后向前（m由M-1--0级）推进一级，则此系数为后级系数中偶数序号的那一半。 看出：第m级的系数为 2.4 FFT算法中一些概念 2点DFT 2点DFT 2点DFT 2点DFT 2点DFT 2点DFT 2点DFT 2点DFT 两个2点 DFT 两个2点 DFT 两个2点 DFT 两个2点 DFT 两个 4点DFT 两个 4点DFT 两个 N/2点 DFT X1(k) …... X2(k) X(k) 由于每一步分解都是基于在每级按输入时间序列的次序是属于偶数还是奇数来分解为两个更短的序列，所以称为“