传递函数辨识(2):脉冲响应两点法和三点法

发布时间:2023-04-15 20:59:10

传递函数辨识(2:脉冲响应两点法和三点法丁锋;徐玲;刘喜梅【摘要】本工作利用系统的脉冲响应观测数据,提出了辨识一阶系统、二阶系统传递函数参数的两点法、三点法等,以及确定传递函数参数的差分方程法和面积法.提出的方法能够避免直接求解超越方程,且原理简单,实现方便.%Bymeansofthesystemimpulseresponsedata,thispaperpresentstwo-pointmethodsandthree-pointmethodsforidentifyingtheparametersoffirst-ordersystemsandsecond-ordersystems,whicharedescribedbytransferfunctions,andpresentsthedifferenceequationmethodandtheareamethodforidentifyingtransferfunctions.Theproposedalgebraicmethodsofdeterminingtheparametersofthetransferfunctionshavesimplemechanismandeasetounderstand,andavoidsolvingsometranscendentalequations.【期刊名称】《青岛科技大学学报(自然科学版)》【年(,期】2018(039002【总页数】15(P1-15【关键词】传递函数;参数估计;系统辨识;阶跃响应;脉冲响应【作者】丁锋;徐玲;刘喜梅【作者单位】青岛科技大学自动化与电子工程学院,山东青岛266042;江南大学联网工程学院,江苏无锡214122;江南大学物联网工程学院,江苏无锡214122;青岛科技大学自动化与电子工程学院,山东青岛266042
【正文语种】【中图分类】TP273传递函数是一种参数模型。传递函数辨识就是利用系统的观测数据确定传递函数模型的参数,为后续控制器的设计奠定模型基础。一些统计辨识方法可用于传递函数的参数估计[1-5],也可用于信号模型的参数估[6-11]。在最近的连载论文中,讨论了基于阶跃响应的传递函数参数辨识的两点法和三点法[12]。一些传递函数参数辨识方法可参见文献[13-20]线性系统的传递函数辨识方法有阶跃响应法、脉冲响应法、频率响应法。在经典辨识方法中,脉冲函数和阶跃函数是辨识实验中常用的输入信号。脉冲响应是指系统在脉冲输入信号作用下所产生的零状态输出响应。阶跃响应是指系统在阶跃输入信号作用下所产生的零状态输出响应。脉冲响应和阶跃响应容易获得,因而被广泛用于辨识实验。本工作研究脉冲响应下的传递函数辨识方法,即在脉冲响应曲线上选取某些时间点tk的观测数据y(tk,利用观测数据构造关于传递函数参数的方程组,研究代数求解传递函数未知参数的两点法、三点法等,求出待辨识的传递函数的参数,如一阶系统的时间常数T和增益K等,典型二阶系统的阻尼比ξ和无阻尼自然振荡角频wn等。脉冲响应实验进行参数辨识的原理如图1所示。设系统的输入为脉冲函数δ(t,输出为y(t。在系统稳定工况下,假设系统的输入输出基准线均为零,施加系统的输入量δ(t,即在系统输入端施加一个单位脉冲信号,然后用函数记录仪或数据采集系统记录系统的输入和输出。假设辨识实验中数据采集时刻为tk,收集到各个采样时刻的输出响应观测数据为y(tkk=123,…,L。根据这些观测数据
y(tk可以绘制出脉冲响应曲线来辨识传递函数模型的参数。1数据采集系统示意图Fig.1Data-takingsystem数学上,单位脉冲函数定义为δ(tdt=1,或者定义为它代表宽度为零,高度为无穷大,面积为1的一个冲击。脉冲函数的拉普拉斯(Laplace变换(简称拉氏变换δ(te-stdt=1实际上理想的脉冲信号是很难实现的,它只是一种数学的抽象。工程上使用的近似脉冲函数为其中ε很小。1一阶传递函数系统本节讨论一阶传递函数系统参数辨识的任意两点法和倍数两点法。1.1一阶系统任意两点法考虑下列一阶系统其中:K为系统增益,T为时间常数。一阶系统传递函数的辨识就是估计出这两个参数。假设系统输入为单位脉冲信号δ(t,其拉氏变换为R(s=1。因此在单位脉冲信号δ(t作用下,系统输出响应的拉氏变换为Y(s=G(sR(s=G(sY(s的拉氏反变换得到系统的单位脉冲响应:y(t=L-1[Y(s]=L-1[G(s]=
(1由于脉冲响应y(t中含有两个参数KT,期望用代数方法求解这两个参数。在代数求解方法中,两个未知数需要两个独立方程求解。在辨识实验中,收集脉冲响应观测数据(tky(tkk=12,…,L。选择t1t2两个时刻的观测数据(t1y(t1(t2y(t2来确定参数KT,见图22单位脉冲响应Fig.2Unitimpulseresponse一般情况下,稳定系统的脉冲响应过程持续的时间较短,两个观测时间点的选择不宜太接近,也不宜太远,通常在其幅值的25%75%之间。将这两个时刻的观测数据代入式(1可以得到两个方程:(2(3(2两边除以式(3两边得到两边取对数可得到时间常数T:获得参数T后,代入式(2(3中,可得到增益:K=y(t1Tet1/T=K=y(t2Tet2/T=
实际中观测数据总是存在误差(即干扰噪声,所以求解得到的增益和时间常数记作其估计和为获得高精度的参数估计,可以多次选择不同时刻的观测数据,代入脉冲响应表达式中,构造出方程组,求得不同观测数据下的参数估计,然后取平均值作为最终的参数估计。令不同的观测时刻为tk(k=12,…,L,对应的观测数据为y(tk,计算得到的参数估计为(4(5(6然后取代数平均值作为增益K和时间常数T的估计:(71.2一阶系统倍数两点法为简化计算,考虑倍数观测时间点t1t2=2t1的观测数据(t1y(t1(2t1y(2t1,将其代入式(1中,得到方程组:a=e-t1/T,那么上述方程组可以简化为
(8(9(9与式(8两边相除,得(10两边取对数得ln[y(2t1/y(t1]=-t1/T由上式可求得系统时间常数(11由式(8和式(10使用式(11,得到系统增益为了提高估计精度,可以用多个倍数两点的数据(t2y(t2(2t2y(2t2(t3y(t3(2t3y(2t3,…,(tLy(tL(2tLy(2tL,重复上面的过程,得到KT的一系列估计和最后取其平均值作为增益K和时间常数T的估计:(12
(131.3一阶系统脉冲响应图解法一阶系统的脉冲响应为t=0时,有这表明,在t=0时刻,脉冲响应指数曲线与纵轴相交于当t=T时,有因此,在纵轴0.36788倍的位置对应横坐标的时间tx即为常数T的估计,通过作图方法即可确定参数T的估计然后利用t=tx的脉冲响应数据y(tx,可计算出增K的估计1.4一阶系统脉冲响应差分方程法考虑下列一阶系统的传递函数其中:b为根轨迹增益,a为负极点。若系统的输入为单位脉冲函数,则输出为脉冲响应函数y(t。因为单位脉冲函数只作用于时刻t=0,而在其他时刻系统的输入为零,所以系统的输出y(t是从t=0始的脉冲响应函数。假设采样周期为h,采样起始时间为t1,利用脉冲响应y(t一阶系统可以用一阶差分方程y(t1+cy(t1+h=0(14来表示,其中c为待定系数。根据上述差分方程可以得到
(15这里的观测数据可以使用任意两个相邻观测时间的数据。一阶系统的脉冲响应函数为y(t=be-at那么y(t1=be-at1(16y(t1+h=be-a(t1+h(17将式(16(17代入到式(14中,有be-at1+cbe-a(t1+h=be-at1(1+ce-ah=0(18由于be-at1≠0,要想等式(18成立,那么有1+ce-ah=0x=e-ah,有1+cx=0进一步求得(19(19两边求对数,得到a的估计(20代入式(16或式(17得到参数b的估计为
(21(22差分方程法仅需要少量的观测数据就能得到一阶系统的参数估计,由于参与估计运算的数据较少,因此这种方法的精度很难保证,为提高精度,一种简单的方法就是利用多个观测数据,多次计算,然后取平均值。使用多次采集到的观测数据(tky(tkk=12,…,L,利用差分方程方法得到的一阶系统参数ab的估计为(23(24对于时间常数增益形式的一阶惯性系统传递函数同样可以采用差分方程法,在计算过程中,先将传递函数转换为b/(s+a的形式,用上面的求解公式得到ab的估计和然后根据确定系统增益K和时间常数T的估计。2二阶传递函数系统二阶系统传递函数通常分为典型二阶系统、二阶等容惯性系统、二阶等容惯性加零
点系统、二阶惯性无零点系统、二阶惯性加零点系统等几种形式。这里针对二阶系统传递函数,利用脉冲响应观测数据来估计系统的特征参数。2.1二阶等容惯性系统的任意两点法考虑下列二阶等容惯性系统传递函数其中:K为系统的增益,T为系统的时间常数,它们为待辨识传递函数的参数。单位脉冲函数的拉氏变换为R(s=1。在单位脉冲输入信号作用下,系统的输出响应为(25选取任意两个观测时刻t1t2的观测数据(t1y(t1(t2y(t2,代入式(25得到下列方程组:(26(27(26两边除以式(27两边,整理得到上式两边取对数得由此可以得到二阶等容惯性系统时间常数T的估计为将得到时间常数T的估计代入式(26或式(27可以得到增益K的估计
(28(292.2二阶等容惯性系统的倍数两点法二阶等容惯性系统通过任意两点观测数据,可以用代数方法求解得到参数KT的估计。也可借鉴一阶传递函数参数估计的倍数两点法思想,旨在提出一种较为简单的确定二阶等容惯性系统参数估计的代数求解方法。将两倍采样时刻观测数据(t1y(t1(2t1y(2t1代入式(25得到下列方程组:(30(31定义a=e-t1/T,式(30(31可以简化为(32(33(33(32两边相除得到(34两边取对数得到时间常数T的估计
(35由式(32及式(34可得(36将式(35代入式(36可进一步得到增益K的估计为获得更高精度的代数解参数估计,采集多个两倍观测时间点的观测数据(t2y(t2(2t2y(2t2(t3y(t3(2t3y(2t3,…,(tLy(tL(2tLy(2tL,取平均值,得到二阶等容惯性传递函数系统增益K和时间常数T的估计。(37(382.3二阶等容惯性加零点系统的任意两点法考虑下列二阶等容惯性加零点传递函数其中:T1T2为系统的时间常数,K为系统的增益。两个时间常数包括两种可能:T1>T2T2>T1。如果G(s的阶跃响应没有超调,或G(s脉冲响应不变符号,说明T1>T2。如果G(s的阶跃响应有超调,或G(s脉冲响应改变符号,说明T2>T1。下面的讨论假定T1>T2。读者可针对T2>T1加以讨论。由于单位脉冲信号的拉氏变换为1,在单位脉冲输入信号作用下,系统G(s的脉冲响应为
y(t=L-1[Y(s]=L-1[G(s]=(39由于二阶等容惯性加零点系统脉冲响应中含有3个未知参数KT1T2,按照代数求解方法,3个未知数可以通过3个方程建立方程组来求得,一般的方法是选择单位脉冲响应曲线上3个不同的观测点,就可以估计出参数KT1T2。任意选3个观测时刻数据(t1y(t1(t2y(t2(t3y(t3代入二阶等容惯性加零点系统单位脉冲响应(39中,得到下列方程组(40(41(42方程组(41(42是一个关于未知参数的超越方程,求其代数解几乎是不可能的,为简化计算,考虑采用倍数时间点的观测数据建立方程组,用代数方法求解系统的参数估计。选取倍数时间点t1=τ,t2=2τ,t3=3τ的观测数据(τ,y(τ,(2τ,y(2τ,(3τ,y(3τ,代入单位脉冲响应表达式(39中,可得到下列方程组
(43(44(45为计算简便,令于是方程组(43(45可以简化为y(τ=Kα(τβ+γ,(46y(2τ=Kα2(2τβ+γ,(47y(3τ=Kα3(3τβ+γ。(48由式(46可得(49将式(49代入式(47可得(50将式(49代入式(48可得(51由式(50可得
Kα2τβ=y(2τ-αy(τ。(52将式(52代入式(51可得y(3τ=2α[y(2τ-αy(τ]+α2y(τ=2αy(2τ-2α2y(τ+α2y(τ=2αy(2τ-α2y(τ。将上式进一步整理可以得到关于α的一元二次方程y(τα2-2y(2τα+y(3τ=0。根据一元二次方程求根公式得到α的解为根据定义α=e-τ/T1>0,所以有等式两边取对数可以求出T1=(53在计算出α的基础上(即已经得到T1e-τ/T1,为求解另外两个参数KT2将式(44与式(43两边相除得到由此可得y(2τ[T1T2+(T1-T2τ]=y(τ[T1T2+(T1-T22τ]e-τ/T1。
整理得y(2ττT1+(T1-τy(2τT2=2y(ττT1e-τ/T1+(T1-2τy(τe-τ/T1T2。由此可以求出(54在计算出T1T2的估计之后,代入式(43可以计算出增益K(55(53(55给出利用单位脉冲响应求解二阶等容惯性加零点系统传递函数参数估计的算法如下:(56(57(58单位脉冲响应二阶等容惯性加零点系统传递函数参数估计代数求解方法的步骤如下。1收集单位脉冲响应观测数据,绘制单位脉冲响应曲线。2确定一个合适的采样时刻τ,采集观测数据(τ,y(τ,(2τ,y(2τ,(3τ,y(3τ,用式(56计算时间常数的估计3用式(57计算时间常数的估计
4用式(58计算增益的估计1倍数三点法确定传递函数的参数需要根据三个倍数时刻的观测数据先计算出时间常数T1的估计然后根据计算出T2的估计最后在此基础上计算出增益K的参数估计在计算出后,可以根据3个观测数据点构成3个方程组中的任意两个方程来计算和为了提高估计精度,可以取多个倍数三点的观测数据计算得到T1T2K的一系列估计,最后取其平均值作为系统参数估计。具体做法如下。取多个不同的采样时τ=τk,k=12,…,L,利用观测数据(τk,y(τk,(2τk,y(2τk和(3τk,y(3τk,得到计算T1T2K的估计算法如下:(59(60(612.4二阶惯性无零点系统的倍数三点法考虑下列二阶惯性无零点系统的传递函数其中K为系统的增益,T1T2为系统的时间常数,它们都是需要辨识的参数。不妨设T1系统在单位脉冲信号作用下的输出响应为
(62将任意3个观测时刻的观测数据(t1y(t1(t2y(t2(t3y(t3代入单位脉冲响应式(62中,得到关于未知参数KT1T2的超越方程,故无法求得参数KT1T2的代数解,原因是存在指数项。为此,同样考虑采用倍数时间点,利用指数函数的特点,简化方程组以求得二阶惯性无零点系统传递函数的参数。选取倍数时间(t1=τ,t2=2τ,t3=3τ观测数据(τ,y(τ,(2τ,y(2τ,(3τ,y(3τ,代入单位脉冲响应表达式(62中,可得到下列方程组(63(64(65为简化方程组求解,作如下定义α=e-τ/T1,β=e-τ/T2。于是方程组(63(65可以简化为(66(67(68
将式(67与式(66两边相除,可得(69将式(68与式(66两边相除,可得利用式(69可得进一步可得(70此时αβ为两个未知数,且它们之和(69与之积(70是已知的,根据方程根与系数的关系可以得到下列方程:x2-(α+βx+αβ=0,(71该方程的两个根即为αβ。将式(69(70代入式(71,根据二次方程的求根公式得αβ的代数解为根据定义α=e-τ/T1,β=e-τ/T2以及αβ的代数解,可以得到时间常数T1T2的估计为(72
(73在得到时间常数估计和后,可将其代入式(63(65中的任意一个方程,即可求得增益K的估计(也可以都代入后取其平均值。将和代入式(63得到增益K的估计(742对于二阶惯性加零点系统T1,T3≠T1,T3≠T2,基于脉冲响应数据,有待读者找到代数求解系统增益K和时间常数T1T2T3的方法。2.5二阶惯性加零点系统的差分方程法考虑下列二阶惯性加零点系统传递函数T1,T3≠T1,T3≠T2,(75其中:K为系统的增益,T1T2T3为系统的时间常数。对于二阶惯性加零点系统(75,基于脉冲响应数据,采用多点法,还难以找到代数方法求解系统增益K和时间常数T1T2T3。下面使用差分方程法间接估计这些参数。将上述传递函数G(s转化为部分分式形式,可得定义
(76那么传递函数G(s可以进一步表示为(77其中:abcd为未知参数。考虑到直接利用单位脉冲响应求得参数KT1T2T3比较困难,这里将传递函数进行部分分式展开,通过间接的方法来估计这些未知参数。给系统施加单位脉冲信号,收集单位脉冲响应观测数据。设采样周期为h,那么该二阶系统可以用下列二阶差分方程来描述:y(t1+c1y(t1+h+c2y(t1+2h=0y(t1=-c1y(t1+h-c2y(t1+2h(78其中:t1为起始观测时刻,c1c2为待定系数。将采集到的4个时刻的观测数据y(t1+khk=0123,代入到二阶差分方(78中,可以得到下列关系:y(t1=-c1y(t1+h-c2y(t1+2h(79y(t1+h=-c1y(t1+2h-c2y(t1+3h(80为求解方程组(79(80,将其写成如下矩阵方程
(81y(t1+hy(t1+3h-y2(t1+2h≠0时,解矩阵方程(81可得于是得到二阶差分方程中的待定系数c1c2的解:(82(83二阶系统传递函数在单位脉冲作用下的响应可以通过传递函数的拉氏反变换求得,be-at+de-ct(84将观测数据(t1+khy(t1+khk=0123,代入单位脉冲响应(84中,可y(t1=be-at1+de-ct1(85y(t1+h=be-a(t1+h+de-c(t1+h(86y(t1+2h=be-a(t1+2h+de-c(t1+2h
(87将式(85(87代入二阶差分方程(78be-at1+de-ct1+c1(be-a(t1+h+de-c(t1+h+c2(be-a(t1+2h+de-c(t1+2h=0将上述方程整理得be-at1[1+c1e-ah+c2(e-ah2]+de-ct1[1+c1e-ch+c2(e-ch2]=0由于be-at1≠0,de-ct≠0,上式方程对任意t1都成立,因此下列条件必须满足:1+c1e-ah+c2(e-ah2=01+c1e-ch+c2(e-ch2=0由于c1c2已经求出,定义x1=e-ahx2=e-ch,那么上述两个方程可以转化为一个一元二次方程c2x2+c1x+1=0该方程的两个根分别为e-ahe-ch,利用二次方程求根公式得根据定义x1=e-ahx2=e-ch取对数可以得到参数ac的估计:(88(89在求得参数ac后,选取两个时刻的观测数据y(t1y(t2代入脉冲响应(84中,
得到下列方程组:y(t1=be-at1+de-ct1y(t2=be-at2+de-ct2写成矩阵方程为求解该矩阵方程可得到参数bd的解:使用参数ac的估计,可求出参数bd的估计:(90(91在得到参数估计之后,可进一步计算参数KT1T2T3的估计。由于直接根abcd的定义式(76求解KT1T2T3比较麻烦,为此采用比较系数的方法,将传递函数(77通分得(92比较式(92和式(75可得(93联立式(82(83(88(91,能够得到二阶惯性加零点系统传递函数参数估计算法:
(94(95(96(97(98(99(100二阶惯性加零点系统传递函数参数估计算法计算参数估计和的步骤如下。1施加单位脉冲信号,设置采样周期h,收集单位脉冲响应观测数据。2选取4个时刻的观测数据y(t1y(t1+hy(t1+2hy(t1+3h3用式(94计算参数估计c1,用式(95计算参数估计c24用式(96计算参数估计用式(97计算参数估计5任意选取两个时刻的观测数据y(t1y(t2,用式(98计算参数估计用式(99算参数估计6用式(100计算参数估计和3计算参数估计和的两个时刻的观测数据可以选择之前计算参数估计和时的3个观测数据中的2个,也可以使用其它时刻的观测数据。这个算法仅仅利用了少量的观测数据,因此不能保证参数估计精度,为此可以利用更多时刻的观测数据,
进行多次计算,然后取平均值。4这里的估计方法是针对二阶惯性加零点系统传递函数提出的,从上述推导可以看出,该方法适用于所有可以通过部分分式分解得到形如式(77的二阶系统传递函数。3典型欠阻尼二阶系统本节针对典型欠阻尼二阶系统传递函数,基于脉冲响应观测数据,研究倍数三点法、特殊两点图解法、面积法辨识传递函数的参数。3.1典型二阶系统倍数三点法实际中有些系统不是由惯性环节构成的,也就是系统包含了复数极点,如典型二阶系统,其传递函数为(101其中:0<ξ<1为阻尼比,wn为无阻尼自然角频率,ξwn为待辨识的参数。在脉冲输入信号作用下,典型二阶系统的输出响应为(102定义(103那么式(102可以简化为y(t=asin(wte-ξwnt。(104
选择脉冲响应倍数观测时刻t1t2=2t1t3=3t1的观测数据(t1y(t1(2t1y(t2(3t1y(t3,定义β=e-ξwnt1,(105将观测数据代入式(104,可以得到下列方程组:y(t1=aβsin(wt1,(106y(t2=aβ2sin(2wt1,(107y(t3=aβ3sin(3wt1。(108将式(107按三角函数展开,可得y(t2=2aβ2sin(wt1cos(wt1=2βy(t1cos(wt1。由上式可进一步得(109由式(106可得(110根据三角函数的性质sin2(wt1+cos2(wt1=1,将式(109(110代入得到由此可得
(111将式(108按三角函数展开,可得y(t3=aβ3[3sin(wt1-4sin3(wt1]=3aβ3sin(wt1-4aβ3sin3(wt1=3β2y(t1-(112整理得到a2y(t3=3a2β2y(t1-4y3(t1使用式(111得到上式两边同乘以4y2(t1得到4a2y2(t1y(t3=3a2y(t1y2(t2+12y5(t1-16y5(t1整理得(113两边开方得到(114由式(111可得使用式(113可得开方得到
(115将式(115代入式(109,可得两边取反余弦,可得l=012,…。对于适当的t1,取l=0,由上式求得(116根据定义式(103(105,以及由式(114(116计算出的参数a,βw,得到下列方程组(117(118e-ξwnt1=β。(119求解上述方程组就可以得到参数ξwn的估计,(117和式(118两边相乘得到两边开方得使用式(114(116得到wn的估计
(120(119两边取对数得-ξwnt1=lnβ。进一步得使用式(115和式(120得到ξ的估计(121(120(121给出了确定典型二阶系统传递函数参数ξwn的倍数三点求解算法。为了获得更高精度的参数估计,可以多次选取3倍数时间观测数据(t2y(t2(2t2y(2t2(3t2y(3t2(t3y(t3(2t3y(2t3(3t3y(3t3,…,(tLy(tL(2tLy(2tL(3tLy(3tL,利用式(120(121多次计算,得到ξwn的一系列估计和最后取其平均值:最后取其平均值:3.2典型二阶系统特殊两点法考虑下列欠阻尼典型二阶系统其中:ξ为阻尼比,且0<ξ<1,wn为无阻尼自然角频率。欠阻尼典型二阶系统的单位脉冲响应为单位脉冲响应曲线如图3所示。
3单位脉冲响应曲线特殊采样数据点Fig.3Particularpointsoftheimpluseresponsecurvey(t关于时间t的一阶导数得(122收集各个时刻的观测数据,可以绘制出单位脉冲响应曲线。从图3响应曲线上读(t1y(t1,此刻对应的脉冲响应一阶导数为零,即由于所以有可得(123由此可得脉冲响应最大峰值时间对应的脉冲响应最大峰值为
(124从脉冲响应曲线上读出(t2y(t2,此时y(t2=0。将(t20代入脉冲响应中得到因为所以上式成立的条件为于是,有l=1,由上式可得(125将式(125代入式(123(126求解得到ξ的估计(127求出ξ后,由式(125可求得参数wn的估计5这里计算欠阻尼典型二阶系统传递函数阻尼比和无阻尼自然角频率的参数估计方法是在脉冲响应曲线上找到一个极值点时刻t1和一个响应为零的时刻t2(3所示来估计参数的。由于仅仅依靠两个观测时间来计算参数估计,因此这种方法的参数估计精度不能保证。但这种方法计算简单,可以用于精度要求不高的场合,也可采用多次平均值代替。
3.3典型二阶系统面积法考虑下列典型二阶系统传递函数(128其中:ξ为阻尼比,且0<ξ<1,wn为无阻尼自然角频率。收集单位脉冲响应观测数据,绘制出脉冲响应曲线y(t如图4所示,从图中可以读t1t2。脉冲响应曲线开始的两个波形与横轴的面积S1S2可以采用梯形方法近似计算。假设采样周期为h,采样时刻为tk=khk=012,…,收集单位脉冲响应观测数据y(tk。当h较小时,面积S1S2可近似计算出:4单位脉冲响应曲线Fig.4Unitimpluseresponsecurve(129(130这里[x]表示不大于x的最大整数。如x=1963[x]=1963;x=3.09[x]=3在利用梯形方法近似计算出面积S1S2后,再推导出面积与系统阻尼比ξ和自然角频率wn的关系,并进行求解。典型二阶系统(128的单位脉冲响应的解析表达式为先求y(t与时间轴交点的坐标tk。令y(t=0,得到所以脉冲响应曲线开始的两个波形与横轴的面积S1S2可以通过下列积分计算:
利用积分公式(131依次可得π2ξ2=ln2(1-S1(1-ξ2,[π2+ln2(1-S1]ξ2=ln2(1-S1由此可求出阻尼比的估计阻尼比也可按下列方式求。计算面积
(132取对数得到由此可求出阻尼比的估计自然角频率wn的估计为此外,设脉冲响应曲线第k+1个波形与横轴的面积为(133使用第1个面积S1与第k+1个面积Sk+1计算出的阻尼比估计为k=123,…,对应的自然角频率wn的估计为也可以取平均值作为ξw的估计
(134(135(136随着k的增加,面积Sk逐渐趋于零,因此L的取值不能太大。4基于单位脉冲响应,针对一阶系统传递函数、二阶惯性系统传递函数、二阶等容惯性系统传递函数、典型二阶系统传递函数等,讨论了两点法、倍数两点法、倍数三点法、面积法估计传递函数的参数。[1]丁锋.系统辨识新论[M].北京:科学出版社,2013.DINGFeng.SystemIdentification-NewTheoryandMethods[M].Beijing:SciencePress,2013.[2]丁锋.系统辨识—辨识方法性能分析[M].北京:科学出版社,2014.DINGFeng.SystemIdentification-PerformanceAnalysisforIdentificationMethods[M].Beijing:SciencePress,2014.[3]丁锋.系统辨识—辅助模型辨识思想与方法[M].北京:科学出版社,2017.DINGFeng.SystemIdentification-AuxiliaryModelIdentificationIdeaandMethods[M].Beijing:SciencePress,2017.[4]丁锋.系统辨识—迭代搜索原理与辨识方法[M].北京:科学出版社,2018.DINGFeng.SystemIdentification-IterativeSearchPrincipleandIdentificationMethods[M].Beijing:SciencePress,2018.
[5]丁锋.系统辨识—多新息辨识理论与方法[M].北京:科学出版社,2016.DINGFeng.SystemIdentification-Multi-InnovationIdentificationTheoryandMethods[M].Beijing:SciencePress,2016.[6]丁锋,徐玲,刘喜梅.信号建模(1:单频率信号[J].青岛科技大学学报(自然科学,2017,38(1:1-13.DINGFeng,XULing,LIUXimei.Signalmodeling-PartA:Single-frequencysignals[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,2017,38(1:1-13.[7]丁锋,徐玲,刘喜梅.信号建模(2:双频率信号[J].青岛科技大学学报(自然科学,2017,38(2:1-17.DINGFeng,XULing,LIUXimei.Signalmodeling-PartB:Dual-frequencysignals[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,2017,38(2:1-17.[8]丁锋,徐玲,刘喜梅.信号建模(3:多频信号模型的递推参数估计[J].青岛科技大学学报(自然科学版,2017,38(3:1-12.DINGFeng,XULing,LIUXimei.Signalmodeling-PartC:Recursiveparameterestimationformulti-frequencysignalmodels[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,2017,38(3:1-12.[9]丁锋,徐玲,刘喜梅.信号建模(4:多频信号模型的迭代参数估计[J].青岛科技大学学报(自然科学版,2017,38(4:1-11.DINGFeng,XULing,LIUXimei.Signalmodeling-PartD:Iterativeparameterestimationformulti-frequencysignalmodels[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,
2017,38(4:1-11.[10]丁锋,徐玲,刘喜梅.信号建模(5:多频信号模型的递阶参数估计[J].青岛科技大学学报(自然科学版,2017,38(5:1-15.DINGFeng,XULing,LIUXimei.Signalmodeling-PartE:Hierarchicalparameterestimationformulti-frequencysignalmodels[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,2017,38(5:1-15.[11]丁锋,徐玲,刘喜梅.信号建模(6:多频信号模型的递阶迭代参数估计[J].青岛科技大学学报(自然科学版,2017,38(6:1-13.DINGFeng,XULing,LIUXimei.Signalmodeling-PartF:Hierarchicaliterativeparameterestimationformulti-frequencysignalmodels[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,2017,38(6:1-13.[12]丁锋,徐玲,刘喜梅.传递函数辨识(1:阶跃响应两点法和三点法[J].青岛科技大学学报(自然科学版,2018,39(1:1-14.DINGFeng,XULing,LIUXimei.Transferfunctionidentification.PartA:Two-pointandthree-pointmethodsbasedonthestepresponses[J].JournalofQingdaoUniversityofScienceandTechnology(NaturalScienceEdition,2018,39(1:1-14.[13]丁锋,谢新民,周荣宝.经典辨识方法及其在电厂温度系统中的应用[C]//东北大学出版社,中国控制与决策学术年会论文集,厦门,1994:426-430.DINGFeng,XIEXinmin,ZHOURongbao.Classicidentificationmethodanditsapplicationinthetemperaturecontrolsystemsinpowersupply[C].NortheastUniversityPress,ProceedingofChineseControlandDecision
Conference,Xiamen,1994:426-430.[14]陈垒.连续系统的传递函数模型辨识[D].无锡:江南大学,2012.CHENLei.IdentificationforTransferFunctionModelsofContinuous-TimeSystems[D].Wuxi:JiangnanUniversity,2012.[15]CHENL,LIJH,DINGR.Identificationofthesecond-ordersystemsbasedonthestepresponse[J].MathematicalandComputerModelling,2011,53(5/6:1074-1083.[16]XUL,CHENL,XIONGWL.ParameterestimationandcontrollerdesignfordynamicsystemsfromthestepresponsesbasedontheNewtoniteration[J].NonlinearDynamics,2015,79(3:2155-2163.[17]XUL.ApplicationoftheNewtoniterationalgorithmtotheparameterestimationfordynamicalsystems[J].JournalofComputationalandAppliedMathematics,2015,288:33-43.[18]XUL,DINGF.Theparameterestimationalgorithmsfordynamicalresponsesignalsbasedonthemulti-innovationtheoryandthehierarchicalprinciple[J].IETSignalProcessing,2017,11(2:228-237.[19]XUL,DINGF.Parameterestimationforcontrolsystemsbasedonimpulseresponses[J].InternationalJournalofControl,AutomationandSystems,2017,15(6:2471-2479.[20]徐玲.基于移动数据窗的传递函数多新息随机梯度辨识方法[J].控制与决策,2017,32(6:1091-1096.XULing.Movingdatawindowbasedmulti-innovationstochasticgradientidentificationmethodfortransferfunctions[J].ControlandDecision,2017,32(6:1091-1096.

传递函数辨识(2):脉冲响应两点法和三点法

相关推荐