IIR滤波器.ppt

6.1 滤波器的基本概念 1、数字滤波器概念 指输入，输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些成分的器件。 现代滤波器种类：维纳滤波器、卡尔曼滤波器、自适应滤波器。 7、典型模拟滤波器 1）巴特沃斯滤波器 2）切比雪夫滤波器 3）椭圆滤波器 例1：N=3阶Butterworth 滤波器振幅平方函数的极点分布如图，求其系统函数。 解： 由前面分析计算出6极点如图 考虑到系统的稳定性， Butterworth 滤波器的系统函数 是由s平面左半部分的极点 （SP3，SP4，SP5）组成的， 它们分别为： 如用 归一化s，即s’=s/Ωc，得归一化的三阶BF： 例2 已知滤波器通带截止频率为5kHz，通带最大衰减为2dB；阻带截止频率为12kHz，阻带最小衰减为30dB，设计巴特沃斯低通滤波器。 解：过程见板书 总结：低通巴特沃斯滤波器的设计步骤 1）根据技术指标，求出滤波器阶数N 2）根据阶数N查表，得到归一化传输函数Ha(p) 3）将Ha(p)去归一化。将p=s/ Ωc代入Ha(p)，得到实际的滤波器传输函数H(S). MATLAB为低通模拟巴特沃斯滤波器的产生提供了函数buttap，其调用格式为： [z,p,k]=buttap(N) [z,p,k]=butter(N,Rp,Wn,’s’) [z,p,k]=butter(N ,Rp,Wn,’ftype’,’s’) 其中N表示滤波器的阶次，Rp为通带最大衰减(dB)，Wn为截止频率，ftype为滤波器的类型，S表示模拟滤波器，Z表示零点，P表示极点，K表示增益。 巴特沃斯滤波器的阶次可以通过buttord函数确定，该函数的调用格式为： [N,Wn]=buttord(Wp,Ws,Rp,Rs,’ s’) N表示巴特沃斯滤波器的最小阶次，Wn为其截止频率。 例 绘制巴特沃斯低通模拟原型滤波器的平方幅频响应曲线，其中滤波器的阶次分布为2，5，10，20.（exam3_1.m） 总 结 巴特沃斯滤波器的频率特性在通带和阻带内部都是随着频率单调变化的，如果在通带的边缘能够满足指标，那么在阻带的内部会超过设计的指标要求，造成滤波器的阶次N比较高。 二、切比雪夫滤波器（chebyshev） 特点：在通带内是 等波纹变化，在阻带内单调衰减，过渡迅速。 切比雪夫I型滤波器的振幅平方函数为： 式中Ωc为有效通带截止频率，ε是与通带波纹有关的参量，ε大,波纹大 ，0 ε1 ；VN（x）为N阶切比雪夫多项式 例 通带截止频率为4KHz，通带最大衰减1dB。阻带频率为8KHz。阻带最小衰减15dB，试设计切比雪夫滤波器。 切比雪夫模拟滤波器matlab实现 切比雪夫I型低通滤波器可通过MATLAB函数cheb1ap实现，其调用格式为： [z,p,k]= cheb1ap(N,Rp) [b,a]=cheb1 (N,Rp,Wn,’s’) [b,a]=cheb1 (N ,Rp,Wn,’ftype’,’s’) 切比雪夫I型滤波器的阶次可以通过cheb1ord函数确定，该函数的调用格式为： [N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,’ s’) N表示巴特沃斯滤波器的最小阶次，Wn为其截止频率。 例： 绘制切比雪夫I型低通模拟滤波器的平方幅频响应曲线，其中滤波器的阶次分布为2，4，6，8.（exam3_2.m） 总结 对于阶次相同（系统函数复杂程度一样）的巴特沃斯滤波器和切比雪夫滤波器相比，切比雪夫的过渡带下降更快，即有更窄的过渡带，性能更好些。 椭圆滤波器（考尔滤波器） 特点：幅值响应在通带和阻带内都是等波纹的，对于给定的阶数和给定的波纹要求，椭圆滤波器能获得较其它滤波器更窄的过渡带宽，就这点而言，椭圆滤波器是最优的。 其振幅平方函数为： 典型的椭园滤波器振幅平方函数如图 3.10 6.3 根据模拟滤波器设计IIR滤波器 利用模拟滤波器设计数字滤波器，就是从已知的模拟滤波器传递函数Ha(s)设计数字滤波器传递函数H(z)，这是一个由S平面到Z平面的映射变换，这种映射变换应遵循两个基本原则： ?? 1) H(z)的频响要能模仿Ha(s)的频响，即S平面的虚轴应映射到Z平面的单位圆ejω上。 ???? 2) Ha(s)的因果稳定性映射到H(z)后保持不变，即S平面从左半平面Re{S}＜0映射到Z平面的单位圆内|z|1。 一. 脉冲响应不变法 脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n) 模仿模拟滤波器的冲击响应ha(t), 使h(n)正好等于ha(t)的采样值，即