dsp的IIR滤波器设计.pptVIP

IIR滤波器设计 切比雪夫I型滤波器 [z,p,k]=cheb1ap(N,Rp) [num,den]=cheby1(N,Rp,Wn,’s’) [num,den]=cheby1(N,Rp,Wn,’type’,’s’) [N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,’s,’) 切比雪夫II型滤波器 [z,p,k]=cheb2ap(N,Rp) [num,den]=cheby2(N,Rp,Wn,’s’) [num,den]=cheby2(N,Rp,Wn,’type’,’s’) [N,Wn]=cheb2ord(Wp,Ws,Rp,Rs,’s,’) 椭圆滤波器 [z,p,k]=ellipap(N,Rp,Rs) [num,den]=ellip(N,Rp,Rs,Wn,’s’) [num,den]=ellip(N,Rp,Rs,Wn,’type’,’s’) [N,Wn]=ellipord(Wp,Ws,Rp,Rs,’s,’) 贝塞尔滤波器滤波器 [z,p,k]=besselap(N) [num,den]=besself(N,Wn) [num,den]=besself(N,Wn,’type’) 将模拟原型传输函数Ha(s)变换成所需的数字IIR传输函数G(z)的基本思路就是要把s域映射到z域，从而使数字滤波器能保持模拟滤波器的基本频率响应特性，因此，这种映射函数必须满足以下要求： s平面的虚轴（jΩ）必须映射到z平面的单位圆上。 稳定的模拟传输函数能变换为稳定的数字传输函数。 IIR滤波器设计的双线性变换法 有很多种变换方法可以将一个模拟传输函数Ha(s)变换成一个数字传输函数G(z)，从而使z域的数字传输函数保留s域的模拟传输函数的基本性质 . 在这些变换中，更多地使用双线性变换法来设计基于模拟原型滤波器变换的IIR数字滤波器。 双线性变换法(Bilinear transformation )：是从频域出发，使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换法。 双线性变换定义： 以上变换是一个一一映射，它将s平面上的一点映射为z平面上的一点，或将z平面上的一点映射为s平面上的一点。数字传输函数G(z)和原型模拟传输函数 Ha(s)之间的关系为 : 双线性变换是通过应用梯形数值积分方法来从Ha(s)的微分方程得到G(z)的差分方程的一种变换。参数T表示数值积分的步长 双线性变换的推导（梯形积分） 双线性变换主要用在: 用反双线性变换来将数字滤波器的性能指标转换为模拟原型滤波器的性能指标 用双线性变换来从模拟传输函数得到所希望的数字滤波器的传输函数G(z) 参数T对 G(z)的表达式没有影响, 可以选择T=2来简化设计的过程 0 ? j? Im(z) Re(z) -1 1 0 s域到z域的映射关系： 令： ? ? 双线性变换扭曲了数字频率和模拟频率的关系。 ? ? ? 稳定性不变，但形状发生变换（频率轴有畸变）； 映射高度非线性，由此引起的频率轴失真称为频率弯折 例: 模拟滤波器传输函数为 -3dB频率为2000弧度/秒(318.31Hz)，求对应的数字滤波器传输函数，采样率为1500Hz 频率响应之间的差异见上图，数字滤波器的滚降比模拟滤波器要陡，且带宽也不同，这种误差是由于双线性变换的扭曲造成的。扭曲误差可以通过对-3dB频率的预扭曲来克服。 例： 模拟传输函数为 求数字滤波器，-3dB频率为318.3Hz，采样率为1500Hz 二阶模拟陷波器的传输函数为 IIR陷波器的设计 令： 陷波频率?0 和3-dB陷波带宽 Bw 与常数? 和 ? 相关 ?/? ?/? 设计一个低通IIR数字滤波器G(z)，性能指标如下: 低通IIR滤波器设计 将数字截止频率预弯折为模拟截止频率 三阶归一化低通Butterworth传输函数: 第一种方法分为以下几个步骤： 高通、带通和带阻IIR数字滤波器设计 Step 1：用式 预畸所求数字滤波器GD(z)的数字频率指标，从而得到一个等价的模拟滤波器HD(s) 的频率指标。 Step 2：选取一种合适的频率变换，将HD(s)的频率指标转换成原型低通滤波器HLP(s)的频率指标。 Step3:设计模拟低通滤波器HLP(s). Step4: 用步骤2中频率变换的逆变换将传输函数HLP(s)转换为HD(s) . Step5:对传输函数HD(s)进行如式(9.14)所示的双线性变换，从而得到所求的数字IIR传输函数GD(z) 第二种方法分为以下几个步骤： 步骤1：用式 预畸所求数字滤波器GD(z)的数字频率指标，从而得到一个等价的模拟滤波器HD(s)的频