答案及评分标准(现代控制理论06试卷B)
一、已知系统的微分方程如下,写出其状态空间表达式(10分)yy4y5y3u
解:(1)选择状态变量yx1,yx2,yx3(3分)
x1x2
(2)写出微分方程组(4分) x2x3
x5x
314x2x33u
yx1
(3)写出状态空间表达式(3分)
x1010x10
x
2001x
20u
x3541x33
x
y[100]x
2
x
3
二、已知系统传递函数,求出系统的约旦标准型的实现
W(s)6(s1)
s(s2)(s3)2
解:(1)把原传递函数展开成部分分式形式(9分)W(s)6(s1)410/331/3
s(s2)(s3)2(s3)2s3s2s
注:分子或者分母的数字算错一个扣1分,全对给9分。
(2)直接写出其并联型对角阵实现形式(6分):
31
30
1
xxu
2
01
1
y(410/331/3)x
1 (15分)
注:A阵元素必须和C阵元素一一对应,错一组扣1分;B阵元素错一个扣1分,扣完为止。
三、求下列状态空间表达式的解:(15分)
010
xxu 231
初始状态x(0),输入u(t)1(t)是单位阶跃函数。 解:(1)计算(t)(10分)
用课文中给出的三种方法求均可。此处仅给出其中一种方法,如果考生采用其它方法,请阅卷老师酌情扣分扣分。
求得特征值,1122(2分)
0
0
0(t)求(4分) 1(t)
0(t)11et11et2ete2t
2tt2t (t)1t
12e12eee
12
11
求(t)(4分)
2ete2t(t)e0(t)I1(t)At2t
2e2e
Atete2t
t2te2e
00
(2)计算x(t),把B,x(0), u(t)1(t)代入求解 (5分)
10
t
x(t)(t)x(0)(t)Bu()d
e(t)e2(t)
(t)d2(t)2e0e
1t12tee22t2t
ee
t
四、判断下列系统的能控性和能观测性并说明理由。(10分)
610004060000x00310x0
00030100002120101yx 30010
200u 10
解:(1) 能控性判断(5分)
61
A阵中约当块最后一行对应的B阵相应行为零向量,故系统不
06
完全能控。
(2) 能观性判断(5分)
6131
A阵中约当块及第一列对应的C阵相应列不为零向
0603
量,且C阵最后一列不为零向量,故系统完全能观。
五、设系统
(A,B,C,D)中:(15分)
0100
0101A000,B1,C,D 10120030
(1)试判断系统的能控性、能观测性;
(2)若不能控或不能观测,试给出其能控性分解或能观测性分解的描述形式。
解:(1) 判断系统的能控性、能观测性:
01001010
2① AB00010,AB00000300003
00
10 00
MB
AB
010
100
A2B
000
∵rankM=2<3∴该系统不完全能控。(4分)
010
010000② CA 000101013003
0100100002
CA000 101003009
01
C
0NCA
0
2CA0
0
100100
010 309
∵rankN=3∴该系统完全能观。(4分)
(2) ∵该系统不完全能控∴作能控性分解。
010
取M中的两个线性无关列10,再任取一个与其线性无关的列0,001010010
1
组成变换阵T100,求得T100T。(3分)
001001
从而求得能控性分解的系数阵:
0001
ˆT1AT100;BˆCT100 ˆT1B0;CA0110030
其中二维能控子系统为:
ˆ1x
001
ˆxu 1
100
101
ˆ1uy(4分) x
012
六、设二阶线性定常系统的状态方程为:(15分)
11
xx
23
试用李雅普诺夫第二法判断该系统原点的稳定性。
0
解:由x0知xe 是系统的一个平衡态。(2分)
0
p11
设P
p21
p12
,QI,其中P为实对称阵 p22
T
代入李雅普诺夫方程AP+PA=-I
p12p11
p22p21
p121110
(4分) p222301
11p11
得:p2321
可得联立方程组为:
T
2p114p121
p114p122p220 2p6p1
1222
74
解得P
58
5
8
(3分)
38
817
=﹥0(3分) 3648
根据希尔维斯特判据知:
7
741﹥0,2
54
8
所以P是正定的,系统的平衡点是大范围渐近稳定的。且李雅普诺夫函
数为:
12
(3分) V(x)=xTPx=(14x1210x1x23x2)
8
七、已知受控系统状态方程为:(20分)
1000
011x0u ,设计状态反馈阵,将极点配置在-3,-3,-3处。x 0161
解:(1)(5分)判断系统的能控性:
0
Qc,A,A20
1
1
17, 635
∵Qc10,∴Qc满秩,系统完全能控,所以通过状态反馈可任意配置极点。
(2)(5分)求出加入状态反馈阵k0
T
k1k2后闭环特征多项式:
k1
k2
100000T0110kf()I(A)0000161003(7k2)2(7k1k2)k0 (3)(5分)确定希望的闭环特征多项式:f*()(3)33922727 (4)(5分)计算其状态反馈阵:
令 f()f*() ,就可得到: kk0
T
k1k227182