状态和状态空间表达式-Read课件

根据系统的输入输出关系建立状态空间模型(2/2)由高阶常微分方程建立状态空间模型由传递函数建立状态空间模型根据状态变量图及系统方块图列写状态空间模型多输入多输出线性系统非线性系统根据系统的输入输出关系建立状态空间模型(2/2)1由高阶常微分方程建立状态空间模型(1/1)1.1由高阶常微分方程建立状态空间模型本节主要讨论由描述系统输入输出关系的常微分方程建立系统的状态空间模型,分别讨论由不含输入量导数项和由含输入量导数项的微分方程建立状态空间模型。由高阶常微分方程建立状态空间模型(1/1)1.1由高阶常2微分方程中不包含输入量的导数项(1/9)1.1.1微分方程中不包含输入量的导数项描述单输入单输出线性系统的输入输出间动态行为,不包含有输入量的导数项时的线性定系数常微分方程为y(n)+a1y(n-1)+…+any=bu其中y和u分别为系统的输出和输入;n为系统的阶次。这里所要研究的是建立上述常微分方程描述的动态系统的如下状态空间数学模型--状态空间模型本节问题的关键是如何选择状态变量。微分方程中不包含输入量的导数项(1/9)1.1.1微分方程3微分方程中不包含输入量的导数项(2/9)由微分方程理论知,若初始时刻t0的初值y(t0),y’(t0),…,y(n-1)(t0)已知,则对给定的输入u(t),微分方程有唯一解,也即系统在tt0的任何瞬时的动态都被唯一确定。因此,选择状态变量为如下变量x1(t)=y(t),x2(t)=y’(t),…,xn(t)=y(n-1)(t)可完全刻划系统的动态特性。取输出y和y的各阶导数(也称相变量)为状态变量,物理意义明确,易于接受。微分方程中不包含输入量的导数项(2/9)由微分方程理论知,若4微分方程中不包含输入量的导数项(3/9)将上述选择的状态变量代入输入输出的常微分方程,有如下状态方程和输出方程y=x1微分方程中不包含输入量的导数项(3/9)将上述选择的状态变量5微分方程中不包含输入量的导数项(4/9)将上述状态方程和输出方程写成矩阵形式有微分方程中不包含输入量的导数项(4/9)将上述状态方程和输出6微分方程中不包含输入量的导数项(5/9)该状态空间模型可简记为:其中微分方程中不包含输入量的导数项(5/9)该状态空间模型可简记7微分方程中不包含输入量的导数项(6/9)上述式子清楚说明了状态空间模型中系统矩阵A与微分方程中的系数a1,a2,…,an之间,输入矩阵B与方程中系数b之间的对应关系。通常将上述取输出y和y的各阶导数为状态变量称为相变量。微分方程中不包含输入量的导数项(6/9)上述式子清楚说明了状8微分方程中不包含输入量的导数项(8/9)-例2-1例

将以下系统输入输出方程变换为状态空间模型y”’+6y”+11y’+6y=6u解

本例中a1=6a2=11a3=6b=6因此,当选择输出y及其1阶与2阶导数等相变量为状态变量时,由课本式(1-4)和(1-7)可得状态空间模型如下P11例1.2微分方程中不包含输入量的导数项(8/9)-例2-1例将以下9微分方程中包含输入量的导数项(1/11)1.1.2.

微分方程中包含输入量的导数项描述单输入单输出线性系统的输入输出间动态行为的微分方程的一般表达式为y(n)+a1y(n-1)+…+any=b0u(n)+…+bnu所要研究的是建立上述常微分方程描述的动态系统的如下状态空间数学模型--状态空间模型建立该状态空间模型的关键是如何选择状态变量?微分方程中包含输入量的导数项(1/11)1.1.2.微分方10微分方程中包含输入量的导数项(2/11)若按照前面的方法那样选取相变量为状态变量,即x1(t)=y(t),x2(t)=y’(t),…,xn(t)=y(n-1)(t)则可得如下状态方程根据微分方程解的存在性和唯一性条件,要求输入u(t)为分段连续,而上述状态方程中输入u的各阶导数可能不连续,从而使微分方程解的存在性和唯一性的条件不成立。因此,状态方程中不应有输入u的导数项出现,即不能直接将输出y的各阶导数项取作状态变量。微分方程中包含输入量的导数项(2/11)若按照前面的方法那样11微分方程中包含输入量的导数项(3/11)为避免状态方程中显示地出现输入的导数,通常,可利用输出y和输入u以及其各阶导数的线性组合来组成状态变量,其原则是:使状态方程中不显含输出u的各阶导数。基于这种思路选择状态变量的方法很多,下面先介绍一种,其他的方法将在后续章节中陆续介绍。微分方程中包含输入量的导数项(3/11)为避免状态方程中显示12微分方程中包含输入量的导数项(4/11)根据上述原则,选择状态变量如下其中i(i=0,1,…,n)为待定系数。用anan-1an-2.....a1分别乘于上式两边，移项后可得：与P13方法有所不同，但本质一样微分方程中包含输入量的导数项(4/11)根据上述原则,选择状13微分方程中包含输入量的导数项(4/11)y(n)=xn+1+β0u(n)+β1u(n-1)+…+βnu以上两式左边和右边分别相加后，左边等于原线性方程的左边，所以右边相加的结果也应等于原线性方程的右边。由此可解得：微分方程中包含输入量的导数项(4/11)y(n)=xn+1+14微分方程中包含输入量的导数项(5/11)因此,有微分方程中包含输入量的导数项(5/11)因此,有15微分方程中包含输入量的导数项(7/11)则该高阶微分方程可转化描述为如下不含有输入导数项的状态空间模型微分方程中包含输入量的导数项(7/11)则该高阶微分方程可转16微分方程中包含输入量的导数项(9/11)-例2-2例

