现代控制理论-第三章 2 能控性.pdf

1、第三章 线性系统的结构特性 1 3.2 系统的能控性 线性控制系统的能控性和能观性，是现代控制理论中两个非常重要的概念。1960 年美籍匈牙利人 R.E.Kalman 发表“控制系统的一般理论”等论文，引入状态空间法分析系统，提出能控性、能观测性、最佳调节器和卡尔曼滤波等概念，奠定了现代控制理论的基础。能控、能观性概念在刚提出来的时候，并未受到应有的重视，但现在已成为控制理论中的基本概念。无论在分析或综合一个现代控制系统时，总要研究一下，它是否能控与能观测。 现代控制理论，其数学模型采用状态方程，立足于状态变量法，包含最优控制、最优估计、系统辨识等理论。最优控制以能控性为基础，一般而言能控才能

2、得到最优解；最优估计以能观性为基础。 以下的讨论，均假定系统为线性定常（LTI）连续系统。 1.系统能控性的含义及例子 系统的能控性指系统的输入 )(tU 对系统的状态 )(tX 的控制能力，即衡量系统在 )(tU 作用下其内部状态转移的能力。 以后会知道，能控性是系统的一种内在性质，是系统的结构性质。 例 3-9 设 SISO 离散系统状态方程为 )(10)()(2001)1()1(2121kukxkxkxkx+=+，=11)0(X 试分析系统的输入 )(ku 对系统状态 )(kX 的影响。 解：用递推法解系统的状态方程，得 +=)1()0(22)1()()(121kuukxkxkkknul

3、l由上式可见，系统的第一个状态变量1x ，无 论 u 取何值，永远不受 u 的控制。所以，系统的状态不能通过它的输入 u 的作用而转移到任意所需的状态上去。 例 3-10 见讲义 P81 2.系统能控性的定义 设线性定常连续系统的状态方程为 矩阵。为矩阵，为，其中 mnBnnAtBUtAXtX += )()()(null系统能控性的定义为 ： 如果有一个控制作用 )(tU ，能在有限的时间 T 内，把系统从任意初始状态 0)0( X ，转移到终了状态 () 0XT= ,则称系统状态完全第三章 线性系统的结构特性 2 能控，简称系统能控。若系统中有一个或一些初始状态不能在有限的时间 T 内，在输

4、入控制 )(tU 的作用下转移到终了状态 () 0XT= ，则称系统状态不完全能控，简称系统不能控。 能控态，不能控态。 有限时间 T，存在 U(t)，状态能控，系统能控。 系统能达性的定义为 ： 若把系统初态规定为状态空间原点，即 0)0( =X ,终了状态规定为状态空间中的任意一个非0有限点 )(tX ， 若存在一个控制作用 )(tU ,能在有限时间 T 内使系统完成从初态到终态的转移，则称系统具有状态能达性。 讨论 ： 能控与能达的几何解释见图 3-11 所示。能控指在有限的时间 T 内，00X ；能达指在有限的时间 T 内，fX0 。 图 3-11 在线性定常连续系统中，能控性与能达性

5、是等价的。 推论： 如果存在控制作用 )(tU ，能在有限时间 T 内，把系统从任意初始状态0X 转移到有限终了状态fX ，则称系统状态完全能控。 有这么多的定义原因，是对线性定常连 续系统而言以上三个定义（甚至更多）等价，但对于其他类型的系统，这些定义则有可能不等价。比如对线性定常时间离散系统而言，系统的能控性与能达性就不等价。所以，有必要有许多概念、定义，以便更深入细致地描述系统的能控性质的差异。 能控性的含义是指控制作用 )(tU 对状态变量 )(tX 的影响程度。在第二章 3.1.4 节已指出 )(tX 的运动只能通过 )(tU 来影响，但未讨论 )(tU 对 )(tX 的影第三章 线

6、性系统的结构特性 3 响能力问题。能控性是控 制系统控制能力的标志， )(tU 对 )(tX 的影响能力由能控性来描述。 例3-8 已知系统的状态模拟图如图 3-12 所示，试判别其能控性。 图 3-12 解： 从图 3-12 可以看出， 系统的输入 u 只能影响系统的状态1x ， 不能影响2x（2x 趋于零的运动是自由运动，有限时间 T 内做不到） ，因此系统是不能控的。 例 3-9 利用定义判断图 3-13 所示电路系统的能控性 图 3-13 第三章 线性系统的结构特性 4 从上面的三个例子可以看出，对于简单系统来说，可以根据系统能控性的定义，从状态方程的解或状态模拟图来判断；对较复杂的系

7、统，求解往往很困难，其状态图一般也较复杂，这时就需要借助能控性判据来判别系统的能控性。 注意：判据皆从定义而来，定义是源、本。当忘记了判据等公式，或遇到特殊问题时，可以从定义入手解决问题，有时可能导致新的发现。这是面对“问题本身”的方法。 3.能控性判据之一 系统的状态能控性指的是系统的输入 )(tU 与系统的状态 )(tX 之间的关系，它与系统的输出 )(ty 无关。也就是说，系统的状态能控性仅与系统的状态方程BUAXX +=null有关，即只与 A、 B 阵有关。为了简单起见，把上面的状态方程简记作 A、 B 或 A、 B , 矩阵。为矩阵，为其中 mnBnnA 第三章 线性系统的结构特性

8、 5 第三章 线性系统的结构特性 6 (1) 定理（能控性判据） ： 线性定常连续系统 A、 B 状态完全能控的充分必要条件是其能控性矩阵 nmnnBABAABBS= 12null 满秩，即 nSrank = 。 第三章 线性系统的结构特性 7 第三章 线性系统的结构特性 8 1.定理（能控性判据） ： 线性定常连续系统 A、 B 状态完全能控的充分必要条件是其能控性矩阵 nmnnBABAABBS= 12null 满秩，即 nSrank = 。 证明： =+=TATTAATdBueXTXdBueXeTX00)()()0(0)(,)()0()( dBuAbAbIbTnn)()()()(01110

9、+= null 凯哈定理 =100)()(njTjjdubBA =10)(njijTBfA （1） 其中 =TjjmdubTf01,)()()( ，为一个确定量（向量） 将（1）写成矩阵形式： 11101)()()()0(=nmnnmnnTfTfTfBAABBXnullnull （2） 这实质上是由几个方程求 nm 个未知量的线性代数方程组。考虑到 )0(X 为任意的，则方程（2）的解110 nfff null、 存在的充分必要条件为： 能控性矩阵 nmnnBAABBS=1null 满秩，即 nSrank = ，或 0TSS . 若方程（2）的解存在，即110 nfff null、 存在，则

10、u(t)必定存在（由于if 是从状态方程的解（1）中得到的， （2）又表明if 存在，又jb 仅与 A 阵有关，故 u 必存在） ，系统A，B能控。矩阵 S 称为系统A，B的能控性矩阵。U(t)非唯一。 上面对 u(t)未加限制， 最优控制是求解某 u(t)使系统的某指标最小， 其时 u(t)应唯一。即最优控制问题是求 )()0( TXXf 中某指标 J 最小的 u(t)（用格拉姆矩阵证明则无 uf 的问题，在那，甚至可以构造出 u） 。 第三章 线性系统的结构特性 9 例 3-9 已知系统状态方程为 UXX+=111112310920231null，试判别系统状态是否完全能控。 解：系统的能

11、控性矩阵为 =4422114422114523122BAABBS 将上式中能控性矩阵的第二行加到第三行，则经过此初等变换后，有 =000000442211452312S 。 所以， nSrank = 32 。因为初等变换不改变矩阵的秩，所以系统状态不完全能控 例 310 试判别能控标准型状态的能控性。 解：在能控标准型状态方程中，有 nnnaaaaA=nullnullnullnullnullnullnullnullnull321100001000010，11000=nB null 所以，其能控性矩阵为 nnnnnaaBABAABBS=nullnullnullnullnullnull11012系

12、统的能控性矩阵 S 为下三角阵。因为 1=S ，所以 nSrank = 。即能控标准型的状态是完全能控的。 (2) 定理 ：若 SISO 线性定常连续系统状态完全能控，则存在一个非奇异变换矩阵Q, 可将该系统的原状态方程变换为能控标准型状态方程。 定理 ：若 SISO 线性定常连续系统的动态方程为： 第三章 线性系统的结构特性 10 CXYtBuAXX=+= )(null若系统能控，即能控性矩阵 nmnnBAABBS=1null 非奇异， 则存在一个非奇异变换矩阵nnQ,使 QXX = 而得到能控标准型的状态方程： XCYuBXAX=+=其中：=naaaQAQA211100010nullnul

13、lnullnull， =100nullQBB ，1=CQC 变换矩阵： =1nqAqAqQnull，11100= Sqnnull ，为行向量。 因此，能控标准型具有一般意义,即只要系统 能控，则可用能控标准型作为其模型而不失一般性。 综上，能控标准型一定能控，若能控也能化为能控标准型。 例 3-11 XYUXX 01,112101=+=null,将其化为能控标准型。 解 ： （1）先判断系统的能控性： = = ,23111 nrankSrank 状态完全能控。所以，可以化为能控标准型。 （2）求非奇异变换 Q： 第三章 线性系统的结构特性 11 = =32101012,102121212110

14、21212123311111111QAQAQqAqQSqS可得：=10QBB 121=CQC . 4.能控性判据之二 判别系统 A、 B 状态能控性的另一种方法是利用线性非奇异变换，将 A 阵对角化或约当化，然后根据变换后的 B 阵来判别系统的能控性。 (1) 定理：线性非奇异变换不改变系统的状态能控性。 证明 ： 变换前系统 A、 B 的能控性矩阵 S 为： BAABBSn 1= null 变换后系统 A、 B 的能控性矩阵 S 为： BABABSn 1= null 111111BPPAPBAPPPBPn = null 11BAABBPn= null SP1= 变换矩阵 P 非奇异，即 nP

15、rank = 。由线性代数有关定理知： 1SrankSrankSrankSPrankSrankSrankSPrankSrank=证毕。 第三章 线性系统的结构特性 12 (2) 定理：设系统 A、 B 具有两两相异的特征值，则其状态能控的充要条件是系统经线性非奇异变换后的对角标准型状态方程 UBXXn+=001nullnull的 B 阵中不包含元素全为零的行。其中： 11,ii imnbB bb bb=nullnull 证明 ：若要系统 A、 B 状态能控，其能控性矩阵 S ： nmnnnnnnnnnnbbbbbbbbbBABABS=12122221111111nullnullnullnull

16、nullnull 应满秩，即 nSrank = ， 当且仅当 nibinull,2,1,0 = 时才能成立，即 B 阵不能含有元素全为零的行。B 中全零行对应的iX 为不可控状态变量，tie为不可控模态。解释其含义：全零行，相应iX 不受 u 控制。 注意，当 A 阵含有重特征值，且仍能对角化时，上述结论不成立。在单输入系统中，更有结论：即使 B 中不含 0 元素， A、 B 也是状态不能控的，因为此时能控性矩阵 S 中至少有两行线性相关，所以 S 不满秩；在 MI 系统中，若重根对应的 B 中各行独力，则系统能控。 在对角标准型状态方程中，因为各状态变量已经解耦，所以各状态变量只能通过输入

17、)(tU 来直接进行控制。这样一来，一旦 B 阵中出现全零行，则相应的状态变量便不能控。 (3)定理：设 n维系统 A、 B 有重特征值为：11m为 重根，kkm为 重根，=kiinm1；且当 ji 时，有ji 。则系统状态完全能控的充要条件是在系统第三章 线性系统的结构特性 13 经线性非奇异变换后的约当标准形 UBXJJJXn+=nullnull 21中，和每个约当块iJ ),2,1( ki null= 的最后一行相对应的 B 阵中的那些相应行，其每行元素不全为零。 若两个约当块有相同特征值，上述结论不成立；若想要上述结论成立，则需要对应的 B 中相应行是线性独立的。 例 312 判别下列

18、各系统的状态能控性。 1 uXX+=312300020001null2 uXX+=021300020001null3 UXX+=312001300020001null4 UXX+=123100300020001null5 uXX+=112002null6 uXX+=210100020012null7 UXX+=300021100020012null解： 1 因为系统的状态方程为对角标准型， 其特征值互不相等， 且 B 阵无全零行。所以，系统状态完全能控。 2 因为系统的状态方程为对角标准型， 其特征值互不相等， 但 B 阵有全零行。所以，系统状态不完全能控。不能控的状态变量是3x 。 3 因为

19、系统的状态方程为对角标准型， 其特征值互不相等， 且 B 阵无全零行。所以，系统状态完全能控。 4 因为系统的状态方程为对角标准型， 其特征值互不相等， 但 B 阵有全零行。所以，系统状态不完全能控。不能控的状态变量是1x 。 5 因为系统的状态方程为对角标准型，其特征值相等，且为 SI 系统。所以，系统状态不完全能控。 第三章 线性系统的结构特性 14 6 因为系统的状态方程为约当标准型，其约当块的特征值互不相等，且 B 阵相应行无全零行。所以，系统状态完全能控。 7 因为系统的状态方程为约当标准型，其约当块的特征值互不相等，但 B 阵相应行有全零行。所以，系统状态不完全能控。不能控的状态变

20、量是2x 。 5.输出能控性 上面讨论了系统的状态能控性。对实际系统，很关心其输出，因此有必要讨论一下控制系统的输出能控性。 (1) 输出能控性的定义 ： 设线性定常连续控制系统的动态方程为 =+=CXYBUAXXnull，其中 U 为 m 维， Y 为 l维。 若存在容许控制 )(tU ， 能在有限时间 T(0T)内， 将系统由任意初始值 Y(0)引导到预先任意指定值 Y(T)，则称系统是输出能控的。 (2) 输出能控性判据 ： 系统 A,B,C输出能控的充要条件是其输出能控性矩阵 CB CAB CA1nBnml满秩，即有 l个输出时，其秩为 l。 (3) 系统的输出能控性与状态能控性不等价,即互不包含 。 例 313 判别下列系统的状态能控性和输出能控性。 1 uXX+=150154null2 uXX+=112001nullXy 11 = XY=2142解： 1 因为 2151255 = nrankABBrankSrank lrankCABCBrank = 1306 所以，系统状态不完全能控，但系统输出能控。 2 因为 nrankABBrankSrank = 22111 第三章 线性系统的结构特性 15 2153106= lrankCABCBrank 所以，系统状态完全能控，但系统输出不能控。