当前位置:文档之家› 《自动控制原理》习题及解答03-Ed

《自动控制原理》习题及解答03-Ed

《自动控制原理》习题及解答03-Ed
《自动控制原理》习题及解答03-Ed

第三章习题及答案

3-1 已知系统脉冲响应如下,试求系统闭环传递函数Φ(s)。

t e t k 25.10125.0)(-=

解 Φ()()./(.)s L k t s ==+00125125

3-2 设某高阶系统可用下列一阶微分方程近似描述

T c t c t r t r t ??

+=+()()()()τ

其中,0<(T-τ)<1。试证系统的动态性能指标为 T T T t d ??

?

??????

??-+=τln 693.0

t T r =22. T T T t s ??

???

?

-+=)ln(

3τ 解 设单位阶跃输入s

s R 1)(= 当初始条件为0时有:

1

1

)()(++=Ts s s R s C τ 1

11

11)(+--

=

?

++=

Ts T s s Ts s s C τ

τ C t h t T T

e t T

()()/==---1τ 1) 当 t t d = 时

h t T T

e t t

d ()./==---051τ

12=--T T e t T d τ/ ; T

t T T d

-??? ??-=-τln 2ln ?????

???? ??-+=∴

T T T t d τln 2ln

2) 求t r (即)(t c 从0.1到0.9所需时间) 当 T

t e T

T t h /219.0)(---

==τ; t T T T 201=--[ln()ln .]τ

当 T

t

e

T

T t h /111.0)(---==τ; t T T T 109=--[ln()ln .]τ 则 t t t T T r =-==21

09

01

22ln ... 3) 求 t s

T

t s s e

T

T t h /195.0)(---==τ ∴=--t T T T s [ln ln .]τ005=-+T T T

[ln ln ]τ20=+-T T T [ln

]3τ

3-3 一阶系统结构图如题3-3图所示。要求系统闭环增益2=ΦK ,调节时间4.0≤s t (s ),试确定参数21,K K 的值。

解 由结构图写出闭环系统传递函数

111)(212211211

+=+=+

=ΦK K s

K K K s K s

K K s K s

令闭环增益21

2

==

ΦK K , 得:5.02=K 令调节时间4.03

32

1≤=

=K K T t s ,得:151≥K 。 3-4 在许多化学过程中,反应槽内的温度要保持恒定, 题3-4图(a )和(b )分别为开环和闭环温度控制系统结构图,两种系统正常的K 值为1。

(1) 若)(1)(t t r =,0)(=t n 两种系统从开始达到稳态温度值的63.2%各需多长时间? (2) 当有阶跃扰动1.0)(=t n 时,求扰动对两种系统的温度的影响。

解 (1)对(a )系统: 1

101

110)(+=

+=

s s K s G a , 时间常数 10=T Θ 632.0)(=T h (a )系统达到稳态温度值的63.2%需要10个单位时间;

对(a )系统:1101

10101100

10110100

)(+=+=Φs s s b , 时间常数 10110=

T Θ 632.0)(=T h (b )系统达到稳态温度值的63.2%需要0.099个单位时间。

(2)对(a )系统: 1)

()

()(==

s N s C s G n 1.0)(=t n 时,该扰动影响将一直保持。

对(b )系统: 101101

101

1010011)

()

()(++=++

==

Φs s s s N s C s n 1.0)(=t n 时,最终扰动影响为001.0101

1

1.0≈?

。 3-5 一种测定直流电机传递函数的方法是给电枢加一定的电压,保持励磁电流不变,测出电机的稳态转速;另外要记录电动机从静止到速度为稳态值的50%或63.2%所需的时间,利用转速时间曲线(如题3-5图)和所测数据,并假设传递函数为

)

