第6章 无限脉冲响应数字滤波器的设计B精要.ppt

第6章 无限脉冲响应数字滤波器的设计 6.1 数字滤波器的基本概念 6.2 模拟滤波器的设计 6.3 用脉冲响应不变法设计IIR数字低通滤波器 6.4 用双线性变换法设计IIR数字低通滤波器 6.5 数字高通、带通和带阻滤波器的设计 经典滤波器从功能上分又可分为： 低通滤波器（LP):Low pass filter 带通滤波器(BP):Bandpass filter 高通滤波器(HP):High pass filter 带阻滤波器(BS):Bandstop filter 四.低通巴特沃斯滤波器的设计步骤 (1)根据技术指标Ωp,αp,Ωs和αs，求出滤波器的阶数N。 (2) 求截止频率Ωc (3)求出归一化极点pk，得到归一化传输函数Ha(p)。 (4)将Ha(p)去归一化。将p=s/Ωc代入Ha(p)，得到实际的滤波器传输函数Ha(s)。 MATLAB信号处理工具箱函数buttap, buttord和butter是巴特沃斯滤波器设计函数。 　(2) 由(6.2.44)和 (6.2.45)式求阶数N和ε： 　　(3) 将极点pk、N和ε代入(6.2.39)式求Ga(p)： 由(6.2.46)式求出N=5时的极点pi，代入上式，得到： 　　(4) 将Ga(p)去归一化，得到： 　　2. 用MATLAB设计切比雪夫滤波器 　　MATLAB信号处理工具箱函数cheb1ap，cheb1ord和cheby1是切比雪夫Ⅰ型滤波器设计函数。其调用格式如下: 　　1） ［z, p, k］= cheb1ap(N, Rp) 　　2） ［N, wpo］= cheb1ord(wp, ws, Rp, As) 　　3） ［N, wpo］= cheb1ord(wp, ws, Rp, As, ′s′) 　　4） ［B, A］= cheby1(N, Rp, wpo, ′ftype′) 　　5） ［B, A］= cheby1(N, Rp, wpo, ′ftype′, ′s′) 　　【例6.2.3】 设计切比雪夫Ⅰ型和切比雪夫Ⅱ型模拟低通滤波器。要求与例6.2.2相同。 　　解 设计程序ep623.m如下: ％例6.2.3设计程序: ep623.m ％设计切比雪夫Ⅰ型模拟低通滤波器 wp=2*pi*3000; ws=2*pi*12000; Rp=0.1; As=60; 　　％设置指标参数 ［N1, wp1］=cheb1ord(wp, ws, Rp, As, ′s′); 　　　　　　　　 　　％计算切比雪夫Ⅰ型模拟低通滤波器阶数和通带边界频率 ［B1, A1］=cheby1(N1, Rp, wp1, ′s′); 　　％计算切比雪夫Ⅰ型模拟低通滤波器系统函数系数 subplot(2, 2, 1); fk=0:12000/512:12000; wk=2*pi*fk; Hk=freqs(B1, A1, wk); plot(fk/1000, 20*log10(abs(Hk))); grid on xlabel(′频率(kHz)′); ylabel(′幅度(dB)′) axis(［0, 12, -70, 5］) 　　运行结果: 　　N=5 　　切比雪夫Ⅰ型模拟低通滤波器通带边界频率: wp1 =1.8850e+004 　　切比雪夫Ⅰ型模拟低通滤波器系统函数分子分母多项式系数: 　　B=1.2187e+011 　　A =［ 1 3.2873e+004 9.8445e+008 1.6053e+013 1.8123e+017 9.7448e+020 ］ 　　滤波器损耗函数如图6.2.9所示。 图6.2.9 五阶切比雪夫Ⅰ型模拟低通滤波器损耗函数 （例6.2.3的设计结果 6.2.4 椭圆滤波器的设计 　　椭圆（Elliptic）滤波器在通带和阻带内都具有等波纹幅频响应特性。由于其极点位置与经典场论中的椭圆函数有关，所以由此取名为椭圆滤波器。又因为在1931年考尔（Cauer）首先对这种滤波器进行了理论证明，所以其另一个通用名字为考尔（Cauer）滤波器。椭圆滤波器的典型幅频响应特性曲线如图6.2.10所示。 图6.2.10 椭圆滤波器幅频响应特性曲线 　　由图6.2.10(a)可见，椭圆滤波器通带和阻带波纹幅度固定时，阶数越高, 过渡带越窄； 由图6.2.10(b)可见，当椭圆滤波器阶数固定时，通带和阻带波纹幅度越小, 过渡带就越宽。所以椭圆滤波器的阶数N由通带边界频率Ωp、阻带边界频率Ωs、通带最大衰减? p和阻带最小衰减 ?s共同决定。后面对五种滤波器的比较将证实，椭圆滤波器可以获得对理想滤波器幅频响应的最好逼近，是一种性能