根据其它数学模型建立状态空间模型

单击此处编辑母版标题样式,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,Ch.2,控制系统的状态空间模型,2.3,根据其它数学模型建立状态空间模型,本节讨论由描述线性定常系统的其它数学模型,通过选择适当的状态变量建立系统的状态空间模型,.,由系统的输入输出关系模型求其状态空间模型的问题称为系统的,实现,问题,本节的内容为：,由高阶常微分方程建立状态空间模型,由传递函数建立状态空间模型,由系统方框图建立状态空间模型,2.3.1,由高阶常微分方程建立状态空间模型,本节主要讨论由描述系统输入输出关系的常微分方程建立系统的状态空间模型,分别讨论由,不含输入量导数项和,含输入量导数项的,微分方程建立状态空间模型,.,本节关键问题:,如何选择状态变量,保持系统的输入输出间的动态和静态关系不变,微分方程中不包含输入量的导数项(1,/9),1.,微分方程中不包含输入量的导数项,描述单输入单输出线性系统的输入输出间动态行为,不包含有输入量的导数项时的线性定系数常微分方程为,y,(,n,),a,1,y,(,n,-1),a,n,y,bu,(2.1),其中,y,和,u,分别为系统的输出和输入,n,为系统的阶次,.,这里所要研究的是建立上述常微分方程描述的动态系统的如下状态空间模型,本节问题的关键是如何选择状态变量,微分方程中不包含输入量的导数项(2,/9),由微分方程理论知, 若初始时刻,t,0,的初值,y,(,t,0,),y,(,t,0,),y,(,n,1),(,t,0,),已知, 则对给定的输入,u,(,t,),微分方程,(2.1),有唯一解,也即系统在,t,t,0,的任何瞬时的动态都被唯一确定,.,因此,选择状态变量如下,x,1,(,t,),y,(,t,),x,2,(,t,),y,(,t,), ,x,n,(,t,),y,(,n,-1),(,t,),可完全刻划系统的动态特性,取输出,y,及其各阶导数为状态变量,物理意义明确,易于接受,微分方程中不包含输入量的导数项(3,/9),将上述选择的状态变量代入输入输出的常微分方程,有如下状态方程,和输出方程,y,x,1,微分方程中不包含输入量的导数项(4,/9),将上述状态方程和输出方程写成矩阵形式有其中,x, ,x,1,x,2,x,n,T,u, ,u,y, ,y,.,微分方程,:,y,(,n,),a,1,y,(,n,-1),a,n,y,bu,状态变量,:,x,1,y,x,2,y,(1), ,x,n,y,(,n,-1),微分方程中不包含输入量的导数项(,5/9),该状态空间模型可简记为,:,其中,微分方程中不包含输入量的导数项(,6/9),上述式子清楚说明了状态空间模型中系统矩阵,A,与微分方程,(2.1),中的系数,a,1,a,2,a,n,之间,输入矩阵,B,与方程,(2.1),中系数,b,之间的对应关系,.,通常将上述取输出,y,及其各阶导数为状态变量称为相变量,.,上述状态空间模型中的系统矩阵具有特别形式,该矩阵的最后一行与其矩阵特征多项式的系数有对应关系,前,n,-1,行为,1,个,n,-1,维的零向量与,(,n,-1),(,n,-1),的单位矩阵,.,微分方程中不包含输入量的导数项(,7/9),上述实现状态空间模型的模拟结构图如下图所示,y,x,1,微分方程中不包含输入量的导数项(,8/9),例,将以下系统输入输出方程变换为状态空间模型,y,”,6,y,”,11,y,6,y,2,u,解,本例中,a,1,6,a,2,11,a,3,6,b,2,因此,当选择输出,y,及其,1,阶与,2,阶导数等相变量为状态变量时,可得状态空间模型如下,微分方程中不包含输入量的导数项(,9/9),其系统结构图如下所示,微分方程中包含输入量的导数项(1,/10),2.,微分方程中包含输入量的导数项,描述单输入单输出线性系统的输入输出间动态行为的微分方程的一般表达式为,y,(,n,),+,a,1,y,(,n,-1),+,a,n,y,=,b,0,u,(,n,),+,b,n,u,本小节所要研究的是建立上述常微分方程描述的动态系统的如下状态空间数学模型-状态空间模型,建立该状态空间模型的关键是如何选择状态变量,微分方程中包含输入量的导数项(2,/10),若按照前面的方法那样选取相变量为状态变量,即,x,1,(,t,),y,(,t,),x,2,(,t,),y,(,t,), ,x,n,(,t,),y,(,n,-1),(,t,),则可得如下状态方程,上述状态方程中输入,u,的各阶导数可能不连续,从而使微分方程解的存在性和唯一性的条件不成立,.,因此,状态方程中不应有输入,u,的导数项出现,即不能直接将输出,y,的各阶导数项取作状态变量,.,微分方程中包含输入量的导数项(3,/10),为避免状态方程中显式地出现输入的导数,通常,可利用输出,y,和输入,u,以及其各阶导数的线性组合来组成状态变量,其原则是:,使状态方程中不显含输入,u,的各阶导数,基于这种思路选择状态变量的方法很多,下面只介绍一种,微分方程中包含输入量的导数项(,4/10),根据上述原则,选择状态变量如下,其中,i,(,i,=0,1,n,),为待定系数。对各式两边求导数得到,:,y,(,n,),+,a,1,y,(,n,-1),+,a,n,y,=,b,0,u,(,n,),+,b,n,u,微分方程中包含输入量的导数项(,5/10),因此,有,微分方程中包含输入量的导数项(,6/10),因此,有,0,0,0,0,n,微分方程中包含输入量的导数项(,7/10),若待定系数,i,(,i,0,1,n,),满足如下关系式,0,b,0,1,b,1,a,1,0,2,b,2,a,1,1,a,2,0,n,1,b,n,1,a,1,n,2,a,n,1,0,则有,n,微分方程中包含输入量的导数项(,8/10),i,(,i,0,1,n,),满足,:,0,b,0,1,b,1,a,1,0,2,b,2,a,1,1,a,2,0,n,1,b,n,1,a,1,n,2,a,n,1,0,n,b,n,a,1,n,1,a,n,1,1,a,n,0,x,i,满足,:,微分方程中包含输入量的导数项(,9/10),待定系数,:,状态方程,:,输出方程,:,微分方程,:,y,(,n,),+,a,1,y,(,n,-1),+,a,n,y,=,b,0,u,(,n,),+,b,n,u,微分方程中包含输入量的导数项(,10/10),上述实现状态空间模型的模拟结构图如右图所示,由传递函数建立状态空间模型,(1/,5),2.3.2,由传递函数建立状态空间模型,下面讨论由描述系统输入输出关系的传递函数建立系统的状态空间模型,关键问题:,1. 如何选择状态变量,2. 保持系统的输入输出间的动态和静态关系不变,由于传递函数与线性常系数常微分方程有直接的对应关系,故前面讨论的由高阶线性微分方程建立状态空间模型的方法同样适用于将传递函数模型变换为状态空间模型,.,类似地,本节讨论的由传递函数建立状态空间模型的方法亦适用于对微分方程建立状态空间模型,.,由传递函数建立状态空间模型,(,2/5),实际物理系统传递函数中分子多项式阶次小于或等于其分母多项式阶次,此时称该传递函数为真有理传递函数,.,而分子多项式阶次小于分母多项式阶次时,则称为严格真有理传递函数,.,本节讨论描述,单输入单输出(,SISO),线性系统的输入输出间动态行为的如下传递函数,由传递函数建立状态空间模型,(,3/5),对上述传递函数,由长除法,有,其中,由传递函数建立状态空间模型,(,4/5),本节所要研究的是建立该传递函数所描述的动态系统的状态空间模型,(,A,B,C,D,),上述常数项,d,即为状态空间模型,(,A,B,C,D,),中的直联矩阵,D;,严格真有理传递函数,G,(,s,),对应可建立(,A,B,C,D,),中的(,A,B,C,),即,由传递函数建立状态空间模型,(,5/5),下面分传递函数,极点互异和,有重极点,两种情况讨论如何建立状态空间模型,传递函数中极点互异时的变换(1/8),1.,传递函数中极点互异时的变换,对于传递函数,G,(,s,),其特征方程为,s,n,+,a,1,s,n,-1,+,a,n,= 0,若其特征方程的,n,个特征根,s,1,s,2,s,n,互异,则用部分分式法可将,G,(,s,)表示为,如下并联分解,其中,k,1,k,2,k,n,为待定系数,其计算公式为,传递函数中极点互异时的变换(2/8),下面以,k,1,计算式的推导过程为例说明的,k,i,的计算式,将,G,(,s,),乘以,s,-,s,1,有,因此,由于特征根,s,1,s,2,s,n,互异,有,下面讨论通过选择状态变量求得相应的状态空间模型,传递函数中极点互异时的变换(3/8),考虑到,输出,y,(,t,),和输入,u,(,t,),的拉氏变换满足,因此,若选择状态变量,x,i,(,t,)使其,拉氏变换满足,则,经反变换可得系统状态方程为,传递函数中极点互异时的变换(4/8),相应地,系统输出,y,(,t,),的拉氏变换为,Y,(,s,) =,k,1,X,1,(,s,)+,k,2,X,2,(,s,)+,k,n,X,n,(,s,),因此,经拉氏反变换可得如下输出方程,y,=,k,1,x,1,+,k,2,x,2,+,k,n,x,n,整理上述状态方程和输出方程可得如下,状态空间模型,传递函数中极点互异时的变换(5/8),上述用部分分式法建立的状态空间模型中的系统矩阵有一个重要特征,即,A,为对角线矩阵,系统矩阵,A,具有上述对角线形式的状态空间模型即为下一节将详细讨论的所谓对角线规范形,事实上,对角线规范形其实是将系统转换为,n,个一阶子系统(惯性环节)的并联, 如右图所示,对角线规范形的结构图,传递函数中极点互异时的变换(6/8),例,用部分分式法将,2.3.1,节例,子,中微分方程对应的下述传递函数变换为状态空间模型,传递函数中极点互异时的变换(7/8),解,由系统特征多项式,s,3,6,s,2,11,s,6,可求得系统极点为,s,1,1,s,2,2,s,3,3,于是有,其中,传递函数中极点互异时的变换(8/8),故当选择状态变量为,G,(,s,),分式并联分解的各个一阶惯性环节的输出,可得如下状态空间模型,将上述结果与,2.3.1,节例子,的结果相比较可知,即使对同一个系统,采用不同的建立状态空间模型的方法,将得到不同的状态空间模型。,即,状态空间模型不具有唯一性。,传递函数中有重极点时的变换(1/13),2.,传递函数中有重极点时的变换,当系统特征方程有重根时,传递函数不能分解成如式,的情况,亦得不到只有单极点情形时的对角线规范形状态方程,不失一般性, 为清楚地叙述变换方法, 以下设系统特征方程有6个根,其值分别为,s,1,s,1,s,1,s,4,s,5,s,5,即,s,1,为3重极点,s,5,为2重极点,.,相应地,用部分分式法可将所对应的传递函数表示为,传递函数中有重极点时的变换(2/13),其中,k,ij,为待定系数,其计算公式为,其中,l,为极点,s,i,的重数,.,传递函数中有重极点时的变换(3/13),下面以系数,k,13,的计算公式的推导为例来说明,k,ij,的计算式,将,G,(,s,),的两端乘以(,s,-,s,1,),3,有,对等式两边求2次导数后,因此，有,传递函数中有重极点时的变换(4/13),下面讨论通过选择状态变量求得相应的状态空间模型,如何选择状态变量,?,考虑到,输出,y,(,t,),和输入,u,(,t,),的拉氏变换满足,传递函数中有重极点时的变换(5/13),选择状态变量,x,i,(,t,)使其,拉氏变换满足,则有,传递函数中有重极点时的变换(6/13),即有,则经反变换可得系统状态方程为,传递函数中有重极点时的变换(7/13),相应地,系统输出,y,(,t,),的拉氏变换为,Y,(,s,),k,11,X,1,(,s,),k,12,X,2,(,s,),k,13,X,3,(,s,),k,41,X,4,(,s,),k,51,X,5,(,s,),k,52,X,6,(,s,),经拉氏反变换可得如下输出方程,y,k,11,x,1,k,12,x,2,k,13,x,3,k,41,x,4,k,51,x,5,k,52,x,6,传递函数中有重极点时的变换(8/13),因此,整理可得如下矩阵描述的状态空间模型,传递函数中有重极点时的变换(9/13),上述用部分分式法建立的状态空间模型中的系统矩阵有一个重要特征,即,A,为,块,对角矩阵,且每个矩阵方块为只有一个重特征值的特定矩阵块(约旦块)。,系统矩阵,A,具有上述特定块对角形式的状态空间模型即为下一节还将介绍的所谓约旦规范形。,事实上, 约旦规范形是将系统转换为多个子系统(惯性环节)的串-并联。,如下图所示。,传递函数中有重极点时的变换(10/13),传递函数中有重极点时的变换(11/13),例,用部分分式法将,下述传递函数变换为状态空间模型,传递函数中有重极点时的变换(12/13),解,由系统特征多项式,s,3,+5,s,2,+8,s,+4,可求得系统有二重极点,s,1,2,和单极点,s,2,1,于是有,其中,传递函数中有重极点时的变换(13/13),故,当选择状态变量为,G,(,s,),分式串,-,并联分解的各个一阶惯性环节的输出,可得如下状态空间模型,由系统方框图建立状态空间模型,(1/5),2.3.3,由系统方框图建立状态空间模型,下面讨论由系统方框图建立系统的状态空间模型,基本思想:,1. 将系统的各个环节分解为积分,惯性和比例环节,2. 选取积分环节和惯性环节的输出为系统状态,由系统方框图建立状态空间模型,(2/5),例,建立,下图示系统的状态空间模型,s,z,s,p,K,s,1,s,a,u,y,z,p,s,p,K,s,1,s,a,u,y,x,1,x,2,x,3,然后选择状态变量,解,:,先转化为积分,比例和惯性环节,由系统方框图建立状态空间模型,(3/5),再列出状态关系和输出关系,:,整理成状态空间模型形式便是,:,由系统方框图建立状态空间模型,(4/5),例,建立,下图示系统的状态空间模型,解,:,已经是积分和比例环节,.,选择状态变量,x,1,x,2,x,3,.,1,s,u,y,1,s,1,s,3,2,6,x,3,x,2,x,1,列出状态关系和输出关系就是,:,由系统方框图建立状态空间模型,(5/5),整理成状态空间模型形式便是,:,习题,习题,2-2,描述系统的微分方程为其中,y(t,),为输出量,u(t,),为输入量,.(1),选状态变量 写出系统的状态空间模型,;(2),选另一组变量 满足,:,写出系统以 为状态变量的状态空间模型,.,习题,2-3,求以下传递函数模型对应的状态空间模型,