()()()(a s s K

s V s s G +=Θ=

可求得K 和a 的值。

若实测结果是:加10伏电压可得每分钟1200转的稳态转速,而达到该值50%的时间为1.2秒,试求电机传递函数。 [提示:注意

)()(s V s Ω=a s K +,其中dt

d t θ

ω=)(,单位是弧度/秒] 解 依题意有: 10)(=t v (伏)

ππ

ω4060

21200)(=?=∞ (弧度/秒) (1)

πωω20)(5.0)2.1(=∞= (弧度/秒) (2)

设系统传递函数 a

s K

s V s s G +=Ω=

)()()(0 应有

πω401010lim )()(lim )(0

00

==+??

=?=∞→→a

K a s K s s s V s G s s s (3)

[][]

at

e a K a s s L a K a s s K L s V s G L t -----=??????+-=??

??

??+=?=1101110)(10)()()(1101ω 由式(2),(3)

[][]

ππω20140110)2.1(2.12.1=-=-=

--a a e e a

K

得 5.012.1=--a

e

解出 5776.02

.15

.0ln =-=

a (4) 将式(4)代入式(3)得 2586.74==a K π

3-6 单位反馈系统的开环传递函数)

5(4

)(+=

s s s G ,求单位阶跃响应)(t h 和调节时间

t s 。

解:依题,系统闭环传递函数

)1)(1(4

)

4)(1(4

454)(2

12T s T s s s s s s ++=

++=++=Φ ???==25.0121T T

)4)(1(4

)()()(++=

Φ=s s s s R s s C =4

1210++++s C s C s C

1)

4)(1(4

lim

)()(lim 00

0=++=Φ=→→s s s R s s C s s

3

4

)4(4lim

)()()1(lim 0

1

1-=+=Φ+=→-→s s s R s s C s s

3

1

)1(4lim

)()()4(lim 0

4

2=+=Φ+=→-→s s s R s s C s s

t t e e t h 43

1

341)(--+-=

Θ 421

=T T , ∴3.33.3111==???

? ??=T T T t t s s 。 3-7 设角速度指示随动系统结构图如题3-7图。若要

求系统单位阶跃响应无超调,且调节时间尽可能短,问开环增益K 应取何值,调节时间s t 是多少?

解 依题意应取 1=ξ,这时可设闭环极点为

02,11T -=λ。

写出系统闭环传递函数

K

s s K

s 101010)(2

++=Φ 闭环特征多项式

2

0022

02

1211010)(???

? ??++=???? ?

?+=++=T s T s T s K s s s D 比较系数有 ????

???=???

? ??=K T T 10110

2

2

00 联立求解得 ?

?

?==5.22

.00K T 因此有 59.075.40''==T t s 1''<

3-8 给定典型二阶系统的设计指标:超调量

σ%5≤%,调节时间 s t 3<(s ),峰值时间1

解 依题

σ%5≤%, ? )45(707.0?≤≥βξ;

35

.3

<=

n

s t ωξ,? 17.1>n ωξ;

n

p t ωξπ

21-=

1<,? n ωξ21-14.3>

综合以上条件可画出满足要求的特征根区域如图解3-8所示。

3-9 电子心律起博器心率控制系统结构图如题3-9图所示,其中模仿心脏的传递函数相当于一纯积分环节,要求:

(1) 若ξ=0.5对应最佳响应,问起博器增益K 应取多大?

(2) 若期望心速为60次/分钟,并突然接通起博器,问1秒钟后实际心速为多少?瞬时最

大心速多大?

解 依题,系统传递函数为

2

222205

.005.0105.0)(n

n n s s K s s K

s ωξωω++=+

+=Φ ???

?????==n n K

ωξω205.0105.0 令 5.0=ξ可解出 ???==2020

n

K ω

将 1=t (秒)代入二阶系统阶跃响应公式 ()

βωξξ

ξω+---

=-t e t h n t n 22

1sin 11)(

可得 000024.1)1(=h (次/秒)=00145.60(次/分)

5.0=ξ时,系统超调量 σ%=1

6.3% ,最大心速为

163.1163.01(=+=)p t h (次/秒)=78.69(次/分)

3-10 机器人控制系统结构图如题3-10图所示。试确定参数21,K K 值,使系统阶跃响应的峰值时间

5.0=p t (s ),超调量σ%2=%。

解 依题,系统传递函数为

22212121211

2.)1()

1()1(1)1()(n n n

s s K K s K K s K s s s K K s s K s ωξωω++=+++=+++

+=ΦΦ

由 ??

???=-=≤=--5.0102.0212n p o

o t e ωξπσξπξ 联立求解得

??

?==1078

.0n

ωξ 比较)(s Φ分母系数得

??

?

??=-===146.012100

1221K K K n n ξωω 3-11 某典型二阶系统的单位阶跃响应如题3-11图所示。试确定系统的闭环传递函数。

解 依题,系统闭环传递函数形式应为

2

2

22.)(n

n n

s s K s ωξωω++=ΦΦ 由阶跃响应曲线有:

21

)(lim )()(lim (0

==?

Φ=Φ=∞Φ→→K s

s s s R s s h s s )

???

????

=-===-=--o o

o o n p e t 25225.221212ξξπσξωπ 联立求解得 ???==717.1404

.0n

ωξ,所以有

95

.239.19

.5717.1717.1404.02717.12)(2222++=+??+?=Φs s s s s

3-12 设单位反馈系统的开环传递函数为 )

12.0(5

.12)(+=

s s s G

试求系统在误差初条件1)0(,10)0(==e e &作用下的时间响应。

解 依题意,系统闭环传递函数为 5

.6255

.62)(1)()()()(2++=+==

Φs s s G s G s R s C s 当0)(=t r 时,系统微分方程为 0)(5.62)(5)(=+'+''t c t c t c 考虑初始条件,对微分方程进行拉氏变换

[]

[]0)(5.62)0()(5)0()0()(2=+-+'--s C c s C s c c s s C s

整理得 ()

())0()0(5)(5.6252

c c s s C s s

'++=++ (1)

对单位反馈系统有 )()()(t c t r t e -=, 所以

1

10)0()0()0(10

1000()0()0(-=-='-'='-=-=-=e r c e r c )

将初始条件代入式(1)得 2

225.7)5.2(26

)5.2(105.6255110)(++++-=++--=

s s s s s s C 2

2225.7)5.2(5

.747.35.7)5.2()5.2(10

++-+++-=s s s

)8.705.7sin(6.105.7sin 47.35.7cos 10)(5.25.25.2?+-=--=---t e t e t e

t c t t t

3-13设题3-13图(a)所示系统的单位阶跃响应如题3-13图(b )所示。试确定系统参数,

1

K

2

K和a。

解由系统阶跃响应曲线有

?

?

?

?

?

=

-

=

=

=

o

o

o

o

p

t

h

3.

33

3

)3

4(

1.0

3

)

(

σ

系统闭环传递函数为

2

2

2

2

1

2

2

1

2

)

(

n

n

n

s

s

K

K

as

s

K

K

s

ω

ξω

ω

+

+

=

+

+

=

Φ(1)

?

?

?

?

?

=

=

=

-

=

-

-

o

o

o

o

n

p

e

t

3.

33

1.0

1

2

1

2

ξ

ξπ

σ

ω

ξ

π

联立求解得

?

?

?

=

=

28

.

33

33

.0

n

ω

ξ

由式(1)

?

?

?

=

=

=

=

22

2

1108

2

1

n

n

a

K

ξω

ω

另外3

lim

1

)

(

lim

)

(

2

1

2

2

1

=

=

+

+

=

?

Φ

=

K

K

as

s

K

K

s

s

s

h

s

s

3-14 题3-14图是电压测量系统,输入电压

)(t

e

t

伏,输出位移)

(t

y厘米,放大器增益

10

=

K,丝杠每转螺距1毫米,电位计滑臂每移

动1厘米电压增量为0.4伏.当对电机加10伏阶

跃电压时(带负载)稳态转速为1000转/分钟,

达到该值63.2%需要0.5秒。画出系统方框图,

求出传递函数)(/)(s E s Y t ,并求系统单位阶跃响应的峰值时间p t 、超调量o o

σ、

调节时间s t

和稳态值

)(∞h 。

解 依题意可列出环节传递函数如下

比较点: )()()(s F s E s E t -= 伏 放大器:

10)

()

(==K s E s U a 电动机: 1

5.03515.060101000

1)()(+=+?=+=Ωs s s T K s U s m m a (转/秒/?伏)

丝杠:

1.0)

()

(1==ΘK s s Y (厘米?转) 电位器:

4.0)

()

(2==K s Y s F (伏?厘米) 画出系统结构图如图解3-14所示

系统传递函数为

3

42310

)()()(2+

+==Φs s s E s Y s t ???????===866.02232n n ωξω

∴ 44.512

''=-=

n

p t ωξπ

o o o o e 433.02

1==--ξξπ

σ

5.35

.3''==

n

s t ξω

5.21

)(lim )(0

=?

Φ=∞→s

s s h s 3-15 已知系统的特征方程,试判别系统的稳定性,并确定在右半s 平面根的个数及纯虚根。

(1)1011422)(2

3

4

5

+++++=s s s s s s D =0 (2)483224123)(2

3

4

5

+++++=s s s s s s D =0 (3)22)(4

5

--+=s s s s D =0

(4)0502548242)(2

3

4

5

=--+++=s s s s s s D

解(1)1011422)(2

3

4

5

+++++=s s s s s s D =0

Routh : S 5 1 2 11

S 4 2 4 10 S 3 ε 6 S 2

εε124- 10

S 6 S 0 10

第一列元素变号两次,有2个正根。

(2)483224123)(2

3

4

5

+++++=s s s s s s D =0

Routh : S 5 1 12 32

S 4 3 24 48

S 3

3122434?-= 32348

316?-= 0 S 2

424316

4

12?-?= 48 S 121644812

0?-?= 0 辅助方程 124802s +=,

S 24 辅助方程求导:24s=0 S 0 48

系统没有正根。对辅助方程求解,得到系统一对虚根 s j 122,=±。 (3)22)(4

5

--+=s s s s D =0

Routh : S 5 1 0 -1

S 4 2 0 -2 辅助方程 0224=-s S 3 8 0 辅助方程求导 083

=s

S 2 ε -2 S ε16

S 0 -2

第一列元素变号一次,有1个正根;由辅助方程0224

=-s 可解出: ))()(1)(1(2224

j s j s s s s -+-+=-

))()(1)(1)(2(22)(4

5

j s j s s s s s s s s D -+-++=--+= (4)0502548242)(2

3

4

5

=--+++=s s s s s s D Routh : S 5 1 24 -25

S 4 2 48 -50 辅助方程 05048224=-+s s S 3 8 96 辅助方程求导 09683

=+s s

S 2 24 -50 S 338/3

S 0 -50

第一列元素变号一次,有1个正根;由辅助方程0504822

4=-+s s 可解出: )5)(5)(1)(1(2504822

4

j s j s s s s s -+-+=-+

)5)(5)(1)(1)(2(502548242)(2345j s j s s s s s s s s s s D -+-++=--+++=

3-16 题3-16图是某垂直起降飞机的高度控制系统结构图,试确定使系统稳定的K 值范围。

解 由结构图,系统开环传递函数为:

)

4()

124()(232++++=s s s s s K s G

??

?==3

4

v K K k 系统型别开环增益

0244)(2

345=+++++=K Ks Ks s s s s D Routh : S 5 1 4 2K S 4 1 4K K

S 3 -4(1-K) K 1

S 2 )

1(4)1615(K K K -- K 067.11516=>?

K

S )

1(41647322

K K K --+- 933.0536.0<

S 0 K 0>?K

∴使系统稳定的K 值范围是: 933.0536.0<

3-17 单位反馈系统的开环传递函数为

)

5)(3()(++=

s s s K

s G

为使系统特征根的实部不大于-1,试确定开环增益的取值范围。

解 系统开环增益 15K K k =。特征方程为: 0158)(2

3

=+++=K s s s s D 做代换 1-'=s s 有:

0)8(25)1(15)1(8)1()(2323=-+'+'+'=+-'+-+-'='K s s s K s s s s D

Routh : S 3 1 2 S 2 5 K-8 S 5

18K - 18

K S 0 K-8 8>?

K

使系统稳定的开环增益范围为:

15

1815158〈K K k =< 。 3-18 单位反馈系统的开环传递函数为

)

12)(1()

1()(+++=

s Ts s s K s G

试在满足 1,0>>K T 的条件下,确定使系统稳定的T 和K 的取值范围,并以T 和K 为坐标画出使系统稳定的参数区域图。

解 特征方程为:

0)1()2(2)(2

3

=+++++=K s K s T Ts s D Routh : S 3 T 2 K +1 0>?T S 2 T +2 K 2->?T S

T

TK K +-+221 1

4

2-+

?

K

综合所得条件,当1>K 时,使系统稳定的参数取值 范围如图解3-18中阴影部所示。

3-19 题3-19图是核反应堆石墨棒位置控制闭环系统,其目的在于获得希望的辐射水

平,增益4.4就是石墨棒位置和辐射水平的变换系数,辐射传感器的时间常数为0.1秒,直流增益为1,设控制器传递函数1)(=s G c 。

(1) 求使系统稳定的功率放大器增益K 的取值范围; (2) 设20=K ,传感器的传递函数1

1

)(+=

s s H τ(τ不一定是0.1),求使系统稳定的τ的取值范围。

解 (1)当控制器传递函数1)(=s G c 时

K

s s s s

K s R s C s 64.2)11.0)(6()

11.0(64.2)()()(++++==

Φ 04.2660164.26)10)(6()(2

3

=+++=+++=K s s s K s s s s D

4.2636.360

164.269604.2616601:0

1

23

>→<→-K K

s K K

s

K s s Routh

∴ 36.360<

(2)20=K ,1

1

)(+=s s H τ时 8

.52)1)(6()

1(8.52)()()(++++==Φs s s s s R s C s ττ

08.526)16(8.52)1)(6()(23=++++=+++=s s s s s s s D τττ

8

.52357

.00

1

68.166167.08.52166:0

1

23

s s

s s Routh <→+-->→+ττττττ

∴ 357.00<<τ

3-20 题3-20图是船舶横摇镇定系统结构图,引入内环速度反馈是为了增加船只的阻尼。

(1) 求海浪扰动力矩对船只倾斜角的传递函数

)

()

(s M s d Θ;

(2) 为保证d M 为单位阶跃时倾斜角θ的值不超过0.1,且系统的阻尼比为0.5,求a K 、

1K 和g K 应满足的方程;

(3) 取a K =1时,确定满足(2)中指标的1K 和g K 值。

解 (1)

)5.01()5.02.0(5.01

2.05.012.05.0112.05.0)()(122

122a g a a g a d K K s K K s s s K K s s s K K s s s M s ++++=++++++++=Θ (2)令: 1.05.015

.0)()(1lim )()()(lim )(100

≤+=Θ??=Θ?

=∞→→a

d s d d s K K s M s s s s M s s M s θ

得 81≥a K K 。 由 )()(s M s d Θ 有: ??

?

??=+=+=5.025.02.05.011n a g

n K

K K K ωξω, 可得 a g a K K K K 15.0125.02.0+=+

(3)1=a K 时,81≥K ,525.02.0≥+g K ,可解出 072.4≥g K 。

3-21 温度计的传递函数为

1

1

+Ts ,用其测量容器内的水温,1min 才能显示出该温度的98%的数值。若加热容器使水温按10oC/min 的速度匀速上升,问温度计的稳态指示误差有多大?

解法一 依题意,温度计闭环传递函数

=

Φ)(s 1

1

+Ts 由一阶系统阶跃响应特性可知:o o T h 98)4(=,因此有 min 14=T ,得出 min 25.0=T 。 视温度计为单位反馈系统,则开环传递函数为

Ts

s s s G 1

)(1)()(=Φ-Φ=

?

?

?==11v T

K 用静态误差系数法,当t t r ?=10)( 时,C T K

e ss ?===

5.21010

解法二 依题意,系统误差定义为 )()()(t c t r t e -=,应有 1

111)()(1)()()(+=+-=-==

ΦTs Ts

Ts s R s C s R s E s e C T s Ts Ts s

s R s s e s e s ss ?==?+=Φ=→→5.21010

1lim )()(lim 2

3-22 系统结构图如题3-22图所示。试求局部反馈加入前后系统的静态位置误差系数、静态速度误差系数和静态加速度误差系数。

解 局部反馈加入前,系统开环传递函数为

)

1()

12(10)(2

++=

s s s s G

∞==→)(lim 0

s sG K s v

10)(lim 20

==→s G s K s a

局部反馈加入后,系统开环传递函数为

)20()12(1012011(10

12)(2

+++=++

+?+=s s s s s s s s s s G )

()

∞==→)(lim 0

s G K s p

5.0)(lim 0

==→s sG K s v

0)(lim 20

==→s G s K s a

3-23 已知单位反馈系统的开环传递函数为

)

22)(4()

1(7)(2++++=

s s s s s s G

试分别求出当输入信号t t t r ),(1)(=和2

t 时系统的稳态误差[)()()(t c t r t e -=]。

==→)(lim 0

s G K s p

解 )22)(

4()

1(7)(2

++++=

s s s s s s G ???==1

87v K 由静态误差系数法

)(1)(t t r =时, 0=ss e

t t r =)(时, 14.17

8

===K A e ss

2)(t t r =时, ∞=ss e

3-24 系统结构图如题3-24图所示。已知)(1)()()(21t t n t n t r ===,试分别计算

)()(),(21t n t n t r 和作用时的稳态误差,并说明积分环节设置位置对减小输入和干扰作用下的

稳态误差的影响。

解 )

1)(1()(21++=

s T s T s K

s G

?

??=1v K

)(1)(t t r =时, 0=ssr e ;

K

s T s T s s T s T s T s K s T s s N s E s en ++++-=+++

+-

==Φ)1)(1()1()

1)(1(1)1(1

)()

()(21121211

)(1)(1t t n =时, K

s s s s N s s e en s en s ssn 11)

(lim )()(lim 1110

10

-=Φ=Φ=→→ K

s T s T s s T s s T s T s K s T s N s E s en ++++-=+++

+-

==Φ)1)(1()1()

1)(1(1)1(1

)()

()(21121222

)(1)(2t t n =时, 01

)

(lim )()(lim 2120

20

=Φ=Φ=→→s

s s s N s s e en s en s ssn 在反馈比较点到干扰作用点之间的前向通道中设置积分环节,可以同时减小由输入和干扰因引起的稳态误差。

3-25 系统结构图如题3-25图所示,要使系统对)(t r 而言是II 型的,试确定参数0K 和

τ的值。

解 )1(

)1)(1()

1()

1)(1()1(1)1)(1()1()(02121021+-+++=

+++-

+++=s K K s T s T s K s T s T s K K s T s T s K s G ττττ )

1()()

1(0021221K K s K K T T s T T s K -+-+++=

ττ

依题意应有:???=-+=-00

1021

0τK K T T K K 联立求解得

?

?

?+==2101T T K

K τ 此时系统开环传递函数为 2

2121)()(s T T K

s T T K s G ++=

考虑系统的稳定性,系统特征方程为

0)()(21221=+++=K s T T K s T T s D

当 1T ,2T ,0>K 时,系统稳定。

3-26 宇航员机动控制系统结构图如题3-26图所示,其中控制器可以用增益2K 来表示。

宇航员及其装备的总转动惯量2

25m kg I ?=。

(1) 当输入为斜坡信号t t r =)((m )时,试确定3K 的取值,使系统稳态误差ss

e 1

=(cm ); (2) 采用(1)中的3K 值,试确定21,K K 的取值,使系统超调量σ%限制在10%以内。

解 (1)系统开环传递函数为

)

()

()()()(3212

132121I

K K K s s I K K K K K s I s K K s E s C s G +=+==

?????==

1

13v K K t t r =)(时,令 01.01

3≤==

K K

e ss , 可取 01.03=K 。 (2)系统闭环传递函数为

I

K K s I K K K s I K K s R s C s 2

132122

1)

()()(++==Φ

???

???

?==I K K K I K K n 2213

21ξω 由 o o o o

e

102

1≤=--ξξπ

σ

,可解出 592.0≥ξ。取 6.0=ξ进行设计。

将25=I ,01.03=K 代入6.022

13==I

K K K ξ表达式,可得

36000021≥K K

3-27 大型天线伺服系统结构图如题3-27图所示,其中ξ=0.707,n ω=15,τ=0.15(s)。

(1) 当干扰)(110)(t t n ?=,输入0)(=t r 时,为保证系统的稳态误差小于0.01o,试确定

a k 的取值;

(2) 当系统开环工作(a k =0),且输入0)(=t r 时,确定由干扰)(110)(t t n ?=引起的系

统响应稳态值。

相关主题
文本预览
相关文档 最新文档