《现代控制理论》课后习题全部答案(最打印版)
(完整word版)现代控制理论习题解答(第三章)
第三章 线性控制系统的能控性和能观性01010（ 1） A10 1B( 2) A 0 0 1 ,B 011024311113 10 1 1（ 3） A0 10 1 0 3 0 , B00 ( 4) AB0 0 11 001211【解】：(1)11U c B AB 1 1, rankU c n 2 ，所以系统完全能控。
c 0 1 c(2)10 0 1 2U c B AB A 2B1 1 11 1 17前三列已经可使 rankU c n 3 ，所以系统完全能控（后续列元素不必计算） 。
（3）A 为约旦标准型， 且第一个约旦块对应的B 阵最后一行元素全为零， 所以系统不完全 能控。
（4）A 阵为约旦标准型的特殊结构特征， 所以不能用常规标准型的判别方法判系统的能控 性。
同一特征值对应着多个约旦块，只要是单输入系统，一定是不完全能控的。
可以求一下能控判别阵。
1213 1223B AB A 2B A 3B2 3 U c1 1 12 13 1 11 12 31111rankU c 2 ，所以系统不完全能控。
3 1110 10 0 x0 3 0x 0 0ux0 01x 0u (1)0 0 12(2)61161101yxy10 0x1 10解】：1）311 已知 A 0 30,B0 001220 0 D CB CAB CA 2B 0 0 前两列已经使 rank D CBCAB110 1 0 00 ， C ，D1 1 0 0 031112CA B m2， 所以系统输出能控。
（2） 系统为能控标准型，所以状态完全能控。
又因输出矩阵 状态维数 n ，所以状态能控则输出必然能控。
C 满秩，且输出维数 m 小于1 0x0 01xx1 1 (1)2 43 ； (2) 1 x 0；011y1 1xyx12 12 1 0 4 0 0x0 20xx4 0x(3)；（4）0 030 1y0 1 1x y11 4x解】：1）已知 A01 00 242-3-3 判断下列系统的能观性。
现代控制理论课后习题答案Word版
绪论为了帮助大家在期末复习中能更全面地掌握书中知识点，并且在以后参加考研考博考试直到工作中，为大家提供一个理论参考依据，我们11级自动化二班的同学们在王整风教授的带领下合力编写了这本《现代控制理论习题集》（刘豹第三版），希望大家好好利用这本辅助工具。
根据老师要求，本次任务分组化，责任到个人。
我们班整体分为五大组，每组负责整理一章习题，每个人的任务由组长具体分配，一个人大概分1~2道题，每个人任务虽然不算多，但也给同学们提出了要求：1.写清题号，抄题，画图（用CAD或word画）。
2.题解详略得当，老师要求的步骤必须写上。
3.遇到一题多解，要尽量写出多种方法。
本习题集贯穿全书，为大家展示了控制理论的基础、性质和控制一个动态系统的四个基本步骤，即建模、系统辨识、信号处理、综合控制输入。
我们紧贴原课本，强调运用统一、联系的方法分析处理每一道题，将各章节的知识点都有机地整合在一起，力争做到了对控制理论概念阐述明确，给每道题的解析赋予了较强的物理概念及工程背景。
在课后题中出现的本章节重难点部分，我们加上了必要的文字和图例说明，让读者感觉每一题都思路清晰，简单明了，由于我们给习题配以多种解法，更有助于发散大家的思维，做到举一反三！这本书是由11级自动化二班《现代控制理论》授课老师王整风教授全程监管，魏琳琳同学负责分组和发布任务书，由五个小组组组长李卓钰、程俊辉、林玉松、王亚楠、张宝峰负责自己章节的初步审核，然后汇总到胡玉皓同学那里，并由他做最后的总审核工作，绪论是段培龙同学和付博同学共同编写的。
本书耗时两周，在同学的共同努力下完成，是二班大家庭里又一份智慧和努力的结晶，望大家能够合理使用，如发现错误请及时通知，欢迎大家的批评指正！2014年6月2日第一章 控制系统的状态空间表达式1-1 试求图1-27系统的模拟结构图，并建立其状态空间表达式图1-27系统方块结构图解：系统的模拟结构图如下：图1-30双输入--双输出系统模拟结构图系统的状态方程如下：u K K x K K x K K x X K x K x x x x J K x J x J K x J K x x J K x x x pp p p n p b1611166131534615141313322211+--=+-==++--===••••••令y s =)(θ，则1x y =所以，系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡••••••6543211654321111111126543210000010000000000000010010000000000010x x x x x x y uK K x x x x x x K K K K K K J K J J K J K J K x x x x x x p p pp n p b1-2有电路如图1-28所示。
现代控制理论课后题及答案
第2章 “控制系统的状态空间描述”习题解答2.1有电路如图P2.1所示，设输入为1u ，输出为2u ，试自选状态变量并列写出其状态空间表达式。
图P2.1解 此题可采样机理分析法，首先根据电路定律列写微分方程，再选择状态变量，求得相应的系统状态空间表达式。
也可以先由电路图求得系统传递函数，再由传递函数求得系统状态空间表达式。
这里采样机理分析法。
设1C 两端电压为1c u ，2C 两端的电压为2c u ，则212221c c c du u C R u u dt++= (1) 112121c c c du u duC C dt R dt+= (2) 选择状态变量为11c x u =，22c x u =，由式(1)和(2)得：1121121121212111c c c du R R C u u u dt R R C R C R C +=--+ 2121222222111c c c du u u u dt R C R C R C =--+ 状态空间表达式为：12111211212121212122222221111111R R C x x x u R R C R C R C x x x u R C R C R C y u u x +⎧=--+⎪⎪⎪=--+⎨⎪⎪==-⎪⎩即: 12121121211112222222211111R R C R C R R C R C x x u x x R C R C R C +⎡⎤⎡⎤-⎢⎥⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦--⎢⎥⎢⎥⎣⎦⎣⎦[]11210x y u x ⎡⎤=-+⎢⎥⎣⎦2.2 建立图P22所示系统的状态空间表达式。
1图P2.2解 这是一个物理系统，采用机理分析法求状态空间表达式会更为方便。
令()f t 为输入量，即u f =，1M ，2M 的位移量1y ，2y 为输出量， 选择状态变量1x =1y ，2x = 2y ，3x =1dy dt，24dyx dt =。
《现代控制理论》课后习题答案1.pdf
《现代控制理论》第一章习题解答1.1 线性定常系统和线性时变系统的区别何在？ 答：线性系统的状态空间模型为：xAx Bu y Cx Du=+=+线性定常系统和线性时变系统的区别在于：对于线性定常系统，上述状态空间模型中的系数矩阵A ，B ，C 和中的各分量均为常数，而对线性时变系统，其系数矩阵D A ，B ，C 和中有时变的元素。
线性定常系统在物理上代表结构和参数都不随时间变化的一类系统，而线性时变系统的参数则随时间的变化而变化。
D 1.2 现代控制理论中的状态空间模型与经典控制理论中的传递函数有什么区别？ 答: 传递函数模型与状态空间模型的主要区别如下:传递函数模型（经典控制理论）状态空间模型（现代控制理论） 仅适用于线性定常系统 适用于线性、非线性和时变系统用于系统的外部描述 用于系统的内部描述基于频域分析基于时域分析1.3 线性系统的状态空间模型有哪几种标准形式？它们分别具有什么特点？答: 线性系统的状态空间模型标准形式有能控标准型、能观标准型和对角线标准型。
对于阶传递函数n 1212101110()n n n n n n n b s b s b s b G s d s a s a s a −−−−−−++++=+++++""，分别有[]012101210100000100000101n n n xx ua a a a yb b b b x du−−−⎧⎡⎤⎪⎢⎥⎪⎢⎥⎪⎢⎥=+⎪⎢⎥⎨⎢⎥⎪⎢⎥⎪−−−−⎣⎦⎪=+⎪⎩"" ###%##"""⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⑴ 能控标准型：[]0011221100010********001n n n b a b a xa x ub a b y xdu −−−⎧−⎡⎤⎡⎤⎪⎢⎥⎢⎥−⎪⎢⎥⎢⎥⎪⎢⎥⎢⎥=−+⎪⎢⎥⎢⎥⎨⎢⎥⎢⎥⎪⎢⎥⎢⎥⎪−⎣⎦⎣⎦⎪=+⎪⎩"" "######""⑵ 能观标准型：[]1212001001001n n p p x x up y c c c x du⎧⎡⎤⎡⎤⎪⎢⎥⎢⎥⎪⎢⎥⎢⎥=+⎪⎢⎥⎢⎥⎨⎢⎥⎢⎥⎪⎣⎦⎣⎦⎪⎪=+⎩"" ##%##""⑶ 对角线标准型： 式中的和可由下式给出，12,,,n p p p "12,,,n c c c "12121012111012()n n n n n n n n nb s b s b s bc c c G sd d s a s a s a s p s p s p −−−−−−++++=+=++++++−−−"""++能控标准型的特点：状态矩阵的最后一行由传递函数的分母多项式系数确定，其余部分具有特定结构，输出矩阵依赖于分子多项式系数，输入矩阵中的元素除了最后一个元素是1外，其余全为0。
现代控制理论课后习题答案
现代控制理论课后习题答案第⼀章习题1.2求下列多项式矩阵()s D 和()s N 的两个不同的gcrd:()2223(),()1232s s s s s s s s s ??++== ? ?+-??D N 解：()()22232321s s s s s s s++ =++ ? ?D S N S ； ()3r 2,1,2E -：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 ，根据拉⽒反变换求解转移矩阵tA e 。
(2) 已知412102113-?? ?= ? ?-??A ，根据C-H 有限项展开法求解转移矩阵t A e 。
解：(1)11()41s s s --??-= ?--??I A1110.50.50.250.2511(3)(1)(3)(1)13131()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 tt t s ------??+-??=-= ??? ?-+?A L I A (2)由2412()12(1)(3)0113λλλλλλ--?? ?=--=--= ? ?--??A I -，得1,233,1λλ== 对1,23λ=，可以计算1,2()2rank λ=A I -，所以该特征值的⼏何重数为1。
《现代控制理论》课后习题全部答案(最完整打印版)
第一章习题答案1-1 试求图1-27系统的模拟结构图，并建立其状态空间表达式。
11K s K K p +sK s K p 1+s J 11sK n 22s J K b -++-+-)(s θ)(s U 图1-27系统方块结构图解：系统的模拟结构图如下：)(s U )(s θ---+++图1-30双输入--双输出系统模拟结构图1K pK K 1pK K 1+++pK n K ⎰⎰⎰11J ⎰2J K b ⎰⎰-1x 2x 3x 4x 5x 6x系统的状态方程如下：u K K x K K x K K x X K x K x x x x J K x J x J K x J K x x J K x x x pp p p n p b1611166131534615141313322211+--=+-==++--===∙∙∙∙∙∙阿令y s =)(θ，则1x y =所以，系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∙∙∙∙∙∙654321165432111111112654321000001000000000000010010000000000010x x x x x x y uK K x x x x x x K K K K K K J K J J K J K J K x x x x x x p p pp npb1-2有电路如图1-28所示。
以电压)(t u 为输入量，求以电感中的电流和电容上的电压作为状态变量的状态方程，和以电阻2R 上的电压作为输出量的输出方程。
R1L1R2L2CU---------Uc---------i1i2图1-28 电路图解：由图，令32211,,x u x i x i c ===，输出量22x R y =有电路原理可知：∙∙∙+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=∙∙∙写成矢量矩阵形式为：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000*********x x x R y u L x x x CCL L R L L R x x x 。
(完整word版)《现代控制理论》第3版课后习题答案
《现代控制理论参考答案》第一章答案1-1 试求图1-27系统的模拟结构图，并建立其状态空间表达式。
图1-27系统方块结构图解：系统的模拟结构图如下：图1-30双输入--双输出系统模拟结构图系统的状态方程如下：u K K x K K x K K x X K x K x x x x J K x J x J K x J K x x J K x x x pp p p n p b1611166131534615141313322211+--=+-==++--===••••••令y s =)(θ，则1x y =所以，系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡••••••654321165432111111112654321000001000000000000010010000000000010x x x x x x y uK K x x x x x x K K K K K K J K J J K J K J K x x x x x x p p pp npb1-2有电路如图1-28所示。
以电压)(t u 为输入量，求以电感中的电流和电容上的电压作为状态变量的状态方程，和以电阻2R 上的电压作为输出量的输出方程。
U图1-28 电路图解：由图，令32211,,x u x i x i c ===，输出量22x R y =有电路原理可知：•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000*********x x x R y u L x x x CCL L R L L R x x x 。
《现代控制理论》课后习题答案
=
3 2
， c2
=
2s + 5 lim s→−3 s + 1
=
1 2
。
从输入通道直接到输出通道上的放大系数 d = 1，由此可得:
⎡ x1
⎢ ⎣
x 2
⎤ ⎥ ⎦
=
⎡− 1
⎢ ⎣
0
0⎤ − 3⎥⎦
⎡ ⎢ ⎣
x1 x2
⎤ ⎥ ⎦
+
⎡1⎤ ⎢⎣1⎥⎦u
y
=
⎡ ⎢⎣
3 2
1 2
⎤ ⎥⎦
⎡ ⎢ ⎣
x1 x2
u
d
d
b2
dt
dt
d
b1
m
dt
b0
因此，两个环节调换后的系统状态变量图为
u
d
d
b2
dt
dt
d
b1
dt
b0
m
−∫
−∫
y −∫
a0
a1
a2
进一步简化，可得系统状态变量图为 u
b0
b1
b2
− ∫ x1
− ∫ x2
− ∫ x3 y
a0
a1
a2
3
取 y = x3 ， y = x2 ， y = x1 ，可以得到两个环节调换后的系统的状态空间模型为
a(s)
1 a(s)
=
s3
+
1 a2s2 +
a1s
+
a0
， b(s)
=
b2 s 2
+ b1s
+ b0
。
2
由于 s−3 y 相当于对 y 作 3 次积分，故 y = 1 可用如下的状态变量图表示： m a(s)
现代控制理论课后题答案(第二章-第六章)
1
R R2C1 1 1 1 1 x1 x2 u1 x R1 R2C1 R2C1 R2C1 1 1 1 2 x1 x2 u1 x R2C2 R2C2 R2C2 y u2 u1 x1
即:
R1 R2C1 1 R1 R2C1 x x 2 1 R2C2
x2
u
R1
R3
y
R2
图 P2.8 RL 电网络
解 采用机理分析法求状态空间表达式。由电路原理可得到如下微分方程
2 x1 x2 R3 R2 x2 L2 x
1 x1 x2 R3 u x1 L1 x / R1
y x1 x2 R3
(2) 解 采用拉氏变换法求取状态空间表达式。对微分方程 (2)在零初试条件 下取拉氏变换得：
2s 3Y ( s ) 3sY ( s ) s 2U ( s ) U ( s ) 1 2 1 s Y (s) s 1 2 2 U ( s ) 2s 3 3s s 3 3 s 2
dy1 dy ， x4 2 。 dt dt
3 Kx1 B1 M1 x
2
d ( x2 x1 ) dt
对 M 2 有：
4 f (t ) B M2x
经整理得：
1
d ( x2 x 1) dx B 2 dt dt
2
状态方程为：
1 x3 x x 2 x4 B B K 3 x1 1 x3 1 x4 x M1 M1 M1 B B B 1 4 1 x3 ( 1 2 ) x4 u x M2 M2 M2 M2
1 1 R2C1 x1 R2C1 u1 x2 1 1 R C R2C2 2 2
现代控制理论第版课后习题答案
现代控制理论参考答案第一章答案1-1 试求图1-27系统的模拟结构图,并建立其状态空间表达式; 解：系统的模拟结构图如下： 系统的状态方程如下： 令y s =)(θ,则1x y =所以,系统的状态空间表达式及输出方程表达式为1-2有电路如图1-28所示;以电压)(t u 为输入量,求以电感中的电流和电容上的电压作为状态变量的状态方程,和以电阻2R 上的电压作为输出量的输出方程; 解：由图,令32211,,x u x i x i c ===,输出量22x R y =有电路原理可知：•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为：1-4 两输入1u ,2u ,两输出1y ,2y 的系统,其模拟结构图如图1-30所示,试求其状态空间表达式和传递函数阵;解：系统的状态空间表达式如下所示： 1-5系统的动态特性由下列微分方程描述列写其相应的状态空间表达式,并画出相应的模拟结构图; 解：令..3.21y x y x y x ===，，,则有 相应的模拟结构图如下： 1-6 2已知系统传递函数2)3)(2()1(6)(+++=s s s s s W ,试求出系统的约旦标准型的实现,并画出相应的模拟结构图解：ss s s s s s s s W 31233310)3(4)3)(2()1(6)(22++++-++-=+++=1-7 给定下列状态空间表达式[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321321321100210311032010x x x y u x x x x x x ‘（1） 画出其模拟结构图（2） 求系统的传递函数 解：2⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+-=-=31103201)()(s s s A sI s W 1-8 求下列矩阵的特征矢量3⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=6712203010A 解：A 的特征方程 061166712230123=+++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+---=-λλλλλλλA I 解之得：3,2,1321-=-=-=λλλ当11-=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---3121113121116712203010p p p p p p 解得： 113121p p p -== 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P 或令111-=p ,得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P 当21-=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---32221232221226712203010p p p p p p解得： 1232122221,2p p p p =-= 令212=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1423222122p p p P或令112=p ,得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=21213222122p p p P 当31-=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---33231333231336712203010p p p p p p 解得： 133313233,3p p p p =-= 令113=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3313323133p p p P 1-9将下列状态空间表达式化成约旦标准型并联分解2⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡32121321321110021357213311201214x x x y y u x x x x x x解：A 的特征方程 0)3)(1(311212142=--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=-λλλλλλA I 当31=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--3121113121113311201214p p p p p p 解之得 113121p p p == 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P 当32=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--1113311201214312111312111p p p p p p解之得 32222212,1p p p p =+= 令112=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0013222122p p p P 当13=λ时,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--332313332313311201214p p p p p p 解之得3323132,0p p p == 令133=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1203323133p p p P约旦标准型1-10 已知两系统的传递函数分别为W 1s 和W 2s试求两子系统串联联结和并联连接时,系统的传递函数阵,并讨论所得结果 解：1串联联结 2并联联结1-11 第3版教材已知如图1-22所示的系统,其中子系统1、2的传递函数阵分别为 求系统的闭环传递函数 解：1-11第2版教材 已知如图1-22所示的系统,其中子系统1、2的传递函数阵分别为 求系统的闭环传递函数 解：1-12 已知差分方程为试将其用离散状态空间表达式表示,并使驱动函数u 的系数b 即控制列阵为 1⎥⎦⎤⎢⎣⎡=11b 解法1： 解法2：求T,使得⎥⎦⎤⎢⎣⎡=-111B T 得⎥⎦⎤⎢⎣⎡=-10111T 所以 ⎥⎦⎤⎢⎣⎡-=1011T 所以,状态空间表达式为第二章习题答案2-4 用三种方法计算以下矩阵指数函数At e ;2 A=1141⎛⎫⎪⎝⎭解：第一种方法： 令0I A λ-=则11041λλ--=-- ,即()2140λ--=;求解得到13λ=,21λ=- 当13λ=时,特征矢量11121p p p ⎡⎤=⎢⎥⎣⎦由 111Ap p λ=,得11112121311341p p p p ⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦即112111112121343p p p p p p +=⎧⎨+=⎩,可令112p ⎡⎤=⎢⎥⎣⎦当21λ=-时,特征矢量12222p p p ⎡⎤=⎢⎥⎣⎦由222Ap p λ=,得121222221141p p p p -⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦即1222121222224p p p p p p +=-⎧⎨+=-⎩ ,可令212p ⎡⎤=⎢⎥-⎣⎦则1122T ⎡⎤=⎢⎥-⎣⎦,111241124T -⎡⎤⎢⎥=⎢⎥⎢⎥-⎢⎥⎣⎦第二种方法,即拉氏反变换法：第三种方法,即凯莱—哈密顿定理 由第一种方法可知13λ=,21λ=-2-5 下列矩阵是否满足状态转移矩阵的条件,如果满足,试求与之对应的A 阵;3()22222222t tt t ttt t e e e e t e e e e --------⎡⎤--Φ=⎢⎥--⎣⎦ 4()()()()3333112412t t t t t tt t e e e e t e e e e ----⎡⎤+-+⎢⎥Φ=⎢⎥⎢⎥-++⎢⎥⎣⎦解：3因为 ()10001I ⎡⎤Φ==⎢⎥⎣⎦,所以该矩阵满足状态转移矩阵的条件 4因为()10001I ⎡⎤Φ==⎢⎥⎣⎦,所以该矩阵满足状态转移矩阵的条件2-6 求下列状态空间表达式的解：初始状态()101x ⎡⎤=⎢⎥⎣⎦,输入()u t 时单位阶跃函数;解： 0100A ⎡⎤=⎢⎥⎣⎦ 因为 01B ⎡⎤=⎢⎥⎣⎦,()()u t I t =2-9 有系统如图所示,试求离散化的状态空间表达式;设采样周期分别为T=和1s,而1u 和2u 为分段常数; 图 系统结构图 解：将此图化成模拟结构图 列出状态方程则离散时间状态空间表达式为 由()At G T e =和()0TAt H T e dtB =⎰得：当T=1时 ()()()()11111001111k e e x k x k u k e ke ----⎡⎤-⎡⎤+=+⎢⎥⎢⎥--⎢⎥⎣⎦⎣⎦当T=时 ()()()()()0.10.10.10.11001110.90.1k e e x k x k u k e k e ----⎡⎤-⎡⎤⎢⎥+=+⎢⎥⎢⎥---⎣⎦⎣⎦第三章习题3-1判断下列系统的状态能控性和能观测性;系统中a,b,c,d 的取值对能控性和能观性是否有关,若有关,其取值条件如何 1系统如图所示： 解：由图可得： 状态空间表达式为：由于•2x 、•3x 、•4x 与u 无关,因而状态不能完全能控,为不能控系统;由于y 只与3x 有关,因而系统为不完全能观的,为不能观系统; 3系统如下式：解：如状态方程与输出方程所示,A 为约旦标准形;要使系统能控,控制矩阵b 中相对于约旦块的最后一行元素不能为0,故有0,0≠≠b a ;要使系统能观,则C 中对应于约旦块的第一列元素不全为0,故有0,0≠≠d c ; 3-2时不变系统试用两种方法判别其能控性和能观性; 解：方法一：方法二：将系统化为约旦标准形;⎥⎦⎤⎢⎣⎡=1-111T ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=21212121T 1- B T -1中有全为零的行,系统不可控;CT 中没有全为0的列,系统可观; 3-3确定使下列系统为状态完全能控和状态完全能观的待定常数i i βα和 解：构造能控阵：要使系统完全能控,则211αα≠+,即0121≠+-αα 构造能观阵：要使系统完全能观,则121αα-≠-,即0121≠+-αα 3-4设系统的传递函数是1当a 取何值时,系统将是不完全能控或不完全能观的 2当a 取上述值时,求使系统的完全能控的状态空间表达式; 3当a 取上述值时,求使系统的完全能观的状态空间表达式; 解：1 方法1 ：)6)(3)(1()()()(++++==s s s as s u s y s W 系统能控且能观的条件为Ws 没有零极点对消;因此当a=1,或a=3或a=6时,系统为不能控或不能观; 方法2：系统能控且能观的条件为矩阵C 不存在全为0的列;因此当a=1,或a=3或a=6时,系统为不能控或不能观;2当a=1, a=3或a=6时,系统可化为能控标准I 型3根据对偶原理,当a=1, a=2或a=4时,系统的能观标准II 型为 3-6已知系统的微分方程为：u y y y y 66116...=+++试写出其对偶系统的状态空间表达式及其传递函数; 解：63611603210=====b a a a a ，，，，系统的状态空间表达式为 传递函数为其对偶系统的状态空间表达式为： 传递函数为61166)(23+--=s s s s W 3-9已知系统的传递函数为 试求其能控标准型和能观标准型;解：345213486)(222++++=++++=s s s s s s s s W系统的能控标准I 型为 能观标准II 型为3-10给定下列状态空间方程,试判别其是否变换为能控和能观标准型;解：[]100210311032010=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=C b A ，， 3-11试将下列系统按能控性进行分解1[]111,100,340010121-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=C b A 解：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--==9310004102b A Ab bM rankM=2<3,系统不是完全能控的; 构造奇异变换阵c R ：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==010*********R Ab R b R ，，,其中3R 是任意的,只要满足c R 满秩;即⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=031100010c R 得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=-010*******c R 3-12 试将下列系统按能观性进行结构分解1 []111,100,340010121-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=C b A 解： 由已知得[]111,100,340010121-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=C b A 则有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=4742321112CA CA C Nrank N=2<3,该系统不能观构造非奇异变换矩阵10R -,有10111232001R --⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦ 则0311210001R --⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦3-13 试将下列系统按能控性和能观性进行结构分解1[]211,221,102322001=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=C b A 解：由已知得211121226202M A Ab Ab ⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥-⎣⎦rank M=3,则系统能控 rank N=3,则系统能观所以此系统为能控并且能观系统取211121226202c T ⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦,则1217344173215344c T -⎡⎤-⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦则002105014A ⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦,12100c B T b -⎡⎤⎢⎥==⎢⎥⎢⎥⎣⎦,[]271323c c cT == 3-14求下列传递函数阵的最小实现; 1 ()111111w s s ⎡⎤=⎢⎥+⎣⎦解： 01α=,01111B ⎡⎤=⎢⎥⎣⎦,1001c A -⎡⎤=⎢⎥-⎣⎦ 1001c B ⎡⎤=⎢⎥⎣⎦,1111c C ⎡⎤=⎢⎥⎣⎦,0000c D ⎡⎤=⎢⎥⎣⎦ 系统能控不能观取101101R -⎡⎤=⎢⎥⎣⎦,则01101R -⎡⎤=⎢⎥⎣⎦所以10010ˆ01A R AR --⎡⎤==⎢⎥-⎣⎦,1011ˆ01c B R B -⎡⎤==⎢⎥⎣⎦ 010ˆ10c C C R ⎡⎤==⎢⎥⎣⎦,00ˆ00D ⎡⎤=⎢⎥⎣⎦所以最小实现为ˆ1m A =,[]ˆ11m B =,1ˆ1m C ⎡⎤=⎢⎥⎣⎦,00ˆ00m D ⎡⎤=⎢⎥⎣⎦ 验证：()()1111ˆˆˆ111m mm C sI A B w s s -⎡⎤-==⎢⎥+⎣⎦3-15设1∑和2∑是两个能控且能观的系统1试分析由1∑和2∑所组成的串联系统的能控性和能观性,并写出其传递函数； 2试分析由1∑和2∑所组成的并联系统的能控性和能观性,并写出其传递函数; 解： 11∑和2∑串联当1∑的输出1y 是2∑的输入2u 时,331222x x x x =-++010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=--+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]001y x = 则rank M=2<3,所以系统不完全能控; 当2∑得输出2y 是1∑的输入1u 时011034100021x x u ⎡⎤⎡⎤⎢⎥⎢⎥=--+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]210y x = 因为 2001016124M bAbA b ⎡⎤⎢⎥⎡⎤==-⎣⎦⎢⎥⎢⎥--⎣⎦rank M=3 则系统能控因为2210321654c N cA cA ⎡⎤⎡⎤⎢⎥⎢⎥==--⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦rank N=2<3 则系统不能观 21∑和2∑并联010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=--+⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦,[]211y x = 因为rank M=3,所以系统完全能控 因为rank N=3,所以系统完全能观现代控制理论第四章习题答案4-1判断下列二次型函数的符号性质：1222123122313()31122Q x x x x x x x x x x =---+-- 2222123122313()4262v x x x x x x x x x x =++--- 解：1由已知得110∆=-<,2112013-∆==>-,31111711302411112--∆=--=-<--- 因此()Q x 是负定的 2由已知得110∆=>,2113014-∆==>-,3111143160131--∆=--=-<--因此()Q x 不是正定的 4-2已知二阶系统的状态方程：试确定系统在平衡状态处大范围渐进稳定的条件;解：方法1：要使系统在平衡状态处大范围渐进稳定,则要求满足A 的特征值均具有负实部;即：有解,且解具有负实部; 即：1122112212210a a a a a a +<>且方法2：系统的原点平衡状态0e x =为大范围渐近稳定,等价于T A P PA Q +=-;取Q I =,令11121222PP P P P ⎡⎤=⎢⎥⎣⎦,则带入TA P PA Q +=-,得到 若 11211211222111221122122112222204()()0022a a a a a a a a a a a a a a +=+-≠,则此方程组有唯一解;即 其中11221221det A A a a a a ==- 要求P 正定,则要求 因此11220a a +<,且det 0A >4-3试用lyapunov 第二法确定下列系统原点的稳定性;11123x x -⎡⎤=⎢⎥-⎣⎦ 21111x x -⎡⎤=⎢⎥--⎣⎦解：1系统唯一的平衡状态是0e x =;选取Lyapunov 函数为2212()0V x x x =+>,则 ()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;2系统唯一的平衡状态是0e x =;选取Lyapunov 函数为2212()0V x x x =+>,则 ()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;4-6设非线性系统状态方程为： 试确定平衡状态的稳定性;解：若采用克拉索夫斯基法,则依题意有： 取P I =很明显,()Q x 的符号无法确定,故改用李雅普诺夫第二法;选取Lyapunov 函数为2212()0V x x x =+>,则()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;4-9设非线性方程：试用克拉索夫斯基法确定系统原点的稳定性; 解：1采用克拉索夫斯基法,依题意有：x →∞,有()V x →∞; 取P I =则2121013()132x Q x x ⎡⎤-+=⎢⎥-+⎣⎦,根据希尔维斯特判据,有： 2221121210310310132x x x -∆=∆==->-+，（）,()Q x 的符号无法判断; 2李雅普诺夫方法：选取Lyapunov 函数为421233()042V x x x =+>,则 ()V x •是负定的;x →∞,有()V x →∞;即系统在原点处大范围渐近稳定;4-12试用变量梯度法构造下列系统的李雅普诺夫函数 解：假设()V x 的梯度为： 计算()V x 的导数为：选择参数,试选112212211,0a a a a ====,于是得：12x V x ⎛⎫∇= ⎪⎝⎭,显然满足旋度方程12122121,0V V x xx x x x ∂∇∂∇∂∂===∂∂∂∂即,表明上述选择的参数是允许的;则有：如果121211202x x x x -><或,则()V x •是负定的,因此,1212x x <是12x x 和的约束条件; 计算得到()V x 为：()V x 是正定的,因此在121211202x x x x -><即范围内,0e x =是渐进稳定的;现代控制理论第五章习题答案5-1已知系统状态方程为：试设计一状态反馈阵使闭环系统极点配置为-1,-2,-3; 解：依题意有：2011012112M bAbA b ⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥⎣⎦3rankM =,系统能控; 系统0(,,)A b C =∑的特征多项式为：则将系统写成能控标准I 型,则有010*********x x u ⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦; 引入状态反馈后,系统的状态方程为：()x A bK x bu =++,其中3K ⨯为1矩阵,设[]012K k k k =,则系统(,,)K A bK C =∑的特征多项式为：根据给定的极点值,得到期望特征多项式为：比较*()()f f λλ与各对应项系数,可解得：012599k k k =-=-=-，则有：[]-5-9-9K =;5-3有系统：（1） 画出模拟结构图;（2） 若动态性能不满足要求,可否任意配置极点 （3） 若指定极点为-3,-3,求状态反馈阵; 解1系统模拟结构图如下：2系统采用状态反馈任意配置极点的充要条件是系统0(,,)A b C =∑完全能控; 对于系统0(,,)A b C =∑有： []0111M bAb ⎡⎤==⎢⎥-⎣⎦2rankM =,系统能控,故若系统动态性能不满足要求,可任意配置极点;3系统0(,,)A b C =∑的特征多项式为：则将系统写成能控标准I 型,则有010231x x u ⎡⎤⎡⎤=+⎢⎥⎢⎥--⎣⎦⎣⎦; 引入状态反馈后,系统的状态方程为：()x A bK x bu =++,设[]01K k k =,则系统(,,)KA bK C =∑的特征多项式为：根据给定的极点值,得到期望特征多项式为：比较*()()f f λλ与各对应项系数,可解得：[]017373k k K =-=-=--，; 5-4设系统传递函数为试问能否利用状态反馈将传递函数变成 若有可能,试求出状态反馈K ,并画出系统结构图;解：6522)3)(2)(1()2)(1()(232--+-+=+-++-=s s s s s s s s s s s W由于传递函数无零极点对消,因此系统为能控且能观; 能控标准I 型为 令[] 210k k k K =为状态反馈阵,则闭环系统的特征多项式为由于状态反馈不改变系统的零点,根据题意,配置极点应为-2,-2,-3,得期望特征多项式为比较 )(λf 与 )(*λf 的对应项系数,可得 即[]52118---=K 系统结构图如下：5-5使判断下列系统通过状态反馈能否镇定;11222 A 011,01011b ---⎡⎤⎡⎤⎢⎥⎢⎥=-=⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦解：系统的能控阵为：2240010115M bAbA b -⎡⎤⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥-⎣⎦3rankM =,系统能控; 由定理 5.2.1可知,采用状态反馈对系统0(,,)A b C =∑任意配置极点的充要条件是(,,)A b C =∑完全能控;又由于3rankM =,系统0(,,)A b C =∑能控,可以采用状态反馈将系统的极点配置在根平面的左侧,使闭环系统镇定; 5-7设计一个前馈补偿器,使系统 解耦,且解耦后的极点为1,1,2,2----; 解：0()()() d W s W s W s = 5-10已知系统：试设计一个状态观测器,使观测器的极点为-r,-2rr>0;解：因为1001c N cA ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎣⎦满秩,系统能观,可构造观测器; 系统特征多项式为[]21det det 0I A λλλλ-⎡⎤-==⎢⎥⎣⎦,所以有10010,0,10a a L ⎡⎤===⎢⎥⎣⎦ 于是11001100x T ATx T bu x u --⎡⎤⎡⎤=+=+⎢⎥⎢⎥⎣⎦⎣⎦ 引入反馈阵12g G g ⎡⎤=⎢⎥⎣⎦,使得观测器特征多项式：根据期望极点得期望特征式：比较()f λ与()*f λ各项系数得：即223r G r ⎡⎤=⎢⎥⎣⎦,反变换到x 状态下2201321023r r G TG r r ⎡⎤⎡⎤⎡⎤===⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ 观测器方程为：。
现代控制理论课后作业答案程鹏王艳东
现代控制理论第一次作业1-1.由图1-1所示，可得：1311322323313112121()331()122x u x s x u x x x u x x x u x s x x x x y x x u s y x x u⎧=-⎪+=--⎧⎪⎪⎪=--=-⎪⎪⇒+⎨⎨=⎪⎪=⎪⎪=++⎩⎪⎪=++⎩ 则状态空间可表示为：()301101112000110x x uy x u--⎛⎫⎛⎫⎪ ⎪=--+ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭=+ 1-4.由101,111A B ⎛⎫⎛⎫==⎪ ⎪⎝⎭⎝⎭。
11210()()1110111(1)1s s sI A s s s s ---⎛⎫Φ=-= ⎪--⎝⎭⎛⎫ ⎪-⎪= ⎪ ⎪--⎝⎭则，110[()]t Attt e e L sI A tee --⎛⎫=-= ⎪⎝⎭， ()010()()1()t ttA t t t e eBu d u d t e e τττττττττ----⎛⎫⎛⎫=⎪ ⎪-⎝⎭⎝⎭⎰⎰，()1u τ= 则，()0()(0)()1010212tAtA t t t t t t t t x t e x e Bu d e e te e te e te τττ-=+⎛⎫⎛⎫-⎛⎫=+ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎛⎫-= ⎪⎝⎭⎰1-5.（1）极点多项式为： 由()2rank G s =， 一阶子式公分母：2(1)s s + 二阶子式公分母：22(1)s s + 极点多项式为：22(1)s s + （2）零点多项式为：二阶子式：2222212(1)()212(1)(1)s s s s s s s s --+-++=++ 零点多项式为：1(1)()2s s -+现代控制理论第二次作业1-7.系统的状态方程为：x Ax bu =+其中，01101001n A a a a -⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥---⎣⎦，001b ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦。
1101111101111011000()011**1001**111n n n n n n n n n s sI A b s s s s s s s s s s s ααααααααα----------⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥-⎢⎥⎢⎥+⎣⎦⎣⎦⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥++++⎢⎥⎢⎥⎣⎦⎣⎦⎡⎤⎢⎥⎢⎥=⎢⎥++++⎢⎥⎣⎦1-8.反证法：设1121[]n rank b AbA b n n -=<，则存在一个非零向量α使得： 11111011[]00n n n bAbA b b Ab A b αααα---=+++=不防设110n α-≠,则，11111201211()n n n n A b b Ab A b αααα----=-+++两边同乘A ，则11111201211()n n n n A b Ab A b A b αααα---=-+++则可看出1n A b 能用12(,,)n b Ab A b -线性表出,以此类推，可得11+1n n n A b A b A b （，）均可由12(,,)n b Ab A b -线性表出，则：121[]n rank b AbA b n n -=≠ 与已知矛盾，假设不成立，所以有111[]n rank b AbA b n -=1-9.（1）解：010110001A ⎡⎤⎢⎥=--⎢⎥⎢⎥⎣⎦,011b ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，[]121C =可控性矩阵2011110111U bAbA b -⎡⎤⎢⎥⎡⎤==-⎣⎦⎢⎥⎢⎥⎣⎦，det 0U ≠,故系统可控。
《现代控制理论》刘豹著(第版)课后习题答案
《现代控制理论》刘豹著（第3版）课后习题答案《现代控制理论》刘豹著(第3版)课后习题答案第一章习题答案1-1试求图1-27系统的模拟结构图，并建立其状态空间表达式。
解：系统的模拟结构图如下：系统的状态方程如下：令，则所以，系统的状态空间表达式及输出方程表达式为1-2有电路如图1-28所示。
以电压为输入量，求以电感中的电流和电容上的电压作为状态变量的状态方程，和以电阻上的电压作为输出量的输出方程。
解：由图，令，输出量有电路原理可知：既得写成矢量矩阵形式为：1-3参考例子1-3（P19）.1-4两输入，，两输出，的系统，其模拟结构图如图1-30所示，试求其状态空间表达式和传递函数阵。
解：系统的状态空间表达式如下所示：1-5系统的动态特性由下列微分方程描述列写其相应的状态空间表达式，并画出相应的模拟结构图。
解：令，则有相应的模拟结构图如下：1-6（2）已知系统传递函数,试求出系统的约旦标准型的实现，并画出相应的模拟结构图解：1-7给定下列状态空间表达式‘（1）画出其模拟结构图（2）求系统的传递函数解：（2）1-8求下列矩阵的特征矢量（3）解：A的特征方程解之得：当时，解得：令得（或令，得）当时，解得：令得（或令，得）当时，解得：令得1-9将下列状态空间表达式化成约旦标准型（并联分解）（2）解：A的特征方程当时，解之得令得当时，解之得令得当时，解之得令得约旦标准型1-10已知两系统的传递函数分别为W1(s)和W2(s)试求两子系统串联联结和并联连接时，系统的传递函数阵，并讨论所得结果解：（1）串联联结（2）并联联结1-11（第3版教材）已知如图1-22所示的系统，其中子系统1、2的传递函数阵分别为求系统的闭环传递函数解：1-11（第2版教材）已知如图1-22所示的系统，其中子系统1、2的传递函数阵分别为求系统的闭环传递函数解：1-12已知差分方程为试将其用离散状态空间表达式表示，并使驱动函数u的系数b(即控制列阵)为（1）解法1：解法2：求T,使得得所以所以，状态空间表达式为第二章习题答案2-4用三种方法计算以下矩阵指数函数。
现代控制理论课后习题答案
前言本书是为了与张嗣瀛院士等编写的教材《现代控制理论》相配套而编写的习题解答。
本书对该教材中的习题给予了详细解答，可帮助同学学习和理解教材的内容。
由于习题数量较多，难易程度不同，虽然主要对象是研究型大学自动化专业本科学生，但同时也可以作使用其它教材的专科、本科、以及研究生的学习参考书。
书中第5、6、8章习题由高立群教授组织编选和解答；第4、7 章由井元伟教授组织编选和解答，第1、2章由郑艳副教授组织编选和解答。
由于时间比较仓促，可能存在错误，请读者批评、指正。
另外有些题目解法和答案并不唯一，这里一般只给出一种解法和答案。
编者 2005年5月第2章 “控制系统的状态空间描述”习题解答2.1有电路如图P2.1所示，设输入为1u ，输出为2u ，试自选状态变量并列写出其状态空间表达式。
图P2.1解 此题可采样机理分析法，首先根据电路定律列写微分方程，再选择状态变量，求得相应的系统状态空间表达式。
也可以先由电路图求得系统传递函数，再由传递函数求得系统状态空间表达式。
这里采样机理分析法。
设1C 两端电压为1c u ，2C 两端的电压为2c u ，则212221c c c du u C R u u dt++= (1) 112121c c c du u duC C dt R dt+= (2) 选择状态变量为11c x u =，22c x u =，由式(1)和(2)得：1121121121212111c c c du R R C u u u dt R R C R C R C +=--+ 2121222222111c c c du u u u dt R C R C R C =--+ 状态空间表达式为：12111211212121212122222221111111R R C x x x u R R C R C R C x x x u R C R C R C y u u x +⎧=--+⎪⎪⎪=--+⎨⎪⎪==-⎪⎩即: 12121121211112222222211111R R C R C R R C R C x x u x x R C R C R C +⎡⎤⎡⎤-⎢⎥⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦--⎢⎥⎢⎥⎣⎦⎣⎦[]11210x y u x ⎡⎤=-+⎢⎥⎣⎦2.2 建立图P22所示系统的状态空间表达式。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一章习题答案1-1 试求图1-27系统的模拟结构图，并建立其状态空间表达式。
图1-27系统方块结构图解：系统的模拟结构图如下：图1-30双输入--双输出系统模拟结构图系统的状态方程如下：u K K x K K x K K x X K x K x x x x J K x J x J K x J K x x J K x x x pp p p n pb1611166131534615141313322211+--=+-==++--===••••••阿 令y s =)(θ，则1x y =所以，系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡••••••654321165432111111112654321000001000000000000010010000000000010x x x x x x y uK K x x x x x x K K K K K K J K J J K J K J K x x x x x x p p pp n p b1-2有电路如图1-28所示。
以电压)(t u 为输入量，求以电感中的电流和电容上的电压作为状态变量的状态方程，和以电阻2R 上的电压作为输出量的输出方程。
L1L2U图1-28 电路图解：由图，令32211,,x u x i x i c ===，输出量22x R y =有电路原理可知：•••+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x C x C x x L x L R x u L x L x L R x =+-=+-=+--=•••写成矢量矩阵形式为：[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000*********x x x R y u L x x x CCL L R L L R x x x 。
1-3 参考例子1-3（P19）.1-4 两输入1u ，2u ，两输出1y ，2y 的系统，其模拟结构图如图1-30所示，试求其状态空间表达式和传递函数阵。
1u 2u 图1-30双输入--双输出系统模拟结构图解：系统的状态空间表达式如下所示：[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡432121432134561243210101000000010*******x x x x y u b b x x x x a a a a a a xx x x⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--+-=-3456121010001)(a a a s a as a sA sI ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--+-=-=--211345612100000001010001)()(b b a a a s a as a sB A sI s W ux []⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--+-=-=--2113456121000000010100010101)()(b b a a a s a as a sB A sIC s W uy1-5系统的动态特性由下列微分方程描述u u u y y y y 23375)2(......++=+++列写其相应的状态空间表达式，并画出相应的模拟结构图。
解：令..3.21y x y x y x ===，，，则有[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡321321321132100573100010x x x y u x x x x x x 。
相应的模拟结构图如下：1-6 （2）已知系统传递函数2)3)(2()1(6)(+++=s s s s s W ,试求出系统的约旦标准型的实现，并画出相应的模拟结构图解：ss s s s s s s s W 31233310)3(4)3)(2()1(6)(22++++-++-=+++= ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡--=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡432143214321313310411100000020*********x x x x y u x x x x x x x x1-7 给定下列状态空间表达式[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡321321321100210311032010x x x y u x x x x x x ‘（1） 画出其模拟结构图 （2） 求系统的传递函数 解：（2）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+-=-=31103201)()(s s s A sI s W )1)(2)(3()3(2)3(2+++=+++=-s s s s s s A sI()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++---++-+++++=--)2)(1(150)3()3(2033)1)(2)(3(1)(21s s s s s s s s s s s s A sI ()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+++++++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++---++-+++++=-=-)3)(12()3()3()1)(2)(3(1210)2)(1(150)3()3(2033)1)(2)(3(1)()(21s s s s s s s s s s s s s s s s s s s s B A sI s W ux[])1)(2()12()1)(2)(3(1)3)(12()3()3(100)()(1+++=+++⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++=-=-s s s s s s s s s s s B A sI C s W uy 1-8 求下列矩阵的特征矢量（3）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=6712203010A 解：A 的特征方程 061166712230123=+++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+---=-λλλλλλλA I 解之得：3,2,1321-=-=-=λλλ当11-=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---3121113121116712203010p p p p p p 解得： 113121p p p -== 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P（或令111-=p ，得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P ）当21-=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---32221232221226712203010p p p p p p 解得： 1232122221,2p p p p =-= 令212=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1423222122p p p P（或令112=p ，得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=21213222122p p p P ）当31-=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---33231333231336712203010p p p p p p 解得： 133313233,3p p p p =-= 令113=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3313323133p p p P1-9将下列状态空间表达式化成约旦标准型（并联分解）⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡32121321321110021357213311201214x x x y y u x x x x x x（2）解：A 的特征方程 0)3)(1(311212142=--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=-λλλλλλA I 1,332,1==λλ当31=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--3121113121113311201214p p p p p p 解之得 113121p p p == 令111=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1113121111p p p P当32=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--1113311201214312111312111p p p p p p 解之得 32222212,1p p p p =+= 令112=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0013222122p p p P当13=λ时，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--332313332313311201214p p p p p p 解之得3323132,0p p p == 令133=p 得 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=1203323133p p p P⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=101201011T ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=-1102112101T⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=-4325183572131102112101B T⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡=302413101201011110021CT 约旦标准型x ~y ux ~x ~⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3024134325181000300131-10 已知两系统的传递函数分别为W 1(s)和W 2(s)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++=2102111)(1s s s s s W ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+++=0114131)(2s s s s W试求两子系统串联联结和并联连接时，系统的传递函数阵，并讨论所得结果 解：（1）串联联结⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++++++++++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+++==)2)(1(1)1(1)4)(3)(2(75)3)(1(121021110114131)()()(2212s s s s s s s s s s s s s s s s s s W s W s W（2）并联联结⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+++±⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++=±=01141312102111)()()(11s s s s s s s s W s W s W1-11 （第3版教材）已知如图1-22所示的系统，其中子系统1、2的传递函数阵分别为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=210111)(1s s s s W ⎥⎦⎤⎢⎣⎡=10012)s (W求系统的闭环传递函数 解：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=2101111001210111)()(211s s s s s s s W s W ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++-++=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-++=+2301121001210111)()(1s s s s s s s s I s W s W I[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++++++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++++=+-320)3(12112012331)()(121s s s s s s s s s s s s s s s W s W I[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++-+=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+-+++++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++++=+=-310)3(1211101)1)(2(3312111112012331)()()()(1121s s s s s s s s s s s s s ss s s s s s s s s s W s W s W I s W1-11（第2版教材） 已知如图1-22所示的系统，其中子系统1、2的传递函数阵分别为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=2121111s s s )s (W ⎥⎦⎤⎢⎣⎡=10012)s (W求系统的闭环传递函数 解：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=212111100121211111s s s s s s )s (W )s (W⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++-++=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+=+232112100121211111s s s s s s s s )s (W )s (W I []⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++-+++++=+-1221232512111s s s s s s s )s (s )s (W )s (W I[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++-++++++++-+++++=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++-++++-++++-++++++=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++-+++++=+=-252)25)(2(66251)25()2()83()1(1121)2(222)2(1)2(32)2(325)1(2112112212325)1()()()()(222322222221111s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s ss s s s s s s s s s s s W s W s W I s W1-12 已知差分方程为)(3)1(2)(2)1(3)2(k u k u k y k y k y ++=++++试将其用离散状态空间表达式表示，并使驱动函数u 的系数b(即控制列阵)为（1）⎥⎦⎤⎢⎣⎡=11b解法1：21112332)(2+++=+++=z z z z z z W )(11)(2001)1(k u k x k x ⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡--=+ [])(11)(k x k y =解法2：)(2)(3)()(3)(2)1()()1(2121221k x k x k y u k x k x k x k x k x +=+--=+=+ [])(23)()(10)(3210)1(k x k y k u k x k x =⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡--=+ 求T,使得⎥⎦⎤⎢⎣⎡=-111B T 得⎥⎦⎤⎢⎣⎡=-10111T 所以 ⎥⎦⎤⎢⎣⎡-=1011T⎥⎦⎤⎢⎣⎡---=⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡=-15041011321010111AT T [][]13101123-=⎥⎦⎤⎢⎣⎡-=CT 所以，状态空间表达式为[])(13)()(11)(1504)1(k z k y k u k z k z -=⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡---=+第二章习题答案2-4 用三种方法计算以下矩阵指数函数At e 。