线性系统的结构分解

将状态空间按照能控性和能观性进行结构分解,在理论上揭示了状态空间的本质特征,也为最小实现问题提供了理论依据。

对于线性定常系统:

{x˙=Ax+Buy=Cx

一、按能控性分解

能控性判别矩阵:

rankM=rank(BABAn1B)=n1<n

则存在非奇异变换 Rc:前 n1列向量取能控性判别矩阵的 n1线性无关的列,另外 nn1 个列可以任意选取(只要保证 Rc 线性无关即可,可以尽量选取的简单)

得到变换后的状态空间表达式:

x=Rcx^{x^˙=A^x^+B^uy=C^x^A^=Rc1ARc=(A^11A^120A^22)B^=Rc1B=(B^10)C^=CRc=(C^1C^2)x^=(x^1x2)u=(u1u2)

能控子空间x^˙1=A^11x^1+B^1u1

二、按能观性分解

能观性判别矩阵:

rankN=rank(CCACAn1)=n1<n

则存在非奇异变换 Ro其逆矩阵 Ro1 的前 n1行向量取能控性判别矩阵的 n1线性无关的行,另外 nn1 个行可以任意选取(只要保证 Ro1 线性无关即可,可以尽量选取的简单)

注意

能观性分解时,首先得到的是非奇异变换的逆矩阵!而不是非奇异变换本身!!!

得到变换后的状态空间表达式:

x=Rox~{x~˙=A~x~+B~uy=C~x~A~=Ro1ARo=(A~110A~21A~22)B~=Ro1B=(B~1B~2)C~=CRo=(C~10)x~=(x~1x~2)

能观子空间

x^~1=A~11x~1+B~1u1y=C~1x~1

三、按能控能观性分解

如果系统不完全能观,不完全能控,则可以分解为四个子空间(实际上是先进行能控性分解,再进行能观性分解)

A=R1AR=(A110A130A21A22A23A2400A33000A34A44)B=R1B=(B1B200)C=CR=(C10C30)x=(xcoxco¯xc¯oxc¯o¯)