将以下系统输入输出方程变换为状态空间模型y”’+5y”+8y’+4y=2u”+14u’+24u解

本例中a1=5a2=8a3=4b0=0b1=2b2=14b3=24因此,有0=b0=01=b1-a10=22=b2-a11-a20=43=b3-a12-a21-a30=-12微分方程中包含输入量的导数项(9/11)-例2-2例将以下17微分方程中包含输入量的导数项(10/11)-例2-2即得系统的状态空间模型为P14例1.3,1.4同一个控制系统可以有不同的状态空间表达式。因为它们选取的状态变量不同！微分方程中包含输入量的导数项(10/11)-例2-2即得系统18由传递函数建立状态空间模型(1/6)1.2由传递函数建立状态空间模型下面讨论由描述系统输入输出关系的传递函数建立系统的状态空间模型。关键问题:1.如何选择状态变量2.保持系统的输入输出间的动态和静态关系不变由传递函数建立状态空间模型(1/6)1.2由传递函数建立状19线性定常微分方程由传递函数建立状态空间模型(2/6)由于传递函数与线性定系数常微分方程有直接的对应关系,故前面讨论的由高阶线性微分方程建立状态空间模型的方法同样适用于将传递函数建立变换为状态空间模型。类似地,本节讨论的由传递函数建立状态空间模型的方法亦适用于对微分方程建立状态空间模型。传递函数第一章第一节方法第一章第二节方法建立状态空间模型方法对线性定常系统拉氏变换线性定常微分方程由传递函数建立状态空间模型(2/6)由于传递20由传递函数建立状态空间模型(3/6)实际物理系统传递函数中分子多项式阶次小于或等于其分母多项式阶次,此时称该传递函数为真有理传递函数。而分子多项式阶次小于分母多项式阶次时,则称为严格真有理传递函数。单输入单输出(SISO)线性系统的输入输出间动态行为的如下传递函数由传递函数建立状态空间模型(3/6)实际物理系统传递函数中分21由传递函数建立状态空间模型(4/6)对上述传递函数,由长除法,有其中由传递函数建立状态空间模型(4/6)对上述传递函数,由长除法22由传递函数建立状态空间模型(5/6)本节所要研究的是建立该传递函数所描述的动态系统的状态空间模型(A,B,C,D)。上述常数项d即为状态空间模型(A,B,C,D)中的直联矩阵D;严格真有理传递函数G(s)对应可建立(A,B,C,D)中的(A,B,C)。即由传递函数建立状态空间模型(5/6)本节所要研究的是建立该传23由传递函数建立状态空间模型(6/6)下面分传递函数极点互异和有重极点两种情况讨论如何建立状态空间模型。由传递函数建立状态空间模型(6/6)下面分传递函数24传递函数中极点互异时的变换(1/8)1.2.1.

传递函数中极点互异时的变换对于传递函数G(s),其特征方程为sn+a1sn-1+…+an=0若其特征方程的n个特征根s1,s2,…,sn互异,则用部分分式法可将G(s)表示为如下并联分解其中k1,k2,…,kn为待定系数,其计算公式为传递函数中极点互异时的变换(1/8)1.2.1.传递函数中25传递函数中极点互异时的变换(3/8)考虑到,输出y(t)和输入u(t)的拉氏变换满足因此,若选择状态变量xi(t)使其拉氏变换满足则,经反变换可得系统状态方程为传递函数中极点互异时的变换(3/8)考虑到,输出y(t)和输26传递函数中极点互异时的变换(4/8)相应地,系统输出y(t)的拉氏变换为Y(s)=k1X1(s)+k2X2(s)+…+knXn(s)因此,经拉氏反变换可得如下输出方程y=k1x1+k2x2+…+knxn整理上述状态方程和输出方程可得如下状态空间模型传递函数中极点互异时的变换(4/8)相应地,系统输出y(t)27传递函数中极点互异时的变换(6/8)-例2-3例

