状态空间分析法

1、第第8 8章章 状态空间分析法状态空间分析法 8.1 8.1 状态空间法的基本概念状态空间法的基本概念8.1.1 8.1.1 状态与状态变量状态与状态变量 “状态”和“状态空间”并不是新概念，长期以来，二者就在质点和刚体的古典力学中得到广泛应用。在经典控制理论中，“相平面”概念（在本书中未涉及）就是特殊的二维状态空间，但是，将“状态”和“状态空间”的概念在古典动力学基础上加以发展，并使之适合于控制过程的描述，这还是在20世纪60年代前后的事。下面通过举例说明状态、状态变量。【例8-1】无源网络如图8-1所示，试分析网络的数学模型。图8-1 RLC无源网络解：由图8-1可得方程ddddCCuCi

2、tiLRiuut 消去中间变量可得 （8-2） 用传递函数表示为 （8-3） 式（8-1）、式（8-3）均是描述该网络的数学模型，均可表示系统的状态。分析式（8-1）的两个一阶方程组 （8-4）用向量矩阵方程表示为 u （8-5）在该网络中，如果i(t0)、uC (t0)的初始值和tt0时输入电压均已知，则tt0时网络状态、uC (t)可完全确定，因此，可以说i(t)、uC (t)是这个二阶系统的一组状态变量，由此可得到如下概念。 22ddddCCCuuLCRCuutt2( )1( )1CUsU sLCsRCsd1dd11dCCCuuitCiRiuiutLLL 10011CCuuCiRiLLL

3、状态：即动力学系统状态，是指完整地和确定地描述系统的时域行为的一组最小变量。如果给定t=t0时刻这组变量的值和tt0时输入的时间函数，那么系统在tt0任何瞬时的行为就完全确定下来了，这样一组变量称为状态变量 8.1.2 8.1.2 状态向量与状态空间状态向量与状态空间（1）状态向量。状态向量是以状态变量为元组成的向量。如x1(t)，x2(t)、x3(t)，xn(t)是系统的一组状态变量，则状态向量就是以这组状态变量为分量的向量，即（2）状态空间。以x1(t)，x2(t)，x3(t)，xn(t)为坐标轴所组成的正交n维空间，称为状态空间，状态空间中的每一点，都代表状态变量的唯一和特定的一组值。由

4、系统的状态变量构成的一阶微分方程组，称为系统的状态方程。在例8-1中式（8-4）描述的就是系统的状态方程。如果将状态变量用一般符号来表示，即式（8-4）中令x1 = uC，x2 = i，并写成矩阵形式，则状态方程为12( )( )( )( )nx tx tX tx t （8-6）或 = AX + bu式中 ， ， 在指定系统输出的情况下，该输出与状态变量间的函数关系式，称为系统的输出方程。【例8-1】中指定x1 = uC作为输出，输出一般用y表示，则有y = uC，即y = x1，系统的输出方程为 （8-7）或 y = CT X式中 ，CT = 1 0 112210011xxCuxRxLLLX

5、12xXx101CARLL 01bL1210 xyx12xXx8.2 8.2 状态空间描述状态空间描述8.2.1 8.2.1 状态变量的选取状态变量的选取 对于状态变量的选择，原则上是任意选择的，但在有些系统的分析中，为了便于问题的分析或为了更进一步地了解系统内部变量及其之间的相互关系，状态变量应根据系统的要求来选取。下面举例说明如何选取状态变量。【例8-2】如图8-2所示的质量弹簧系统，当外力F(t)作用时，系统产生运动，质量及弹簧弹性系数见图，如不计算摩擦，试求：（1）以质量块m2的位移y(t)为输出量、外力F(t)为输入量，列写系统的微分方程； 图8-2 质量弹簧系统（2）自选一定数目的

6、状态变量，建立上述系统的状态方程和输出方程。 解：（1）设质量块m1的位移为z，根据牛顿定律有F(t) k1 (z y) = m1 （8-8）同理，对质量块m2有k1(z y) k2 y = m2 （8-9）zz联立式（8-8）、式（8-9）消去中间变量z，得出系统的微分方程为m1m2 y(4) + (k1 + k2)m1 + k1m2 + k1k2 y = k1F(t)（8-10）式（8-10）即为所求的系统的微分方程。（2）设状态变量x1 = z，x2 = = ，x3 = y，x4 = = 由式（8-8）可得由式（8-9）可得由此，状态方程和输出方程分别为 （8-11） （8-12） 即为

7、所求。yz1xy3x112131111( )kkxzxxF tmmm 12143122kkkxyxxmm 11112112133411242201000001( )00010000 xxkkmmxxmF txxxkkkxmm12340010 xxyxx8.2.2 8.2.2 线性系统的状态空间描述线性系统的状态空间描述 状态方程和输出方程综合起来，构成对一个系统完整的动态描述，称为系统的状态空间表达式，如式（8-6）和式（8-7）所示。而式（8-11）和式（8-12）就是图8-2所示系统的空间表达式。下面介绍列写状态空间描述的一般方法。设单输入单输出系统，其状态变量为x1，x2，xn，则状态方

