[工学]1111数字信号处理第四章2011.ppt

数字信号处理Digital Signal Processing 3.一次复数运算换算成实数运算量 （1）利用DFT运算的系数 的固有特性，改善DFT的运算效率。 例：求 N=23=8点FFT变换（1）先按N=8--N/2=4，做4点的DFT： 将N=8点分解成2个4点的DFT的信号流图  蝶形流图 (4)一个完整N=8的按时间抽取FFT的运算流图  (2) 因子的分布 三、按时间抽取FFT算法的特点 1.原位运算（in-place) 2.码位倒读规则 定义：当数据输入到存储器以后，每一组运算的结果，仍然存放在这同一组存储器中直到最后输出。 原位运算的结构，可以节省存储单元，降低设备成本。 设输入序列长度为N=2M(M为正整数)，将该序列的频域的输出序列X(k)(也是M点序列)，按其频域序号的奇偶分解为越来越短的子序列，称为基2按频率抽取的FFT算法。也称为Sander-Tukey算法。 三、蝶形流图表示 （1）N=8--N/2=4： 将N=8点分解成2个4点的DIF的信号流图 （b) 蝶形流图 (4)一个完整N=8的按频率抽取FFT的运算流图 （1）DIF与DIT两种算法结构倒过来。 DIF为输入顺序，输出倒序。最后进行变址运算 DIT为输入倒序，输出顺序。先进行变址运算。 （2）DIF与DIT根本区别：在于蝶形结不同。 DIT的蝶形为先复数相乘后加减。 DIF的蝶形为先加减后复数相乘。 以上所讨论的FFT的运算方法同样可用于IDFT的运算，简称为IFFT。即快速傅里叶反变换。 利用FFT计算IFFT的思路 将下列两式进行比较: 前面的IFFT算法，要改变FFT的程序和参数才能实现。 另一种IFFT算法，则可以完全不必改动FFT程序。 DFT变换： 频域上X(k)按k的奇偶分为两组，在时域上x(n)按n的顺序分前后两部分，现将输入x(n)按n的顺序分前后两部分: 前半子序列 x(n), 0≤n≤N/2-1; 后半子序列 x(n+N/2), 0≤n≤N/2-1; 例：N=8时，前半序列为：x(0),x(1),x(2),x(3); 后半序列为： x(4),x(5),x(6),x(7); (1).分组 二、算法步骤 (2).代入DFT中 （3）按 k 的奇偶将 X(k) 分成两部分，进行变量置换 当 k 为偶数时： 当 k 为奇数时： 即： 蝶形单元： x(n) x(n+N/2) 与时间抽取法的推演过程一样，由于N=2M, N/2仍为偶数，所以可以将N/2点DFT的输出X(k)再分为偶数组和奇数组，这样就将一个N/2点的DFT分成两个N/4点DFT的输入，也是将N/2点的DFT的输入上、下对半分后通过蝶形运算而形成，直至最后为2点DFT。 先将N=8 DFT分解成2个4点DFT： 可知：时域上：x(0),x(1),x(2),x(3)为前半序列 x(4),x(5),x(6),x(7)为后半序列 产生新的子序列： x1(n)、 x2(n) 频域上：X(0),X(2),X(4),X(6)由x1(n)的DFT给出 X(1),X(3),X(5),X(7),由x2(n)的DFT给出 例子：求 N=23=8点DIF 4点 DFT x(0) x(1) x(2) x(3) 4点 DFT x(4) x(5) x(6) x(7) X(0) X(2) X(4) X(6) X(1) X(3) X(5) X(7) X1(k) 前半部分序列 后半部分序列 x1(n) x2(n) X2(k) 因为4点DIF还是比较麻烦，所以再继续分解。 (2)N/2(4点)--N/4(2点)FFT 2点DFT 2点DFT x1(0) x1(1) x1(2) x1(3) x3(0) x3(1) x4(0) x4(1) X(0) X(4) X(2) X(6) 2点DFT 2点DFT x2(0) x2(1) x2(2) x2(3) x5(0) x5(1) x6(0) x6(1) X(1) X(5) X(3) X(7) 两点DFT可用一个蝶形结表示。取x3(0)、x3(1)为例。 (3) 2点DFT 蝶形流图 X(0) X(4) x3(0) x3(1) x(0) x(1) x(2) x(3) x(4) x(5) x(6) x(7) X(0) X(4) X(2) X(6) X(1) X(5) X(3) X(7) L=1 L=2 L=3 X(0) =1 X(1) =2 X