数字信号处理4 FFT.pptxVIP

第9次 @ 6.3回顾：DFS 和 DFT周期卷积，圆周位移，圆周卷积快速傅里叶变换（Fast Fourier Transform，FFT）： 按时间抽取基-2快速傅立叶变换（DIT-FFT）DIT-FFT 原理DIT-FFT 结构DIT-FFT 计算量回顾：DFS和DFT周期卷积圆周位移圆周卷积第五章，快速傅里叶变换(FFT)直接算DFT：计算量问题2. FFT类型I（DIT-FFT）3. FFT类型II（DIF-FFT）4. FFT的应用重难点： 计算量分析，FFT运算流图绘制，快速卷积法直接DFT任务已知：x(n), N 求：X(k)公式：程序长度N序列DFT的计算量复数乘法：N2 次复数加法：N(N-1) 次.FFTFFT的含义：FFT是DFT的算法实现形式，是对DFT直接计算过程进行简化而得到的快速算法。FFT算法有“按时间抽取(DIT)”和“按频率抽取(DIF)”两大类基础：旋转因子( )特性对称性周期性可约性周期内对称性旋转因子空间分布DIT-FFT：按时间抽取基2快速傅立叶变换公式推导DIT-FFT公式推导前半部分DIT-FFT公式推导后半部分DIT-FFT图形表达：DIT的蝶形2个长N/2的序列的DFT 合并为 1个长N的序列的DFTDIT-FFT计算实例求成分组合DIT-FFT构造运算流图DIT-FFT进一步分解DIT-FFT2级分解后的DFT运算DIT-FFT2点DFT的计算及图形表达x(0)X(0)0WNx(4)X(1)- -10WNx(2)X(2)--102WWNNx(6)X(3)- -1--10WNx(1)X(4)--101WWNNx(5)X(5)--1--102WWNNx(3)X(6)--1--1023WWWNNNx(7)X(7)--1--1--1DIT-FFT8点DFT的完整运算流图DIT-FFTFFT的结构23个点的FFT共3级运算每级23/2个蝶形长度N序列DFT的计算量每个蝶形：1次复数乘法，2次复数加法蝶形个数：log2N×N/2乘法总量：log2N×N/2加法总量：log2N×NDIT-FFT倒位序的求法自然顺序（I） 二进制数 倒位序二进制数 倒位序（J） 0123456700000101001110010111011100010001011000110101111104261537DIT-FFTDIT运算流图的绘制Step1：N条平行线Step2：输出原位序排列（同址运算）； 输入倒位序排列（变址运算）Step3：绘制蝶形蝶形的层数各层蝶形的形状蝶形的参数第10次 @ 6.8回顾： DFT直接计算：计算量第1种FFT：DIT-FFT原理， DIT-FFT蝶形图DIT-FFT运算流图：结构，计算量本次：第2种FFT： DIF-FFT原理，DIF-FFT蝶形图DIF-FFT运算流图：结构，计算量FFT应用：谱分析、快速卷积DIT-FFT：按时间抽取基2快速傅立叶变换问题：N点序列的DFT求解：方法计算量直接计算：奇偶分解：完全分解：?合并?DIF-FFT：按频率抽取基2快速傅立叶变换公式推导乘法结合律降低计算量DIF-FFT公式推导DIF-FFT公式推导偶半部分奇半部分DIF-FFT图形表示: DIF的蝶形1个长N的序列分解成2个长N/2的序列DIF-FFTDIF与DIT对比蝶形1：X1(k)、X2(k)合并成X(k) 蝶形2 ：x(n)分解成x1(n)、x2(n)抽取=奇偶分列DIT：x(n)奇偶X(k)前后DIF：x(n)前后X(k)奇偶回顾：FFT原理基本思想：改变计算顺序，降低计算量直接DFT：FFT：x(0)X(0)0WNx(4)X(1)- -10WNx(2)X(2)--102WWNNx(6)X(3)- -1--10WNx(1)X(4)--101WWNNx(5)X(5)--1--102WWNNx(3)X(6)--1--1023WWWNNNx(7)X(7)--1--1--1回顾： DIT-FFT8点DFT的完整运算流图DIF-FFT运算流图DIF流图与DIT流图对称DIF-FFTDIF运算流图的绘制Step1：N条平行线Step2：输入原位序排列（同址运算）； 输出倒位序排列（变址运算）Step3：绘制蝶形蝶形的层数各层蝶形的形状蝶形的参数IFFT含义用于IDFT的FFT运算，称为快速傅里叶反变换，简称IFFT。算法1 思考：IFFT的运算流图与FFT的有什么差别？IFFT算法2原理 步骤 IFFT算法1，需要对FFT内部代码进行修改；IFFT算法2，只需要在FFT前后补充代码。FFT应用谱分析低通滤波AD转换(采样)加窗处理(截取)DFTFFT应用谱分析Step1.算参数Step2.求V(k)Step3.分析V(k)幅度谱相位谱功率谱例5.6.1 有一频谱分析用的