数字信号处理原理与实现（第3版） 课件汇总 刘泉 5.1 直接计算DFT的问题及改进的基本途径---8.5频率采样法设计matlab .pptx

2快速傅利叶变换

5.1直接计算DFT的问题及改进的基本途径3N点有限长序列x(n)正变换：5.1.1直接计算DFT的运算量反变换:

5.1直接计算DFT的问题及改进的基本途径4?运算量?

5.1直接计算DFT的问题及改进的基本途径5??????????

5.1直接计算DFT的问题及改进的基本途径6????算法的基本思想：利用??????系数的特性，合并??????运算中的某些项，把长序列??????→短序列??????，从而减少其运算量。FFT算法分类:时间抽选法DIT:Decimation-In-Time频率抽选法DIF:Decimation-In-Frequency

5.1直接计算DFT的问题及改进的基本途径7?

5.1直接计算DFT的问题及改进的基本途径8所以带入DFT中：??展开：

5.1直接计算DFT的问题及改进的基本途径9?有：合并有些项????

5.1直接计算DFT的问题及改进的基本途径10结论由此找出其它各项的类似归并方法：乘法次数可以减少一半。例：?

5.1直接计算DFT的问题及改进的基本途径112、将长序列DFT利用对称性和周期性分解为短序列DFT--思路因为DFT的运算量与N2成正比的如果一个大点数N的DFT能分解为若干小点数DFT的组合，则显然可以达到减少运算工作量的效果。

5.1直接计算DFT的问题及改进的基本途径122、将长序列DFT利用对称性和周期性分解为短序列DFT--方法把N点数据分成二半：其运算量为：再分二半：这样一直分下去，剩下两点的变换。

5.1直接计算DFT的问题及改进的基本途径132、将长序列DFT利用对称性和周期性分解为短序列DFT--结论快速付里时变换(FFT)就是在此特性基础上发展起来的，并产生了多种FFT算法，其基本上可分成两大类：按抽取方法分：时间抽取法(DIT);频率抽取法(DIF)按“基数”分：基-2FFT算法；基-4FFT算法；混合基FFT算法；分裂基FFT算法

5.2按时间抽选的基-2FFT算法(DIT-FFT)165.2.1算法原理设序列点数N=2L，L为整数。若不满足，则补零N为2的整数幂的FFT算法称基-2FFT算法。将序列x(n)按n的奇偶分成两组：

5.2按时间抽选的基-2FFT算法(DIT-FFT)17则x(n)的DFT:??????

5.2按时间抽选的基-2FFT算法(DIT-FFT)18再利用周期性求X(k)的后半部分和是以N/2为周期的

5.2按时间抽选的基-2FFT算法(DIT-FFT)19蝶式计算结构也即为蝶式信号流图上面频域中前/后半部分表示式可以用蝶形信号流图表示。X1(k)X2(k)作图要素：(1)左边两路为输入(2)右边两路为输出(3)中间以一个小圆表示加、减运算（右上路为相加输出、右下路为相减输出)(4)如果在某一支路上信号需要进行相乘运算，则在该支路上标以箭头，将相乘的系数标在箭头旁。(5)当支路上没有箭头及系数时，则该支路的传输比为1。

5.2按时间抽选的基-2FFT算法(DIT-FFT)20

5.2按时间抽选的基-2FFT算法(DIT-FFT)21分解后的运算量：复数乘法复数加法一个N/2点DFT(N/2)2N/2(N/2–1)两个N/2点DFTN2/2N(N/2–1)一个蝶形12N/2个蝶形N/2N总计运算量减少了近一半

5.2按时间抽选的基-2FFT算法(DIT-FFT)22N/2仍为偶数，进一步分解：N/2N/4

5.2按时间抽选的基-2FFT算法(DIT-FFT)23同理：其中：

5.2按时间抽选的基-2FFT算法(DIT-FFT)24图5.4DIT-FF二次分解

5.2按时间抽选的基-2FFT算法(DIT-FFT)25这样逐级分解，直到2点DFT当N=8时，即分解到X3(k)，X4(k)，X5(k)，X6(k)，k=0,1

5.2按时间抽选的基-2FFT算法(DIT-FFT)26

5.2按时间抽选的基-2FFT算法(DIT-FFT)275.2.2DIT-FFT的运算量当N=2L时，共有L级蝶形，每级N/2个蝶形，每个蝶形有1次复数乘法2次复数加法。复数乘法：复加次数：比较DFT:

5.2按时间抽选的基-2FFT算法(DIT-FFT)305.2.3DIT-FFT算法的特点1.原位计算m表示第m级迭代，k，j表示数据所在的行数5.6按时间抽取蝶形运算结构

5.2按时间抽选的基-2FFT算法(DIT-FFT)31例：N=8FFT运算,输入:看出：用原位运算结构运算后，A(0)