数字信号处理丁玉美版教案第4章.ppt

第四章 快速傅立叶变换 （FFT） Fast Fourier Transform 4.1 引言 ? DFT 是数字信号处理的基础，是 一种重要变换。 ? 快速算法的引出，使数字信号处 理技术得到广泛应用。 ? 各种快速算法不断被采用 4.2 基 2 抽取FFT 算法 （the Decimation-In-Time Radix-2 FFT Algorithm） FFT分为两大类： 时域抽取FFT （Decimation-In-Time FFT，简称DIT-FFT） 频域抽取FFT （Decimation-In-Frequency FFT， 简称DIF-FFT） 4.2.1 直接计算DFT的问题及改进的途径 k=0, 1, …, N-1 可见，DFT 与 IDFT 的计算成本基本相同。 直接计算N点DFT 时： 对应一个k需要N次复数乘和（N-1）次 复数加；对所有N个k值，则需要 N2复数乘和N (N-1)次复数加 。 其中: 一次复数乘需要4次实数乘和2次实数加方能 完成。当N较大时，运算量很大。 例如：当N=8时，DFT需要64次复数乘； 当N=1024时，DFT所需复数乘为1048576 次，即一百多万次复数乘。 为了减少运算次数，改进计算的方法： 1）利用旋转因子的对称性、周期性和可约性； 旋转因子 （twiddle factor） 2）使长序列变短。 4.2.2 基2时域抽取法原理 （库利-图基算法） 设序列x（n）的长度为N，且M为自然数 N-point DFT,N is even 将其一分为二，分成偶数和奇数序列项 （the even-indexed and odd-indexed terms） 则N/2点的序列为： even: x1(r)=x(2r) , r=0, 1, …, N/2-1 odd: x2(r)=x(2r+1) , r=0, 1, …, N/2-1 偶数序列 the range: 0?2n?N-2 (N is even) 0?n?(N/2)-1 奇数序列 the range: 1?2n+1?N-1 (N is even) 0?2n?N-2 0?n?(N/2)-1  由于 所以 上式说明,按n的奇偶性将 分解为两个N/2长的序列 和 ,则N点DFT可分解为两个N/2点的DFT来计算. 其中 N/2-point DFT: k=0,1,…,N/2-1 因此，可写出两个N/2点的方程: 而 同理 而 所以  表示上述算法可用蝶形结（ butterfly） Example : 如N=4 x(n)=x0, x1, x2, x3 even: x0, x2 even: x0 odd: x2 odd: x1, x3 even: x1 odd: x3 bit reversal/shuffling process  the butterfly（蝶形运算） 如N4，则 将 x1(r) 再按r的奇偶进一步分解成两个 N/4点长的子序列：  其中 由X3（k）和X4（k）的周期性和 的对称性，得 同理，得 其中 8点DIT-FFT运算流图 4.2.3 IDFT的高效算法 IDFT的运算规律 求IFFT，也可用DIT-FFT的流程来实现。 Example: Determine the x(n) by IDFT. X=[5, -1, 1, -1] Solution: n=0,1,2,3 所以 x(n)=[1, 1, 2, 1] %the program in MATLAB: X=input(‘Please input X=’)