快速离散傅里叶变换.pptVIP

基2时间抽取的FFT算法 重庆师范大学 运算规律及编程思想 2 3、旋转因子变化规律 N点DIT―FFT运算流图中，每个蝶形都要乘以旋转因子WNp ，p称为旋转因子的指数。N＝8 ＝23 时各级的旋转因子 第一级：M=1，有1个旋转因子： WNp =WN/4J =W2MJ J=0 第二级：M=2，有2个旋转因子： WNp =WN/2J =W2MJ J=0,1 第三级：M=3，有4个旋转因子： WNp =WNJ =W2MJ J=0,1,2,3 基2时间抽取的FFT算法 重庆师范大学 运算规律及编程思想 2 3、旋转因子变化规律 对于N = 2L的一般情况，第M级共有2M-1个不同的旋转因子： WNp =W2MJ J=0,1,2,… ,2L-1－1 2M =2M?L?2L = N?2M?L 故： 按照上面两式可以确定第M级运算的旋转因子。 J J?2L-M WNp =W2MJ =WN?2M-L = WN p=J×2L-M， J=0,1,2,… ,2M-1－1 基2时间抽取的FFT算法 重庆师范大学 运算规律及编程思想 2 3、旋转因子变化规律 同一级中，同一旋转因子对应蝶形数目 第M级FFT运算中，同一旋转因子用在2L-M个蝶形中； 同一级中，蝶形运算使用相同旋转因子之间相隔的“距离” 第M级中，蝶距：D=2M。 同一蝶形运算两输入数据的距离 在输入倒序，输出原序的FFT变换中，第M级的每一个蝶形的2个输入数据相距：B=2M-1。 基2时间抽取的FFT算法 重庆师范大学 运算规律及编程思想 2 4、蝶形运算规律 序列经过时域抽选后，存入数组中，如果蝶形运算的两个输入数据相距B个点，应用原位计算，蝶形运算可表示成如下形式： XM-1(J) X M-1 (J+B) XM (J)= XM-1(J)+ WNp? X M-1 (J+B) XM (J) = XM-1(J)?WNp? X M-1 (J+B) WNp p=J×2L-M， J=0,1,2,… ,2L-1－1 基2时间抽取的FFT算法 重庆师范大学 运算规律及编程思想 2 5、DIT-FFT程序框图 根据DIT-FFT原理和过程，DIT-FFT的完整程序框图包括以下几部分： (1)倒序:输入自然顺序序列x(n)，根据倒序规律，进行倒序处理； (2)循环层1:确定运算的级数，M=1?L (N=2L)；确定一蝶形两输入数据距离B=2M-1； (3)循环层2:确定M级的2M-1个旋转因子；旋转因子指数p=J2L-M，J=0?B-1； (4)循环层3:对于同一旋转因子，用于同一级2L-M个蝶形运算中：k的取值从J到N-1，步长为2M (使用同一旋转因子的蝶形相距的距离); (5)完成一个蝶形运算。 基2时间抽取的FFT算法 重庆师范大学 运算规律及编程思想 2 5、DIT-FFT程序框图 基2时间抽取的FFT算法 重庆师范大学 例 2 已知x(n)={1,2,3,4}，利用基2FFT算法流图计算 6 -2 2 j 10 -2 -2+2j -2-2j DFT{x(n)}= {10, -2+2j, -2, -2-2j} 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 已知x(n)={1,2,3,4}，补零，插零序列DFT x1(n)={1,2,3,4} x2(n)={1,2,3,4,0,0,0,0} x3(n)={1,0,2,0,3,0,4,0} DFT{x1(n)}={10, -2+2j, -2, -2-2j} DFT{x2(n)}={10, -0.4142-7.2426j, -2+2j, 2.4142-1.2426j, -2, 2.4142+1.2426j , -2-2j, -0.4142-7.2426j} DFT{x3(n)}={10, -2+2j, -2, -2-2j, 10, -2+2j, -2, -2-2j} 基2频率抽取的FFT算法 重庆师范大学 算法原理 3 基2频率抽取的FFT算法 重庆师范大学 算法原理 3 将X(k)分解成偶数组与奇数组， k取偶数(k=2r,r=0,1,…,N/2-1)时