现代控制理论习题解答

《现代控制理论》第1章习题解答

1.1 线性定常系统和线性时变系统的区别何在？ 答：线性系统的状态空间模型为：

x Ax Bu

y Cx Du

=+=+

线性定常系统和线性时变系统的区别在于：对于线性定常系统，上述状态空间模型中的系数矩阵A ，B ，C 和D 中的各分量均为常数，而对线性时变系统，其系数矩阵A ，B ，C 和

D 中有时变的元素。线性定常系统在物理上代表结构和参数都不随时间变化的一类系统，

而线性时变系统的参数则随时间的变化而变化。

1.2 现代控制理论中的状态空间模型与经典控制理论中的传递函数有什么区别？ 答: 传递函数模型与状态空间模型的主要区别如下:

1.3 线性系统的状态空间模型有哪几种标准形式？它们分别具有什么特点？

答: 线性系统的状态空间模型标准形式有能控标准型、能观标准型和对角线标准型。对于n 阶传递函数

121210

1

110

()n n n n n n n b s b s b s b G s d s a s a s a ------++++=+++++， 分别有

⑴ 能控标准型： []012

101210100000100000101n n n x x u a a a a y b b b b x du

---?????

????????????

???=+??

????????

?

?????----????

?

=+?? ⑵ 能观标准型： []001122110

001000

100010

1n n n b a b a x a x u b a b y x du

---?-??

?????

??-?????

?????=-+????

?

?????

??????-????

?=+??

⑶ 对角线标准型： []1212

001001001n n p p x x u p y c c c x du ???????????????=+??????

???????????=+? 式中的12,,,n p p p 和12,,,n c c c 可由下式给出，

12121012

111012

()n n n n n

n n n n

b s b s b s b

c c c G s

d d s a s a s a s p s p s p ------++++=+=++

+

+++++--- 能控标准型的特点：状态矩阵的最后一行由传递函数的分母多项式系数确定，其余部分具有特定结构，输出矩阵依赖于分子多项式系数，输入矩阵中的元素除了最后一个元素是1外，其余全为0。

能观标准型的特点：能控标准型的对偶形式。 对角线标准型的特点：状态矩阵是对角型矩阵。 1.4 对于同一个系统，状态变量的选择是否惟一？

答：对于同一个系统，状态变量的选择不是惟一的，状态变量的不同选择导致不同的状态空间模型。

1.5 单输入单输出系统的传递函数在什么情况下，其状态空间实现中的直接转移项D 不等

于零，其参数如何确定？

