1
2
0
2
1
1
1 0 2 3 1 1 0 0
分别记C阵中的第1,3,4,5,8,列 c111 ,c112 ,c113 ;c121 ,c122 。 这里c下标中的第一个数字为1,表示每个约当块的第1列,
第二个数字表示是第几个特征值,第三个数字表示第几 个约当块。按照如上规定可知c121表示第二个特征值的第
0 0 0 1 0 0 2 0 0 0
y 1 0 0 0 x
>> A=[0,1,0,0; 3,0,0,2; 0,0,0,1;0,-2,0,0]; >> B=[0; 1; 0; 0]; >> C=[1,0,0,0]; >> D=0; >> Uo=[C, C*A, C*A^2, C*A^3]; >> rank(Uo) ans =
x1 x2
2
0
0 x1
3
x2
y 1 0 x
x1 1
x2
x3 x4
x5
y 5 0
1 1
2 1 2
2 0 0 x
x1
x2
1
x3 x4
,
2 x5
利用Matlab判定系统可观性
0 1 0 0 0 x 3 0 0 2 x 1 u
一个约当块的第一列所对应的C中的列向量 0 0 1。T
1 2 0
c111 c112 c113 1 1 2
1 2 3
0 0
c121 c122 0 1
1 0
由于上述矩阵分别都是列满秩的,系统可观。
定理2-6 若线性定常连续系统 x Ax Bu ,系统可观的充分必要条件为下列矩阵的列复数 域上线性无关,即