基时域频域FFT算法.ppt

第1页，共33页，星期日，2025年，2月5日4点序列{2，3，3，2}DFT的计算复杂度复数加法N(N-1)复数乘法N2如何提高DFT的运算效率?问题的提出第2页，共33页，星期日，2025年，2月5日一般性DFT:直接计算的计算量：N-1复数乘法N对N个不同X[m]复数加法N(N-1)复数乘法N2对一固定的m复数加法如计算1024点DFT：复数乘法次数：N2=10242=220=1048576第3页，共33页，星期日，2025年，2月5日解决问题的思路1.将长序列DFT分解为短序列的DFT2.利用旋转因子的周期性、对称性、可约性。第4页，共33页，星期日，2025年，2月5日1)周期性(periodicity)2)复共轭对称性(complexconjugate)3)当N是偶数时的性质第5页，共33页，星期日，2025年，2月5日解决问题的方法将时域序列逐次分解为一组子序列，利用旋转因子的特性，由子序列的DFT来实现整个序列的DFT。基2时间抽取(Decimationintime)FFT算法基2频率抽取(Decimationinfrequency)FFT算法第6页，共33页，星期日，2025年，2月5日3.1基2时域抽取FFT算法算法推导:N=2M时域抽取(Decimationintime)]a][1kx][2lx][][21lxkxa+[][21lxkxa--1(一)算法原理基本蝶形(butterfly)第7页，共33页，星期日，2025年，2月5日算法推导:N=2M时域抽取第8页，共33页，星期日，2025年，2月5日例：]0[X]1[X-11]0[x]1[x2点FFT流图N=2基2时域FFT流图的推导第9页，共33页，星期日，2025年，2月5日N=4时域基2FFT流图的推导-1X[2]X[3]-12点DFT2点DFTx[0]x[3]x[1]x[2]X[0]X[1]2点DFT04W14Wx[0]x[3]x[1]x[2]X[3]X[1]X[2]X[0]-1-1-1-1第10页，共33页，星期日，2025年，2月5日N=8基2时域FFT流图的推导08W18W28W38W4点DFT4点DFTX[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第11页，共33页，星期日，2025年，2月5日08W18W28W38WX[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]04W14W04W14W-1-1-1-1-1-1-1-1-1-1-1-18点FFT流图第三级第一级第二级第12页，共33页，星期日，2025年，2月5日(二)算法工作量a][1kx][2lx][][21lxkxa+][][21lxkxa--1一个蝶形计算量：复数乘法1，复数加法2复数乘法:N，NNM2log22=复数加法：N=2M时，分解级数为M，每级蝶形数N/2计算1024点复数乘法次数：DFT：N2=10242=220=1048576；FFT：第13页，共33页，星期日，2025年，2月5日(三)流图规律1)原位运算(In-placeComputations)a][1kx][2lx][][21lxkxa+][][21lxkxa--1优点：节省存储空间第14页，共33页，星期日，2025年，2月5日(三)流图规律第二级的蝶形系数为，蝶形节点的距离为2。第一级的蝶形系数均为，蝶形节点的距离为1。第三级的蝶形系数为，蝶形节点的距离为4。第M级的蝶形系数为，蝶形节点的距离为N/2。2)FFT算法流图旋转因子规律第15页，共33页