现代控制理论课后习题答案
- 格式:pdf
- 大小:194.00 KB
- 文档页数:25
现代控制理论课后习题答案
第⼀章习题1.2求下列多项式矩阵()s D 和()s N 的两个不同的gcrd:
()2223(),()1232s s s s s s s s s ??++== ? ?
+-??D N 解:()()22
232321s s s s s s s
++ =++ ? ?
D S N S ; ()3r 2,1,2
E -:223381s s s s s s ??++ ?-- ? ???;()3r 2,3,3E :223051s s s s s ??++ ?
- ? ???;
()3r 1,3,2E s --:01051s s ?? ?- ? ;()3r 2,1,5E s -:01001s ?? ?
;()3r 3,1,1E -:01000s ?? ? ? ???;()1r 2,3E :01000s ?? ? ? ???;()1r 1,2E :00100s ?? ?
;
所以⼀个gcrd 为001s ??
;取任⼀单模矩阵预制相乘即可得另⼀个gcrd 。1.9 求转移矩阵t A e (1)已知1141??=
A ,根据拉⽒反变换求解转移矩阵t
A e 。
(2) 已知412102113-?? ?
= ? ?-??
A ,根据C-H 有限项展开法求解转移矩阵t A e 。
解:(1)11()41s s s --??-= ?--??
I A
1110.50.5
0.250.2511(3)(1)(3)(1)1
3131()4141110.50.5(3)(1)(3)(1)(3)(1)3131s s s s s s s s s s s s s s s s s s s s s s s --
+-
--+-+??-+-+ ? ?
-=== ? ?---+ ?
-+ ? ?-+-+-+-+?
I A 3311330.5e 0.5e 0.25e 0.25e e ()e e 0.5e 0.5e t t t t t t t
t t s ------??+-??=-= ??? ?-+?
A L I A (2)由241
2()1
2(1)(3)0113λλλλλλ--?? ?
=--=--= ? ?--??
A I -,得1,233,1λλ== 对1,23λ=,可以计算1,2()2rank λ=A I -,所以该特征值的⼏何重数为1。⽽31λ=的⼏
何重数显然是1。由此可写出其Jordan 形为310030001?? ?= ? ???
J 。
将特征值代⼊⽅程2012()()()e t t t t λααλαλ++=,并计及第⼀个根是⼏何重数为1的重
根,可得如下⽅程组330301233312
1012332
953()e e e 442()3()9()e 33()6()e ()e e 2e 22()()()e 11()e e e 442t t t t t
t t t t
t t t t t t t t t t t t t t t t t t ααααααααααα?
=-+??++=+=?=---++=???
=-+??
于是,2012e ()()()t t t t ααα?
=+=
A A A I +
1.14给出下列有理函数的Smith-Mcmillan 规范型
222
21(1)(2)(2)()12
(2)s s s s s s s s s ??+
+++ = - ++
G 解:()s G 的最⼩公分母为2()(1)(2)d s s s =++,于是()s G 可表⽰为
22211(1)()()()(1)(2)(1)(1)(2)s s s s d s s s s s s s ??
+== ?++-+++??
G N 对()s N 进⾏⾏列变换,导出Smith 规范型 3r 3r 1c 2c 2c 2r 23r 2r 2
(1,2,1):(2,1,(25)):(1,2):(2,3):(2,1,):(1,13):
(2,1,2(1)(21)):
10(2,12):()0(1)(20.5)s s s s s s s s s s +++-??
-= ?+++??
E E E E E E E E Λ
基此得到Smith-Mcmillan 规范型为22210
(1)(2)
()(20.5)0(2)s s s s s s s ?
++
=
++ ?
+
M ,显然它满⾜有
关互质性与整除性的属性。
第⼆章习题2.1求⽆源电路⽹络的模型。
下⾯是RC 电路,请根据原理画出其⽅框图,并求其传递函数和状态空间描述。R
i
i
u o
u C
图P2.1 ⼀阶RC ⽹络
解:利⽤基尔霍夫电压定律及电容元件特性可得i o
u u i R
-=,o idt u c =?
对其进⾏拉⽒变换得i o ()()()U s U s I s R -=
o ()
()I s U s sC
= 因此图P2.1即可转换为下图运算电路形式。
R
)
(s U i sC
1)
(s I )
(s U o
图SP2.1-1 ⼀阶RC
⽹络运算电路
由此分别得到下图a 和b ,将图a 和b 组合起来即得到图c ,图c 为该⼀阶RC ⽹络的⽅框
图。)(s U i )
(s U o -
+R
1
)
(s I (
I U sC
1-+
1(I 1 -+)(s I )(s U o sC 1-+(s I 1
-+)(s U i -
+R
1)
(s I )
(s U o sC
1)
(s U o
(a) (b) (c)图SP2.1-2 ⼀阶RC ⽹络的⽅框图
显然系统是⼀阶的,只有⼀个状态变量,令o x u =,i u u =则状态空间描述为 11x x uRC RC y x =-+=
2.2求⽆源电路⽹络的模型。
已知R -L -C ⽹络如图所⽰,试求以i u 为输⼊,o u 为输出的传递函数和状态空间描述。i
u o
u +
-
+
-
1
R 2
R 3
R L
C
1i 3
i 2
i
图P2.2 R-L-C ⽹络
解:电感⽅程:3o 11i d d i L
u R i u t ++=...(1) 电容⽅程:c 2d d u
C i t
= (2)
有6个变量,列出微分⽅程模型时保留2个,因此要消掉4个变量,还需要列出3个⽅程:
由KVL :11c 22i R i u R i u ++=...(3) 由KCL :123i i i =+...(4) 在输出端:33o R i u = (5)
将(5)代⼊(1)、(4)可消去3i ,然后将(4)代⼊(1)、(3)消去1i 得到: oo o 12i 33
o 12c 22i 3d ()(6)d ()(7)u u L u R i u R t R u R i u R i u R ?+++=
+++=
⽤⽅程(2)消去c u :将(7)整理为1o c 122i 3
()R u u R R i u R +++=后取时间的导数,再将(2)代⼊,得到:o 122i 123d d d ()d d d u R i i u
R R R t C t t
+++= o 2i 1323123d d d ()d d d u i u
R C
R i R R R C R C t t t
+++=...(8) 最后,将(6)整理为o 1323i 31o d [()]d u
R R i R u L R R u t
=-++,代⼊1(8)R ?得到
2o o 13i 31o 2o o i i
12331132d d {[()]}
d d d d d d (){[()]}d d d d u u
R C R u L R R u t t
u u u u R R C R L R R R R C
t t t t +-++++-++=
经整理,可得到系统的微分⽅程模型为2o 121323o
23312i o
i 2
13131313d d ()d d d d u R R C R R C R R C L u R R C R R R LC u u u R R R R t R R t R R t ++++++=+++++
所以其传递函数阵为233
1212212132312()()()1
()R R C R s R R LC R R LC
W s R R C R R C R R C L s s R R LC
+
++=++++++
这⾥写出能控标准I 型的状态状态空间描述11121323221213232121201011(),()()x x u R R C R R C R R C L x x R R LC x R R R C y x R R LC R R LC ??
?=++++ ? ? ? ?-- ?+= ? ?++??
2.3求有源电路⽹络的模型。
已知有源⽹络如下图,求传递函数与状态空间模型。-+
i 1
i 2KR 2
R 2
R 1R 1C 1
C 2
u i
u o
图P2.3 有源电路⽹络
解:1122o ,,c c x u x u y u ===令i 11
1121111111122222222221111211121222o 12
2222222222111210220()22u x x C x R C R C x x R R R C R u s R x x R x x R C R W s R R C R R x u x x R C R C C x y u RR x R R C --=+ ? ? ? ??=++ ? ? ? ??- ? ?=-??=-????-??-=+?=- ? ? ??1111C s R C +
2.4求弹簧阻尼器系统的模型。
已知弹簧阻尼器系统(忽略质量)如下图,f 是阻尼系数数,x 、x B 、y 分别是A 、B 、C 三点的位置变量,以x 为输⼊,y 为输出,求传递函数。并思考:传递函数求出的阶次与系统中储能元件的个数⼀致吗?为什么?x
y
f B A
C
B
x 1
k 2
k
图P2.4 弹簧阻尼器系统
解:设B 点位移为B x ,根据B 、C 点⼒平衡关系列写⽅程 对于B 点1()
()B B d x y k x x f
dt
--= 对于C 点
2()B d x y f k y dt
-=
上⾯两个⽅程两边同时进⾏拉⽒变换(初始条件为0),有[][]1()()()()B B k X s X s sf X s Y s -=-
[]2()()()B sf X s Y s k Y s -=
解上述⽅程组,得11212()