上海交大杜秀华老师《现代控制理论》第四章 线性系统的能控性和能观性6

三、 多输入多输出系统的标准形旺纳姆（Wonham ）标准形和龙伯格（Luenberger ）标准形。

1.多输入多输出系统的能控标准形考虑线性定常系统:Σx Ax Bu y Cx=+=x 为n 维状态向量，u 为p 输入向量，y 为q 维输出向量如果系统能控，则系统的能控性矩阵的秩为n ，即cQ 中有n 个线性无关列。

111121212[]c Q b b b Ab Ab Ab A b A b A b n n n p p p ---=对多输入系统，1p >，c Q 中有np 列，所以，在c Q 中可以找出很多种n 个线性无关列的情况。

这里介绍两种寻找n 个线性无关列的方法，以构成状态变换阵，将状态空间描述形式变换为旺纳姆能控标准形和龙伯格能控标准形。

定理 [旺纳姆能控标准形]对完全能控的线性定常系统，存在线性非奇异变换1x Px Q x -==使状态空间表达式转化为旺纳姆能控标准形:Σx A x B u y C xcW c c c =+=式中111211222A A A 0A A A Q AQ 00A m m c mm -⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎣⎦11(11(1)11010011,2,,0001A ,i i ii i mννννααα⨯-⎡⎤⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥---⎢⎥⎣⎦）112()00,1,,00A i j ij ij ij ij j i m νννγγγ⨯⎡⎤⎢⎥⎢⎥==+⎢⎥⎢⎥⎢⎥⎣⎦1(1)11(1)001001B Q B m p c n m np ββββ+-+⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦C CQ c =（无特殊形式）证明：见书 例 求如下系统的旺纳姆能控标准形121100*********A B -⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦计算系统的能控性矩阵2101204010101001042BABA B c Q ⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥⎣⎦3c Q rank =，系统完全能控。

按2111{,,} b A b A b 和2222{,,} b Ab A b 搜索，可知112,,b Ab b 是线性独立的。

得1,221==νν。

21111121211111121100110,0,0,000014410A b Ab b b Ab A b αααα=--⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥====--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦解得 441211=αα，＝－111111121e b Ab e b α=+=得111212121131040000101100e e b e b -⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=-+=====⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎣⎦由22111211112212-Ab b e e αγγ=++，即21121122203111000010αγγ-⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=-++⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦解得 20,121==-=122121 , γγα由此可得1111122201000-44210000101A A A Q AQ B Q B A c c --⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥=====⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦式中，变换矩阵为 310001100Q -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦系统的旺纳姆能控标准形为011004421000101x x u ⎡⎤⎡⎤⎢⎥⎢⎥=-+⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦龙伯格能控标准形相对于旺纳姆能控标准形在系统极点配置中应用更广，这是由于在构成状态变换阵时，相对均衡地考虑了各个输入的作用。

考虑线性定常系统:Σx Ax Bu y Cx=+=x 为n 维状态向量，u 为p 输入向量，y 为q 维输出向量。

如果系统完全能控，设矩阵B 的列线性无关，这在实际上通常总能成立，或可进行处理获得，b i 为B 的第i 列，那么，系统的能控性矩阵c Q 可表示为111121212[]c Q b b b Ab Ab Ab A b A b A b n n n p p p ---=对上式从左到右依次找出n 个线性无关列，即若某个列不能表示成其左边各线性独立列的线性组合就为线性无关列，否则为线性相关列。

考虑到c Q 的组成特点，如果,1,2,,A b i m m p = 和其左边线性独立列线性相关，则1Ab i m+也和其左边线性独立列线性相关，又由于假设B 为满秩，则可将此搜索结果重新排列，并组成非奇异矩阵P ：121111111222P b Ab A b b Ab A b b Ab A b pp p p μμμ----⎡⎤=⎣⎦式中，n p =+++μμμ 21， },,{21p μμμ 为系统能控性指数集。

对定义1P -矩阵求逆，并表示成分块矩阵形式：1111111()e e P P e e pTT T p T p μμ--⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦式中，块的行数为,1,2,,i i p μ= 。

再在矩阵P 中，取出各个块矩阵的末行，即为1212,,,e e e pT T T p μμμ ，并按如下方式构造变换矩阵P cL：111111111e e A e A P e e A e A p pppTTT cLT p Tp Tp μμμμμμμμ--⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦定理[龙伯格能控标准形]对完全能控的线性定常系统，存在线性非奇异变换1ˆxP x Q x cL cL -== 使状态空间表达式转化为龙伯格能控标准形ˆˆˆˆ:ˆˆΣx A x B u y Cx cL c cc=+= 1112112122212ˆˆˆˆˆˆˆˆˆˆA A A A A A A Q AQ A AA pp c cL cL p p pp -⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎣⎦(01001ˆ1,2,,0001***A ,i i ii i p μμ⨯⎡⎤⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦）()00ˆ,,,1,2,,000***A i j iji j i j p μμ⨯⎡⎤⎢⎥⎢⎥=≠=⎢⎥⎢⎥⎣⎦1001**ˆ001B Q B c cL-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ˆC CQ c cL =（无特殊形式）*表示可能的零元或非零元。

证明：从略。

例 求如下系统的龙伯格能控标准形1000303110114101101000A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥==⎢⎥⎢⎥--⎢⎥⎢⎥-⎣⎦⎣⎦计算系统的能控性矩阵c Q223312121212000114413100331010300114413135100011339c Q b b Ab Ab A b A b A b A b =⎡⎤⎣⎦⎡⎤⎢⎥-----⎢⎥=⎢⎥⎢⎥-----⎣⎦ 式中，21212,,,b b Ab A b 线性无关，所以，系统完全能控，且11μ=，23μ=。

由此构成矩阵1P -，并求其逆阵P 。

121212001410310014130013P b b Ab A b -⎡⎤⎢⎥--⎢⎥==⎡⎤⎣⎦⎢⎥⎢⎥--⎣⎦1102101330041001P -⎡⎤⎢⎥-⎢⎥=⎢⎥--⎢⎥⎣⎦将P 分为两个块阵，第一个块的行数为1，第二个块的行数为3。

取出P 中的第一行11e T 和第四行23e T，组成变换矩阵P cL ，并求其逆阵Q cL 。

1123232231102100110000010e eP e A e A T T cL T T -⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦10010123000010110Q P cL cL -⎡⎤⎢⎥-⎢⎥==⎢⎥⎢⎥-⎣⎦由此可求出系统的龙伯格能控标准形为1101010001000ˆˆ000100113401A Q AQ B Q B=c cL cLc cL--⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦2.多输入多输出系统的能观标准形利用能观性和能控性的对偶关系，可以方便地得出完全能观系统的旺纳姆能观标准形和龙伯格能观标准形。

定理 [旺纳姆能观标准形]对完全能观的线性定常系统，存在线性非奇异变换1xP x Q x o o -== 使状态空间表达式转化为旺纳姆能控标准形:Σx A x B u y Cx oW o oo=+= 式中111212212A A A A P AP A AA o o ol l ll -⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎣⎦1101001,1,2,,0001A i ii i i lγαα-⎡⎤⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦1200,1,,,1,2,,00A i ij ij ij iji j l j l νςςς⎡⎤⎢⎥⎢⎥==+=⎢⎥⎢⎥⎢⎥⎣⎦1B P B o o-= （无特殊形式） (1)1(1)1010101C CP o o l l n p pn ββββ++⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦定理 [龙伯格能观标准形]对完全能观的线性定常系统，存在线性非奇异变换1x P x Q x o o -==使状态空间表达式转化为龙伯格能控标准形:Σx A x B u y C xoL o o o =+=式中1112121221112A A A A A A A P AP A A A q q o oL oL q q qq -⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎣⎦0*1001,1,2,,00001*A ii i q⎡⎤⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦00*00*,00*A ij i j ⎡⎤⎢⎥⎢⎥=≠⎢⎥⎢⎥⎣⎦1B P B o oL -= （无特殊形式）01**01C CP o oL ⎡⎤⎢⎥==⎢⎥⎢⎥⎣⎦能观标准形在系统观测器设计中有重要应用。

应该指出，单变量系统的传递函数与能控标准形和能观标准形的参数有明显的关系。

已知传递函数可直接写出系统的能控标准形和能观标准形，不必求出变换矩阵，反之，由标准形也可以直接写出传递函数。

但在多变量系统中，不存在这样的对应关系，必须求出变换矩阵，才能求出标准形。

4.7 线性系统的结构分解介绍如何通过非奇异变换，将系统的状态空间描述按能控性和能观性进行结构分解。

通过状态空间描述的结构分解可以建立状态空间描述和系统传递函数之间的关系，为最小实现问题的提出提供理论依据。

状态空间描述的结构分解在系统综合状态反馈和系统镇定问题中也有重要作用。

一、 按能控性的系统结构分解按能控性分解是通过线性非奇异变换，将系统的完全能控状态和不能控状态，以明显的方式区分开来表示。

考虑线性定常系统:Σx Ax Bu y Cx=+=式中，x 为n 维状态向量，u 为p 输入向量，y 为q 维输出向量。