按时间抽取快速傅里叶变换DIT.pdf

第4章 快速傅里叶变换 (FFT) Fast Fourier Transform 4.1 引言 (1) 快速傅立叶变换FFT • 有限长序列 通过离散傅里叶变换(D F T) 将其频域 离散化成有限长序列。但其 计算量太大, 很难实时地处理问题 , 因 此引出了快速 傅里叶变换(FFT) 。 • FFT 并不是一种新的变换形式,它只是 DFT 的一种快速算法。并且根据对序列 分解与选取方法的不同而产生了FFT 的多种算法。 • FFT 在离散傅里叶反变换 、线性卷积 和线性相关等方面也有重要应用。 (2) FFT产生故事 当时加文(Garwin)在研究中极需要一个计算 傅里叶变换的快速方法。他注意到图基 (J.W.Turkey)正在写有关傅里叶变换的文章，因 此详细询问了图基关于计算傅里叶变换的技术 知识。图基概括地对加文介绍了一种方法，它 实质上就是后来的著名的库利(J.W. Cooley) 图基 算法。在加文的迫切要求下，库利很快设计出 一个计算机程序。1965年库利-- 图基在 Mathematic of Computation 杂志上发表了著名的 “机器计算傅里叶级数的一种算法”文章，提出一 种快速计算DFT 的方法和计算机程序--揭开了 FFT发展史上的第一页。 (3) 主要内容 • 直接计算DFT算法存在的问题及改进途 径。 • 多种FFT算法(时间抽取算法，即DIT算 法；频率抽取,即DIF法) • FFT 的应用 4.2 直接计算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 W ( ) ( )  ( ) N N k 0 2 其中x(n)为复数， kn j N kn 也为复数 W e N 所以DFT与IDFT二者计算量相同。 2. 以DFT为例，计算DFT复数运算量 • 计算一个X(k)(一个频率成分)值，运算量为 N 1 例，k=1则 n X (1) x n W( )