线卷积的fft算法.pptxVIP

线卷积的FFT算法第1页/共70页点击进入§4-1 引言§4-2按时间抽取(DIT)的FFT算法§4-3 DIF的FFT算法目录§4-4 IFFT算法§4-5线性卷积的FFT算法第2页/共70页§4-1引言一.DFT的计算工作量 两者的差别仅在指数的符号和因子1/N. 第3页/共70页一个X(k)的值的工作量,如X(1) 通常x(n)和 都是复数,所以计算一个 X(k)的值需要N次复数乘法运算,和 次 复数加法运算.那么,所有的X(k)就要N2次复 数乘法运算,N(N-1)次复数加法运算.当N很 大时,运算量将是惊人的,如N=1024,则要完 成1048576 次(一百多万次)运算.这样,难以做到实时处理.第4页/共70页二.改进的途径 1.的对称性和周期性对称性:周期性:得:第5页/共70页 利用上述特性,可以将有些项合并,并将DFT分解为短序列,从而降低运算次数,提高运算速度.1965年,库利(cooley)和图基(Tukey)首先提出FFT算法.对于N点DFT,仅需(N/2)log2N 次复数乘法运算.例如N=1024=210 时，需要(1024/2)log2 210 =512*10=5120次。5120/1048576=4.88% ,速度提高20倍第6页/共70页§4-2 按时间抽取(DIT)的FFT算法 —库利-图基算法一.算法原理(基2FFT)(一)N/2点DFT1.先将 按n的奇偶分为两组作DFT,设N=2L ,不足时,可补些零。这样有: n为偶数时: n为奇数时:因此，第7页/共70页(n为偶数) (n为奇数)由于: 所以,上式可表示为: 第8页/共70页其中,2.两点结论: (1) X (k),X (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出 X(k)的k=个；即前一半的结果。1 2k1 2N第9页/共70页3.X(k)的后一半的确定由于 （周期性）,所以: 同理, 这就是说,X1(k),X2(k)的后一半,分别 等于其前一半的值。第10页/共70页又由于 ,所以 可见,X(k)的后一半，也完全由X1(k), X2 (k)的前一半所确定。 *N点的DFT可由两个N/2点的DFT来计算。第11页/共70页4.蝶形运算由X1(k)、X 2(k)表示X(k)的运算是一种特殊的运算-碟形运算 前一半 后一半实现上式运算的流图称作蝶形运算（N/2个蝶形)1 1(前一半)11-1(后一半)第12页/共70页5.计算工作量分析按奇、偶分组后的计算量：(1)N/2点的DFT运算量:复乘次数: 复加次数:(2)两个N/2点的DFT运算量:复乘次数: 复加次数: (3)N/2个蝶形运算的运算量:复乘次数: 复加次数: 复乘:总共运算量:复加: *但是,N点DFT的复乘为N2 ;复加N(N-1);与分解后相比可知,计算工作点差不多减少 一半。第13页/共70页 例如 N=8 时的DFT,可以分解为两个 N/2=4点的DFT.具体方法如下: (1)n为偶数时,即 分别记作: 第14页/共70页 (2) n为奇数时,分别记作:第15页/共70页(3)对X (k)和X (k)进行蝶形运算,前半部为 X(0) X(3),后半部分为X(4) X(7) 整个过程如下图所示:1 2~~X1(0)X(0) x1(0)=x(0) x1(1)=x(2) N/2点 x1(2)=x(4) DFT x1(3)=x(6) x2(0)=x(1) x2(1)=x(3) N/2点 x2(2)=x(5) DFT x2(3)=x(7) X1(1)X(1)X1(2)X(2)X1(3)X(3)X2(0)X(4)0-1WX2(1)NX(5)1-1WX2(2)NX(6)2W-1NX2(3)X(7)-13WN第16页/共70页(二) N/4点DFT 由于N=2 L ,所以 N/2仍为偶数,可以进 一步把每个N/2点的序列再按其奇偶部分 分解为两个N/4的子序列。例如，n为偶数时的 N/2点分解为:进行N/4点的DFT,得到(偶中偶)(偶中奇)第17页/共70页从而可得到前N/4的X1(k)后N/4的X1(k)为第18页/共70页同样对n为奇数时，N/2点分为两个N/4点 的序列进行DFT，则有:(奇中偶)(奇中奇)第19页/共70页 例如,N=8时的DFT可分解为四个N/4的DFT, 具体步骤如下:(1) 将原序列x(n)的“偶中偶”部分:构成N/4点DFT，从而得到X3(0), X3(1)。第20页/共70页