数字信号处理及应用 工业和信息化普通高等教育“十二五”规划教材立项项目 新教学课件 作者 卢光跃 黄庆东 包志强 第4章 快速傅里叶变换FFT).pptVIP

第4章 快速傅里叶变换 通信与信息工程学院 数字信号处理教学团队 第4章 快速傅里叶变换 4.1 引言 4.2 直接计算DFT的问题及改进的途径 4.3 按时间抽取（DIT）的基2-FFT算法 4.4 按频率抽取（DIF）的基2-FFT算法 4.5 IDFT的高效算法 4.6 FFT的其他应用—快速卷积 学习目标 理解按时间抽取的基-2FFT算法的原理、运算流图、所需计算量和算法特点； 理解按频率抽取的基-2FFT算法的原理、运算流图、所需计算量和算法特点； 理解IFFT算法； §4.1 引 言 快速傅里叶变换（FFT）并不是一种新的变换，而是离散傅里叶变换（DFT）的一种快速算法。  在信号的频谱分析、 系统的分析、 设计和实现中都会用到DFT的计算。 DFT的计算量太大，很难对问题进行实时处理，难以实用。 1965年首次发现了DFT运算的一种快速算法以后，人们开始认识到DFT运算的一些内在规律，从而很快地发展和完善了一套高速有效的运算方法， 这就是快速傅里叶变换（FFT）的算法。 FFT出现后使DFT的运算大大简化，运算时间一般可缩短一二个数量级之多，使DFT的运算在实际中真正得到了广泛的应用。 §4.2 直接计算DFT的问题及改进的途径 一、 直接计算DFT的运算量问题设x(n)为N点有限长序列，其DFT为 k=0, 1, …, N-1 反变换（IDFT）为 n=0, 1, …, N-1 运算量 一个X(k) 复数乘法 复数加法 N N-1 N2 N(N-1) N个X(k) 一次复数乘法需用四次实数乘法和二次实数加法； 一次复数加法需二次实数加法。 二、 改进的途径 x(n) 长序列————短序列 WNkn §4.3 按时间抽取（DIT）的基 -2 FFT算法 一、 算法原理 1. 一次分解 N→N/2 设序列x(n)长度为N，且满足N=2M，M为正整数。按n的奇偶把x(n)分解为两个N/2点的子序列： k=0,1,…..,N/2-1 利用周期性求X(k)的后半部分 x(0) x(1) x(2) x(3) x(4) x(5) x(6) x(7) -1 -1 -1 -1 WN0 WN1 WN2 WN3 N/2点 DFT N/2点 DFT 2. 二次分解 N/2→N/4 X2(k)也可进行同样的分解： -1 -1 -1 -1 WN0 WN1 WN2 WN3 N/2点DFT分解为两个N/4点DFT 一个N=8点DFT分解为四个N/4点DFT N=8，四个N/4点DFT即四个2点DFT，其输出分别为： X3(k)，X4(k)，X5(k)，X6(k)，k=0, 1 二、 按时间抽取的FFT算法与直接计算DFT运算量的比较 N=2M时，共有M级蝶形， 每级都由N/2个蝶形运算组成； 每个蝶形需要一次复乘、 二次复加； 每级运算都需N/2次复乘和N次复加； M级运算总共需要： 复乘数 复加数 以乘法为例，直接DFT复数乘法次数是N2，FFT复数乘法次数是(N/2) log2N； 直接计算DFT与FFT算法的复数乘法计算量之比为 ： 当N=2048时，这一比值为372.4，即直接计算DFT的运算量是FFT运算量的372.4倍； 当点数N越大时，FFT的优点更为明显。 对一个连续时间信号xa(t)采样1秒得到一个4096个采样点的序列，若计算采样信号的4096点DFT。 假定仅对200≤f≤300Hz频率范围所对应的DFT采样点感兴趣。 （1）若直接用DFT，要计算这些值需要多少次复乘？ MN=101×4096=413696 （2）若用基2的DIT-FFT计算，则需要多少次复乘？ N/2log2N=4096/2log24096=24576 思考题 N点x(n)序列，构造新序列 y(n)=x(n/2),n为偶数时， y(n)=0，n为奇数时， n=0,1,2,…2N-1， 求：Y(K)与X(K)的关系，K=0,1,2,…. 2N-1 例 题： 解：由DIT-FFT可得 y(2n)= y1(n)=x(n), Y1(K)=X(K) y(2n+1)= y2(n)= 0, Y2(K)=0 n=0,1,….,N-1 已知x(n)长度为N（N为偶数）。定义两个长度为N/2的序列如下: 其中，G(k)=DFT[g(n)], H(k)=DFT[h(n)],分别为N/2长。 试利用G(k)和H(k)表示出X(k)=DFT[x(n)]，k=0,1,……N-1。 补充题 提示 提示 三、 按时间抽取的FFT算法