05数字滤波器结构..pptVIP

数字信号处理 5 数字滤波器结构 滤波器表示 分类 IIR滤波器 FIR滤波器 引言 引言 数字滤波 —— 输入、输出均为时域离散信号； 利用离散系统特性，对输入信号进行加工和变换，改变输入序列的频谱或信号波形，保留有用频率的信号分量，抑制无用的信号分量输入； 换句话说，通过一定运算关系，以改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的算法。 因此，从概念上说，数字滤波与模拟滤波相同，只是信号的类型和实现方法不同。 5.1 滤波器表示 5.1 滤波器表示 5.1 滤波器表示 5.1 滤波器表示 不同的信号流图代表不同的运算方法，对于同一个系统函数，可以有很多种信号流图与之对应； 运算结构是重要的，不同结构所需的存储单元及乘法次数是不同的，前者影响复杂性，后者影响运算速度； 另外，在有限精度（有限字长）情况下，不同运算结构的误差稳定性等都是不同的，因此根据不同的使用要求，选择不同的计算结构是有实际意义的。 5.2 分类 滤波器的种类很多，分类方法也不同； 从功能上分：低通、高通、带通、带阻； 从实现方法上分：FIR、IIR； 从设计方法上来分：切比雪夫、 巴特沃斯； 从处理信号分：经典滤波器、现代滤波器。 5.2 分类 1. 经典滤波器 假定输入信号x(n)中的有用成分和希望去除的成分，各自占有不同的频带。当x(n)经过一个线性系统（即滤波器）后即可将欲去除的成分有效地去除。 若信号和噪声的频谱相互重叠，那么经典滤波器将无能为力。 5.2 分类 5.2 分类 IIR滤波器： —— 系统的单位抽样响应h(n)延伸至无限长。 系统函数 至少存在一个不能被零点抵消的极点，即在有限长Z平面（0|Z|∞)存在极点，换句话说，分母多项式不能被完全化简掉。不难推出：因果稳定的IIR滤波器其全部极点一定在单位园内； 从时域差分方程来看，至少有一个ak不为0，即输出不但与输入有关，而且还与输出的位移及其加权值有关。换句话说，存在着输出对输入的反馈，所以，这种结构有时也被称为“递归结构”。 5.2 分类 IIR系统按照其有、无零点，可分为： 全极点系统：H(z)分子只有b0一项，在有限Z平面上无零点，又称为自回归系统（AR系统）； 零极点系统：既有零点又有极点，又称为自回归滑动平均系统（ARMA系统）。 5.2 分类 FIR滤波器： —— 系统的单位抽样响应h(n)有限长。 无极点系统：所有的分母多项式系数ak均为0，此时系统即为FIR系统。由于没有极点，故系统又称为全零点系统或滑动平均系统（MA系统）。 FIR系统既可采用没有反馈的“非递归型结构”，也可采用“递归型结构”(利用零极点相互抵消的办法实现) 5.3 IIR滤波器 5.3 IIR滤波器 5.3 IIR滤波器 特点： 两个网络级联：第一个横向结构M节延时网络实现零点，第二个有反馈的N节延时网络实现极点； 共需(N+M)级延时单元； 系数ai、bi不是直接决定单个零极点，因而不能很好地进行滤波器性能控制； 极点对系数的变化过于灵敏，从而使系统频率响应对系统变化过于灵敏，也就是对有限精度（有限字长）运算过于灵敏，容易出现不稳定或产生较大误差。 直接型 —— 函数 P130 filter函数 格式： 例 1 —— 应用函数 例 2 —— P129 【例5.3.1】 设IIR数字滤波器的系统函数H(z)为： 例 2 —— P129 【例5.3.1】 5.3 IIR滤波器 例 3 —— P131 【例5.3.2】 设系统函数H(z)如下式： 级联型 —— 函数 1. 设计函数：直接型→级联型 function [b0,B,A] = dir2cas(b,a); b0 = b(1); b = b/b0; a0 = a(1); a = a/a0; b0 = b0/a0; M = length(b); N = length(a); if N M b = [b zeros(1,N-M)]; % 补0 elseif M N a = [a zeros(1,M-N)]; % 补0，使a、b等长 N = M; else NM = 0; end 级联型 —— 函数 K = floor(N/2); B = zeros(K,3); A = zeros(K,3); if K*2 == N; b = [b 0]; a = [a 0]; end broots = cplxpair(roots(b)); % 共轭复根对 aroots = cplxpair(roots(a)); for i=1:2:2*K Brow =