用部分分式法将例2-1中微分方程对应的下述传递函数变换为状态空间模型传递函数中极点互异时的变换(6/8)-例2-3例用部分分式28传递函数中极点互异时的变换(7/8)解

由系统特征多项式s3+6s2+11s+6可求得系统极点为s1=-1s2=-2s3=-3于是有其中传递函数中极点互异时的变换(7/8)解由系统特征多项式其中29传递函数中极点互异时的变换(8/8)故当选择状态变量为G(s)分式并联分解的各个一阶惯性环节的输出,可得如下状态空间模型传递函数中极点互异时的变换(8/8)故当选择状态变量为G(s30传递函数中有重极点时的变换(1/13)1.2.2.传递函数中有重极点时的变换当系统特征方程有重根时,传递函数不能分解成如式的情况,亦得不到如式(2-26)所示的状态方程。不失一般性,为清楚地叙述变换方法,以下设系统特征方程有6个根,其值分别为s1,s1,s1,s4,s5,s5,即s1为3重极点,s2为2重极点。相应地,用部分分式法可将所对应的传递函数表示为传递函数中有重极点时的变换(1/13)1.2.2.传递函数中31传递函数中有重极点时的变换(2/13)其中kij为待定系数,其计算公式为其中l为极点si的重数。传递函数中有重极点时的变换(2/13)其中kij为待定系数,32传递函数中有重极点时的变换(4/13)下面讨论通过选择状态变量求得相应的状态空间模型。如何选择状态变量?考虑到,输出y(t)和输入u(t)的拉氏变换满足传递函数中有重极点时的变换(4/13)下面讨论通过选择状态变33传递函数中有重极点时的变换(5/13)选择状态变量xi(t)使其拉氏变换满足则有传递函数中有重极点时的变换(5/13)选择状态变量xi(t)34传递函数中有重极点时的变换(6/13)即有则经反变换可得系统状态方程为传递函数中有重极点时的变换(6/13)即有则经反变换可得系统35传递函数中有重极点时的变换(7/13)相应地,系统输出y(t)的拉氏变换为Y(s)=k11X1(s)+k12X2(s)+k13X3(s)+k41X4(s)+k51X5(s)+k52X6(s)经拉氏反变换可得如下输出方程y=k11x1+k12x2+k13x3+k41x4+k51x5+k52x6传递函数中有重极点时的变换(7/13)相应地,系统输出y(t36传递函数中有重极点时的变换(8/13)因此,整理可得如下矩阵描述的状态空间模型传递函数中有重极点时的变换(8/13)因此,整理可得如下矩阵37传递函数中有重极点时的变换(11/13)-例2-4例

用部分分式法将下述传递函数变换为状态空间模型传递函数中有重极点时的变换(11/13)-例2-4例用部分38传递函数中有重极点时的变换(12/13)解由系统特征多项式s3+5s2+8s+4可求得系统有二重极点s1=-2和单极点s2=-1,于是有其中传递函数中有重极点时的变换(12/13)解由系统特征多项式39传递函数中有重极点时的变换(13/13)故当选择状态变量为G(s)分式串-并联分解的各个一阶惯性环节的输出,可得如下状态空间模型传递函数中有重极点时的变换(13/13)故当选择状态变量为G40多输入多输出线性系统(1/5)1.3根据状态变量图及方块图列写状态空间表达式a.状态变量图一般由积分器、放大器和加法器组成.b.一般取积分环节的输出作为状态变量.遇到的问题：a.如何根据已有的状态变量图列写状态空间表达式？b.如何根据已有的微分方程或传递函数，转化为状态变量图，从而列写状态空间表达式？多输入多输出线性系统(1/5)1.3根据状态变量图及方块图41多输入多输出线性系统(1/5)多输入多输出线性系统(1/5)42多输入多输出线性系统(1/5)等价于输出项无导数项的微分方程多输入多输出线性系统(1/5)等价于输出项无导数项的微分方程43多输入多输出线性系统(1/5)设描述系统的微分方程为此系统为MIMS系统，用前述基于微分方程或传递函数的方法显然难于建立状态空间模型，同SISO系统一样,该系统的实现也是非唯一的。下面采用状态变量图方法来建立状态空间模型。多输入多输出线性系统(1/5)设描述系统的微分方程为此系统44多输入多输出线性系统(2/5)该系统的方程也可表示为对每一个方程积分,直至消除导数符号为止。为此,有多输入多输出线性系统(2/5)该系统的方程也可表示为对每一个45多输入多输出线性系统(3/5)故可得模拟结构图,如图所示。系统模拟结构图多输入多输出线性系统(3/5)故可得模拟结构图,如图所示。系46多输入多输出线性系统(4/5)取每个积分器的输出为一个状态变量,如图2-13所示。则式(2-33)的一种状态空间实现为相应地输出方程为多输入多输出线性系统(4/5)取每个积分器的输出为一个状态变47多输入多输出线性系统(5/5)因此,该双输入双输出系统的矩阵形式状态空间模型为多输入多输出线性系统(5/5)因此,该双输入双输出系统的矩阵48多输入多输出线性系统(1/5)设描述系统的传递函数为：如何把之转换为状态变量图，并列写出状态空间模型？？多输入多输出线性系统(1/5)如何把之转换为状态变量图，并列49多输入多输出线性系统(1/5)

分子分母同除于s的n次方:

引入中间变量Z(s):

引入中间变量Z(s):多输入多输出线性系统(1/5)分子分母同除于s的n次方:50多输入多输出线性系统(1/5)

可得：参看课本P20例1.6多输入多输出线性系统(1/5)可得：参看课本P20例1.651微分方程中包含输入量的导数项(7/11)则状态空间模型为：参看课本P15-16,比较转换过程和转换结果微分方程中包含输入量的导数项(7/11)则状态空间模型为：参52多输入多输出线性系统(1/5)小结：多输入多输出线性系统(1/5)小结：53多输入多输出线性系统(1/5)1.4多输入多输出线性系统下面,以双输入双输出的三阶系统为例介绍由描述MIMO系统的高阶微分方程组如何建立状态空间模型。设描述系统的微分方程为同SISO系统一样,该系统的实现也是非唯一的。下面采用模拟结构图的方法,按高阶导数项求解方法来建立状态空间模型。多输入多输出线性系统(1/5)1.4多输入多输出线性系统同54多输入多输出线性系统(2/5)因此,该系统的方程也可表示为对每一个方程积分,直至消除导数符号为止。为此,有多输入多输出线性系统(2/5)因此,该系统的方程也可表示为对55多输入多输出线性系统(3/5)故可得模拟结构图,如图所示。系统模拟结构图多输入多输出线性系统(3/5)故可得模拟结构图,如图所示。系56多输入多输出线性系统(4/5)取每个积分器的输出为一个状态变量,如图2-13所示。则式(2-33)的一种状态空间实现为相应地输出方程为多输入多输出线性系统(4/5)取每个积分器的输出为一个状态变57多输入多输出线性系统(5/5)因此,该双输入双输出系统的矩阵形式状态空间模型为P17-21例1.5已知系统方块图的MIMO系统的状态空间转换方法。多输入多输出线性系统(5/5)因此,该双输入双输出系统的矩阵58非线性系统(1/10)1.5非线性系统倒立摆系统是一个多变量、存在严重非线性的非自治不稳定性系统,经常被用来研究和比较各种控制方法的性能。其结构和飞机着陆、火箭飞行及机器人的关节运动等有很多相似之处,因而对倒立摆系统平衡的控制方法在航空及机器人等领域有着广泛的用途,人们对倒立摆控制的研究也越来越感兴趣。下面通过一个一级倒立摆的例子,来简述对非线性系统来说,如何通过描述其动力学模型的常微分方程建立状态空间模型。非线性系统(1/10)1.5非线性系统59非线性系统(2/10)图为某一级倒立摆结构示意图。图

一级倒立摆示意图非线性系统(2/10)图为某一级倒立摆结构示意图。图一级倒60非线性系统(3/10)图中所示的带轮小车可以前后移动来平衡一根杆,此杆由其底部的一个支点来支撑。该系统中还有一个电机,一根连接电机与小车的皮带和一些滑轮。还有一些传感器,用来测量小车的速度、位置、杆底部与铅垂线所成的角度及其微分。其控制任务是由电机通过皮带施加合适的力f给小车从而使杆不倒,并使小车不超过左右边界。一级倒立摆有两个运动自由度,一个沿水平方向运动,另一个绕轴转动。非线性系统(3/10)图中所示的带轮小车可以前后移动来平衡一61非线性系统(4/10)解通过对滑轮小车和摆竿的受力分析和推导,且忽略交流电机的动特性并且假设交流电机由u到f的静态增益为1,得到倒立摆系统的动力学描述如下:其中c是小车与导轨的摩擦系数;f为施加在小车水平方向上的外力;u为作用在驱动电机上的电压,其为控制变量;非线性系统(4/10)解通过对滑轮小车和摆竿的受力分析和推62非线性系统(5/10)J为转动惯量,x为小车的水平位移,由与电机相连的电位计测得;

连接的电位计测得的信号经微分而得;为杆与垂线的夹角,并取顺时针方向为正方向,由安装在小车上并与杆的基座相连的电位计测得;

为小车的水平速度,由与电机

为杆转动的角速度,由安装在小车上并与杆的基座相连的电位计测得的信号经微分而得。非线性系统(5/10)J为转动惯量,连接的电位计测得63非线性系统(6/10)整理上式,得到:其中非线性系统(6/10)整理上式,得到:其中64非线性系统(7/10)对该倒立摆系统,选取状态变量:由上式得到该倒立摆系统的状态空间模型为可以明显的看到，状态方程是非线性化的非线性系统(7/10)对该倒立摆系统,选取状态变量:由上式65非线性系统(8/10)由于数学方法的局限以及工程系统实现的困难,在进行系统分析与控制时,复杂的非线性模型将导致难于分析求解及控制。因此,常将非线性模型在其平衡点(工作点)附近对其进行Taylor级数展开至一阶线性方程,以获得简化的数学模型,实现系统分析与控制。这种处理也是工程中的常用方法,如若摆杆相对于垂直线的角度保持足够小(如),则常有如下线性展开近似非线性系统(8/10)由于数学方法的局限以及工程系统实现的困66非线性系统(9/10)因此,对本例来说,在平衡点

附近,非线性状态方程的近似线性化状态方程为

其中非线性系统(9/10)因此,对本例来说,在平衡点附近,非线67非线性系统(10/10)和相应的输出方程:

至此,得到一级倒立摆系统的状态空间形式的线性化数学模型。此例为本征非线性，而课本P28例1.13为本质非线性非线性系统(10/10)和相应的输出方程:至此,得到一级倒68根据系统的输入输出关系建立状态空间模型(2/2)由高阶常微分方程建立状态空间模型由传递函数建立状态空间模型根据状态变量图及系统方块图列写状态空间模型多输入多输出线性系统非线性系统根据系统的输入输出关系建立状态空间模型(2/2)69由高阶常微分方程建立状态空间模型(1/1)1.1由高阶常微分方程建立状态空间模型本节主要讨论由描述系统输入输出关系的常微分方程建立系统的状态空间模型,分别讨论由不含输入量导数项和由含输入量导数项的微分方程建立状态空间模型。由高阶常微分方程建立状态空间模型(1/1)1.1由高阶常70微分方程中不包含输入量的导数项(1/9)1.1.1微分方程中不包含输入量的导数项描述单输入单输出线性系统的输入输出间动态行为,不包含有输入量的导数项时的线性定系数常微分方程为y(n)+a1y(n-1)+…+any=bu其中y和u分别为系统的输出和输入;n为系统的阶次。这里所要研究的是建立上述常微分方程描述的动态系统的如下状态空间数学模型--状态空间模型本节问题的关键是如何选择状态变量。微分方程中不包含输入量的导数项(1/9)1.1.1微分方程71微分方程中不包含输入量的导数项(2/9)由微分方程理论知,若初始时刻t0的初值y(t0),y’(t0),…,y(n-1)(t0)已知,则对给定的输入u(t),微分方程有唯一解,也即系统在tt0的任何瞬时的动态都被唯一确定。因此,选择状态变量为如下变量x1(t)=y(t),x2(t)=y’(t),…,xn(t)=y(n-1)(t)可完全刻划系统的动态特性。取输出y和y的各阶导数(也称相变量)为状态变量,物理意义明确,易于接受。微分方程中不包含输入量的导数项(2/9)由微分方程理论知,若72微分方程中不包含输入量的导数项(3/9)将上述选择的状态变量代入输入输出的常微分方程,有如下状态方程和输出方程y=x1微分方程中不包含输入量的导数项(3/9)将上述选择的状态变量73微分方程中不包含输入量的导数项(4/9)将上述状态方程和输出方程写成矩阵形式有微分方程中不包含输入量的导数项(4/9)将上述状态方程和输出74微分方程中不包含输入量的导数项(5/9)该状态空间模型可简记为:其中微分方程中不包含输入量的导数项(5/9)该状态空间模型可简记75微分方程中不包含输入量的导数项(6/9)上述式子清楚说明了状态空间模型中系统矩阵A与微分方程中的系数a1,a2,…,an之间,输入矩阵B与方程中系数b之间的对应关系。通常将上述取输出y和y的各阶导数为状态变量称为相变量。微分方程中不包含输入量的导数项(6/9)上述式子清楚说明了状76微分方程中不包含输入量的导数项(8/9)-例2-1例

将以下系统输入输出方程变换为状态空间模型y”’+6y”+11y’+6y=6u解

本例中a1=6a2=11a3=6b=6因此,当选择输出y及其1阶与2阶导数等相变量为状态变量时,由课本式(1-4)和(1-7)可得状态空间模型如下P11例1.2微分方程中不包含输入量的导数项(8/9)-例2-1例将以下77微分方程中包含输入量的导数项(1/11)1.1.2.

微分方程中包含输入量的导数项描述单输入单输出线性系统的输入输出间动态行为的微分方程的一般表达式为y(n)+a1y(n-1)+…+any=b0u(n)+…+bnu所要研究的是建立上述常微分方程描述的动态系统的如下状态空间数学模型--状态空间模型建立该状态空间模型的关键是如何选择状态变量?微分方程中包含输入量的导数项(1/11)1.1.2.微分方78微分方程中包含输入量的导数项(2/11)若按照前面的方法那样选取相变量为状态变量,即x1(t)=y(t),x2(t)=y’(t),…,xn(t)=y(n-1)(t)则可得如下状态方程根据微分方程解的存在性和唯一性条件,要求输入u(t)为分段连续,而上述状态方程中输入u的各阶导数可能不连续,从而使微分方程解的存在性和唯一性的条件不成立。因此,状态方程中不应有输入u的导数项出现,即不能直接将输出y的各阶导数项取作状态变量。微分方程中包含输入量的导数项(2/11)若按照前面的方法那样79微分方程中包含输入量的导数项(3/11)为避免状态方程中显示地出现输入的导数,通常,可利用输出y和输入u以及其各阶导数的线性组合来组成状态变量,其原则是:使状态方程中不显含输出u的各阶导数。基于这种思路选择状态变量的方法很多,下面先介绍一种,其他的方法将在后续章节中陆续介绍。微分方程中包含输入量的导数项(3/11)为避免状态方程中显示80微分方程中包含输入量的导数项(4/11)根据上述原则,选择状态变量如下其中i(i=0,1,…,n)为待定系数。用anan-1an-2.....a1分别乘于上式两边，移项后可得：与P13方法有所不同，但本质一样微分方程中包含输入量的导数项(4/11)根据上述原则,选择状81微分方程中包含输入量的导数项(4/11)y(n)=xn+1+β0u(n)+β1u(n-1)+…+βnu以上两式左边和右边分别相加后，左边等于原线性方程的左边，所以右边相加的结果也应等于原线性方程的右边。由此可解得：微分方程中包含输入量的导数项(4/11)y(n)=xn+1+82微分方程中包含输入量的导数项(5/11)因此,有微分方程中包含输入量的导数项(5/11)因此,有83微分方程中包含输入量的导数项(7/11)则该高阶微分方程可转化描述为如下不含有输入导数项的状态空间模型微分方程中包含输入量的导数项(7/11)则该高阶微分方程可转84微分方程中包含输入量的导数项(9/11)-例2-2例

将以下系统输入输出方程变换为状态空间模型y”’+5y”+8y’+4y=2u”+14u’+24u解

本例中a1=5a2=8a3=4b0=0b1=2b2=14b3=24因此,有0=b0=01=b1-a10=22=b2-a11-a20=43=b3-a12-a21-a30=-12微分方程中包含输入量的导数项(9/11)-例2-2例将以下85微分方程中包含输入量的导数项(10/11)-例2-2即得系统的状态空间模型为P14例1.3,1.4同一个控制系统可以有不同的状态空间表达式。因为它们选取的状态变量不同！微分方程中包含输入量的导数项(10/11)-例2-2即得系统86由传递函数建立状态空间模型(1/6)1.2由传递函数建立状态空间模型下面讨论由描述系统输入输出关系的传递函数建立系统的状态空间模型。关键问题:1.如何选择状态变量2.保持系统的输入输出间的动态和静态关系不变由传递函数建立状态空间模型(1/6)1.2由传递函数建立状87线性定常微分方程由传递函数建立状态空间模型(2/6)由于传递函数与线性定系数常微分方程有直接的对应关系,故前面讨论的由高阶线性微分方程建立状态空间模型的方法同样适用于将传递函数建立变换为状态空间模型。类似地,本节讨论的由传递函数建立状态空间模型的方法亦适用于对微分方程建立状态空间模型。传递函数第一章第一节方法第一章第二节方法建立状态空间模型方法对线性定常系统拉氏变换线性定常微分方程由传递函数建立状态空间模型(2/6)由于传递88由传递函数建立状态空间模型(3/6)实际物理系统传递函数中分子多项式阶次小于或等于其分母多项式阶次,此时称该传递函数为真有理传递函数。而分子多项式阶次小于分母多项式阶次时,则称为严格真有理传递函数。单输入单输出(SISO)线性系统的输入输出间动态行为的如下传递函数由传递函数建立状态空间模型(3/6)实际物理系统传递函数中分89由传递函数建立状态空间模型(4/6)对上述传递函数,由长除法,有其中由传递函数建立状态空间模型(4/6)对上述传递函数,由长除法90由传递函数建立状态空间模型(5/6)本节所要研究的是建立该传递函数所描述的动态系统的状态空间模型(A,B,C,D)。上述常数项d即为状态空间模型(A,B,C,D)中的直联矩阵D;严格真有理传递函数G(s)对应可建立(A,B,C,D)中的(A,B,C)。即由传递函数建立状态空间模型(5/6)本节所要研究的是建立该传91由传递函数建立状态空间模型(6/6)下面分传递函数极点互异和有重极点两种情况讨论如何建立状态空间模型。由传递函数建立状态空间模型(6/6)下面分传递函数92传递函数中极点互异时的变换(1/8)1.2.1.

传递函数中极点互异时的变换对于传递函数G(s),其特征方程为sn+a1sn-1+…+an=0若其特征方程的n个特征根s1,s2,…,sn互异,则用部分分式法可将G(s)表示为如下并联分解其中k1,k2,…,kn为待定系数,其计算公式为传递函数中极点互异时的变换(1/8)1.2.1.传递函数中93传递函数中极点互异时的变换(3/8)考虑到,输出y(t)和输入u(t)的拉氏变换满足因此,若选择状态变量xi(t)使其拉氏变换满足则,经反变换可得系统状态方程为传递函数中极点互异时的变换(3/8)考虑到,输出y(t)和输94传递函数中极点互异时的变换(4/8)相应地,系统输出y(t)的拉氏变换为Y(s)=k1X1(s)+k2X2(s)+…+knXn(s)因此,经拉氏反变换可得如下输出方程y=k1x1+k2x2+…+knxn整理上述状态方程和输出方程可得如下状态空间模型传递函数中极点互异时的变换(4/8)相应地,系统输出y(t)95传递函数中极点互异时的变换(6/8)-例2-3例

用部分分式法将例2-1中微分方程对应的下述传递函数变换为状态空间模型传递函数中极点互异时的变换(6/8)-例2-3例用部分分式96传递函数中极点互异时的变换(7/8)解

由系统特征多项式s3+6s2+11s+6可求得系统极点为s1=-1s2=-2s3=-3于是有其中传递函数中极点互异时的变换(7/8)解由系统特征多项式其中97传递函数中极点互异时的变换(8/8)故当选择状态变量为G(s)分式并联分解的各个一阶惯性环节的输出,可得如下状态空间模型传递函数中极点互异时的变换(8/8)故当选择状态变量为G(s98传递函数中有重极点时的变换(1/13)1.2.2.传递函数中有重极点时的变换当系统特征方程有重根时,传递函数不能分解成如式的情况,亦得不到如式(2-26)所示的状态方程。不失一般性,为清楚地叙述变换方法,以下设系统特征方程有6个根,其值分别为s1,s1,s1,s4,s5,s5,即s1为3重极点,s2为2重极点。相应地,用部分分式法可将所对应的传递函数表示为传递函数中有重极点时的变换(1/13)1.2.2.传递函数中99传递函数中有重极点时的变换(2/13)其中kij为待定系数,其计算公式为其中l为极点si的重数。传递函数中有重极点时的变换(2/13)其中kij为待定系数,100传递函数中有重极点时的变换(4/13)下面讨论通过选择状态变量求得相应的状态空间模型。如何选择状态变量?考虑到,输出y(t)和输入u(t)的拉氏变换满足传递函数中有重极点时的变换(4/13)下面讨论通过选择状态变101传递函数中有重极点时的变换(5/13)选择状态变量xi(t)使其拉氏变换满足则有传递函数中有重极点时的变换(5/13)选择状态变量xi(t)102传递函数中有重极点时的变换(6/13)即有则经反变换可得系统状态方程为传递函数中有重极点时的变换(6/13)即有则经反变换可得系统103传递函数中有重极点时的变换(7/13)相应地,系统输出y(t)的拉氏变换为Y(s)=k11X1(s)+k12X2(s)+k13X3(s)+k41X4(s)+k51X5(s)+k52X6(s)经拉氏反变换可得如下输出方程y=k11x1+k12x2+k13x3+k41x4+k51x5+k52x6传递函数中有重极点时的变换(7/13)相应地,系统输出y(t104传递函数中有重极点时的变换(8/13)因此,整理可得如下矩阵描述的状态空间模型传递函数中有重极点时的变换(8/13)因此,整理可得如下矩阵105传递函数中有重极点时的变换(11/13)-例2-4例

用部分分式法将下述传递函数变换为状态空间模型传递函数中有重极点时的变换(11/13)-例2-4例用部分106传递函数中有重极点时的变换(12/13)解由系统特征多项式s3+5s2+8s+4可求得系统有二重极点s1=-2和单极点s2=-1,于是有其中传递函数中有重极点时的变换(12/13)解由系统特征多项式107传递函数中有重极点时的变换(13/13)故当选择状态变量为G(s)分式串-并联分解的各个一阶惯性环节的输出,可得如下状态空间模型传递函数中有重极点时的变换(13/13)故当选择状态变量为G108多输入多输出线性系统(1/5)1.3根据状态变量图及方块图列写状态空间表达式a.状态变量图一般由积分器、放大器和加法器组成.b.一般取积分环节的输出作为状态变量.遇到的问题：a.如何根据已有的状态变量图列写状态空间表达式？b.如何根据已有的微分方程或传递函数，转化为状态变量图，从而列写状态空间表达式？多输入多输出线性系统(1/5)1.3根据状态变量图及方块图109多输入多输出线性系统(1/5)多输入多输出线性系统(1/5)110多输入多输出线性系统(1/5)等价于输出项无导数项的微分方程多输入多输出线性系统(1/5)等价于输出项无导数项的微分方程111多输入多输出线性系统(1/5)设描述系统的微分方程为此系统为MIMS系统，用前述基于微分方程或传递函数的方法显然难于建立状态空间模型，同SISO系统一样,该系统的实现也是非唯一的。下面采用状态变量图方法来建立状态空间模型。多输入多输出线性系统(1/5)设描述系统的微分方程为此系统112多输入多输出线性系统(2/5)该系统的方程也可表示为对每一个方程积分,直至消除导数符号为止。为此,有多输入多输出线性系统(2/5)该系统的方程也可表示为对每一个113多输入多输出线性系统(3/5)故可得模拟结构图,如图所示。系统模拟结构图多输入多输出线性系统(3/5)故可得模拟结构图,如图所示。系114多输入多输出线性系统(4/5)取每个积分器的输出为一个状态变量,如图2-13所示。则式(2-33)的一种状态空间实现为相应地输出方程为多输入多输出线性系统(4/5)取每个积分器的输出为一个状态变115多输入多输出线性系统(5/5)因此,该双输入双输出系统的矩阵形式状态空间模型为多输入多输出线性系统(5/5)因此,该双输入双输出系统的矩阵116多输入多输出线性系统(1/5)设描述系统的传递函数为：如何把之转换为状态变量图，并列写出状态空间模型？？多输入多输出线性系统(1/5)如何把之转换为状态变量图，并列117多输入多输出线性系统(1/5)

分子分母同除于s的n次方:

引入中间变量Z(s):

引入中间变量Z(s):多输入多输出线性系统(1/5)分子分母同除于s的n次方:118多输入多输出线性系统(1/5)

可得：参看课本P20例1.6多输入多输出线性系统(1/5)可得：参看课本P20例1.6119微分方程中包含输入量的导数项(7/11)则状态空间模型为：参看课本P15-16,比较转换过程和转换结果微分方程中包含输入量的导数项(7/11)则状态空间模型为：参120多输入多输出线性系统(1/5)小结：多输入多输出线性系统(1/5)小结：121多输入多输出线性系统(1/5)1.4多输入多输出线性系统下面,以双输入双输出的三阶系统为例介绍由描述MIMO系统的高阶微分方程组如何建立状态空间模型。设描述系统的微分方程为同SISO系统一样,该系统的实现也是非唯一的。下面采用模拟结构图的方法,按高阶导数项求解方法来建立状态空间模型。多输入多输出线性系统(1/5)1.4多输入多输出线性系统同122多输入多输出线性系统(2/5)因此,该系统的方程也可表示为对每一个方程积分,直至消除导数符号为止。为此,有多输入多输出线性系统(2/5)因此,该系统的方