答: 当传递函数)(s G 的分母与分子的阶次相同时，其状态空间实现中的直接转移项D 不等于零。

转移项D 的确定：化简下述分母与分子阶次相同的传递函数

11

10

111)(a s a s a s b s b s b s b s G n n n

n n n n ++++++++=---- 可得：

d a s a s a s c s c s c s G n n n n n ++++++++=----0

11

10

111)( 由此得到的d 就是状态空间实现中的直接转移项D 。

1.6 在例1.

2.2处理一般传递函数的状态空间实现过程中，采用了如图1.12的串联分解，

试问：若将图1.12中的两个环节前后调换，则对结果有何影响？ 答: 将图1.12中的两个环节调换后的系统方块图为： 图中，

32

210

11()a s s a s a s a =+++，2

210()b s b s b s b =++。

由于3s y -相当于对y 作3次积分，故1

()

y m a s =

可用如下的状态变量图表示：

因为2

s b 相当于对

b 作2次微分，故()m

b s u

=可用如下的状态变量图表示:

因此，两个环节调换后的系统状态变量图为

进一步简化，可得系统状态变量图为

取3y x =，2y x =，1y x =，可以得到两个环节调换后的系统的状态空间模型为

001122001

001[0

01]a b x a x b u

a b y x

-????????=-+????-????????=

两个环节调换前的状态空间模型是：

012012010000101[]x x u a a a y b b b x

????

????=+????????---????=

显然，调换前后的状态空间实现是互为对偶的。 1.7 已知系统的传递函数

2()6

()56

Y s s U s s s +=++ 试求其状态空间实现的能控标准形和能观标准形。 答： 系统的能控标准形为:

[]01065161x x u y x

?????=+?????--?????

?=?

系统的能观标准形为:

[]06615101x x u y x

?-????=+?????-?????

?=?

1.8 考虑由下图描述的二阶水槽装置，

图1.18 二阶水槽装置图

该装置可以看成是由两个环节串联构成的系统，它的方块图是：

图1.19 二阶水槽系统的方块图

试确定其状态空间模型。

答：图1.19中两个环节的状态空间模型分别为:

??

?=+-=22

22222x y u b x a x

和 ?

?

?=+-=11111x y u b x a x

又因为21x u u +=，所以

1121111u b x b x a x

++-= 22222u b x a x

+-= 1x y =

进一步将其写成向量矩阵的形式，可得：

??

?

?????????+???????????

?--=??????21212121121000u u b b x x a b a x x []??

?

???=2101x x y

1.9 考虑以下单输入单输出系统:

61166+++=y y y y u

试求该系统状态空间模型的对角线标准形。 答： 由微分方程可得:

3

21)3)(2)(1(6

61166)(32123+++++=+++=+++=

s c s c s c s s s s s s s G

其中，

3)

3)(2(6

lim

11=++=-→s s c s

6)

3)(1(6

lim 22-=++=-→s s c s

3)2)(1(6

lim 33=++=-→s s c s

故该系统状态空间模型的对角线标准形为:

u x x x x x x ????

??????+????????????????????---=??????????111300020001321321 []??

??

??????-=321363x x x y

1.10 已知单输入单输出时不变系统的微分方程为:

()4()3()()6()8()++=++y t y t y t u t u t u t

试求:（1）建立此系统状态空间模型的对角线标准形；

（2）根据所建立的对角线标准形求系统的传递函数。

答: （1）由微分方程可得:

3

45

213486)(222++++=++++=s s s s s s s s G

记

12122525

()(1)(3)13

43c c s s G s s s s s s s ++=

==+++++++，

其中，

23352lim

11=++=-→s s c s ，2

1

152lim 32=++=-→s s c s 。

从输入通道直接到输出通道上的放大系数1=d ，由此可得:

u x x x x

??

????+???????

?????--=??????1130012121 u x x y +????????????=212123 （2） 由于??????--=3001A ，??

????=11B ，???

???=2123C ，1=D ，因此 1()()30113

11011(1)(2)221.50.5

113

G s C sI A B D

s s s s s s -=-++??????=+??????+++??????

=

++++ 1.11 已知系统的传递函数为

25

()(3)(5)

+=

++s G s s s

（1） 采用串联分解方式，给出其状态空间模型，并画出对应的状态变量图； （2） 采用并联分解方式，给出其状态空间模型，并画出对应的状态变量图。 答:（1）将()G s 重新写成下述形式：

125()35s G s s s +=

?++ 每一个环节的状态空间模型分别为：

??

?=+-=11

113x y u x x

和 ??

?+-=+-=1

2122255u x y u x x

又因为11y u =， 所以

??

?-=+-=2121153x x x u x x

2152x x y -=

因此，若采用串联分解方式，则系统的状态空间模型为：

u x x x x ??????+???????????

?--=??????0151032121 []??

?

???-=2152x x y

对应的状态变量图为：

（2）将()G s 重新写成下述形式：

0.5 2.5

()35

S G s s -=

+++ 每一个环节的状态空间模型分别为：

1111

30.5x x u

y x =--??

=? ??

?=+-=2

2225.25x y u x x

又由于

??

?+-=--=u x x u x x

5.255.032211 1212y y y x x =+=+

因此，若采用并联分解方式，则系统的状态空间模型为：

u x x x x ??????-+???????????

?--=??????5.25.050032121 []??

?

???=2111x x y

对应的状态变量图为：

-

-

?

1

x ?

2

x x ?2x -

1

x ?

1

x -

1.12 已知系统的状态空间模型为,x Ax Bu y Cx =+=，写出该系统的特征多项式和传递

函数矩阵。

答: 系统的特征多项式为det()sI A -， 传递函数为1()()G s C sI A B -=-。

1.13 一个传递函数的状态空间实现是否惟一？由状态空间模型导出的传递函数是否惟

一？

答: 一个传递函数的状态空间实现不惟一；而由状态空间模型导出的传递函数是惟一的。

1.14 已知系统的状态空间模型为Cx y Bu Ax x

=+=, ，写出其对偶状态空间模型。 答: 其对偶状态空间模型为:

???=+=x

B y u

C x A x T T T ~

~~ 1.15 两个对偶状态空间模型之间的特征多项式和传递函数有什么关系？

答: 对于互为对偶的 ?

??=+=Cx y Bu Ax x

与

?

??=+=x B y u

C x A x T T T ~~~ ，它们对应的特征多项式分别为det()sI A -和det()T sI A -。由于一个矩阵和其装置的特征多项式是相同的，故互为对偶的两

个状态空间模型具有相同的特征多项式。 它们对应的传递函数分别为

1

1()()()det()

C sI A B

G s C sI A B sI A *--=-=-

1

2()()()det()

T T T

T

T T

T B sI A C G s B sI A C sI A *--=-=-

由于det()det()T

sI A sI A -=-，(

)

()()T

T T T C sI A B

B sI A

C *

*-=-，故对偶状态空间模

型之间的传递函数关系为12()()T G s G s =，即互为转置。 1.16 考虑由以下状态空间模型描述的系统:

010651[11]x x u y x

????

=+????--????=

试求其传递函数。 答: 由于

11()()()G s C sI A B D C sI A B --=-+=-

??

????-+++=--s s s s A sI 6156)5(1)(1

故

[]651

)1(6

5110615116)5(1)(2

2+++=+++=????????????-+++=

s s s s s s s s s s s G 1.17 给定系统的状态空间模型:

10000

4310112011000

01x x u y u ????

????=-+????????---??????

=?

???

求系统的传递函数矩阵。

答: 系统的传递函数为1()()G s C sI A B -=-。由于

????

?

?????+----+-++++++=??

??

??????+-+-=---s s s s s s s s s s s s s s s s A sI 414323321163116121134001

)(2

2

22

31

1

因此，

B A sI

C s G 1)()(--=

?????

?

??????????????+----+-+++??????+++=100100414323321161000013116122

223s s s s s s s s s s s s s ??

????+--++++=s s s s s s s 413231161

2

23 1.18 试用MATLAB 软件求出下列传递函数的状态空间实现:

232

1047160

()1456160

s s G s s s s ++=+++ 答: 执行以下的m-文件：

num=[0 10 47 160]; den=[1 14 56 160]; [A,B,C,D]=tf2ss(num,den)

得到：

??????????---=0100011605614A ， ??

??

?

?????=001B ， []1604710=C ， 0=D

由此可知：

u x x x x x x

????

??????+????????????????????---=??????????0010100011605614321321

[]??

??

??????=3211604710x x x y

1.19 试用MATLAB 软件求以下系统的传递函数:

112233123010011011000[100]x x x x u x x x y x x ????????

????????=--+???????????????

?????????????

=??

????

答: 执行以下m-文件：

A=[0 1 0;-1 -1 0;1 0 0];

B=[0;1;0];

C=[1 0 0]; D=0;

[num,den]=ss2tf(A,B,C,D)

可得：

num = 0 0 1.0000 0 den = 1.0000 1.0000 1.0000 0

因此，系统的传递函数为

s

s s s

s G ++=

2

3)( 1.20 试用MATLAB 软件求以下系统的传递函数:

11122233123210010201001301[001]x x u x x u x x x y x x ??????????????????=+??

?????????????????

?????????????

=??

????

答: 执行以下的m-文件：

A=[2 1 0;0 2 0;0 1 3]; B=[0 1;1 0;0 1]; C=[0 0 1]; D=[0 0];

[num,den]=ss2tf(A,B,C,D,1)

[num,den]=ss2tf(A,B,C,D,2) 可得要求的两个传递函数是

3

21()2

()71612

Y s s U s s s s -=-+- 23

2

2()44

()71612

Y s s s U s s s s -+=-+- 1.21 已知系统的状态空间模型为,x Ax Bu y Cx =+=，取线性变换阵为P ，且x Px =，

写出线性变换后的状态空间模型。 答: 把x Px =代入,x Ax Bu y Cx =+=，得

Px APx Bu y CPx

=+=

因此，线性变换后的等价状态空间模型为：

11x P APx P Bu y CPx

--=+=

1.22 线性变换是否改变系统的特征多项式和极点？简单证明之。 答: 假设系统的状态空间模型为

Du

Cx y Bu Ax x

+=+=

经过线性变换Tx x =后，系统的状态模型变为：

u

D x C y u B x A x

+=+=

其中，

1-=TAT A ，TB B =，1-=CT C ，D D =

由于

1111det()det()det()det()det()det()det()

sI A sI TAT sTT TAT T sI A T sI A -----=-=-=-=-

故线性变换不会改变系统的特征多项式和极点。 1.23 已知以下微分方程描述了系统的动态特性:

32y y y u ++=

（1） 选择状态变量12,x y x y ==，写出系统的状态方程；

（2） 根据（1）的结果，由以下的状态变换:

112212

2x x x x x x =+=--

确定新的状态变量12,x x ，试写出关于新状态变量12,x x 的状态空间模型。

答: (1) 由12,x y x y ==可得

??

?

??=+--==11222123x

y u x x x x x 写成矩阵向量形式，可得

[]???

???

???????=??

????+????????????--=??????21212101103210x x y u x x x x (2) 由于112212,2x x x x x x =+=--，即

11221112x x x x ??????=??????--?

?????

容易验证这是一个等价线性变换，故可得

[]???

???

????

???=??????-+????????????--=??????21212111112001x x y u x x x x

1.24 给定系统

[]010100x x u a b d y x

????=+????--????= 试确定参数b a ,和d 的值，以使得该系统模型能等价地转换成以下的对角型

[]30101155z z u y z

-????=+????-????=- 答: 由对角型状态空间模型可知

2551010

()31(3)(1)43

G s s s s s s s -=

+==++++++ 而从原状态空间模型则可得传递函数

210()d

G s s bs a

=

++

由于等价的状态空间模型具有相同的传递函数，故经比较系数可得：

341

a b d ===