第3章 快速傅里叶变换(FFT)NEW.ppt

第3章 快速傅里叶变换(FFT) 3.1 引言 3.2 直接计算DFT及其改进途径 3.3 DIT基-2 FFT算法 3.4 DIF FFT算法 3.5 分裂基FFT算法 3.6 几个问题 3.7 Small-Radix Cooley-Tukey Algorithms 3.1 引言 DFT是信号分析与处理中的一种重要变换。因直接计算DFT的计算量与变换区间长度N的平方成正比，当N较大时，计算量太大，所以在快速傅里叶变换(简称FFT)出现以前，直接用DFT算法进行谱分析和信号的实时处理是不切实际的。直到1965年发现了DFT的一种快速算法以后，情况才发生了根本的变化。 3.2 直接计算DFT及其改进途径 3.3 按时间抽取(DIT)的基-2 FFT算法 FFT算法基本上分为两大类：时域抽取法FFT(Decimation In Time FFT,简称DIT-FFT)和频域抽取法FFT(Decimation In Frequency FFT,简称DIF―FFT)。 设序列x(n)的长度为N，且满足 则x(n)的DFT为 其中X1(k)和X2(k)分别为x1(r)和x2(r)的N/2点DFT。由于X1(k)和X2(k)均以N/2为周期，且 ，所以X(k) 又可表示为 与第一次分解相同，将x1(r)按奇偶分解成两个N/4长的子序列x3(l)和x4(l)，即 用同样的方法可计算出 每一级运算都需要N/2次复数乘和N次复数加(每个蝶形需要两次复数加法)。所以，M级运算总共需要的复数乘次数为 1.原位(同址)运算 由DIT-FFT运算流程图可以看出，DIT-FFT的运算过程很有规律。 点的FFT共进行M级运算，每级由N/2个蝶形运算组成。 3. 序列的倒序 DIT―FFT算法的输入序列的排序看起来似乎很乱，但仔细分析就会发现这种倒序是很有规律的。由于 ，所以顺序数可用M位二进制数(nM-1nM-2…n1n0)表示。 3. 蝶形运算规律 设序列x(n)经时域抽选(倒序)后，存入数组X中。如果蝶形运算的两个输入(数据)节点“ 距离”为 ，应用原位计算，则蝶形运算可表示成如下形式： 4. 编程思想及程序框图 3.4 频域抽取法FFT(DIF―FFT) 3.5 分裂基FFT算法 当n=pq，且p=N/4，q=4时，n可表示为 则上式可写成如下更简明的形式： 3.6 几个问题 1.奈奎斯特采样定律：确定不产生频谱混叠失真的条件。 2.栅栏效应：只能得到离散点频谱，sinc插指。 3.频谱泄漏：有限截断造成频谱能量的散布。（P142） 4.谱间干扰：因频谱泄漏导致主谱线两边出现旁瓣，引起不同频率分量间的干扰。 5.非2的幂次方分裂算法 3.7 Small-Radix Cooley-Tukey Algorithms 再将上式中的k表示为 可得 对k0=0，1,2,3,并用k表示k1，用n表示n0，可以 写出 令 对偶序号使用基-2FFT算法，对奇序号使用基-4FFT算法，称分裂基FFT算法 针对 的算法中具有最少乘法次数，且同址运算。 将 分 成三个序列。 3.5.2 分裂基FFT算法原理二 偶序号的 点DFT 奇序号的 点DFT 利用周期性 分成四段： 的第一级分解得4个分裂基 同样对 、 、 作进一步分解。 和 的第二级分解分别是基-4的4点DFT。 的第二级分解得2个分裂基。 一个基-4的4点DFT和2个基-2的4点DFT。 基-2，基-4等基本碟形结都没有乘法，只有每个分裂基有两次复乘。 运算量： 分裂基碟形数： 3.5.3 运算量 * 为自然数 按n的奇偶把x(n)分解为两个N/2点的子序列 3.3.1算法原理 蝶形运算流图符号 N点DFT的一次时域抽取分解图(N=8) 那么，X1(k)又可表示为 式中 同理，由X3(k)和X4(k)的周期性和Wm N/2的对称 性 Wk+N/4 N/2=-Wk N/2 最后得到： 其中 N点DFT的第二次时域抽取分解图(N=8) N点DIT―FFT