第4章fft算法.pptVIP

第四章 快速傅里叶变换FFT 问题的提出 解决问题的思路与方法（减少运算量的途径） 基2按时间抽取FFT算法 基2按频率抽取FFT算法 FFT算法的实际应用 数字信号处理中DFT运算的用处 问题的提出 减少运算量的途径 FFT 算法简介 按时间抽取(DIT)的基-2FFT算法 基-2FFT算法基本上分为两大类： 按时间抽取法FFT（DIT-FFT） 按频率抽取法FFT（DIF-FFT） 则 x(n) 的DFT为 其中 X1(k) 和 X2(k) 分别为 x1(r) 和 x2(r) 的 N/2 点DFT，即 与第一次分解相似，将 x1(r) 按奇、偶序号分解成两个 N/4 点子序列 x3(s) 和 x4(s) ，即 用同样的方法可以分解x2(n)的DFT，即 DIT-FFT的运算规律： 第二级运算仍可采用这种原位的方式，但是进入蝶形结的 组合关系不同，首先进入蝶形结的是A(0)、A(2)存储单元中的 数据，运算结果仍可送回A(0)、A(2)保存，然后进入蝶形结的 是A(1)、A(3)…，依此类推，每一级运算均可在原位进行，这 种位运算结构可节省存储单元，降低设备成本，还可节省找原 地址的时间。 3. 旋转因子的变化规律及蝶形节点距离 4. 序列的倒位序 DIT-FFT算法的输入序列的排序为“倒位序”，是很有规律的，输出序列为自然顺序。由于N=2L，所以顺序数可用L位二进制数(nL-1nL-2…n1n0)表示，则倒位序为： n0n1…nL-2nL-1 。以N=8为例： 4.4 按频率抽选的基-2FFT(DIF-FFT) DIF-FFT的运算规律： 3. 旋转因子的变化规律及蝶形节点距离 4. 序列的倒位序 DIF-FFT算法的输入序列的排序为自然顺序，输出序列为“倒位序” 。 计算IDFT的快速算法：（两种） 一、IFFT 二、直接调用FFT子程序计算IFFT 基-4 FFT算法 FFT 算法应用 实数序列的FFT 前面讨论的FFT算法都是复数运算，包括序列x(n)也认 为是复数，但大多数场合，信号是实数序列，任何实数都可 看成虚部为零的复数，例如求某实信号y(n)的复谱，可认为 是将实信号加上数值为零的虚部变成复信号(x(n)+ j0)，再用 FFT求其离散付里叶变换。这种作法很不经济,因为把实序列 变成复序列，存储器要增加一倍，且计算机运行时，即使虚 部为零,也要进行涉及虚部的运算,浪费了运算量。合理的解 决方法是利用复数序列FFT对实数序列进行有效计算。 N点实序列的FFT算法 利用N点复序列的FFT计算两个N点实序列FFT 利用FFT计算线性卷积——快速卷积法 重叠相加法——由分段卷积的各段相加构成总的卷积输出 重叠保留法： 数字信号处理的实现 1 采用计算机或微机 适用于各种应用场合，但不能实时处理 2 采用专用的信号处理器 只适用于特定场合，能做到实时处理 3 采用通用的信号处理器 具有适应数字信号处理算法基本运算的指令 具有适应信号处理数据结构的寻址机构 能充分利用算法的并行性，即具有可编程能力，灵活性大，适应性强，且可进行高速信号处理 当 k 取奇数(k =2r+1 , r =0,1 ,…, N/2-1)时 DIF-FFT 蝶形运算流图符号 -1 DIF-FFT 一次分解运算流图（N=8） -1 -1 -1 -1 DIF-FFT 分解运算流图（N=8） -1 -1 -1 -1 -1 -1 -1 -1 x3(0) x3(1) x4(0) x4(1) x5(0) x5(1) x6(0) x6(1) -1 -1 -1 -1 X(0) X(4) X(2) X(6) X(1) X(5) X(3) X(7) 1. 原位运算 N=2L点的FFT共进行L级运算，每级由N/2个蝶形运算组成。 2. DIF-FFT算法与直接计算DFT运算量的比较 复数乘次数为 复数加次数为 -1 第二级的蝶形系数为 （2 组），蝶形节 点的距离为 N /4 第一级的蝶形系数为 （1 组），蝶形节 点的距离为 N /2 第三级的蝶形系数为 (4 组) ，蝶形 节点的距离为 N /8 第L级的蝶形系数均为 （N/2 组），蝶形节点的 距离为 1 …