《数字信号处理及应用第2版》第4章数字滤波网络g.pptVIP

第 4 章 数字滤波网络 数字信号处理及应用 主要内容 4.1 数字滤波器的表示方法 4.2 IIR数字滤波器的基本网络结构 4.3 FIR数字滤波器的基本网络结构 4.4 数字滤波器基本结构的 MATLAB实现 4.1 数字滤波器的表示方法 数字滤波器的功能本质上是将一组输入数字序列通过一定的运算后转变为另一组输出数字序列。其结构表示需要采用加法器，乘法器，延时器等基本运算单元。这些基本单元一般用方框图或信号流图来表示，分别见图4.1所示。 例如，一阶数字滤波器的差分方程可表示为 用方框图和信号流图分别表示为图4.2和图4.3。 图中，1、2、3、4、5、6称为网络节点；x(n)称为输入节点或源节点，表示外部输入或信号源；y(n)称为输出节点或阱节点。节点之间用有向支路连接，每个节点可以有多条输入支路和多条输出支路，任一节点的节点值等于其所有输入支路的信号和。输入支路的信号值等于该支路起点处的节点信号值乘以支路上的传输系数（支路增益），如果支路不标传输系数值，则认为传输系数为1。延迟支路用延迟算子z-1表示，它表示单位延迟。由此可得图4.3的各节点值为 4.2 IIR数字滤波器的基本网络结构 IIR 数字滤波器的结构分类 4.2.1 直接型 对于一个N阶无限长单位脉冲响应（IIR）滤波器，其传递函数可表示为 用差分方程表示为 直接画出差分方程对应系统的信号流图结构称为直接 I型结构。将直接 I型结构中的延时单元z-1尽可能减少的一种流图结构，称为直接Ⅱ型结构。 对于Ⅰ型结构，x(n)和y(n)的延时器相互独立，共需要2N个延时器，其结构如图4.4所示。 对于Ⅱ型结构，交换x(n)和y(n)的级联次序，并将延时器相互合并，只需要N个延时器，其结构如图4.5所示。 4.2.2 级联型 一个N阶的传递函数也可以用它的零、极点来表示 ∵ an 、bn 为实数，∴cn 、dn 均以共轭成对出现，故可将共轭成对的零、极点配成二次有理分式形式： 式中，Hk(z)称为滤波器的二阶基本节，其传递函数的一般形式为 整个滤波器可以看作是二阶基本节的级联，其结构如图4.6所示。 4.2.3 并联型 将传递函数展开成部分分式就可以用并联方式构成滤波器。 对于其中的共轭复根部分，将其合并为二阶实系数分式 式中，L=N+2M。 4.3 FIR数字滤波器的基本网络结构 FIR数字滤波器 有限长单位脉冲响应FIR滤波器的单位脉冲响应是一个有限长序列。其传递函数一般具有如下形式 其差分方程可描述为 由差分方程可以看出，FIR滤波器的阶数为N-1，长度为N。系统输出只是输入的函数，而与过去的输出无直接关系，不含有反馈支路。 4.3.1 直接型 直接型FIR数字滤波器是以给定的形式直接实现差分方程式（4.12），其结构如图4.8所示，稍加改变也可以形成图4.9的结构。 4.3.2 级联型 级联型FIR数字滤波器是把式（4.9）表示的系统函数分解成二阶子系统的乘积，各子系统以直接形式实现，而整个滤波器用二阶子系统的级联形式实现。 将式（4.9）分解为二阶实系数因子 式中，M=N/2，a0i、a1i、a2i表示二阶子系统的实系数。用级联形式实现式（4.13）的传递函数见图4.10。 4.3.3 线性相位型 线性相位型FIR数字滤波器的结构根据其单位脉冲响应的形式可分为偶对称和奇对称两种，其单位冲激响应h(n)具有如下特性 因此，当N为偶数时 当N为奇数时 式中，“+”表示为偶对称；“-”表示为奇对称。 上两式意味着线性相位型比直接型少用N/2个乘法器或乘法运算。其网络结构可用图4.11和图4.12的信号流图来表示。 4.4 数字滤波器基本结构的MATLAB实现 4.4.1 IIR滤波器的MATLAB实现 【例4.1】 用直接型实现系统函数为 的IIR数字滤波器，求单位脉冲响应和单位阶跃信号。 【解】 程序如下: b=[1,-4,10,22,18]; a=[14,11,2,-3,-1]; N=40; h=impz(b,a,N); %直接型单位脉冲相应 n=[-1:30];x=[n=0]; %单位阶跃信号 y=filter(b,a,x);