第6讲 FFT算法设计.ppt

第6讲 FFT算法设计 傅立叶变换将信号从时域转换为频域，可以进行模拟信号的频率分析 离散傅立叶变换(DFT)将信号从频域转换为数字(频)域，可以进行数字信号（模拟信号数字化）的频率分析 为了实现DFT在计算机上的快速实现，提出了快速离散傅立叶变换(FFT) 如何有傅氏变换-DFT-FFT? 欧拉公式： = 令 ,称为旋转因子 = 上式中，k对应数字域，n对应时域 另有推导时需用到的公式： 1） ,l N为l个周期 2） ,N-m为加上一个周期 3) ,其中 4） 推导分析 若序列x(n)的长度为N，且满足N=2M,(M为自然数) 按n的奇偶性把x(n)分解为两个N/2的子序列： x1(r)=x(2r), r=0,1,…,N/2 – 1 x2(r)=x(2r+1), r=0,1,…,N/2 – 1 则x(n)的DFT为： = ，k=0,1,…,N/2 - 1 其中X1(k)和X2(k)均以N/2为周期 = = ，k=0,1,…,N/2 – 1 其中公式 同理，可推出： ， k=0,1,…,N/4 - 1 ， k=0,1,…,N/4 – 1 …… 分到最后，k=0，进行蝶形运算的两个输入就是最初输入序列x(n)的其中两个。 蝶形分解图示 N=8点FFT运算图示 N=16点FFT运算图示 蝶形运算规律 设序列x(n)已经经过时域抽选(倒序)后，存入数组X中。如果蝶形运算的两个输入相距B个点，用原位计算（即计算结果还放在数组的原来位置），则蝶形运算可表示成如下形式： = 其中：p=J*2M-L; J=0,1,…,2L-1-1 L=1,2,…,M 下标L表示第L级运算，XL(J)则表示第L级运算后数组元素X(J)的值。 如果要用实数运算完成上述蝶形运算，可由下面的算法进行： 设： 下标R表示实部 下标I表示虚部 X’R(J)代表上一次的实数值 = = = = 公式(7)(8)(9)主要用于FFT的软件编程 可以看出，流程图总共有3个循环 外循环：次数为级数L的变换范围 中循环：为根据当前L求出各个不同的p，循环次 数为p的个数2L-1 内循环：为每级中每个p对应的蝶形运算个数，循环次数为2M-L 内循环中k值每次变换范围(增量)为2L,这是同一级中每个相同的p对应不同蝶形运算的间隔。 看图推导软件编程规则：方法一 1. 当N=8时，第L级共有2L-1个不同的旋转因子。 因为N=2M，所以有L=1,2,…,M，即L的最大值为M 当L=1时,p=0; (p称为旋转因子指数) 当L=2时,p=0,2; k=2 (k为p的增量) 当L=3时,p=0,1,2,3; k=1 当N=16时 当L=1时,p=0; 当L=2时,p=0,4;