现代控制论基础课件1

《现代控制理论》讲稿贺廉云第1章控制系统的状态空间模型要点：1理解状态空间表示法概念；2掌握状态空间图示法；3掌握连续系统的数学模型转换；4了解多变量系统的传递函数阵及其求法难点：连续系统的数学模型转换一状态空间表示法1.基本术语状态：完全能描述系统时域行为的一个最少变量组。状态变量：是能构成系统状态的变量，能完全描述系统时域行为的一个最少变量组中的每一个变量。状态空间：状态向量X（t）的所有可能值的集合在几何学上叫状态空间。或说由x1轴、x2轴…xn轴所组成的n维空间称为状态空间。状态空间中的每一个点，对应于系统的某一特定状态。反过来，系统在任意时刻的状态都可用状态空间中的一个点来表示。显然，系统在不同时刻下的状态，可用状态空间中的一条轨迹表示。轨迹的形状，完全由系统在0t时刻的初态)(0tx0tt时的输入函数，以系统本身的动力学特性所决定。2状态空间模型的一般形式在现代控制理论中，状态空间模型所能描述的系统可以是单输入单输出的，也可以是多输入多输出的。状态空间表示式是一种采用状态描述系统动态行为（动态特性）的时域描述的数学模型。它包含状态方程输出方程。状态方程是一个一阶向量微分方程，输出方程是一个代数变换方程。U1状态变量y1U2y2……Umx1x2…x4yp图1-1系统表示描述某一动态的一个状态向量x（t）=[x1x2x3…xn]T（这里T为矩阵的转置），如图1-1所示。显然，该系统是n阶系统，若系统有m个输入u1，u2，u3，…，um，有p个输出y1，y2，y3，…，yp，且分别记u（t）=[u1u2u3…un]T和y(t)=[y1y2y3…yp]T位输入和输出向量。则系统的状态空间模型的一般形式为)),(),(()(ttutxftx（1-1）)),(),(()(ttutxty（1-2）式中，f=[f1f2f3…fn]T是n维函数向量；Φ是向量函数。式（1-1）是一阶向量微分方程，也可以看作由n个一阶微分方程所构成的方程组，称其为系统的状态方程；式（1-2）是一个代数方程，表示系统的输出量和输入量以及状态变量之间的关系，称之为系统的输出方程，或称为观测方程。这两个方程总称为系统的状态空间表达式。二状态空间模型的建立要建立状态空间表达式，必须先选取状态变量，状态变量一定要是系统中相互独立的变量。对于同一系统，状态变量选取的不同，所建立的状态空间表达式也不同，通常选取状态变量采取以下三种途径：（1）选择系统中贮能元件的输出物理量作为状态变量，然后根据系统的结构用物理定律列写出状态方程。（2）选择系统的输出及其各阶导数作为状态变量。（3）选择能使状态方程成为某种标准形式的变量作为状态变量。例1：如图1-2所示的电路，试以电压u为输出，以电容C上的电压uC为输出变量，列写其状态空间表达式。CR1R2L1L2图1-2例1电路解：图2电路的贮能元件有电感L1，L2和电容C。根据基尔霍夫定律列写电路方程：2221121222111110idtduCuiRiRiRdtdiLuiRiRdtdiLcc考虑到i1、i2、uc这三个变量是独立的，故可确定为系统的状态变量，经整理上式变为22222122121211111111iCdtduLuiLRRiLRdtdiuLiLRiLRdtdicc现在令x1=i1，x2=i2，x3=uc，将上式写成矩阵形式即为状态方程。uLxxxCLLRRLRLRLRxxx0010101013212221211111321由于前面已指出电容上的电压uc为输出变量，故系统的输出方程为321100xxxy由此可见，该电路的系统矩阵、控制矩阵、输出矩阵分别为uLBCLLRRLRLRLRA001,0101012221211111C=[001]三状态空间模型的图示法1.基本元件)(tx)(tx)(1tx)(3txxKx)(2tx(a)(b)(c)图1-3状态结构基本元件a—积分器b—加法器c—比例器2．一阶标量微分方程BuAxx的一阶系统状态结构图uxx图1-4一阶系统状态结构3.多输入多输出状态方程uxxy图1-5MIMO系统状态结构图∫Kba∫BA∫DC状态方程表达式为CxyBuAxx4.单输入单输出（SISO）线性定常系统微分方程的标准形式为ububububayayaymmmmnnnn1)1(1)(01)1(1)(（1-3）5能控标准形nnxxxx121=1110aaannnInnxxxx121+100u（1-4）][011bbbbymmnnxxxx121（1-5）只要系统状态方程的系数阵A和输出阵B具有式（1-4）的形式，C阵的形式可以任意，则称之为能控标准形，结构图如下nxunx1nx1xy∫∫b0-a1-a2b1-an∫bm…图1-6能控标准形状态结构图6能观标准形ubbbbxxxxaaaaxxxxmmnnnnnnnI011121121112100（1-6）][011bbbbymmnnxxxx121（1-7）只要系统状态空间表达式的A阵和C阵具有式（1-6）和（1-7）的形式，B阵的形式可以任意，则称之为能观标准形，结构图如图1-7所示。…U1x1x2xnxy…图1-7能观标准形状态结构图7对角标准形bm-1b0D=0-a1-an-1bm∫-an∫∫对角标准形的结构图和空间表达式分别如下：1x1xu2x2xynxnx图1-8对角标准形状态结构图uxxxxxxnnn111212121（1-8）nnxxxcccy2121][8约当标准形约当标准形的状态方程式和结构图分别如下C2∫∫λ2C1∫Cnλnλ1uxxxxxxxxxxnjjnjnjj11100111211111121（1-9）njjnjjxxxxxcccccy121111211][1x1x2x2xyujxjx图1-9约当标准形状态结构图四连续系统的数学模型转换对于动态系统，高阶微分方程、传递函数、状态方程表示的数学模型实际上是对系统动态过程的三种不同的形式的描述。微∫∫λ1C12C11∫C1jλ1λ1分方程和传递函数是古典控制理论中描述系统的数学模型，它们对事物外部特征进行描述，只反映系统输入输出之间的关系。而现代控制理论采用的状态空间表示，可以深入反映系统内容状态之间的关系。因此，两者之间存在着内在联系，可以通过适当的手段进行相互转换。1由状态空间模型转换成传递函数系统的状态方程DuCxyBuAxxLG（s）=DBAsIC1][=DAsIBAsIC*][(1-10)AsI是A阵的特征多项式*表示伴随矩阵例2已知某一单一输入输出系统的状态空间表达式为321xxx234100010321xxxu631]001[y321xxx试求其传递函数阵。解：根据式（1-10），可得G（s）=AsIBAsIC*][=63123410012341001]001[ssssss=4326311)2()32(232ssssssT=43235232sssss2传递函数阵的状态空间模型的实现（1）可控标准形的实现对于单输入单输出（SISO）系统，传递函数阵退化成传递函数。要把SISO系统式G(s)=nnnnmmmmasasasbsbsbsb1111110的传递函数形式转换成能控标准性的状态空间模型，即CxybuAxx（1-11）A=11100aaaInnnb=100(1-12)上述A阵是nn方阵，它的维数正好是传递函数的阶数，它的最后一行元素正还是传递函数分母（即系统的特征方程）所对应的稀疏，只不过均相差一个负号,其次对角线的元素均为1，其余为零，而b阵是一个列向量，最后一个元素为1，其余为零。正是b阵中的唯一的1对应友阵A的形式，是的输入信号u能对系统的每一个状态进行控制，因此称其为能控标准行。为了得到A阵和b阵的这种形式，应按下列规律选择状态变量：)1(21,,,nnyxyxyx,于是有uxaxaxaxxxxxnnnn12113221（1-13）现在的问题就是设法求出满足上述关系的输出矩阵c。为此设中间变量Z(s)，对于式G(s)=nnnnmmmmasasasbsbsbsb1111110分子分母同乘以Z(s),则有，设nnnnasasassUsZ111)()(则由nnnnmmmmasasasbsbsbsbsUsY1111110)()(有)()()(1110111mmmmnnnnbsbsbsbasasassUsY即)()()()()()(1111110sZasasassUsZbsbsbsbsYnnnnmmmm（1-14）对上式取拉氏变换zazazazuzbzbzbzbynnnnmmmm1)1(1)(1)1(1)(0（1-15）由式（1-13）有nmmmxxxzxxzxxzxz1)(43)3(322(1-16)将式（1-16）代入式（1-15），并令zx1，则121110xbxbxbxbymmnn（1-17）写成矩阵形式cxxxxbbbymmm2101（1-18）故所求的c阵满足式（1-17）的形式。同理，从式（1-15）的下式可以推出121121xaxaxaxaxunnnnn或写成uxaxaxaxaxnnnnn112211（1-19）可见，式（1-19）正好就是式（1-13）的最后一个表达式。式（1-13）和式（1-17）完全地描述了由传递函数G(s)=nnnnmmmmasasasbsbsbsb1111110所表示的系统。这一组能反映系统状态之间关系和输出变量的方程，也可以直接从能控标准形的结构图写出。对于输入函数中不含u的异数的特殊情况，或传递函数式G(s)=nnnnmmmm