DSP第四章2第五章1.ppt

FFT算法分类: 时间抽选法 DIT: Decimation-In-Time 频率抽选法 DIF: Decimation-In-Frequency 二 、按时间抽选的基-2FFT算法 1、算法原理 设序列点数 N = 2L，L 为整数。 若不满足，则补零 则x(n)的DFT: 再利用周期性求X(k)的后半部分 分解后的运算量： 这样逐级分解，直到2点DFT 当N = 8时，即分解到X3(k)，X4(k)，X5(k)，X6(k)，k = 0, 1 2、运算量 当N = 2L时，共有L级蝶形，每级N / 2个蝶形，每个蝶形有1次复数乘法2次复数加法。 3、算法特点 1）原位计算 2）倒位序 3）蝶形运算 对N = 2L点FFT，输入倒位序，输出自然序， 第m级运算每个蝶形的两节点距离为 2m–1 第m级运算： 蝶形运算两节点的第一个节点为k值，表示成L位二进制数，左移L – m位，把右边空出的位置补零，结果为r的二进制数。 4）存储单元 输入序列x(n) : N个存储单元 4、DIT算法的其他形式流图 输入倒位序输出自然序 输入自然序输出倒位序 输入输出均自然序 相同几何形状 输入倒位序输出自然序 输入自然序输出倒位序 三 、按频率抽选的基-2FFT算法 1、算法原理 按k的奇偶将X(k)分成两部分： 令 N /2仍为偶数，进一步分解：N /2 N /4 逐级分解，直到2点DFT 2、算法特点 1）原位计算 2）蝶形运算 对N=2L点FFT，输入自然序，输出倒位序， 两节点距离：2L-m=N / 2m 蝶形运算两节点的第一个节点为k值，表示成L位二进制数，左移m-1位，把右边空出的位置补零，结果为r的二进制数。 3、DIT与DIF的异同 基本蝶形不同 四 、IFFT算法 比较： 九、线性卷积和线性相关的FFT算法 1、线性卷积的FFT算法 FFT法：以圆周卷积代替线性卷积 比较直接计算和FFT法计算的运算量 1）重叠相加法 2）重叠保留法 2、线性相关的FFT算法 若L点x(n)，M点y(n)，计算线性相关： 思考题： 1）频谱分析 第五章学习目标 理解数字滤波器结构的表示方法 掌握IIR滤波器的基本结构 掌握FIR滤波器的直接型、级联型、线性相位结构，理解频率抽样型结构 了解数字滤波器的格型结构 本章作业练习 P226: 1 2 3 4 6 7 8 (1) 第五章 数字滤波器的基本结构 一、数字滤波器结构的表示方法 数字滤波器的系统函数： 例：二阶数字滤波器 流图结构 节点 二、IIR数字滤波器的基本结构 1、直接Ⅰ型 差分方程: 2、直接Ⅱ型（典范型） 直接型的共同缺点： 3、级联型 将系统函数按零极点因式分解: 将共轭成对的复数组合成二阶多项式，系数即为实数。 为采用相同结构的子网络，也将两个实零点/极点组合成二阶多项式 级联型的特点： 4、并联型 将因式分解的H(z)展成部分分式： 并联型的特点： 转置定理： 原网络中所有支路方向倒转，并将输入x(n)和输出y(n)相互交换，则其系统函数H(z)不改变。 例：设IIR数字滤波器差分方程为： 试用四种基本结构实现此差分方程。 1）系统的单位抽样相应h(n)无限长 IIR数字滤波器的特点： 3）存在输出到输入的反馈，递归型结构 2）系统函数H(z)在有限z平面（ ）上有极点存在 IIR数字滤波器的基本结构： 直接Ⅰ型 直接Ⅱ型（典范型） 级联型 并联型 需N+M个 延时单元 只需实现N阶滤波器所需的最少的N个延时单元，故称典范型。（ ） 系数 ， 对滤波器的性能控制作用不明显 极点对系数的变化过于灵敏，易出现不稳定或较大误差 运算的累积误差较大 当零点为奇数时： 有一个 当极点为奇数时： 有一个 各二阶基本节的排列次序有 种 当M=N时，二阶因子配对方式有 种 调整系数 ， 能单独调整滤波器的第k对零点，而不影响其它零极点 运算的累积误差较小 具有最少的存储器 便于调整滤波器频率响应性能 调整系数 ， 能单独调整滤波器的第k对极点，而不影响其它零极点 当N为奇数时，有一个 组合成实系数二阶多项式： 通过调整系数 ， 可单独调整一对极点位置，但不能单独调整零点位置 各并联基本节的误差互相不影响，故运算误差最小 可同时对输入信号进行运算，故运算速度最高 x3(0) x3(1) -1 -1 x4(0)