8、程的一般形式为 = a11x1 + a12x2 +a1n + b1u = a21x1 + a22x2 +a2n + b2u = an1x1 + an2x2 +ann + bnu输出方程则有如下形式 y = c1x1 + c2x2 +cnxn用向量矩阵表示时的状态空间表达式为 = Ax + bu y = CTx （8-13）式中， n维状态变量1x2xnxx12nxxxx12nxxxx 系统内部状态的联系，称为系统矩阵，为矩阵； 输入对状态的作用，称为输入矩阵或控制矩阵，在此为n1的矩阵。111212122212nnnnnnaaaaaaAaaa12nbbbb8.3 8.3 能控性和能观测性能控性

9、和能观测性能控性和能观测性是系统的一种特性，这两个概念是卡耳曼在20世纪60年代提出的，是现代控制理论的两个基本概念。能控性检查每一状态分量能否被u(t)控制，是指控制作用对系统的影响能力；能观测性表示由观测量y能否判断状态X，它反映系统输出量确定系统状态的可能性。因此，能控性和能观测性从状态的控制能力和状态的识别能力两个方面反映系统本身的内在特性。实际上，现代控制理论中研究的许多问题，如最优控制、最佳估计等，都以能控性和能观测性作为其解存在的条件。实现能控性与能观测性面临的一个问题是，控制作用是否可使系统在有限的时间内，从起始状态导引到要求的状态；另一个问题是，能否通过观测有限时间内的输出量

10、而识别出系统的起始状态，从而识别系统的状态。能控性与能观测性之所以成为现代控制理论中的基本问题，是因为它着眼于对状态的控制。而经典控制理论是着眼于输出控制，受控过程可表示为一个复杂的高阶微分方程，即 y(n)+a1y (n1)+an1 +an y =b0u(n) +b1u(n1)+bn1 + b nu因为被控量y与控制作用u之间存在着明显的依赖关系，所以理论及实践上并不面临能否控制、能否观测的问题。但就系统的状态而言，这个问题仍客观存在，只是由于着眼于输出控制而掩盖了被控量与控制作用之间的依赖关系。在现代控制理论中，用状态空间方程来描述系统，通过对系统的状态方程及输出方程的分析，可以判断系统的

11、能控性和能观测性。也就是说，能控性、能观测性的条件是由系统的状态方程和输出方程的系数矩阵来确定的。8.3.1 8.3.1 线性系统的能控性线性系统的能控性能控性定义：线性系统 =A(t)X + B(t)U，在t0时刻的任意初始值X(t0)=X0，对tat0，taJ（J为系统的时间定义域），可找到容许控制U其元在t0,ta上平方可积），使X(ta)=0，则称系统在t0,ta上是状态能控的。yuX从该定义出发，可以加深对能控性的理解。 （1）系统的初始状态X0，是状态空间中任意非零的有限点，目标状态X(ta)为状态空间的原点。 （2）把系统从初始状态引向目际状态的控制作用，必须满足状态方程解存在唯

12、一性的条件。 （3）把系统由初始状态引向目标状态的时间定义为一个有限区间t0,ta。8.3.2 8.3.2 线性系统的能控性判断依据线性系统的能控性判断依据 在此只讨论线性单输入单输出系统的能控性问题。线性定常系统，即 =AX + BU Y = CX状态完全能控的充分必要条件是其能控性矩阵Q=B AB An1B满秩，即 rankB AB An1B = n对于单输入系统，能控性矩阵为方阵，所以其具体判断依据为Q的行列式不为零，即 |Q| 0该判断依据较简单，但当系统状态不完全能控时，它不能指明哪些状态不能控 X 【例8-3】考察如下系统的能控性。解：由状态方程可知， ， 所以 ， ，rankQ

13、= 3 故系统的状态完全可控。 8.3.3 8.3.3 线性系统的能观测性线性系统的能观测性 前面介绍了系统的状态能控性，下面介绍系统的状态能观测性。能观测定义 在tat0，t0J（J为系统的时间定义域）， 112233122201101011xxxxuxx 201B 122240110110111AB 2005A B 2240 010105QBABA B2400100105( )( )( )XA t XB t UYC t X 根据在t0,ta上的观测值Y(t)，tt0,ta区间能够唯一地确定系统在t0时刻的任意初始状态X0，则称系统在t0,ta上是状态可观测的。能观测性是研究状态和输出量的关

14、系，即通过对输出量在有限时间内的测量，能否把系统的状态识别出来，实质上，可归结为对初始状态识别的问题。8.3.4 8.3.4 线性系统的能观测性判断依据线性系统的能观测性判断依据线性定常系统，即 = AX Y = CX状态完全能观测的充分必要条件是能观测矩阵Qg = CT ATCT (AT )n1CT 满秩，即 满秩【例8-4】系统的动力学方程为Xg1nCCAQCA1122330100001061161xxxxuxx 123451 xyxx试判断该系统的能观测性。解：由系统的状态方程和输出方程可以得到C = 4 5 1 能观测矩阵为，所以rankQg = 23，因此系统不是状态完全能观测的。010451 0016716116CA g2451 671651CQCACA g4516710651Q 8.4 8.4 本本 章章 习习 题题1设系统微分方程为式中，u为输入量，x为输出量。设状态变量x1 = x， ，试列写系统的状态方程与输出方程。2判断下列系统的状