1、第6章 无限冲激响应(IIR） 滤波器设计,6.1 滤波器的基本概念；6.2 模拟低通滤波器设计；6.3 模拟高通、带通及带阻滤波器设计；6.4 冲激响应不变法；6.5 双线性Z变换法；6.6 数字高通、带通及带阻滤波器设计；,1. 滤波原理,6.1 滤波器的基本概念,若 中的有用成分 和希望去除的成分 各自占有不同的频带, 通过一个线性系统可将 有效去除.,分类: 低通(LP), 高通(HP),带通(BP), 带阻(BS),2. 滤波器的分类,加法性噪声,经典滤波器,每一种又有模拟(AF)、数字(DF)两种滤波器.对数字滤波器, 从实现方法上, 有IIR滤波器和FIR滤波器之分, 转移函数分
2、别为:,FIR DF:,IIR DF:,现代滤波器,维纳滤波器是这类滤波器的典型代表,种类：维纳滤波器、卡尔曼滤波器、线性预 测、自适应滤波器,乘法性噪声,卷积性噪声,信号的频谱和噪声道频谱混迭在一起，靠经典的滤波方法难以去除噪声。目标：从含有噪声的数据记录(又称时间序列)中估计出信号的某些特征或信号本身。,3. 滤波器的技术要求,低通：,单位 (dB),若幅度下降到 0.707, 则幅平方下降 0.5 （半功率点）：,若幅度下降到 0.01：,高通：,带通：,带阻：,1. 给定所设计的滤波器的技术指标：,LP, HP,BP, BS,2. 设计出一个 ，使,数字滤波器设计的一般步骤:,给定数字
3、滤波器的技术指标 （更多）,数字IIR滤波器设计的具体步骤：,6.2 模拟低通滤波器的设计,一、概述,上面两式将要求的衰减和模拟滤波器的幅平方特性联系了起来。注意，由于衰减指标只有两个，因此也只能和幅平方特性的两个特殊频率相联系。,我们的目的是由幅平方特性 求出模拟滤波器的转移函数 ，如何将二者联系起来,？,因为,所以，由,注意，由于由于衰减指标只有两个，只能和幅平方特性的两个特殊频率相联系，因此，幅平方特性的表达式必须简化。,将 按不同的原则简化，可得到不同形式的滤波器，即不同的 表达式：,1. 巴特沃思(Butterworth)滤波器,2. 切比雪夫I型(Chebyshev-I)滤波器,3
4、. 切比雪夫II型滤波器,1. 将实际频率归一化，得归一化幅平方特性,2. 求 和,二、Butterworth滤波器的设计,对Butterworth滤波器，通常 ，所以,3. 确定,则：,反映了实际频率,4.巴特沃思滤波器幅频响应的特点:,“最平”幅频 响应滤波器.,例：给定 如下技术指标，设计模拟低通 Butterworth滤波器,Step1.,Step2.,Step3.,Step4.,Step5.,切比雪夫多项式的特点：,三、切比雪夫I型模拟低通滤波器设计,的确是 的多项式,首项系数是,设计步骤：,2. 求 和,为求滤波器的阶次，还要利用另外的条件：,注意：,因此： 时，切比雪夫多项式要重
5、新定义，采用双曲函数：,利用另外的条件：,必须不大于1,则,最后导出：,切比雪夫滤波器的极点分布,反映了实际频率,最后：,模拟高通, 带通, 带阻滤波器设计流程,6.3 模拟高通、带通及带阻滤波器的设计,给定高通、带通或带阻的技术指标,想办法实现高通到低通的转换：,给定高通滤波器的技术指标：,先作频率归一化：,一、模拟高通滤波器的设计,或,由：,得：,于是可得到模拟高通滤波器的转移函数,后面带通、带阻滤波器和低通滤 波器的转换过程大体相同。,对带通滤波器，如何实现频率的归一化,？,定义：有 ：,用带宽归一化,二、模拟带通滤波器的设计,定义 中心频率：,由 于,所以有,N阶低通滤波器转换到带通后
6、, 阶次变为2N.,带阻滤波器频率归一化方法同带通滤波器,定义：有 ：,？,三、模拟带阻滤波器的设计,得到：,及带阻滤波器的转移函数：,以上讨论的是模拟低通、高通、带通及带阻滤波器的设计，然而这并不是我们的目的。我们的目的是设计数字滤波器。首要的问题是如何将数字滤波器的技术指标转换为模拟滤波器的技术指标，最后再实现模拟滤波器到数字滤波器的转换。,6.4 用冲激响应不变法设计 IIR DF,给定数字滤波器的技术指标 （更多）,转换成模拟滤波器的技术指标 （更多）,转换成模拟低通滤波器的技术指标,设计模拟低通滤波器,得到模拟高通、带通、带阻滤波器,得到数字高通、带通、带阻滤波器,利用上一节的方法，
7、可设计出模拟滤波器,最直接到方法，将：,利用：,？,但这样做， 将不再是 的有理多项式，给极零分析带来困难。,令：,基本转换单元：,步骤:,缺点:,？,线性转换关系,6.5 用双线性Z变换法设计 IIR DF,放弃上一节的线性转换关系，找新的关系：,令：,双线性z变换,非线性关系，但是一对一的转换,非线性关系,设计的 AF 并不是按给定的技术指标，但再由 变回 后，保证了 DF的技术要求。,又称为频率的预变形（Freq. Warping) 。例如 ：,Step1.,Step2. 频率转换：,Step3.,Step4.,数字低通滤波器的设计步骤：,所以：,这样：系数 可以省略，因此，双线性Z变换
8、可定义为 ：,这一组定义和前面的定义，对最后的 DF 而言，结果是一样的，差别是中间设计的 AF， 由于缺少了频率定标，将无法给出符合要求的幅频响应。,可求出：,求出：,幅频响应曲线见书 P284,由：,数字高通滤波器设计步骤,6.6 数字高通, 带通及带阻滤波器的设计,对 带通（BP）、带阻（BS）数字滤波器的设计，只需改变图中 Step2 和 Step4：,要求：,按上述转换办法，可以求出：,1buttord.m 确定 LP DF、或 LP AF的阶次；（1） N, Wn = buttord(Wp, Ws, Rp, Rs)；,与本章内容有关的MATLAB文件,对应 数字滤波器。其中 Wp,
9、 Ws分别是通带和阻带的截止频率，其值在 01 之间，1对应抽样频率的一半(归一化频率)。对低通和高通，Wp, Ws都是标量，对带通和带阻，Wp, Ws是12的向量。Rp, Rs 分别是通带和阻带的衰减(dB)。N是求出的相应低通滤波器的阶次，Wn是求出的3dB频率，它和Wp稍有不同。,（2）N, Wn = buttord(Wp, Ws, Rp, Rs,s)： 对应模拟滤波器，式中各个变量的含意和格式（1）相同，但Wp, Ws及Wn的单位为弧度/秒，因此，它们实际上是频率。,2buttap.m 设计模拟低通(Butt)原型滤波器。 z, p, k=buttap(N): N是欲设计的低通原型滤波
10、器的阶次，z, p, k是设计出的极点、零点及增益。,B, A=lp2lp(b, a, Wo), B, A=lp2hp(b, a, Wo) B, A=lp2bp(b, a, Wo, Bw), B, A=lp2bs(b, a, Wo, Bw) b, a 是AF LP 的分子、分母的系数向量，B, A是转换后的的分子、分母的系数向量；在（1）中，Wo是低通或高通滤波器的截止频率；在（2）中，Wo是带通或带阻滤波器的中心频率，Bw是其带宽。,3lp2lp.m、lp2hp.m、lp2bp.m, lp2bs.m,将模拟低通原型转换为实际的低通、高通、带通及带阻滤波器。,（1）,（2）,4bilinear
11、.m ：双线性变换，由模拟滤波器 得到数字滤波器。,Bz, Az=bilinear(B, A, Fs) 式中B, A分别是G（s）的分子、分母多项式的系数向量，Bz, Az分别是H(z)的分子、分母多项式的系数向量，Fs是抽样频率。,5butter.m 本文件可用来直接设计Butterworth数字滤波器，实际上它把 Buttord buttap、 lP2lp Bilinear等文件都包含了进去，从而使设计过程更简捷。,格式（1）（3）用来设计数字滤波器，B，A分别是H（z）的分子、分母多项式的系数向量，Wn是通带截止频率，范围在01之间。若Wn是标量，（1）用来设计低通数字滤波器，若Wn是1
12、2的向量，则（1）用来设计数字带通滤波器；（2）用来设计数字高通滤波器；（3）用来设计数字带阻滤波器，显然，这时的Wn是12的向量；格式（4）用来设计模拟滤波器。,（1） B,A=butter(N,Wn); （2） B,A=butter(N,Wn,high); （3） B,A=butter(N,Wn,stop); （4）B,A=butter(N,Wn,s)。,6cheb1ord.m 求Cheb-型滤波器的阶；7.cheb1ap.m 设计原型低Cheb-I型模拟滤波器；8cheby1.m 直接设计数字Cheb-滤波器。,以上三个文件的调用格式和对应的Butterworth滤波器的文件类似。,9cheb2ord.m; 10. ellipord.m; 11.cheb2ap.m; 12. ellipap.m; 13.besselap.m; 14. cheby2.m; 15. ellip.m; 16.besself.m,16.impinvar.m 用冲激响应不变法实现频率转换；,对应 Cheby-II、椭圆 IIR 滤波器,
