第三章 线性系统的时域分析与校正
习题及答案
3-1 已知系统脉冲响应
t e t k 25.10125.0)(-=
试求系统闭环传递函数)(s Φ。
解 Φ()()./(.)s L k t s ==+00125125 3-2 设某高阶系统可用下列一阶微分方程
T c t c t r t r t ??
+=+()()()()τ
近似描述,其中,1)(0<-<τT 。试证系统的动态性能指标为 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 从1.0到9.0所需时间)
当 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)(---==τ ]ln 3[]20ln [ln ]05.0ln [ln T
T T T T T T T T t s τ
ττ-+=+-=--=∴
3-3 一阶系统结构图如图3-45所示。要求系统闭环增益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-46(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 )系统:1
101
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-47)和所测数据,并假设传递函数为
)
()()()(a s s K
s V s s G +=Θ=
可求得K 和a 的值。
若实测结果是:加10V 电压可得1200min r 的稳态转速,而达到该值50%的时间为1.2s ,试求电机传递函数。
提示:注意
a
s K
s V s +=Ω)()(,其中dt d t θω=
)(,单位是s rad
解 依题意有: 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
12
T s T s s s s s s ++
=++=++=
Φ ??
?==25
.01
21T T
4
1)4)(1(4
)()()(210++++=++=
Φ=s C s C s C s s s s R s s C
1)
4)(1(4
lim
)()(lim 00
0=++=Φ=→→s s s R s s C s s
34
)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-48所示。若要求系统单位阶跃响应无超调,且
调节时间尽可能短,问开环增益K 应取何值,调节时间s t 是多少?
解 依题意应取 1=ξ,这时可设闭环极点为
02,11-=λ。
写出系统闭环传递函数
K
s s K
s 101010)(2
++=Φ 闭环特征多项式
2
002
2
021211010)(???
?
??++=????
??+=++=T s T s T s K s s s D
比较系数有 ????
???=???
? ??=K T T 10110
2
2
00 联立求解得 ???==5.22.00K T 因此有 159.075.40''<''==T t s
3-8 给定典型二阶系统的设计指标:超调量
%5%≤σ,调节时间 s t s 3<,峰值时间s t p 1<,
试确定系统极点配置的区域,
以获得预期的响应特性。
解
依题
%5%≤σ, )45(707.0?≤≥?βξ;
35
.3<=
n
s t ωξ, 17.1>?n ωξ;
n
p t ωξπ
21-=
1<, 14.312>-?n ωξ
综合以上条件可画出满足要求的特征根区域如图解3-8所示。
3-9 电子心脏起博器心律控制系统结构图如题3-49图所示,其中模仿心脏的传递函数相当于一纯积分环节。
(1) 若5.0=ξ对应最佳响应,问起博器增益K 应取多大?
(2) 若期望心速为60次/min ,并突然接通起博器,问1s 钟后实际心速为多少?瞬时最大心
速多大?
解 依题,系统传递函数为
2
2
22205
.005.0105.0)(n
n n s s K s s K
s ωξωω++=+
+=Φ ???
?????==n n K
ωξω205.0105.0 令 5
.0=ξ可解
??
?==20
20
n K ω
将 s t 1=代入二阶系统阶跃响应公式
()
βωξξ
ξω+---
=-t e t h n t n 22
1sin 11)(
可得 min 00145.60000024.1)1(次次==s h
5.0=ξ时,系统超调量 %3.16%=σ,最大心速为
min 78.69163.1163.01(次次)==+=s t h p
3-10 机器人控制系统结构图如图3-50所示。试确定参数21,K K 值,使系统阶跃响应的峰值时间5.0=p t s ,超调量%2%=σ。
解 依题,系统传递函数为
2
22
12121211
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.0121001221K K K n n ξωω 3-11 某典型二阶系统的单位阶跃响应如图3-51所示。试确定系统的闭环传递函数。
解 依题,系统闭环传递函数形式应为
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-52(a )所示系统的单位阶跃响应如图3-52(b )所示。试确定系统参数,1K 2
K 和a 。
解 由系统阶跃响应曲线有
???
??=-===∞o
o o o
p t h 3.333)34(1.03
)(σ
系统闭环传递函数为
2
2
2
2122
12)(n
n n s s K K as s K K s ωξωω++=++=Φ (1)
由 ?????
===-=--o o o
o n
p e t 3.331.012
12
ξξπσωξπ 联立求解得 ???==28.3333.0n ωξ 由式(1)???====22
21108
2
1n n a K ξωω
另外 3lim 1
)(lim )(21
22100
==++=?
Φ=∞→→K K as s K K s s s h s s 3-14 图3-53所示是电压测量系统,输入电压)(t e t 伏,输出位移)(t y 厘米,放大器增益10=K ,丝杠每转螺距1mm ,电位计滑臂每移动1厘米电压增量为0.4V 。当对电机加10V 阶跃电压时(带负载)稳态转速为1000min r ,达到该值63.2%需要0.5s 。画出系统方框图,求出传递函数)(/)(s E s Y ,并求系统单位阶跃响应的峰值时间p t 、超调量o
o σ
、调节时间s t 和稳态值)(∞h 。
解 依题意可列出环节传递函数如下 比较点: )()()(s F s E s E t -= V 放大器:
10)
()
(==K s E s U a 电动机: 1
5.03515.060101000
1)()(+=+?=+=Ωs s s T K s U s m m a r/s/V
丝杠:
1.0)
()
(1==ΘK s s Y cm/r 电位器:
4.0)
()
(2==K s Y s F V/cm 画出系统结构图如图解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)01011422)(2
3
4
5
=+++++=s s s s s s D (2)0483224123)(2
3
4
5
=+++++=s s s s s s D (3)022)(45=--+=s s s s D
(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)(2345+++++=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 辅助方程求导:024=s S 0 48
系统没有正根。对辅助方程求解,得到系统一对虚根 s j 122,=±。 (3)022)(45=--+=s s s s D
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(2224j s j s s s s -+-+=-
))()(1)(1)(2(22)(4
5j s j s s s s s s s s D -+-++=--+= (4)0502548242)(2
345=--+++=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-54是某垂直起降飞机的高度控制系统结构图,试确定使系统稳定的K 值范围。 解 由结构图,系统开环传递函数为:
)
4()
124()(232++++=s s s s s K s G
?
?
?==34
v K K k 系统型别开环增益 0244)(2345=+++++=K Ks Ks s s s s D Routh : S 5 1 4 2K S 4 1 4K K
S 3 )1(4K -- 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 8-K 8>? K 使系统稳定的开环增益范围为: 15 1815158<= ) 12)(1() 1()(+++= s Ts s s K s G 试在满足 1,0>>K T 的条件下,确定使系统稳定的T 和K 的取值范围,并以T 和K 为坐标画出使系统稳定的参数区域图。 解 特征方程为: 0)1()2(2)(23=+++++=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-55是核反应堆石墨棒位置控制闭环系统,其目的在于获得希望的辐射水平,增益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()(23=+++=+++=K s s s K s s s s D 4.2636.360 164.269604.2616601:0 1 2 3 >→<→-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 2 3 s s s s Routh <→+-->→+ττττττ ∴ 357.00<<τ 3-20 图3-56是船舶横摇镇定系统结构图,引入内环速度反馈是为了增加船只的阻尼。 (1) 求海浪扰动力矩对船只倾斜角的传递函数 ) () (s M s N Θ; (2) 为保证N M 为单位阶跃时倾斜角θ的值不超过0.1,且系统的阻尼比为0.5,求2K 、1 K 和3K 应满足的方程; (3) 取2K =1时,确定满足(2)中指标的1K 和3K 值。 解 (1) )5.01()5.02.0(5.01 2.05.012.05.0112.05 .0)()(21322 2 12322K K s K K s s s K K s s s K K s s s M s a N ++++=++++++++=Θ (2)令: 1.05.015 .0)()(1lim )() ()(lim )(2 100 ≤+=??=? =∞→→K K s M s s s s M s s M s N s N N s ΘΘθ 得 821≥K K 。 由 )() (s M s N Θ 有: ?? ? ??=+=+=5.025.02.05.013 23 1n n K K K K ωξω, 可得 21325.0125.02.0K K K K +=+ (3)12=K 时,81≥K ,525.02.03≥+K ,可解出 072.43≥K 。 3-21 温度计的传递函数为 1 1 +Ts ,用其测量容器内的水温,1min 才能显示出该温度的98%的数值。若加热容器使水温按10oC/min 的速度匀速上升,问温度计的稳态指示误差有多大? 解法一 依题意,温度计闭环传递函数 1 1 )(+= ΦTs s 由一阶系统阶跃响应特性可知:o o T h 98)4(=,因此有 min 14=T ,得出 min 25.0=T 。 视温度计为单位反馈系统,则开环传递函数为 Ts s s s G 1 )(1)()(=Φ-Φ= ? ??==11v 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 20 3-22 系统结构图如图3-57所示。试求局部反馈加入前、后系统的静态位置误差系数、静态速度误差系数和静态加速度误差系数。 解 局部反馈加入前,系统开环传递函数为 ) 1() 12(10)(2++= s s s s G ∞==∞ →)(lim s G K s p ∞==→)(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 -=]。 解 )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-58所示。已知 )(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 ???=1 v 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-59所示,要使系统对)(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-60所示。其中控制器可以用增益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 , 可取 30.01K =。 (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 22132 1ξω 由 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-61所示,其中ξ=0.707,n ω=15,τ=0.15s 。 (1) 当干扰)(110)(t t n ?=,输入0)(=t r 时,为保证系统的稳态误差小于0.01o,试确定a K 的取值; (2) 当系统开环工作(a K =0),且输入0)(=t r 时,确定由干扰)(110)(t t n ?=引起的系统 响应稳态值。 解 (1)干扰作用下系统的误差传递函数为 2 222 )2)(1()1()() ()(n a n n n en K s s s s s s N s E s ωωξωττω+++++-==Φ )(110)(t t n ?=时, 令 a en s en s ssn K s s s s s N s e 10 )(10lim )()(lim 0 = Φ?? =Φ??=→→01.0≤ 得: 1000≥a K (2)此时有 ) 2(10)()2()()(22 22222n n n n n n s s s s N s s s s C s E ωξωωωξωω++-=?++-=-= -∞==∞=→)(lim )(0 s sE e e s ss 3-28 单位反馈系统的开环传递函数为 ) 5(25 )(+= s s s G (1) 求各静态误差系数和25.021)(t t t r ++=时的稳态误差ss e ; (2) 当输入作用10s 时的动态误差是多少? 解 (1))5(25 )(+= s s s G ???==1 5v K ∞=+==→→) 5(25 lim )(lim 0 s s s G K s s p 55 25 lim )(lim 00=+==→→s s G s K s s v 05 25lim )(lim 020=+==→→s s s G s K s s a )(1)(1t t r =时, 011 1=+= p ss K e t t r 2)(2=时, 4.052 2=== v ss K A e 235.0)(t t r =时,∞=== 1 3a ss K A e 由叠加原理 ∞=++=321ss ss ss ss e e e e (2) 题意有 25 5) 5()(11)(2+++=+= Φs s s s s G s e 用长除法可得 ++=++++=Φ3 3 32 210008.02.0)(s s s C s C s C C s e 第三章 线性系统的时域分析法 3.1 引言 分析控制系统的第一步是建立模型,数学模型一旦建立,第二步 分析控制性能,分析有多种方法,主要有时域分析法,频域分析法,根轨迹法等。每种方法,各有千秋。均有他们的适用范围和对象。本章先讨论时域法。 实际上,控制系统的输入信号常常是不知的,而是随机的。很难用解析的方法表示。只有在一些特殊的情况下是预先知道的,可以用解析的方法或者曲线表示。例如,切削机床的自动控制的例子。 在分析和设计控制系统时,对各种控制系统性能得有评判、比较的依据。这个依据也许可以通过对这些系统加上各种输入信号比较它们对特定的输入信号的响应来建立。 许多设计准则就建立在这些信号的基础上,或者建立在系统对初始条件变化(无任何试验信号)的基础上,因为系统对典型试验信号的响应特性,与系统对实际输入信号的响应特性之间,存在着一定的关系;所以采用试验信号来评价系统性能是合理的。 3.1.1 典型试验信号 经常采用的试验输入信号: ① 实际系统的输入信号不可知性; ② 典型试验信号的响应与系统的实际响应,存在某种关系; ③ 电压试验信号是时间的简单函数,便于分析。 突然受到恒定输入作用或突然的扰动。如果控制系统的输入量是随时间逐步变化的函数,则斜坡时间函数是比较合适的。 (单位)阶跃函数(Step function ) 0,)(1≥t t 室温调节系统和水位调节系统 (单位)斜坡函数(Ramp function ) 速度 0,≥t t ∝ (单位)加速度函数(Acceleration function )抛物线 0,2 12 ≥t t (单位)脉冲函数(Impulse function ) 0,)(=t t δ 正弦函数(Simusoidal function )Asinut ,当输入作用具有周期性变化时。 通常运用阶跃函数作为典型输入作用信号,这样可在一个统一的基础上对各种控制系统的特性进行比较和研究。本章讨论系统非周期信号(Step 、Ramp 、对正弦试验信号相应,将在第五章频域分析法,第六章校正方法中讨论)作用下系统的响应。 3.1.2 动态过程和稳态过程 电路分析基础练习题 @ 复刻回忆 1-1 在图题1-1 所示电路中。元件 A 吸收功率30W ,元件 B 吸收功率15W ,元件 C 产生功率30W ,分别求出三个元件中的电流I 1、I 2、I 3。 5V A I 15V B I 2 5V C I 3 图题1-1 解I 1 6 A,I 2 3 A ,I 3 6 A 1-5 在图题1-5 所示电路中,求电流I 和电压U AB 。 解I 4 1 2 1 A,U AB 3 10 2 4 4 39 V 1-6 在图题1-6 所示电路中,求电压U。 30V 5 U 2A 50V 1 I 2 5V 3 I 1 2 4V 图题1-6 图题1-7解50 30 5 2 U ,即有U30V 1-8 在图题1-8 所示电路中,求各元件的功率。 解电阻功率:P 3 P22 2 3 42 / 2 12 W, 3 8 W 2A 电流源功率:P2A 2(10 4 6) 0 ,2 P1 A 4 1 4 W 10V 4V 1A 电压源功率:P 10V 10 2 20 W, P 4V 4(1 2 2) 4 W 2-7 电路如图题2-7 所示。求电路中的未知量。 解U S 2 6 I 12 4 A 2 9 3 12 V I 0 2A I 2 I 3 I 3 P3/ U S12 / 12 1 A U S R eq 6 9 R3 I 0 2 R 4 / 3 1 12 12 13 / 3 A P312W 1 U S R eq I 12 36 13/ 3 13 图题2-7 2-9 电路如图题2-9 所示。求电路中的电流解从图中可知, 2 与3 并联,I 1 。 1 2 由分流公式,得 I 2 I 33 5 I1 5 1 1 A 1 3I 1 I 3 I 2 1V I 1 5I 1 3 所以,有 I 1 I 2I 3 3I 1 1 图题2-9 解得I 1 0.5 A 2-8 电路如图题2-8 所示。已知I1 3I 2 ,求电路中的电阻R 。 解KCL :I 1 I 2 60 I1 3I 260mA I 1 2.2k 解得I 1 R 为45 mA, I 215 mA. I 2R R 2.2 45 15 6.6 k图题2-8 解(a) 由于有短路线, (b) 等效电阻为 R AB 6 , R AB 1// 1 (1 1// 1) // 1 0.5 1.5 2.5 1.1 2-12 电路如图题2-12 所示。求电路AB 间的等效电阻R AB 。 3 第三章线性系统的时域分析方法 教学目的:通过本章学习,熟悉控制系统动态性能指标定义,掌握线性系统稳定的充要条件和劳斯判椐的应用,以及稳态误差计算方法,掌握一阶、 二阶系统的时域分析方法。 教学重点:掌握系统的动态性能指标,能熟练地应用劳斯判椐判断系统稳定性,二阶系统的动态响应特性分析。 教学难点:高阶系统的的动态响应特性分析。 本章知识结构图: 系统结构图闭环传递函数 一阶标准式 二阶标准式 特征方程稳定性、稳定域 代数判据 误差传递函数误差象函数终值定理稳态误差开环传递函数系统型别、开环增益 公式 静态误差系数 第九讲 3.1 系统时间响应的性能指标 一、基本概念 1、时域分析方法:根据系统的数学模型求出系统的时间响应来直接分析和评价系统的方法。 (1)响应函数分析方法:建立数学模型→确定输入信号→求出输出响应→ 根据输出响应→系统分析。 (2)系统测试分析方法:系统加入扰动信号→测试输出变化曲线→系统分析。 系统举例分析:举例:原料气加热炉闭环控制系统 2、分析系统的三大要点 (1)动态性能(快、稳) (2)稳态性能(准) (3)稳定性(稳) 二、动态性能及稳态性能 1、动态过程(过渡过程):在 典型信号作用下,系统输出从初始状态到最终状态的响应过程。(衰减、发散、等幅振荡) 2、稳态过程:在典型信号作 用下,当t → ∞ 系统输出量表现的方式。表征输出量最终复现输入量的程度。(稳态误差描述) 3、动态稳态性能指标 图3-1温度控制系统原理图 (1)上升时间tr :从稳态值的10%上升到稳态值的90%所需要的时间。 (2)峰值时间tp :从零时刻到达第一个峰值h(tp)所用的时间。 (3)超调量δ%:最大峰值与稳态值的差与稳态值之比的百分数。(稳) (3-1) %100)(()(%?∞∞-= h h t h p ) δ 线性系统时域分析 理论基础 求解零状态响应 1 2 ?→0 =-∞ 连续时间信号 f (t ) 和 f (t ) 的卷积运算可用信号的分段求和来实现,即: ∞ ∞ f (t ) = f 1 (t )* f 2 (t ) = ?-∞ f 1 (τ ) f 2 (t -τ )d τ = lim ∑ f 1 (k ?) f 2 (t - k ?) ? ? k 如果只求当t = n ?(n 为整数)时 f (t ) 的值 f (n ?) ,则上式可得: ∞ ∞ f (n ?) = ∑ f 1 (k ?) f 2 (n ? - k ?) ? ? = ?∑ f 1 (k ?) f 2[(n - k )?] (2-1) k =-∞ ∞ k =-∞ 式(2-1)中的 ∑ f 1 (k ?) f 2[(n - k )?] 实际上就是连续时间信号 f 1 (t ) 和 f 2 (t ) 经等时间间隔? k =-∞ 均匀抽样的离散序列 f 1 (k ?) 和 f 2 (k ?) 的卷积和。当? 足够小时, f (n ?) 就是卷积积分的结果——连续时间信号 f (t ) 的较好数值近似。 因此,用 MA TL A B 实现连续信号 f 1 (t ) 和 f 2 (t ) 卷积的过程如下: 1、将连续信号 f 1 (t ) 和 f 2 (t ) 以时间间隔? 进行取样,得到离散序列 f 1 (k ?) 和 f 2 (k ?) ; 2、构造与 f 1 (k ?) 和 f 2 (k ?) 相应的时间向量k 1 和k 2(注意,k 1 和k 2 的元素不是整数,而是取样间隔? 的整数倍的时间间隔点); 3、调用 MATLAB 命令 conv()函数计算积分 f (t ) 的近似向量 f (n ?) ; 4、构造 f (n ?) 对应的时间向量k 。 第二章 连续系统的时域分析习题解答 2-1 图题2-1所示各电路中,激励为f (t ),响应为i 0(t )和u 0(t )。试列写各响 解: . 1)p ( ; )1(1)p ( , 111 , 1 111)( )b (; 105.7)625(3 102 ; )(375)()6253(4) ()()61002.041( )a (0202200 204006000f i p f p u p f p p p u i f p p p p p f t u pf i p pu i t f t u p t f t u p =+++=++?++=+=+++= ++= ?=+??==+?=++-- 2-2 求图题2-1各电路中响应i 0(t )和u 0(t )对激励f (t )的传输算子H (p )。 解:. 1 )()()( ; 11)()()( )b (; 625 3105.7)()()( ; 6253375)()()( )a (22 0 20 40 0 +++==+++== +?==+== -p p p p t f t i p H p p p t f t u p H p p t f t i p H p t f t u p H f i f u f i f u 2-3 给定如下传输算子H (p ),试写出它们对应的微分方程。 . ) 2)(1() 3()( )4( ; 323)( )3(; 3 3)( )2( ; 3)( )1( +++=++=++=+= p p p p p H p p p H p p p H p p p H 解:; 3d d 3d d )2( ; d d 3d d )1( f t f y t y t f y t y +=+=+ . d d 3d d 2d d 3d d )4( ; 3d d 3d d 2 )3( 2222t f t f y t y t y f t f y t y +=+++=+ 2-4 已知连续系统的输入输出算子方程及0– 初始条件为: . 4)(0y ,0)(0y )y(0 ),()2(1 3)( )3(; 0)(0y ,1)(0y ,0)y(0 ),()84() 12()( )2(; 1)(0y ,2)y(0 ),()3)(1(4 2)( )1(---2 ---2 --=''='=++==''='=+++-=='=+++= t f p p p t y t f p p p p t y t f p p p t y 1 f u 0(t ) (b) 1 f (t ) 4k 6k 2F u 0(t ) (a) 图题2-1 第三章 线性系统的时域分析与校正 习题及答案 3-1 已知系统脉冲响应 t e t k 25.10125.0)(-= 试求系统闭环传递函数)(s Φ。 解 Φ()()./(.)s L k t s ==+00125125 3-2 设某高阶系统可用下列一阶微分方程 T c t c t r t r t ?? +=+()()()()τ 近似描述,其中,1)(0<-<τT 。试证系统的动态性能指标为 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 从1.0到9.0所需时间) 当 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 =-==2109 01 22ln ... 3) 求 t s T t s s e T T t h /195.0)(---==τ ]ln 3[]20ln [ln ]05.0ln [ln T T T T T T T T T t s τ ττ-+=+-=--=∴ 3-3 一阶系统结构图如图3-45所示。要求系统闭环增益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-46(a )和(b )分别为开环和闭环温度控制系统结构图,两种系统正常的K 值为1。 第三章时域分析法 3-1两相交流电动机,使用在简单位置控制系统中见习题3-1图。假设作误差检测器用的差动放大器增益为10,且它供电给控制磁场。 ω和阻尼系数ζ等于什么? 试问a)无阻尼自然频率n b)相对超调量和由单位阶跃输入引起峰值的时间等于什么? c)写出关于单位阶跃输入下的误差时间函数。 习题3-1图 3-2 差动放大器的增益增加至20,重做习题3-1。并问从你的结果中能得出什么结论? 3-3 两相交流感应电动机采用齿轮传动和负载链接,使用在简单位置系统中,见习题3-1图。假设作误差检测器用的差动放大器增益为20,且由它供电给控制磁场。试问 ω和阻尼系数ζ等于什么? a)无阻尼自然频率n b)相对超调量和由单位阶跃输入下的峰值的时间等于什么? c)写出关于单位阶跃输入下的误差时间函数。 3-4 差动放大器的增益增至40,重做习题3-3.并问从你的结果中能得出什么结论? 3-5 差动放大器的增益减至10,重做习题3-3.并从你的结果中可得出什么结论? 3-6 综合典型有翼可控导弹控制系统,可使用转矩作用于导弹弹体的方法。这些转矩由作用在离重心很远的控制翼面的偏斜来产生。这样做的结果可以用相对小的翼面负载,就能引起较大的转矩。对这一类型控制系统的设计,为使输入指令响应时间最小,就要求控制回路具有高增益。而又必须限制增益在不引起高频不稳定范围内。习题3-6图表示导弹加速度控制操纵系统。给定加速度与加速度计输出量比较,发出驱动控制系统的节本误差信号。由速度陀螺仪的输出作为阻尼。 试求出下列各式: a)确定这一系统的传递函数C(s)/R(s)。 b)对应一下的一组参数: 放大器增益= A k=16,飞行器增益系数=q=4,R k=4, ω和阻尼系数ζ。 确定该系统无阻尼自然频率n c)确定相对超调量和从加速度单位阶跃输入指令所引起的峰值时间。 实验一线性控制系统时域分析 1、设控制系统如图1 所示,已知K=100,试绘制当H分别取H=0.1 ,0.2 0.5,1, 2,5,10 时,系统的阶跃响应曲线。讨论反馈强度对一阶 系统性能有何影响? 图1 答: A、绘制系统曲线程序如下: s=tf('s'); p1=(1/(0.1*s+1)); p2=(1/(0.05*s+1)); p3=(1/(0.02*s+1)); p4=(1/(0.01*s+1)); p5=(1/(0.005*s+1)); p6=(1/(0.002*s+1)); p7=(1/(0.001*s+1)); step(p1);hold on; step(p2);hold on; step(p3);hold on; step(p5);hold on; step(p6);hold on; step(p7);hold on; B 、绘制改变H 系统阶跃响应图如下: 00.050.10.150.20.250.30.350.40.450.5 0.2 0.4 0.6 0.8 1 1.2 1.4 Step Response Time (seconds) A m p l i t u d e 结论: H 的值依次为0.1、0.2、0.5、1、2、5、10做响应曲线。matlab 曲线默认从第一条到第七条颜色依次为蓝、黄、紫、绿、红、青、黑,图中可知随着H 值得增大系统上升时间减小,调整时间减小,有更高的快速性。 2、 二阶系统闭环传函的标准形式为 22 2()2n n n s s s ωψξωω=++,设已知 n ω=4,试绘制当阻尼比ξ分别取0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 5 等值时,系统的单位阶跃响应曲线。求出ξ取值 0.2 ,0.5 ,0.8时的超调量,并求出ξ取值 0.2 ,0.5 ,0.8,1.5,5时的调节时间。讨论阻尼比变化对系统性能的影响。 — P2-1 — 第二章 连续系统的时域分析习题解答 2-1 图题2-1所示各电路中,激励为f (t ),响应为i 0(t )和u 0(t )。试列写各响应关于激励微分算子方程。 解: . 1)p ( ; )1(1)p ( , 111 , 1 111)( )b (; 105.7)625(3 102 ; )(375)()6253(4) ()()61002.041( )a (0202200 204006000f i p f p u p f p p p u i f p p p p p f t u pf i p pu i t f t u p t f t u p =+++=++?++=+=+++= ++= ?=+??==+?=++-- 2-2 求图题2-1各电路中响应i 0(t )和u 0(t )对激励f (t )的传输算子H (p )。 解:. 1 )()()( ; 11)()()( )b (; 625 3105.7)()()( ; 6253375)()()( )a (22 0 20 40 0 +++==+++== +?==+== -p p p p t f t i p H p p p t f t u p H p p t f t i p H p t f t u p H f i f u f i f u 2-3 给定如下传输算子H (p ),试写出它们对应的微分方程。 . ) 2)(1() 3()( )4( ; 323)( )3(; 3 3)( )2( ; 3)( )1( +++=++=++=+= p p p p p H p p p H p p p H p p p H 解:; 3d d 3d d )2( ; d d 3d d )1( f t f y t y t f y t y +=+=+ . d d 3d d 2d d 3d d )4( ; 3d d 3d d 2 )3( 2222t f t f y t y t y f t f y t y +=+++=+ 2-4 已知连续系统的输入输出算子方程及0– 初始条件为: . 4)(0y ,0)(0y )y(0 ),()2(1 3)( )3(; 0)(0y ,1)(0y ,0)y(0 ),()84() 12()( )2(; 1)(0y ,2)y(0 ),()3)(1(4 2)( )1(---2 ---2 --=''='=++==''='=+++-=='=+++= t f p p p t y t f p p p p t y t f p p p t y f (u 0(t ) (b) u 0(t ) (a) 图题2-1 第二章 连续系统的时域分析习题解答 2-1 图题2-1所示各电路中,激励为f (t ),响应为i 0(t )和u 0(t )。试列写各响应关于 激励微分算子方程。 解: . 1)p ( ; )1(1)p ( , 111 , 1 111)( )b (; 105.7)625(3 102 ; )(375)()6253(4) ()()61002.041( )a (0202200 204006000f i p f p u p f p p p u i f p p p p p f t u pf i p pu i t f t u p t f t u p =+++=++?++=+=+++= ++= ?=+??==+?=++-- 2-2 求图题2-1各电路中响应i 0(t )和u 0(t )对激励f (t )的传输算子H (p )。 } 解:. 1 )()()( ; 11)()()( )b (; 625 3105.7)()()( ; 6253375)()()( )a (220 20 40 0 +++==+++==+?==+== -p p p p t f t i p H p p p t f t u p H p p t f t i p H p t f t u p H f i f u f i f u 2-3 给定如下传输算子H (p ),试写出它们对应的微分方程。 . ) 2)(1() 3()( )4( ; 323)( )3(; 3 3)( )2( ; 3)( )1( +++=++=++=+= p p p p p H p p p H p p p H p p p H 解:; 3d d 3d d )2( ; d d 3d d )1( f t f y t y t f y t y +=+=+ . d d 3d d 2d d 3d d )4( ; 3d d 3d d 2 )3( 2222t f t f y t y t y f t f y t y +=+++=+ 2-4 已知连续系统的输入输出算子方程及0– 初始条件为: . 4)(0y ,0)(0y )y(0 ),()2(1 3)( )3(; 0)(0y ,1)(0y ,0)y(0 ),()84() 12()( )2(; 1)(0y ,2)y(0 ),()3)(1(4 2)( )1(---2 ---2 --=''='=++==''='=+++-=='=+++= t f p p p t y t f p p p p t y t f p p p t y 1 f (u 0(t ) (b) @ f (t ) 4k 6k 2F } u 0(t ) (a) 图题2-1 实验二 线性系统时域响应分析 一、实验目的 1.熟练掌握step( )函数和impulse( )函数的使用方法,研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2.通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 二、基础知识及MATLAB 函数 (一)基础知识 时域分析法直接在时间域中对系统进行分析,可以提供系统时间响应的全部信息,具有直观、准确的特点。为了研究控制系统的时域特性,经常采用瞬态响应(如阶跃响应、脉冲响应和斜坡响应)。本次实验从分析系统的性能指标出发,给出了在MATLAB 环境下获取系统时域响应和分析系统的动态性能和稳态性能的方法。 用MATLAB 求系统的瞬态响应时,将传递函数的分子、分母多项式的系数分别以s 的降幂排列写为两个数组num 、den 。由于控制系统分子的阶次m 一般小于其分母的阶次n ,所以num 中的数组元素与分子多项式系数之间自右向左逐次对齐,不足部分用零补齐,缺项系数也用零补上。 1.用MATLAB 求控制系统的瞬态响应 1)阶跃响应 求系统阶跃响应的指令有: step(num,den) 时间向量t 的范围由软件自动设定,阶跃响应曲线 随即绘出 step(num,den,t) 时间向量t 的范围可以由人工给定(例如 t=0:0.1:10) [y ,x]=step(num,den) 返回变量y 为输出向量,x 为状态向量 在MATLAB 程序中,先定义num,den 数组,并调用上述指令,即可生成单位阶跃输入信号下的阶跃响应曲线图。 考虑下列系统: 25 425 )()(2++=s s s R s C 该系统可以表示为两个数组,每一个数组由相应的多项式系数组成,并且以s 第三章 线性系统的时域分析与校正习题及答案 3-1 已知系统脉冲响应t 25.1e 0125.0)t (k -=,试求系统闭环传递函数)s (Φ。 解 [])25.1s /(0125.0)t (k L )s (+==Φ 3-2 设某高阶系统可用下列一阶微分方程)t (r )t (r )t (c )t (c T +τ=+? ? 近似描述,其中,1)T (0<τ-<。试求系统的动态性能指标s r d t ,t ,t 。 解 设单位阶跃输入s s R 1)(= 当初始条件为0时有: 1 Ts 1 s )s (R )s (C ++τ= 1Ts T s 1s 11Ts 1s )s (C +τ--=?++τ= ∴ T /t e T T 1)t (h )t (c -τ--== T )0(h τ=,1)(h =∞,20T T )]0(h )(h [05.0τ -=-∞=? 1) 当 d t t = 时 2T T e T T 1)]0(h )(h [5.0)0(h )t (h t /t d τ += τ--=-∞+=- T /t d e 2 1 -= ; 693T .0t d = 2) 求r t (即)t (c 从1.0)(h ∞到9.0)(h ∞所需时间) 当T /t 2e T T 1)0(h )]0(h )(h [9.0)t (h -τ-- =+-∞=; 当T /t 1e T T 1)0(h )]0(h )(h [1.0)t (h -τ--=+-∞=; )T 1(.0T ln T t 2τ+τ-=, τ +τ -=)T 9(.0T ln T t 1 则 2T .29ln T t t t 12r ==-= 3) 求 s t T /t s s e T T 1)0(h )]0(h )(h [95.0)t (h -τ-- =+-∞= 3T 05.ln0T t s ==∴ 3-3 一阶系统结构如图所示。要求系统闭环增益2k =Φ,调节时间4.0t s ≤s ,试确定参数21k ,k 的值。 解 由结构图写出闭环系统传递函数 1k k s k 1k k s k s k k 1s k )s (212211211 +=+=+ =Φ 第3章 线性系统的时域分析 学习要点 1控制系统时域响应的基本概念,典型输入信号及意义; 2控制系统稳定性的概念、代数稳定判据及应用; 3控制系统的时域指标,一阶二阶系统的阶跃响应特性与时域指标计算; 4高阶系统时域分析中主导极点和主导极点法; 5 控制系统稳态误差概念、计算方法与误差系数,减小稳态误差的方法。 思考与习题祥解 题 思考与总结下述问题。 (1)画出二阶系统特征根在复平面上分布的几种情况,归纳ξ值对二阶系统特征根的影响规律。 (2)总结ξ和n ω对二阶系统阶跃响应特性的影响规律。 (3)总结增加一个零点对二阶系统阶跃响应特性的影响规律。 (4)分析增加一个极点可能对二阶系统阶跃响应特性有何影响 (5)系统误差与哪些因素有关试归纳减小或消除系统稳态误差的措施与方法。 (6)为减小或消除系统扰动误差,可采取在系统开环传递函数中增加积分环节的措施。请问,该积分环节应在系统结构图中如何配置,抗扰效果是否与扰动点相关 答:(1)二阶系统特征根在复平面上分布情况如图所示。 图 二阶系统特征根在复平面上的分布 当0ξ=,二阶系统特征根是一对共轭纯虚根,如图中情况①。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,变化轨迹是 以n ω为半径的圆弧,如图中情况②。 当1ξ=,二阶系统特征根是一对相同的负实根,如图中情况③。 当1ξ>,二阶系统特征根是一对不等的负实根,如图中情况④。 (2)ξ和n ω是二阶系统的两个特征参量。 ξ是系统阻尼比,描述了系统的平稳性。 当0ξ=,二阶系统特征根是一对共轭纯虚根,二阶系统阶跃响应为等幅振荡特性,系统临界稳定。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,二阶系统阶跃响应为衰减振荡特性,系统稳定。ξ越小,二阶系统振荡性越强,平稳性越差; ξ越大,二阶系统振荡性越弱,平稳性越好。因此,二阶系统的时域性能指标超 调量由ξ值唯一确定,即001_ 100%2 ?=-π ξξ σe 。在工程设计中,对于恒值控制系 统,一般取 ξ=~;对于随动控制系统ξ=~。 n ω是系统无阻尼自然振荡频率,反映系统的快速性。当ξ一定,二阶系统的 时域性能指标调节时间与n ω值成反比,即34 s n t ξω≈:。 (3)二阶系统增加一个零点后,增加了系统的振荡性,将使系统阶跃响应的超调量增大,上升时间和峰值时间减小。 所增加的零点越靠近虚轴,则上述影响就越大;反之,若零点距离虚轴越远,则其影响越小。 (4)二阶系统增加一个极点后,减弱了系统的振荡性,将使系统阶跃响应的超调量减小,上升时间和峰值时间减小; 所增加的极点越靠近虚轴,则上述影响就越大;反之,若极点距离虚轴越远,则其影响越小。 (5)系统误差与系统的误差度(开环传递函数所含纯积分环节的个数或系统型别)、开环放大系数,以及作用于系统的外部输入信号有关。如果是扰动误差还与扰动作用点有关。 因此,减小或消除系统稳态误差的措施与方法有:增大开环放大系数,增加系统开环传递函数中的积分环节,引入按给定或按扰动补偿的复合控制结构。 无论采用何种措施与方法减小或消除系统稳态误差,都要注意系统须满足稳定的条件。 (6)采取在系统开环传递函数中增加积分环节的措施来减小或消除系统扰动误差时,所增加的积分环节须加在扰动作用点之前。若所增加的积分环节加在扰动作用点之后,则该积分环节无改善抗扰效果作用。这一点可以通过误差表达式分析得到。 题系统特征方程如下,试判断其稳定性。 (a )0203.002.023=+++s s s ; (b )014844122345=+++++s s s s s ; (c )025266.225.11.0234=++++s s s s 解:(a )稳定; (b )稳定; (c )不稳定。 题 系统结构如题图所示。控制器)1 1()(s T K s G i p c + =,为使该系统稳定,控制器参数p K 、i T 应满足什么关系 第3章 线性系统的时域分析 学习要点 1控制系统时域响应的基本概念,典型输入信号及意义; 2控制系统稳定性的概念、代数稳定判据及应用; 3控制系统的时域指标,一阶二阶系统的阶跃响应特性与时域指标计算; 4高阶系统时域分析中主导极点和主导极点法; 5 控制系统稳态误差概念、计算方法与误差系数,减小稳态误差的方法。 思考与习题祥解 题 思考与总结下述问题。 (1)画出二阶系统特征根在复平面上分布的几种情况,归纳ξ值对二阶系统特征根的影响规律。 【 (2)总结ξ和n ω对二阶系统阶跃响应特性的影响规律。 (3)总结增加一个零点对二阶系统阶跃响应特性的影响规律。 (4)分析增加一个极点可能对二阶系统阶跃响应特性有何影响 (5)系统误差与哪些因素有关试归纳减小或消除系统稳态误差的措施与方法。 (6)为减小或消除系统扰动误差,可采取在系统开环传递函数中增加积分环节的措施。请问,该积分环节应在系统结构图中如何配置,抗扰效果是否与扰动点相关 答:(1)二阶系统特征根在复平面上分布情况如图所示。 图 二阶系统特征根在复平面上的分布 当0ξ=,二阶系统特征根是一对共轭纯虚根,如图中情况①。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,变化轨迹是 以n ω为半径的圆弧,如图中情况②。 @ 当1ξ=,二阶系统特征根是一对相同的负实根,如图中情况③。 当1ξ>,二阶系统特征根是一对不等的负实根,如图中情况④。 (2)ξ和n ω是二阶系统的两个特征参量。 ξ是系统阻尼比,描述了系统的平稳性。 当0ξ=,二阶系统特征根是一对共轭纯虚根,二阶系统阶跃响应为等幅振荡特性,系统临界稳定。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,二阶系统阶跃响应为衰减振荡特性,系统稳定。ξ越小,二阶系统振荡性越强,平稳性越差; ξ越大,二阶系统振荡性越弱,平稳性越好。因此,二阶系统的时域性能指标超 调量由ξ值唯一确定,即001_ 100%2 ?=-π ξξ σe 。在工程设计中,对于恒值控制系 统,一般取 ξ=~;对于随动控制系统ξ=~。 n ω是系统无阻尼自然振荡频率,反映系统的快速性。当ξ一定,二阶系统的 时域性能指标调节时间与n ω值成反比,即34 s n t ξω≈。 (3)二阶系统增加一个零点后,增加了系统的振荡性,将使系统阶跃响应的超调量增大,上升时间和峰值时间减小。 所增加的零点越靠近虚轴,则上述影响就越大;反之,若零点距离虚轴越远,则其影响越小。 (4)二阶系统增加一个极点后,减弱了系统的振荡性,将使系统阶跃响应的超调量减小,上升时间和峰值时间减小; 所增加的极点越靠近虚轴,则上述影响就越大;反之,若极点距离虚轴越远,则其影响越小。 & (5)系统误差与系统的误差度(开环传递函数所含纯积分环节的个数或系统型别)、开环放大系数,以及作用于系统的外部输入信号有关。如果是扰动误差还与扰动作用点有关。 因此,减小或消除系统稳态误差的措施与方法有:增大开环放大系数,增加系统开环传递函数中的积分环节,引入按给定或按扰动补偿的复合控制结构。 无论采用何种措施与方法减小或消除系统稳态误差,都要注意系统须满足稳定的条件。 (6)采取在系统开环传递函数中增加积分环节的措施来减小或消除系统扰动误差时,所增加的积分环节须加在扰动作用点之前。若所增加的积分环节加在扰动作用点之后,则该积分环节无改善抗扰效果作用。这一点可以通过误差表达式分析得到。 题系统特征方程如下,试判断其稳定性。 (a )0203.002.023=+++s s s ; (b )014844122345=+++++s s s s s ; (c )025266.225.11.0234=++++s s s s ! 解:(a )稳定; (b )稳定; (c )不稳定。 第3章 线性系统的时域分析 3.1 学习要点 1控制系统时域响应的基本概念,典型输入信号及意义; 2控制系统稳定性的概念、代数稳定判据及应用; 3控制系统的时域指标,一阶二阶系统的阶跃响应特性与时域指标计算; 4高阶系统时域分析中主导极点和主导极点法; 5 控制系统稳态误差概念、计算方法与误差系数,减小稳态误差的方法。 3.2 思考与习题祥解 题3.1 思考与总结下述问题。 (1)画出二阶系统特征根在复平面上分布的几种情况,归纳ξ值对二阶系统特征根的影响规律。 (2)总结ξ和n ω对二阶系统阶跃响应特性的影响规律。 (3)总结增加一个零点对二阶系统阶跃响应特性的影响规律。 (4)分析增加一个极点可能对二阶系统阶跃响应特性有何影响? (5)系统误差与哪些因素有关?试归纳减小或消除系统稳态误差的措施与方法。 (6)为减小或消除系统扰动误差,可采取在系统开环传递函数中增加积分环节的措施。请问,该积分环节应在系统结构图中如何配置,抗扰效果是否与扰动点相关? 答:(1)二阶系统特征根在复平面上分布情况如图3.1所示。 图3.1 二阶系统特征根在复平面上的分布 当0ξ=,二阶系统特征根是一对共轭纯虚根,如图中情况①。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,变化轨迹是 以n ω为半径的圆弧,如图中情况②。 当1ξ=,二阶系统特征根是一对相同的负实根,如图中情况③。 当1ξ>,二阶系统特征根是一对不等的负实根,如图中情况④。 (2)ξ和n ω是二阶系统的两个特征参量。 ξ是系统阻尼比,描述了系统的平稳性。 当0ξ=,二阶系统特征根是一对共轭纯虚根,二阶系统阶跃响应为等幅振荡特性,系统临界稳定。 当01ξ<<,二阶系统特征根是一对具有负实部的共轭复数根,二阶系统阶跃响应为衰减振荡特性,系统稳定。ξ越小,二阶系统振荡性越强,平稳性越差; ξ越大,二阶系统振荡性越弱,平稳性越好。因此,二阶系统的时域性能指标超 调量由ξ值唯一确定,即001_ 100%2 ?=-π ξξσe 。在工程设计中,对于恒值控制系 统,一般取 ξ=0.2~0.4;对于随动控制系统ξ=0.6~0.8。 n ω是系统无阻尼自然振荡频率,反映系统的快速性。当ξ一定,二阶系统的 时域性能指标调节时间与n ω值成反比,即34 s n t ξω≈ 。 (3)二阶系统增加一个零点后,增加了系统的振荡性,将使系统阶跃响应的超调量增大,上升时间和峰值时间减小。 所增加的零点越靠近虚轴,则上述影响就越大;反之,若零点距离虚轴越远,则其影响越小。 (4)二阶系统增加一个极点后,减弱了系统的振荡性,将使系统阶跃响应的超调量减小,上升时间和峰值时间减小; 所增加的极点越靠近虚轴,则上述影响就越大;反之,若极点距离虚轴越远,则其影响越小。 (5)系统误差与系统的误差度(开环传递函数所含纯积分环节的个数或系统型别)、开环放大系数,以及作用于系统的外部输入信号有关。如果是扰动误差还与扰动作用点有关。 因此,减小或消除系统稳态误差的措施与方法有:增大开环放大系数,增加系统开环传递函数中的积分环节,引入按给定或按扰动补偿的复合控制结构。 无论采用何种措施与方法减小或消除系统稳态误差,都要注意系统须满足稳定的条件。 (6)采取在系统开环传递函数中增加积分环节的措施来减小或消除系统扰动误差时,所增加的积分环节须加在扰动作用点之前。若所增加的积分环节加在扰动作用点之后,则该积分环节无改善抗扰效果作用。这一点可以通过误差表达式分析得到。 题3.2系统特征方程如下,试判断其稳定性。 (a )0203.002.023=+++s s s ; (b )014844122345=+++++s s s s s ; (c )025266.225.11.0234=++++s s s s 解:(a )稳定; (b )稳定; (c )不稳定。 题3.3 系统结构如题3.3图所示。控制器)1 1()(s T K s G i p c + =,为使该系统稳定,控制器参数p K 、i T 应满足什么关系? 第三章 时域分析法习题 3-1设温度计需要在一分钟内指示出响应值的98%,并且假定温度计为一阶系统,试 求时间常数T 。如果将温度计放在澡盆内,澡盆的温度以min /C 10?的速度线性变化,求温度计的误差。 解:min 25.0min 14=?=T T 传递函数1 25.01)(+= Φs s t t r 10)(= 4 25.025.01) 125.0(10)()()(2 2 ++ - =+=?Φ=s s s s s s R s s C t e t t c 45.25.210)(-+-=C t c t r t e e s s ss ?=-==∞ →∞ →5.2)()(lim )(lim 3-4 单位负反馈系统的开环传递函数为) 1(1)(+= s s s G ,求该系统的上升时间r t 、峰 值时间p t 、超调量%σ和调整时间s t 。 解:1 1)(2 ++= Φs s s 12=n ξω,12 =n ω.3 arccos ,5.0,1π ξβξω= ===?n 。 %16%100%,63.31,42.212 1/ 2 2 =?==-= =--= --ξ πξσξ ωπξ ωβπe t t n p n r %)2(84 %),5(63 =?== =?== n s n s t t ξωξω 3-6 系统的单位阶跃响应为t t e e t c 10602.12.01)(---+=,试求: (1) 系统的闭环传递函数; (2) 系统的阻尼比ξ和无阻尼自然振荡频率n ω。 解:(1)s s R s s s s C 1)(,10 2.160 2.01)(= +- ++ = 600 70600 )(2 ++= Φs s s 第3章 线性系统的时域分析 本章讨论线性系统的运动分析。主要介绍连续系统与离散系统的状态空间模型的求解、状态转移矩阵的性质和计算以及连续系统状态方程的离散化。本章最后介绍基于Matlab的状态空间模型求解与控制系统的运动仿真问题的程序设计与仿真计算。 建立了系统的数学描述之后,接下来要对系统作定量和定性分析。定量分析主要研究系统对给定输入信号的响应问题,也就是对描述系统的状态方程和输出方程的求解问题。定性分析主要研究系统的结构性质,如能控性、能观性、稳定性等。本章先讨论用状态空间模型描述的线性系统的定量分析问题,即状态空间模型的求解问题。根据常微分方程理论求解一个一阶定系数线性微分方程组是很容易的,可是求解一个一阶变系数线性微分方程组却非易事。状态转移矩阵的引入,使得定常系统和时变系统的求解公式具有一个统一的形式。为此,本章将重点讨论状态转移矩阵的定义、性质和计算方法,并在此基础上导出状态方程的求解公式。本章讨论的另一个中心问题是连续系统状态方程的离散化,即建立连续系统的离散系统状态方程。随着计算机在控制系统分析、设计和实时控制中的广泛应用,这个问题显得越来越重要。在离散系统状态方程建立的基础上,本章也将讨论相应的状态方程求解问题,并将导出在形式上与连续系统状态方程的解一致的离散系统状态方程的解。 3.1 线性定常连续系统状态方程的解 在讨论一般线性定常连续系统状态方程的解之前,我们先讨论线性定常齐次状态方程的解,以便引入矩阵指数函数和状态转移矩阵的概念。所谓齐次状态方程,就是指状态方程中不考虑输入项的作用,满足方程解的齐次性的一类状态方程。研究齐次状态方程的解,就是研究系统本身在无外力作用下的自由运动。 3.2 状态转移矩阵及其计算 在状态方程求解过程中,关键是状态转移矩阵Φ(t)的计算。对于线性定常连续系统,该问题又归结到矩阵指数函数e At的计算。上一节已经介绍了基于拉氏反变换技术的矩阵指数函数e At的计算方法,下面讲述计算矩阵指数函数的其他3种常用方法。 3.2.1级数求和法 第 三章 时域分析法习题及解答 3-1. 假设温度计可用 1 1 +Ts 传递函数描述其特性,现在用温度计测量盛在容器内的水温。发现需要min 1时间才能指示出实际水温的98%的数值,试问该温度计指示出实际水温从10%变化到90%所需的时间是多少? 解: 41min, =0.25min T T = 3-2. 系统在静止平衡状态下,加入输入信号t t t r +=)(1)(,测得响应为 试求系统的传递函数。 解:2210.90.910(s+1)()=10s (s+10) C s s s s =+-+ 3-3. 某惯性环节在单位阶跃作用下各时刻的输出值如下表所示。试求环节的传递函数。 0 1 2 3 4 5 6 7 h (t ) 解: 设()1 s Ts φ=+ 3-4. 已知系统结构图如图3-49所示。试分析参数a 对输出阶跃响应的影响。 解:1()()111 K K Ts s Kas T Ka s Ts φ+==+++ + 当a>0时,系统响应速度变慢; 0T a K -<<时,系统响应速度变快。 3-5. 设控制系统闭环传递函数为 试在[s ]平面上绘出满足下列各要求的系统特征方程式根的可能分布的区域。 1. 707.01>>ξ, 2≥n ω 2.05.0>>ξ, 24≥≥n ω 3. 5.0707.0>>ξ, 2≤n ω 解:①0.707<<1, 2n ξω≥ ②0<0.5, 24n ξω≤≤≤ ③0.50.707, 2n ξω≤≤≤ 3-6. 已知某前向通路的传递函数(如图3-50所示) 今欲采用负反馈的办法将阶跃响应的调节时间s t 减小为原来的1.0倍,并保证总放大系数不变。试选择H K 和0K 的值。 解: 解得:00.9 =10H K K = 3-7. 设一单位反馈控制系统的开环传递函数 为 试分别求出当1 10-=s K 和1 20-=s K 时系统的阻尼比ξ,无阻尼自然频率n ω,单位阶跃响应的超调量%σ及峰值时间p t ,并讨论K 的大小对系统性能指标的影响。 解: 22()10()1()0.11010G s K K s G s s s K s s K φ= ==+++++ K 增大使%,p t σ↑↓,但不影响调节时间。 3-8. 设二阶控制系统的单位阶跃响应曲线如图3-51所示。如果该系统属于单位反馈控制系统,试确定其开环传递函数。 解: 222 ()2n n n s s s ωφξωω=++ 3-9. 设系统闭环传递函数 试求1.2.0=ξ;s T 08.0=;4.0=ξ;s T 08.0=;0.8ξ=;s T 08.0=时单位阶跃响应的超调量%σ、调节时间s t 及峰值时间p t 。 2.4.0=ξ;s T 04.0=和4.0=ξ;s T 16.0=时单位阶跃响应的超调量%σ、调节时间s t 和峰值时间p t 。 3.根据计算结果,讨论参数ξ、T 对阶跃响应的影响。 题解3-5(1) 题解3-5(2) 题 解线性系统的时域分析法(第七讲)
电路分析基础习题集与答案解析
线性系统的时域分析方法
线性系统时域分析
_第二章连续系统的时域分析习题解答
第三章 线性系统的时域分析与校正 习题及答案
第三章 时域分析法 习题
自动控制原理实验报告《线性控制系统时域分析》
第二章__连续系统的时域分析习题解答
_第二章连续系统的时域分析习题解答
自动控制原理》实验2(线性系统时域响应分析
线性系统的时域分析与校正习题及答案
线性系统的时域分析习题答案
第3章线性系统的时域分析习题答案
第3章线性系统的时域分析习题答案
时域分析法习题与解答
线性系统的时域分析
时域分析法习题及解答