《无限长单位脉冲响应(IIR)滤波器设计(2).ppt》由会员分享,可在线阅读,更多相关《无限长单位脉冲响应(IIR)滤波器设计(2).ppt(36页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、无限长单位脉冲响应无限长单位脉冲响应(IIR)滤波器设计滤波器设计(2)s平面s1平面z平面双线性变换法的映射关系为了将S平面的j轴压缩到S1平面j1轴上的/T到/T一段上,可通过以下的正切变换实现:这里C是待定常数,下面会讲到用不同的方法确定C,可使模拟滤坡器的频率特性与数字源波器的频率特性在不同频率点有对应关系。经过这样的频率变换,当1由 时,由-/T经过变化到/T,即S平面的整个j轴被压缩到S1平面的。将这一关系解析扩展至整个S平面,则得到S平面到S1平面的映射关系:再将 S1 平面通过标准变换关系映射到Z平面,即令通常取C=2/T,考虑z=ej,最后得S平面与Z平面的单值映射关系:双线
2、性换法的主要优点是S平面与Z平面一一单值对应,S平面的虚轴(整个j)对应于Z平面单位圆的一周,S平面的=0处对应于Z平面的=0处,对应即数字滤波器的频率响应终止于折迭频率处,所以双线性变换不存在混迭效应。现在我们看看,这一变换是否符合我们一开始提出的由模拟滤波器设计数字滤波器时,从S平面到Z平面映射变换的二个基本要求:当时,代入(1)式,得:对单位圆,即S平面的虚轴映射到Z平面正好是单位圆。代入(2)式当时,即s左半平面映射在单位圆内,s右半平面映射在单位圆外,因此稳定的模拟滤波器通过双线性变换后,所得到的数字滤波器也是稳定的。如图1。图 双线性变换的频率非线性关系小结 1)与脉冲响应不变法相
3、比,双线性变换的主要优点:S平面与Z平面是单值的一一对应关系(靠频率的严重非线性关系得到的),即整个j轴单值的对应于单位圆一周,关系式为:可见,和为非线性关系,如图2。图2 双线性变换的频率非线性关系由图中看到,在零频率附近,接近于线性关系,进一步增加时,增长变得缓慢,(终止于折叠频率处),所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。2)双线性变换缺点:与成非线性关系,导致:a.数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变,(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。例如,一个模拟微分器,它的幅度与频率是直线关系,但通过双线性变换后,就不可
4、能得到数字微分器 b.线性相位模拟滤波器经双线性变换后,得到的数字滤波器为非线性相位。c.要求模拟滤波器的幅频响应必须是分段恒定的,故双线性变换只能用于设计低通、高通、带通、带阻等选频滤波器。虽然双线性变换有这样的缺点,但它目前仍是使用得最普遍、最有成效的一种设计工具。这是因为大多数滤波器都具有分段常数的频响特性,如低通、高通、带通和带阻等,它们在通带内要求逼近一个衰减为零的常数特性,在阻带部分要求逼近一个衰减为的常数特性,这种特性的滤波器通过双线性变换后,虽然频率发生了非线性变化,但其幅频特性仍保持分段常数的特性。预畸变:即将模拟滤波器的临界频率事先加以畸变,然后通过双线性变换后正好映射到所
5、需要的频率上。利用关系式:将所要设计的数字滤波器临界频率点,变换成对应的模拟域频率,利用此设计模拟滤波器,再通过双线性变换,即可得到所需的数字滤波器,其临界频率正是。如图所示。双线性变换时频率的预畸)计算H(Z)双线性变换比脉冲响应法的设计计算更直接和简单。由于s与z之间的简单代数关系,所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。置换过程:频响:这些都比脉冲响应不变法的部分分式分解便捷得多,一般,当着眼于滤波器的时域瞬态响应时,采用脉冲响应不变法较好,而其他情况下,对于IIR的设计,大多采用双线性变换。3.2 3.2 常用模拟低通滤波器特性常用模拟低通滤波器特性 为了方便学习
6、数字滤波器,先讨论几种常用的模拟低通滤波器设计方法,高通、带通、带阻等模拟滤波器可利用变量变换方法,由低通滤波器变换得到。模拟滤波器的设计就是根据一组设计规范设计模拟系统函数Ha(s),使其逼近某个理想滤波器特性。因果系统中式中ha(t)为系统的冲激响应,是实函数。不难看出定义振幅平方函数式中Ha(s)模拟滤波器系统函数Ha(j)滤波器的频率响应|Ha(j)|滤波器的幅频响应又S=j,2=-S2A(2)=A(-S2)|S=j问题:由A(-S2)Ha(S)对于给定的A(-S2),先在S复平面上标出A(-S2)的极点和零点,由(1)式知,A(-S2)的极点和零点总是“成对出现”,且对称于S平面的实
7、轴和虚轴,选用A(-S2)的对称极、零点的任一半作为Ha(s)的极、零点,则可得到Ha(s)。为了保证Ha(s)的稳定性,应选用A(-S2)在S左半平面的极点作为Ha(s)的极点,零点可选用任一半。三种模拟低通滤波器的设计:1)巴特沃兹滤波器(Butterworth滤波器)(巴特沃兹逼近)特点:具有通带内最大平坦的振幅特性,且随f,幅频特性单调。其幅度平方函数:N为滤波器阶数,如图1图1巴特沃兹滤波器振幅平方函数通带:使信号通过的频带阻带:抑制噪声通过的频带过渡带:通带到阻带间过渡的频率范围c:截止频率。过渡带为零,阻带|H(j)|=0通带内幅度|H(j)|=1,H(j)的相位是线性的。理想滤
8、波器图1中,N增加,通带和阻带的近似性越好,过渡带越陡。通带内,分母/c1,(/c)2N 1,增加,A(2)快速减小。=c,,,幅度衰减 ,相当于3db衰减点。振幅平方函数的极点:令分母为零,得 可见,Butter worth滤波器 的振幅平方函数有2N个极 点,它们均匀对称地分布在|S|=c的圆周上。例:为N=3阶BF振幅平方函数的极点分布,如图。图2三阶A(-S2)的极点分布考虑到系统的稳定性,知DF的系统函数是由S平面左半部分的极点(SP3,SP4,SP5)组成的,它们分别为:系统函数为:令 ,得归一化的三阶BF:如果要还原的话,则有2)切比雪夫(chebyshev)滤波器(切比雪夫多项
9、式逼近)特点:误差值在规定的频段上等幅变化。巴特沃兹滤波器在通带内幅度特性是单调下降的,如果阶次一定,则在靠近截止频率处,幅度下降很多,或者说,为了使通常内的衰减足够小,需要的阶次(N)很高,为了克服这一缺点,采用切比雪夫多项式逼近所希望的 。切比雪夫滤波器的 在通带范围内是等幅起伏的,所以同样的通带衰减,其阶数较巴特沃兹滤波器要小。可根据需要对通带内允许的衰减量(波动范围)提出要求,如要求波动范围小于1db。振幅平方函数为 有效通带截止频率 与通带波纹有关的参量,大,波纹大。0 1 VN(x)N阶切比雪夫多项式,定义为如图1,通带内变化范围1c,随/c,0(迅速趋于零)当 =0时,N为偶数,
10、min,N为奇数,max,切比雪夫滤波器的振幅平方特性有关参数的确定:a、通带截止频率,预先给定b、由通带波纹表为 给定通带波纹值分贝数 后,可求 。c、阶数N由阻带的边界条件确定。(、A事先给定)3、椭圆滤波器(考尔滤波器)特点:幅值响应在通带和阻带内都是等波纹的,对于给定的阶数和给定的波纹要求,椭圆滤波器能获得较其它滤波器更窄的过渡带宽,就这点而言,椭圆滤波器是最优的。其振幅平方函数为RN(,L)雅可比椭圆函数L表示波纹性质的参量 N=5,的特性曲线可见,在归一化通带内(-11),在(0,1)间振荡,而超过L后,在间振荡。这一特点使滤波器同时在通常和阻带具有任意衰减量。下图为典型的椭园滤波器振幅平方函数 椭圆滤波器的振幅平方函数图中和A的定义同切比雪夫滤波器rr当c、s、和A确定后,阶次N的确定方法为:式中为第一类完全椭圆积分查阅有关模拟滤波器手册,确定。上面讨论了三种最常用的模拟低通滤波器的特性和设计方法,设计时按照指标要求,合理选用。一般,相同指标下,椭圆滤波器阶次最低,切比雪夫次之,巴特沃兹最高。以 上 讨 论 了 由 A(2)Ha(s),下 面 讨 论 由Ha(s)H(Z)的变换设计法。