DSP第四章1快速傅里叶变换FFT.pdf

第四章 快速付里叶变换 (FFT) Fast Fourier Transforming 第一节 引言 DFT 的计算复杂性： N −1 DFT : X [k ] ∑x[n]Wkn , 0 ≤k ≤N −1 N n 0 • 对每个特定k , X(k)的DFT的计算：有(N －1) 个复加和N 个 复乘； • 计算整个DFT 共有N (N-1)个复加和N 2 个复乘； • 通过利用W N k n 的对称性和周期性，整个DFT 的复乘计算 量可减至为N log 2 N ⇒ Fast Fourier Transform (FFT). 一、快速付里叶变换FFT • 有限长序列通过离散傅里叶变换(DFT)将其频 域离散化成有限长序列.但其计算量太大(与N 的平方成正比）, 很难实时地处理问题 , 因 此引出了快速 傅里叶变换(FFT) . • FFT 并不是一种新的变换形式,它只是 DFT 的一种快速算法. 并且根据对序列 分解与选取方法的不同而产生了FFT 的多种算法. • FFT 在离散傅里叶反变换 、线性卷积 和线性相关等方面也有重要应用. 。 二、FFT产生故事 当时加文(Garwin)在自已的研究中极需要一个计算 付里叶变换的快速方法。他注意到图基(J.W.Turkey) 正在写有关付里叶变换的文章，因此详细询问了图基 关于计算付里叶变换的技术知识。图基概括地对加文 介绍了一种方法，它实质上就是后来的著名的库利 (Cooley J.W) 图基算法。在加文的迫切要求下，库利很 快设计出一个计算机程序。1965年库利-- 图基在计算 数学 、Mathematic of Computation 杂志上发表了著 名的“机器计算付里级数的一种算法”文章，提出一种 快速计算DFT 的方法和计算机程序--揭开了FFT发展史 上的第一页，促使FFT算法产生原因还有1967年至 1968年间FFT 的数字硬件制成，电子数字计算机的条 件，使DFT 的运算大简化了。 三、本章主要内容 • 1.直接计算DFT算法存在的问题及改进途径。 • 2.多种DFT算法(时间抽取算法DIT算法，频 率抽取算法DIF算法，线性调频Z变换即CZT 法) • 3.FFT 的应用 第二节 直接计算DFT算法 存在的问题及改 进途径 一、直接计算DFT计算量 • 问题提出：设有限长序列x(n),非零 值长度为N,计算对x(n)进行一次DFT 运算，共需多大的运算工作量? 1.比较DFT与IDFT之间的运算量 N −1 DFT kn k 0,1, N −1 x(n) ⎯⎯→⎯ X (k) ∑x(n)W N n 0 IDFT 1 N −1 −kn n 0,1, N −1 ( ) ( ) ( ) X k ⎯⎯⎯→x n ∑X k WN N k 0 2π