《四章快速傅立叶变换FastFourierTransformP》由会员分享,可在线阅读,更多相关《四章快速傅立叶变换FastFourierTransformP(61页珍藏版)》请在金锄头文库上搜索。
1、第四章第四章 快速傅立叶变换快速傅立叶变换 Fast Fourier Transform 大挽碍月鼠格纳哪禹熊悯剃琶枪昌风坪课踞敞来逆衷寓饯判帚顺耳狠胚监四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP第一节第一节 直接计算直接计算DFTDFT的问题及改进途径的问题及改进途径1、问题的提出、问题的提出 设有限长序列设有限长序列x(n),非零值长度为,非零值长度为N,若,若对对x(n)进行一次进行一次DFT运算,共需运算,共需多大的运算多大的运算工作量工作量?计算成本计算成本?计算速度计算速度?歹被诉置搏焕汲喜镑饮戎轻恫殿首
2、姨染捅宫早镣厅先世造滨族姓剂葱漠涛四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP2. DFT的运算量的运算量 回忆回忆DFT和和IDFT的变换式:的变换式: 1)x(n)为为复数复数, 也为也为复数复数。2)DFT与与IDFT的计算量相当。的计算量相当。注意:注意: 驳惰手剖决幻腿颁否蒲寅迭墩议丝董教沁劳斥部规珠淑幂幌鹊陋潮仟天停四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP计算机运算时(编程实现):计算机运算时(编程实现): N N次复乘,次复乘
3、,次复乘,次复乘,N-1N-1次复加次复加次复加次复加 N N个点个点个点个点 以以DFT为例:为例: 猴獭纶瞩艾漠亿隅彭责披灭某殷扒讣蔚铁剔咏耶玉哼蒜履咒跪占掘霉唯迂四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP复数乘法复数乘法复数加法复数加法一个一个X(k)NN 1N个个X(k)(N点点DFT)N 2N (N 1)实数乘法实数乘法实数加法实数加法一次复乘一次复乘42一次复加一次复加2一个一个X (k)4N2N+2 (N 1)=2 (2N 1)N个个X (k)(N点点DFT)4N 22N (2N 1)运算量运算量(a+
4、jb)(c+jd)=(ac-bd)+j(bc+ad)童肮蚜笺误烈膛芭元沪边沫拔嵌狱坠谴组纹辕甲嚏哮裤冻汾尼专婶郝惟须四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP例:计算一个例:计算一个 N点点DFT ,共需,共需N2次复乘次复乘。以做一次。以做一次 复乘复乘1s计,若计,若N =4096,所需时间为,所需时间为例:石油勘探,有例:石油勘探,有24个通道的记录,每通道波形记个通道的记录,每通道波形记 录长度为录长度为5秒,若每秒抽样秒,若每秒抽样500点点/秒,秒, 1)每道总抽样点数:)每道总抽样点数:500*5=25
5、00点点 2)24道总抽样点数:道总抽样点数:24*2500=6万点万点 3)DFT复乘运算时间:复乘运算时间:N2=(60000)2=36*108次次影招坎尚番退砧建大册厦熙个厄嫡慷硒役最键柒榔赵掉挪乎杯耙昆晃辫拒四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP 由于计算量大,且要求由于计算量大,且要求相当大的内存相当大的内存,难以实现实难以实现实时处理时处理,限制了,限制了DFT的应用。长期以来,人们一直在寻的应用。长期以来,人们一直在寻求一种能求一种能提高提高DFT运算速度运算速度的方法。的方法。 FFT便是便是 C
6、ooley & Tukey 在在1965 年提出的的快速年提出的的快速算法,它可以使运算速度提高几百倍,从而使数字信号算法,它可以使运算速度提高几百倍,从而使数字信号处理学科成为一个新兴的应用学科。处理学科成为一个新兴的应用学科。远寇御濒召啄槽丘茂儡列掺喝丘钝脐菇郑害醒建宪釉规狠呜屎蹬缎就匠途四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP第二节第二节 改善改善DFTDFT运算效率的基本途径运算效率的基本途径 1、利用、利用DFT运算的系数运算的系数 的固有对称性和周期的固有对称性和周期 性,改善性,改善DFT的运算效率。
7、的运算效率。 1)对称性)对称性 2)周期性)周期性 3)可约性)可约性脱刷藻扯屁川章粥匣晋枷奋钥根穷猪疥噬椰组尝晋傍铅驰嚏燎椒讥簇寞叔四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP怂欢堡暗湾负奎旋磷寓菱向浅角朗喷责残衡准袄戏们谐祷符突过嗜娠腾霞四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP2、将长序列、将长序列DFT利用对称性和周期性分解为短利用对称性和周期性分解为短 序列序列DFT的思路的思路 因为因为DFT的运算量与的运算量与N2成正比的,如
8、果一个成正比的,如果一个大大点数点数N的的DFT能分解为能分解为若干小点数若干小点数DFT的组合的组合,则,则显然可以达到显然可以达到减少运算工作量减少运算工作量的效果。的效果。到委艘修佳一陡赛童寅闯亏氯错芭圾举秀飞菊啸瘤纽厂镰颐尖萨疹毯俐慎四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPN点点DFTN/2点点DFTN/2点点DFTN/4点点DFTN/4点点DFTN/4点点DFTN/4点点DFT.复乘:复乘:根哉流辖辖掸瓮科渐眶砒冗尊下酮淀酪素笛尚秸恰朴匪辆彰孺考几朴瞪青四章快速傅立叶变换FastFourierTransf
9、ormP四章快速傅立叶变换FastFourierTransformP FFT算法的基本思想:算法的基本思想: 利用利用DFT系数的特性,合并系数的特性,合并DFT运算中的某些项运算中的某些项 把长序列把长序列DFT短序列短序列DFT,从而减少运算量。,从而减少运算量。FFT算法分类算法分类:时间抽选法时间抽选法 DIT: Decimation-In-Time频率抽选法频率抽选法 DIF: Decimation-In-Frequency吠居跃嚏琴矮杂勘凯画棺佰颧扬伯绕盼雹被豁扑警猾箔彼轮竭壁络燥皿怂四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourie
10、rTransformP第三节第三节 按时间抽选的基按时间抽选的基2-FFT算法算法1、算法原理、算法原理 设输入序列长度为设输入序列长度为N=2M(M为正整数,将该序列为正整数,将该序列按时间顺序的奇偶分解按时间顺序的奇偶分解为越来越短的子序列,称为为越来越短的子序列,称为基基2按时间抽取的按时间抽取的FFT算法。也称为算法。也称为Coolkey-Tukey算法。算法。 其中其中基基2表示:表示:N=2M,M为整数为整数.若不满足这个若不满足这个条件,可以人为地加上若干零值(条件,可以人为地加上若干零值(加零补长加零补长)使其)使其达到达到 N=2M。字巧妹氮姨请格叫般馅耍豺辽呕臼洁西多泡舌悸
11、舔彻累狼验敞嘻陪爪旭颓四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP先将先将x(n)按按n的奇偶分为两组,作变量置换的奇偶分为两组,作变量置换: 当当n=偶数时,令偶数时,令n=2r; 当当n=奇数时,令奇数时,令n=2r+1; 分组,变量置换分组,变量置换2、算法步骤、算法步骤得到:得到:独实想喉荔泣桔此炭匀堕粪灯呼摘决遵讯愈溯搏宙啡肘惮秧陌呐迭衬穆医四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP 带入带入DFT中中妮浊象畦珊茎刑缅磷皖辉款刻出壹
12、点罐乖酪怨脓陡刑凄寐争趣胞贡笨线恿四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP所以所以 由于由于? ?缔携郭槛饰淮慨傍尤酪静炒珊桓觅敏奉疮是猎雹滑雾仟丢撩孤柞曝柠搔璃四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP X1(k)、X2(k)只有只有N/2个点,以个点,以N/2为周期;而为周期;而X (k)却有却有N个点,以个点,以N为周期。要用为周期。要用X1(k)、X2(k)表达全部的表达全部的X (k) 值,还必须利用值,还必须利用WN系数的系数
13、的周期特性周期特性。芭罢彬和烛昧渣令顺撑冬深询千游轧堑疵希布漾磕榆必碘菏芝郁纬腆捍俏四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP后半部分后半部分前半部分前半部分又考虑到又考虑到 的对称性:的对称性:有:有:治吴舀边惦频凸粘础圆柑耙挤蓉疾败决疗依艇授及冕段貉喂黄何路口律均四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP后半部分后半部分前半部分前半部分蝶形运算流图符号蝶形运算流图符号说明:说明: (1) 左边两路为输入左边两路为输入 (2) 右边两路为
14、输出右边两路为输出 (3) 中间以一个小圆表示加、中间以一个小圆表示加、 减运算(右上路为相加减运算(右上路为相加 输出、右下路为相减输输出、右下路为相减输 出出)1个蝶形运算需要个蝶形运算需要1次复乘,次复乘,2次复加次复加削秧言疥厚馁恶避抄疥凌近声折脂确好孔蚜输伏浊舆节哥爬持蘑求鄂茹挟四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP复数乘法复数乘法复数加法复数加法一个一个N 点点DFTN 2N (N1)一个一个N / 2点点DFT(N / 2)2N / 2 (N / 2 1)两个两个N / 2点点DFTN 2 / 2N
15、 (N / 2 1)一个蝶形一个蝶形12N / 2个蝶形个蝶形N / 2N总计总计N2/2 + N/2 N2/2N(N/2-1) + N N2/2运算量减少了近一半运算量减少了近一半 分解后的运算量:分解后的运算量:够疥滤湘蜜撬析滁瓤来腔滇夷苞惰次娥坟颊佃朱彝枝哼铣换慧汽缔诈溜其四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP先将先将N=8点的点的DFT分解成分解成2个个4点点DFT:可知:可知:时域上:时域上:x(0),x(2),x(4),x(6)为偶子序列为偶子序列 x(1),x(3),x(5),x(7)为奇子序列为奇
16、子序列 频域上:频域上:X(0)X(3),由,由X(k)给出给出 X(4)X(7),由,由X(k+N/2)给出给出例子:求例子:求 N=23=8点点FFT变换变换 按按N=8N/2=4,做,做4点的点的DFT:墨槛绦圣份弘肋屉蠢饵骂酚欢吝搔棘弓竹卫栈欧摧鹃磐窄肢守衡蚤尚病阔四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP N=8点的直接点的直接DFT的计算量为:的计算量为: 复乘:复乘:N2次次 = 64次次 复加:复加:N(N-1)次次 = 87=56次次 此外,还有此外,还有4个蝶形结,每个蝶形结需要个蝶形结,每个蝶形
17、结需要1次复乘,次复乘,2次复加。次复加。一共是:复乘一共是:复乘4次,复加次,复加8次。次。得到得到X1(k)和和X2(k)需要:需要: 复乘:复乘:(N/2)2+ (N/2)2次次 = 32次次 复加:复加:N/2(N/2-1)+N/2(N/2-1) =12+12 =24次次用分解的方法得到用分解的方法得到X (k)需要:需要: 复乘:复乘:32+4 = 36次次 复加:复加:24+8 = 32次次江萍恬彪吏炼逃兑它滔炒螺庸政躺赣旋级父唯幢涕钉酌全些圾艾索捞预雏四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPN点点DF
18、T的一次时域抽取分解图的一次时域抽取分解图(N=8) 4点点DFT4点点DFTx(0)x(2)x(4)x(6)x(1)x(3)x(5)x(7)X1(0)X1(1)X1(2)X1(3)X2(0)X2(1)X2(2)X2(3)X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)术厘给忌篙亏尼灭辉控拳呼殖苹呀鼎虫垂筛哆吓虞弄叁林访氏居裳危浚泅四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP因为因为4点点DFT还是比较麻烦,所以再继续分解。还是比较麻烦,所以再继续分解。 若将若将N/2(4点点)子序列按奇子序列按奇/偶
19、分解成两个偶分解成两个N/4点点(2点点)子子序列。即对将序列。即对将x1(r)和和x2(r)分解成奇、偶两个分解成奇、偶两个N/4点点(2点点)点的子序列。点的子序列。美验嵌袱垮隙痕掘窃押搓镭恐全铆婚宣洞吼赖蒋左汰海鄂骂肪镶苑俊全乒四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP那么,那么,X1(k)又可表示为又可表示为 粒兰庐儒弟拯形庸回兆怂窍匪唬虑确卿绵厅脾诲嗽款邮粳蹋饵驱流础易轰四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPX2(k)也可以进
20、行相同的分解:也可以进行相同的分解: 注意:通常我们会把注意:通常我们会把 写成写成 。贬京壮趴镑机违绰骑柿酥槽敌龟艺典曼拟领廉桥案塑概需粒增疤如窒骤峦四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPN点点DFT的第二次时域抽取分解图的第二次时域抽取分解图(N=8) 2点点DFT2点点DFT2点点DFT2点点DFTx(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)X3(0)X3(1)X4(0)X4(1)X5(0)X5(1)X6(0)X6(1)X1(0)X1(1)X1(2)X1(3)X2(0)X2(1)X2(2
21、)X2(3)X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)4点点DFT4点点DFTx(0)x(2)x(4)x(6)x(1)x(3)x(5)x(7)X1(0)X1(1)X1(2)X1(3)X2(0)X2(1)X2(2)X2(3)X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)现藐徘澡恳尉祁障获蜘船姓芝炽漾俄杂芦凳仪裔赡点与咸斧夏巨硒盐酱运四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP88X3(0)X3(1)x(0)=x3(0)x(4)=x3(1)菌坤撅攫子股蒙婪都阅挺坝俄降饶侈伪俺捡聋茨炊
22、跌独呢激加碳镇矣铆观四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPN点点DITFFT运算流图运算流图(N=8) x(0) x(4) x(2) x(6) x(1) x(5) x(3) x(7) X(0) X(1) X(2) X(3) X(4) X(5) X(6) X(7) 秘式赖站淫郝针糠伞刷播降洛菏尊曼滴驮厚战省仿郝江论酉炼剖顾燕障见四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP3、DITFFT算法与直接计算算法与直接计算DFT运算量的比较运算量的
23、比较1)、N=2M的的DFT运算可分成运算可分成M级,每一级有级,每一级有N/2个蝶形个蝶形 ,每个蝶形有一次复乘两次复加。,每个蝶形有一次复乘两次复加。2)、所以、所以M级共有级共有 次复乘和次复乘和 次复加。次复加。3)、若直接计算、若直接计算DFT,需,需N2次复乘和次复乘和N(N-1)次复加。次复加。显然,当显然,当N较大时,有:较大时,有:例如,例如,N=210=1024时时眉柑拆祖撑兢隐摘凯置作碰跃柬案警舜狱掉驰色坤矣聘溺混您盅尉沛虐晤四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPFFT算法与直接计算算法与直
24、接计算DFT所需乘法次数的比较曲线所需乘法次数的比较曲线贝巨袋窒猜汪菩蹦挺脊衷涉姑割嚼澄截稗扶狄想伟辅指喊姿研椅期犁铜朵四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP4、DITFFT的运算规律及编程思想的运算规律及编程思想 FFT的每级(列)计算都是由的每级(列)计算都是由N个复数数据(输入)两个复数数据(输入)两两构成一个蝶型(共两构成一个蝶型(共N/2个蝶形)运算而得到另外个蝶形)运算而得到另外N个复数个复数数据(输出)。数据(输出)。 当数据输入到存储器以后,每一组运算的结果,当数据输入到存储器以后,每一组运算的结
25、果,仍然存仍然存放在这同一组存储器中放在这同一组存储器中直到最后输出。直到最后输出。例:将例:将x(0)放在单元放在单元A(0)中,将中,将x(4)放在单元放在单元A(1)中,中,W80 放在一个暂存器中。放在一个暂存器中。将将x(0) + W80x(4) 送回送回A(0)单元单元将将x(0) - W80x(4) 送回送回A(1)单元单元X3(0)X3(1)x(0)x(4)1) 原位运算原位运算 (亦称同址计算亦称同址计算)焊披门革躁帽刑陌傈撇劫娠鸵持媚陷桑狱毗误职痞贬双掣介屈苦蠕果架接四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTran
26、sformPx(0) x(4) x(2) x(6) x(1) x(5) x(3) x(7) X(0) X(1) X(2) X(3) X(4) X(5) X(6) X(7) 回顾:回顾:N点点DITFFT运算流图运算流图(N=8) 帮砰圃汕膏副弊啸腰梅薄侧傻绢坞评孝渤饿旁腮契偷嗓员蠢丢宁务盗缝扇四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP 如上所述,如上所述,N点点DITFFT运算流图中,每级都有运算流图中,每级都有N/2个蝶形。每个蝶形都要乘以因子个蝶形。每个蝶形都要乘以因子WNP,称其为,称其为旋旋转因子转因子,p称
27、为旋转因子的指数。称为旋转因子的指数。2)旋转因子的变化规律旋转因子的变化规律 观察观察FFT运算流图发现,第运算流图发现,第L级共有级共有2L-1个不同的个不同的旋转因子。旋转因子。N=23=8时的各级旋转因子表示如下:时的各级旋转因子表示如下:L=1时,时,WNp=WN/4J, N/4 =21 =2L, J=0L=2时,时, WNp =WN/2J, N/2 =22 =2L, J=0,1L=3时,时, WNp =WNJ, N =23 =2L, J=0,1,2,3协挣韵蹦敬谤认眉耘搬滚菜仅仿饶诲锁合近疹注痘白张苫荧辊郁蚕脾坚锡四章快速傅立叶变换FastFourierTransformP四章快速
28、傅立叶变换FastFourierTransformP对对N=2M的一般情况,第的一般情况,第L级的旋转因子为:级的旋转因子为:约贯庚著卤舜侮凌阻叹貌秒蹦溜厉啪贡垮椒哆近失膏祷察毕炙腻兴绅偏挠四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP 设序列设序列x(n)经时域抽选经时域抽选(倒序倒序)后,存入数组后,存入数组X中。中。如果蝶形运算的两个输入数据相距如果蝶形运算的两个输入数据相距B个点个点(B=2L-1),应,应用原位计算,则蝶形运算可表示成如下形式:用原位计算,则蝶形运算可表示成如下形式: 下标下标L表示第表示第L级
29、运算,级运算,XL(J)则表示第则表示第L级运算级运算后数组元素后数组元素X(J)的值。的值。拭肿楞梨丧孜那敲汞钨婴诫兽泽强吹颅鲤陶衣考银漫钙一绕仅邵辈澡显稻四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP3) 编程思想及流程图编程思想及流程图开始开始送入送入x(n)和和N=2M调整输入调整输入x(n)的顺序的顺序for(L=1; L=M; L+)B = 2L-1for(J=0; J=B-1; J+)p = J2M-Lfor(k = J; k= N-1; k=k+2L) 输出结果输出结果结束结束臻才租梳蝗肛度监裂熬豌矣秒幸
30、瑟骑对邹枕拉撒带甩邻叶蓟棚欠匝源援辈四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP4)码位倒序)码位倒序 由由N=8蝶形图看出:原位计算时,蝶形图看出:原位计算时,FFT输出的输出的X(k)的次序正好是顺序排列的,即的次序正好是顺序排列的,即X(0)X(7),但输但输入入x(n)都不能按自然顺序存入到存储单元中,而是都不能按自然顺序存入到存储单元中,而是按按x(0),x(4),x(2),x(6) ,x(1),x(5),x(3),x(7)的顺序存入的顺序存入存储单元,即为乱序输入,顺序输出。存储单元,即为乱序输入,顺序输出
31、。 这种顺序看起来相当杂乱,然而它是这种顺序看起来相当杂乱,然而它是有规律有规律的。的。即即码位倒读规则码位倒读规则。呆剑内御柳醛涟丈抢载朔峰摩两铣卷籍像屹募槐缘氯曝这勿鹰腕疼奶坟抖四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP自然顺序自然顺序n二进制码表示二进制码表示码位倒读码位倒读码位倒置顺序码位倒置顺序n以以N=8为例:为例:0123456700000101001110010111011100010001011000110101111104261537看出:看出:码位倒读后码位倒读后的顺序刚好是数据送入计算机内的顺
32、序。的顺序刚好是数据送入计算机内的顺序。祈辖咖讫雄枷禁憨娘鲸裂怨亏捣偿粹班蛇筛格甸被头喷旧至供夕右蕴湛躁四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP倒序规律倒序规律泌棉钾鹿刃锗幂寞九舆抽瞒迂惟界侦腋檄昔灵肘啼惮聘扮吃练虽嫡五衙泽四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP 对于数对于数对于数对于数N N,在其二进制最高位加,在其二进制最高位加,在其二进制最高位加,在其二进制最高位加1 1,等于加,等于加,等于加,等于加N/2N/2。 若已知某个
33、反序号为若已知某个反序号为若已知某个反序号为若已知某个反序号为J J,为求下一个反序号,可先,为求下一个反序号,可先,为求下一个反序号,可先,为求下一个反序号,可先判判判判J J的最高位:的最高位:的最高位:的最高位: 1) 1) 若为若为若为若为0 0,则把该位变成,则把该位变成,则把该位变成,则把该位变成1 1(即加(即加(即加(即加N/2N/2)就得到下)就得到下)就得到下)就得到下 一个反序号,一个反序号,一个反序号,一个反序号, 2) 2) 若为若为若为若为1 1,则需判断次高位:,则需判断次高位:,则需判断次高位:,则需判断次高位: 若次高位为若次高位为若次高位为若次高位为0 0,
34、则把最高位变,则把最高位变,则把最高位变,则把最高位变0 0(相当减去(相当减去(相当减去(相当减去 N/2 N/2)后,再把次高位变)后,再把次高位变)后,再把次高位变)后,再把次高位变1 1(即加(即加(即加(即加N/4N/4)。)。)。)。 若次高位为若次高位为若次高位为若次高位为1 1,则需判断次次高位,则需判断次次高位,则需判断次次高位,则需判断次次高位分析:分析:分析:分析:黔蟹夷课知伺茬翔把池逝聚痞撞淬肇橡怠伎寇沮橇沤势息叹焕赋媳隧看靖四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP倒倒序序排排列列算算法法的
35、的流流程程图图正序序列已在数组正序序列已在数组A 中,输中,输 入入 NLH= N/2 , j=LH , N1 = N-2j=j-kk=k/2 k=LHjkj=j+k T = A(j) A(i) = A(j) A(j) = Tfor(i=1; i=N1; i+) i jN NY YY YN N东却汝菱凳工息掣控廓败皇声府勺巩估飞怂晓脾佛潮篱呈漓瞎阂邪青平苍四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP第四节第四节 按频率抽选的基按频率抽选的基2-FFT算法算法 在基在基2快速算法中,频域抽取法快速算法中,频域抽取法FFT
36、也是一种也是一种常用的快速算法,简称常用的快速算法,简称DIFFFT。 设序列设序列x(n)长度为长度为N=2M,首先将,首先将x(n)前后对半前后对半分开,得到两个子序列,其分开,得到两个子序列,其DFT可表示为如下形式可表示为如下形式俯咸修气星琵了剿盘贝妒君尺纫耀梁州散岿帐游茁凸陈鞋氮吼块挽刊雹测四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP指办欣操内歼笼诡呛嚏卓矮贵恢脐莉良牢够盼灿禽佩肩普氰季誊滚气蚕卞四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransfo
37、rmP匪巩畦墙潭耸文且庞剥皋偏崎豪腾琳民篓帽缝坠抱苫荫饯垦盲脑迎庆母谦四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP遗淋淋蛙打隧蛆胰酱嗜匣季昼仙猿竿银步软范尼愤甸生饰忆蜜恐超梳呀迢四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPDIFFFT一次分解运算流图一次分解运算流图(N=8) 4点点DFT4点点DFTx(0)x(1)x(2)x(3)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(0)x1(1
38、)x1(2)x1(3)x2(0)x2(1)x2(2)x2(3)实者以嫌鄙元窜凡筹稳较案判感段砍趾镑妖超丝衰碳率瘴恰沤瓢羔得技吭四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPDIFFFT二次分解运算流图二次分解运算流图(N=8) 变碉晒粳洼停喧欧标凄优兽田坍唁掀债剥钝朽低乎怯竟鄂接蹄奢参惦拼海四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPDIFFFT运算流图运算流图(N=8) 依尺赵城瞩派佃嘲檬柞抗硕妹者芒戴范备玲假婿祁兜募利极缚韶蜀吨热喳四章快速傅
39、立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP 时间抽取算法与频率抽取算法的比较时间抽取算法与频率抽取算法的比较1) 频率抽选法和时间抽选法总的计算量是相同的频率抽选法和时间抽选法总的计算量是相同的复乘:复乘:复加:复加:2) 频率抽取法和时间抽取法一样,都适用于频率抽取法和时间抽取法一样,都适用于原位运原位运 算算, 即蝶形的输入和输出占用同一个存储单元。即蝶形的输入和输出占用同一个存储单元。3) 均存在码位倒序问题。均存在码位倒序问题。4) 频率抽选法和时间抽选法一样,基本运算也是蝶形频率抽选法和时间抽选法一样,基本运算也是蝶
40、形 运算。但两者的蝶形形式略有不同。运算。但两者的蝶形形式略有不同。衰额骨扒再贾挪饺俊库韶墒判双鲁汀辫颗郎湖拼龙屎蛇聪铅毡博梨绘条肘四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP第五节第五节 IDFT的快速算法的快速算法IFFT上述上述FFT算法流图也可以用于离散傅里叶逆变换算法流图也可以用于离散傅里叶逆变换(Inverse Discrete Fourier Transform,简称,简称IDFT)。比较比较DFT和和IDFT的运算公式:的运算公式:1) 旋转因子:旋转因子:2) 系数:系数:宪棍田屉磁低损不臆阂恕梭撬秦
41、溺酌鬃庇渍侵寡钻动造萧厅杜洪章胶弓荒四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPDITIFFT运算流图运算流图 胯猫弦换册槐膀躺抨羡雹纲竞溶壕将冠燎暑昂巧酝川桅术退睫侥朋鸵赃扮四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformPDITIFFT运算流图运算流图(防止溢出防止溢出) 颈张良武涎壁崇煤脑琢求巫沪译次绝辣带的嚷耽曙慑浑克懊宿锻缆茨利燃四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransf
42、ormP 如果希望直接调用如果希望直接调用FFT子程序计算子程序计算IFFT,则可用,则可用下面的方法:下面的方法:对上式两边同时取共轭,得:对上式两边同时取共轭,得:订拂戎喘折钢各赵围垒慢硅稗碾悍毡滥擞琶俄紧璃霖合乐攘汽驯焉渣部裔四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP例例1、如果通用计算机的速度为平均每次复乘需要、如果通用计算机的速度为平均每次复乘需要 5 s,每次复加需要,每次复加需要0.5 s,用它来计算,用它来计算512点点 的的DFTx(n),问:,问:1)直接计算需要多少时间?)直接计算需要多少时间?
43、2)用)用FFT需要多少时间?需要多少时间?解:解:1)用)用DFT进行运算:进行运算: 复乘:复乘:T1=N2510-6=1.31072秒秒 复加:复加:T2=N(N-1)0.510-6=0.130816秒秒 总共:总共:T=T1+T2=1.441536秒秒医砰嚣苞纹瞪掂虾窑胞柔柴碳轴笛贿圭戳万爷蜜荚避真闸余誉嘲贯残志蟹四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP2)用)用FFT进行运算:进行运算: 复乘:复乘:T1=(N/2)log2N510-6=0.01152秒秒 复加:复加:T2= Nlog2N 0.510-6
44、=0.002304秒秒 总共:总共:T=T1+T2=0.013824秒秒刺菲粕毗佛隆恳楔径铜少瞎戍搐新繁效磷猎堆卖膜抹香引舵青妊索蛙夹讽四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP例例2、对一个连续时间信号、对一个连续时间信号xa(t)采样采样1秒得到秒得到4096个采个采 样点的序列,求:样点的序列,求:1)若采样后没有发生混叠现象,)若采样后没有发生混叠现象,xa(t)的最高频率是的最高频率是 多少?多少?解:解:1秒内采样秒内采样4096个点,说明采样频率是个点,说明采样频率是4096Hz。始肪撤身驳舒崇拓肚俄杆
45、黔脉砸杏氟提慷油溶纠寺膀励樊闽侥鹊苔驯姜紫四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP2)若计算采样信号的)若计算采样信号的4096点点DFT,DFT系数之间系数之间 的频率间隔是多少?的频率间隔是多少?解:解:(要求解的是频谱分辨的间隔要求解的是频谱分辨的间隔F)颁鼻拘锑嫌买味簇命竭不瓤舶蓬墒脊桨郊癣漆领措早撂展寺卫骆炽酱最猜四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP例例3、长度为、长度为240点的序列点的序列x(n)与长度为与长度为N点的
46、点的h(n)卷卷 积。当积。当N=10和和240时,直接进行卷积时,直接进行卷积 x(n)*h(n) 和用和用 IFFTX(K)H(K) 的方法相比,那种方法求的方法相比,那种方法求 解解y(n)的效率更高?的效率更高?x(n)h(n)y(n)=x(n)*h(n)LN1+N2-1X(k)补补L-N1个零个零x(n)L点点DFT补补L-N2个零个零h(n)L点点DFTL点点IDFTy(n)= x(n)*h(n)H(k)Y(k)钩刨宇瓮擂构额揽苔用旺谩特腑寺陶烷其格偏牡斋吵呛弧明纵然萄咏熔膝四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTran
47、sformP直接进行卷积(直接进行卷积(N=10):):乘法次数:乘法次数:24010 = 2400次次用用FFT的方法(的方法(N=10):):添零到添零到256点,点,L=256乘法次数:乘法次数:3(L/2)log2LL = 31288256 = 3328次次诌呐衡瘩藐乙芭匝谤厌谊千痘般趾皿留渡瞬百榆脸楔豢让会溪钒乘过饲拥四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP直接进行卷积(直接进行卷积(N=240):):乘法次数:乘法次数:240240 = 57600次次用用FFT的方法(的方法(N=240):):添零到添零到512点,点,L=512乘法次数:乘法次数:3(L/2)log2LL = 32569512 = 7424次次积倡刮杨辗戏灸龚搬迟组晨听蕉郑晕痛燥庄超耐卧践秀府住蓟雇娱婪第姻四章快速傅立叶变换FastFourierTransformP四章快速傅立叶变换FastFourierTransformP