数字信号处理(中澳双语8-1)-第二十讲.ppt

FFT运算量与运算特点 1． N=2L时，共有L=log2N级运算；每一级有N/2个蝶形结。 2．每一级有N个数据中间数据），且每级只用到本级的转入中间数据，适合于迭代运算。 3．计算量： 每级N/2次复乘法，N次复加。（每蝶形只乘一次，加减各一次）。共有L*N/2=N/2log2N 次复乘法；复加法L*N=Nlog2N 次。与直接DFT定义式运算量相比(倍数) N2/(Nlog2N) 。当 N大时，此倍数很大。 郑州大学信息工程学院 * 郑州大学信息工程学院 * 比较DFT 可以直观看出，当点数N越大时，FFT的优点更突出。 DIT FFT中最主要的蝶形运算实现 （1）参与蝶形运算的两类结点（信号）间“距离”（码地址）与其所处的第几级蝶形有关；第m级的“结距离”为 (即原位计算迭代） （2）每级迭形结构为 郑州大学信息工程学院 * 郑州大学信息工程学院 * 蝶形运算两节点的第一个节点为k值，表示成L位二进制数，左移L – m位，把右边空出的位置补零，结果为r的二进制数。 （3） 的确定： 第m级的r取值： 四、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）“组”概念 郑州大学信息工程学院 * 每一级都有N/2个蝶形单元，例如：N=8,则每级都有4个蝶形单元。每一级的N/2个蝶形单元可以分成若干组，每一组具有相同的结构，相同的 因子分布，第m级的组数为： 例：N=8=23，分3级。 m=0级,分成四组,每组系数为 m=1级,分成二组,每组系数为 m=2级,分成一组,每组系数为 (3) 因子的分布 郑州大学信息工程学院 * 结论：每由后向前（m由M-1--0级）推进一级，则此系数为后级系数中偶数序号的那一半。 郑州大学信息工程学院 * 例 用FFT算法处理一幅N×N点的二维图像，如用每秒可做10万次复数乘法的计算机，当N=1024时，问需要多少时间（不考虑加法运算时间）？ 解 当N=1024点时，FFT算法处理一幅二维图像所需复数乘法约为 次，仅为直接计算DFT所需时间的10万分之一。 即原需要3000小时，现在只需要2 分钟。 按频率抽取（DIF）的FFT算法 与DIT-FFT算法类似分解，但是抽取的是X(k)。即分解X(k)成奇数与偶数序号的两个序列。 设： N = 2L，L 为整数。将X(k)按k的奇偶分组前，先将输入x(n)按n的顺序分成前后两半： * （Decimation In Frequency) 一、算法原理 郑州大学信息工程学院 * 下面讨论 郑州大学信息工程学院 * 按k的奇偶将X(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/2仍为偶数，进一步分解：N/2 → N/4 郑州大学信息工程学院 * x3(0) x3(1) -1 -1 x4(0) x4(1) N/4点 DFT N/4点 DFT x1(0) x1(1) x1(2) x1(3) X3(0)=X1(0)=X(0) X4(0)=X1(1)=X(2) X3(1)=X1(2)=X(4) X4(1)=X1(3)=X(6) 按照以上思路继续分解，即一个N/2的DFT分解成两个N/4点DFT，直到只计算2点的DFT，这就是DIF-FFT算法。 郑州大学信息工程学院 * 2个1点的DFT蝶形流图 进一步简化为蝶形流图： 1点DFT x3(0) 1点DFT x3(1) X(0) X(4) X(0) X(4) x3(0) x3(1) 郑州大学信息工程学院 * 二、按频率抽取FFT蝶形运算特点 1）原位计算 郑州大学信息工程学院 * -1 L级蝶形运算，每级N/2个蝶形，每个蝶形结构： m表示第m级迭代，k，j表示数据所在的行数 郑州大学信息工程学院 * 2）蝶形运算 对N=2L点FFT，输入自然序，输出倒