山大数字信号处理课件04快速傅立叶变换.pdfVIP

数字信号处理 Digital Signal Processing(DSP) 主讲: 江铭炎 教授/博导 联系方式: 山东大学 中心校区 信息科学与工程学院南楼203 Email: jiangmingyan@sdu.edu.cn 课程网址: 信息学院网页精品课程数字信号处理 JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 1 / 140 第4章 快速傅里叶变换(FFT) JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 2 / 140 快速傅里叶变换FFT • 有限长序列通过离散傅里叶变换 (DFT) 将其 频域离散化成有限长序列. 但其计算量太大, 很 难实时地处理问题, 因此引出了快速傅里叶变 换 (FFT) . • FFT 并不是一种新的变换形式, 只是DFT的一 种快速算法. 并且根据对序列分解与选取方法 的不同而产生了 FFT 的多种算法 . • FFT 在离散傅里叶反变换、线性卷积和线性 相关等方面也有重要应用。 JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 3 / 140 FFT产生故事 当时加文(Garwin)在自已的研究中极需要一个计算 傅里叶变换的快速方法。他注意到图基(J.W.Turkey)正 在写有关傅里叶变换的文章，图基概括地对加文介绍 了一种方法，它实质上就是后来的著名的库利(Cooley J.W) 图基算法。1965年库利-- 图基在计算数学、 Mathematic of Computation 杂志上发表了著名的“机器 计算傅里级数的一种算法”文章，提出一种快速计算 DFT的方法和计算机程序--揭开了FFT发展史上的第一 页，促使FFT算法产生原因还有1967年至1968年间 FFT的数字硬件制成，电子数字计算机的条件，使DFT 的运算大大简化。 JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 4 / 140 本章主要内容 (1). 直接计算DFT算法存在的问题及改进途径。 (2). 多种DFT算法(时间抽取算法DIT算法，频 率抽取算法DIF算法，线性调频Z变换即CZT 法) (3). FFT的应用 JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 5 / 140 4.1 直接计算DFT的运算 量，减少运算量的途径 JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 6 / 140 4.1.1. 直接计算DFT计算量 • 问题提出： 设有限长序列x(n),非零 值长度为N,计算对x(n)进行一次DFT 运算，共需多大的运算工作量? JMY © Copyrights Reserved, SDU 山东大学 数字信号处理精品课程 7 / 140 (1). 比较DFT与IDFT之间的运算量 N 1 DFT kn x n( ) X k( ) xn(W)   N0,1, k 1 N 