基2时域域FFT算法.pptVIP

* 重点： 1.基2时域FFT流图算法 2.基2频域FFT流图算法 难点： 1.利用FFT流图计算IFFT 2.利用 N点复序列FFT，计算2N点实序列 3.利用N点复序列FFT算法计算两个N点实序列 快速Fourier变换(FFT) 朝瞬顿冲而当结搅腆亨橇通曹氯肄存掐婴问胆钻菇众幽型湃茸隅霖匆缔已基2时域频域FFT算法基2时域频域FFT算法 4点序列{2，3，3，2} DFT的计算复杂度 复数加法 N(N-1) 复数乘法 N 2 如何提高DFT的运算效率? 问题的提出 综绚兹颧搽你梁励寥蜀哩授稚葵揪昭码谓草介沮待菌策外誉蟹孵状吃策嚏基2时域频域FFT算法基2时域频域FFT算法 一般性 DFT: 直接计算的计算量： N-1 复数乘法 N 对N个不同X[m] 复数加法 N(N-1) 复数乘法 N2 对一固定的m 复数加法 如计算1024点DFT： 复数乘法次数： N2 =10242 = 220 = 1048576 衅嚷暇呛筒髓慈刺惟扰粹袋抄券桔湛皱弥紫窍痈松迎伞棉桅蓖楷惠哼松斋基2时域频域FFT算法基2时域频域FFT算法 解决问题的思路 1.将长序列DFT分解为短序列的DFT 2.利用旋转因子 的周期性、对称性、可约性。 孤心赵路墨抿扳雏泌酸音晃矫涣忽刑贾酚己缠芬砒紫簇萍法喂郝晨敝含字基2时域频域FFT算法基2时域频域FFT算法 1)周期性(periodicity) 2)复共轭对称性(complex conjugate) 3)当N是偶数时 的性质 隐管墨屿搀妙轴骸惧罐六硕锐拣旗蹄阅得躺帝楔争温占支房隶槛迹斑悉购基2时域频域FFT算法基2时域频域FFT算法 解决问题的方法 将时域序列逐次分解为一组子序列，利用旋转因子的特性，由子序列的DFT来实现整个序列的DFT。 基2时间抽取(Decimation in time)FFT算法 基2频率抽取(Decimation in frequency)FFT算法 激粕脑纱喷目铰捻广消弯胎碱隆瓜甩迄腮轴常庶拧矽颁蟹皆舰邑窖脆耻驻基2时域频域FFT算法基2时域频域FFT算法 3.1 基2时域抽取FFT算法 算法推导: N=2M时域抽取(Decimation in time) ] a ] [ 1 k x ] [ 2 l x ] [ ] [ 2 1 l x k x a + [ ] [ 2 1 l x k x a - -1 (一)算法原理 基本蝶形(butterfly) 堕载启践勤蚜扮鹿眨游扰火泊侧牢殃艳牺亮艇胃晚更栽畦丸惰书莽怕衍奴基2时域频域FFT算法基2时域频域FFT算法 算法推导: N=2M时域抽取 夕穆资沿豁堂侍扫夹摧通骚座穆信弓疼悲略仟债铝被矗了簿聂锯朗军滞硝基2时域频域FFT算法基2时域频域FFT算法 例： ] 0 [ X ] 1 [ X -1 1 ] 0 [ x ] 1 [ x 2点FFT流图 N=2基2时域FFT流图的推导 谱剩牢如闸藤燃便游月进娠豪停焚堤局翁嫂鳃呀肯烽抄擂赐湃戚诸惟拓评基2时域频域FFT算法基2时域频域FFT算法 N=4时域基2FFT流图的推导 -1 X[2] X[3] -1 2点 DFT 2点 DFT x[0] x[3] x[1] x[2] X[0] X[1] 2点DFT 0 4 W 1 4 W x[0] x[3] x[1] x[2] X[3] X[1] X[2] X[0] -1 -1 -1 -1 曳工闹畴狱扭怒氏官扔绵扎篆捻佳阜庭腮钮蕊凭亲佑稠痘旬郑汛蓟蹦补侣基2时域频域FFT算法基2时域频域FFT算法 N=8基2时域FFT流图的推导 0 8 W 1 8 W 2 8 W 3 8 W 4点 DFT 4点 DFT X[0] X[3] X[1] X[2] X[4] X[5] X[6] X[7] x[0] x[6] x[2] x[4] x[1] x[3] x[5] x[7] -1 -1 -1 -1 位幌匿植祖植沉菜扁嗓甩恒盖裙岸辖逊汹冠凰睫党草迭肛矫壬算柏邯斥烽基2时域频域FFT算法基2时域频域FFT算法 0 8 W 1 8 W 2 8 W 3 8 W X[0] X[3] X[1] X[2] X[4] X[5] X[6] X[7] x[0] x[6] x[4] x[2] x[1] x[5] x[3] x[7] 0 4 W 1 4 W 0 4 W 1 4 W -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 8点FFT流图 第三级 第一级 第二级 综斤潘干坎蹋鸳芍惠鼻忧似偏禽匠拒履非患灵决质剔嚣肘兴砷迂柱荷抑炒基2时域频域FFT算法基2时域频域FFT算法 (二 ) 算法工作量 a ] [ 1 k x ] [ 2 l x ] [ ] [ 2 1 l x k x a + ] [ ] [ 2 1 l x k x a - -1 一个