DSP第四章讲解.ppt

x3(0) x3(1) -1 -1 x4(0) x4(1) N/4点 DFT N/4点 DFT x1(0) x1(1) x1(2) x1(3) X3(0)=X1(0)=X(0) X4(0)=X1(1)=X(2) X3(1)=X1(2)=X(4) X4(1)=X1(3)=X(6) 同理： 其中： 逐级分解，直到2点DFT 当N=8时，即分解到x3(n)，x4(n)，x5(n)，x6(n)，n=0,1 2、算法特点 1）原位计算 -1 L级蝶形运算，每级N/2个蝶形，每个蝶形结构： m表示第m级迭代，k，j表示数据所在的行数 2）蝶形运算 对N=2L点FFT，输入自然序，输出倒位序， 两节点距离：2L-m=N / 2m 第m级运算： 蝶形运算两节点的第一个节点为k值，表示成L位二进制数，左移m-1位，把右边空出的位置补零，结果为r的二进制数。 3、DIT与DIF的异同 基本蝶形不同 DIT: 先复乘后加减 DIF: 先减后复乘 运算量相同 都可原位运算 DIT和DIF的基本蝶形互为转置 四 、IFFT算法 比较： IDFT: DFT: 共轭 FFT 共轭 乘1/ N 直接调用FFT子程序计算IFFT的方法： 五、线性卷积的FFT算法 1、线性卷积的FFT算法 需运算量： 若系统满足线性相位，即： 则需运算量： 若L点x(n)，M点h(n)， 则直接计算其线性卷积y(n) FFT法：以圆周卷积代替线性卷积 1) H(k) = FFT [h(n)] N /2*log2N 4) y(n) = IFFT [Y(k)] N /2*log2N 3) Y(k) = H(k)X(k) N 2) X(k) =FFT [x(n)] N /2*log2N N 比较直接计算和FFT法计算的运算量 讨论： 1）当 2）当 1）重叠相加法 N 2）重叠保留法 舍弃yi(n)的前M-1个点，再将yi(n)顺次连接， 即得y(n)。 分段 右移序列 卷积 N 本章作业练习 P194: 1 2 * 离散时间信号—序列 * 常系数线性差分方程 第一级蝶形系数为WN0，两节点距离为1 第二级蝶形系数为WN0、WNN/4，两节点距离为2 第三级蝶形系数为WN0、WNN/8、WN3N/8，两节点距离为4 第M级蝶形系数为WN0、WN1、……、WN(N/2-1)，两节点距离为N/2 * * 常系数线性差分方程 * 常系数线性差分方程 * 常系数线性差分方程 第四章学习目标 掌握按时间抽选的基-2FFT算法的算法原理、运算流图、所需计算量和算法特点 掌握按频率抽选的基-2FFT算法的算法原理、运算流图、所需计算量和算法特点 理解IFFT算法 理解线性卷积的FFT算法及分段卷积方法 第四章 快速傅里叶变换 FFT: Fast Fourier Transform 1965年，Cooley, Tukey 《机器计算傅里叶级数的一种算法》 一、直接计算DFT的问题及改进途径 运算量 复数乘法 复数加法 一个X(k) N N – 1 N个X(k) (N点DFT) N 2 N (N – 1) 实数乘法 实数加法 一次复乘 4 2 一次复加 2 一个X (k) 4N 2N+2 (N – 1)=2 (2N – 1) N个X (k) (N点DFT) 4N 2 2N (2N – 1) FFT算法分类: 时间抽选法 DIT: Decimation-In-Time 频率抽选法 DIF: Decimation-In-Frequency 二 、按时间抽选的基-2FFT算法 1、算法原理 设序列点数 N = 2L，L 为整数。 若不满足，则补零 将序列x(n)按n的奇偶分成两组： N为2的整数幂的FFT算法称基-2FFT算法。 则x(n)的DFT: 再利用周期性求X(k)的后半部分 分解后的运算量： 复数乘法 复数加法 一个N / 2点DFT (N / 2)2 N / 2 (N / 2 –1) 两个N / 2点DFT N 2 / 2 N (N / 2 –1) 一个蝶形 1 2 N / 2个蝶形 N / 2 N 总计 运算量减少了近一半 N / 2仍为偶数，进一步分解：N / 2 N / 4 同理： 其中： 这样逐级分解，直到2点DFT 当N = 8时，即分解到X3(k)，X4(k)，X5(k)，X6(k)，k = 0, 1 2、运算量 当N = 2L时，共有L级蝶形，每级N / 2个蝶形，每个蝶形有1次复数乘法2次复数加法。 复数乘法： 复数加法： 比较DFT 3、算法特点 1）原位计算 m表示第m级迭代，k，j表示数据所在的行数 2）倒位序