第三章 线性系统的时域分析与校正
习题及答案
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。
(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 电压可得1200m in 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
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-48所示。若要求系统单位阶跃响应无超调,且调节时间尽可能短,问开环增益K 应取何值,调节时间s t 是多少?
解 依题意应取 1=ξ,这时可设闭环极点为
02,11-=λ。
写出系统闭环传递函数
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 101102
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
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 ω
将 s t 1=代入二阶系统阶跃响应公式
()
βωξξ
ξω+---
=-t e t h n t n 22
1sin 11)(
可得 m in 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.012100
1221K 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 2K 和a 。
解 由系统阶跃响应曲线有
???
??=-===∞o
o o o
p t h 3.33)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 阶跃电压时(带负载)稳态转速为1000m in 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)(4
5
=--+=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)(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 4243164
12?-?= 48 S 121644812
0?-?= 0 辅助方程 124802s +=,
S 24 辅助方程求导:024=s
S 0 48
系统没有正根。对辅助方程求解,得到系统一对虚根 s j 122,=±。 (3)022)(4
5
=--+=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(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-54是某垂直起降飞机的高度控制系统结构图,试确定使系统稳定的K 值范围。
解 由结构图,系统开环传递函数为:
)
4()
124()(232++++=s s s s s K s G
??
?==3
4
v K K k 系统型别开环增益 0244)(2
3
4
5
=+++++=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 使系统稳定的开环增益范围为: 151815158<= ) 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-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()(2 3=+++=+++=K s s s K s s s s D 4.2636.360 164.269604.2616601:0 123 >→<→-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-56是船舶横摇镇定系统结构图,引入内环速度反馈是为了增加船只的阻尼。 (1) 求海浪扰动力矩对船只倾斜角的传递函数 ) () (s M s N Θ; (2) 为保证N M 为单位阶跃时倾斜角θ的值不超过0.1,且系统的阻尼比为0.5,求2K 、 1K 和3K 应满足的方程; (3) 取2K =1时,确定满足(2)中指标的1K 和3K 值。 解 (1) )5.01()5.02.0(5.01 2.05.012.05.0112.05.0)()(213222 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 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 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 7v 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 ? ??=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-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(00212 21K K s K K T T s T T s K -+-+++= ττ 依题意应有:?? ?=-+=-0 10210τK K T T K K 联立求解得 ?? ?+==2 101T 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 2213 2 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 ?=时, 令 第三章 线性系统的时域分析法 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 动态过程和稳态过程 武汉工程大学实验报告 专业班号 组别01 教师 姓名同组者(个人) 2222)(n n n s s s G ωζωω++= (1)分别绘出)/(2s rad n =ω,ζ分别取0,,,和时的单位阶跃响应曲线,分析参数ζ对系统的影响,并计算ζ=时的时域性能指标ss s p r p e t t t ,,,,σ。 (2)绘制出当ζ=, n ω分别取1,2,4,6时单位阶跃响应曲线,分析参数n ω对系统的影响。 (3)系统的特征方程式为010532234=++++s s s s ,试用二种判稳方式判别该系统的稳定性。 (4)单位负反馈系统的开环模型为 )256)(4)(2()(2++++= s s s s K s G 试分别用劳斯稳定判据和赫尔维茨稳定判据判断系统的稳定性,并求出使得闭环系统稳定的K 值范围。 三、 实验结果及分析 1.可以用两种方法绘制系统的阶跃响应曲线。 (1)用函数step( )绘制 MATLAB 语言程序: >> num=[ 0 0 1 3 7]; >> den=[1 4 6 4 1 ]; >>step(num,den); >> grid; >>xlabel('t/s');ylabel('c(t)');title('step response'); MATLAB运算结果: (2)用函数impulse( )绘制 MATLAB语言程序: >> num=[0 0 0 1 3 7]; >> den=[1 4 6 4 1 0]; >> impulse(num,den); >> grid; >> xlabel('t/s');ylabel('c(t)');title('step response');MATLAB运算结果: 第三章线性系统的时域分析方法 教学目的:通过本章学习,熟悉控制系统动态性能指标定义,掌握线性系统稳定的充要条件和劳斯判椐的应用,以及稳态误差计算方法,掌握一阶、 二阶系统的时域分析方法。 教学重点:掌握系统的动态性能指标,能熟练地应用劳斯判椐判断系统稳定性,二阶系统的动态响应特性分析。 教学难点:高阶系统的的动态响应特性分析。 本章知识结构图: 系统结构图闭环传递函数 一阶标准式 二阶标准式 特征方程稳定性、稳定域 代数判据 误差传递函数误差象函数终值定理稳态误差开环传递函数系统型别、开环增益 公式 静态误差系数 第九讲 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 ) δ 《机电系统控制基础》大作业一 基于MATLAB的机电控制系统响应分析 哈尔滨工业大学 2013年11月4日 1 作业题目 1. 用MATLAB 绘制系统2 ()25()() 425 C s s R s s s Φ== ++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2 ()25 ()()425 C s s R s s s Φ==++的单位阶跃响应性能指标:上升时间、峰值时间、调节时间和超调量。 3. 数控直线运动工作平台位置控制示意图如下: X i 伺服电机原理图如下: L R (1)假定电动机转子轴上的转动惯量为J 1,减速器输出轴上的转动惯量为J 2,减速器减速比为i ,滚珠丝杠的螺距为P ,试计算折算到电机主轴上的总的转动惯量J ; (2)假定工作台质量m ,给定环节的传递函数为K a ,放大环节的传递函数为K b ,包括检测装置在内的反馈环节传递函数为K c ,电动机的反电势常数为K d ,电动机的电磁力矩常数为K m ,试建立该数控直线工作平台的数学模型,画出其控制系统框图; (3)忽略电感L 时,令参数K a =K c =K d =R=J=1,K m =10,P/i =4π,利用MATLAB 分析kb 的取值对于系统的性能的影响。 2 题目1 单位脉冲响应曲线 单位阶跃响应曲线 源代码 t=[0:0.01:1.6]; %仿真时间区段和输入 nC=[25]; dR=[1,4,25]; fi=tf(nC,dR); %求系统模型 [y1,T]=impulse(fi,t); [y2,T]=step(fi,t); %系统响应 plot(T,y1); xlabel('t(sec)'),ylabel('x(t)'); grid on; plot(T,y2); xlabel('t(sec)'),ylabel('x(t)'); grid on; %生成图形 3 题目2 借助Matlab,可得: ans = 0.4330 0.6860 25.3826 1.0000 即 线性系统时域分析 理论基础 求解零状态响应 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 。 武汉工程大学 实验报告 专业 班号 组别 指导教师 姓名 学号 实验名称 线性系统时域响应分析 一、实验目的 1.熟练掌握step( )函数和impulse( )函数的使用方法,研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2.通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 3.熟练掌握系统的稳定性的判断方法。 二、实验内容 1.观察函数step( )和impulse( )的调用格式,假设系统的传递函数模型为 1 4647 3)(2 342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线试分别绘制。 2.对典型二阶系统 2 22 2)(n n n s s s G ωζωω++= 1)分别绘出)/(2s rad n =ω,ζ分别取0,,,和时的单位阶跃响应曲线,分析参数ζ对系统的影响,并计算ζ=时的时域性能指标ss s p r p e t t t ,,,,σ。 2)绘制出当ζ=, n ω分别取1,2,4,6时单位阶跃响应曲线,分析参数n ω对系统的影响。 3.系统的特征方程式为010532234=++++s s s s ,试用两种判稳方式判别该系统的稳定性。 4.单位负反馈系统的开环模型为 ) 256)(4)(2()(2++++= s s s s K s G 试用劳斯稳定判据判断系统的稳定性,并求出使得闭环系统稳定的K 值范围。 三、实验结果及分析 1.观察函数step( )和impulse( )的调用格式,假设系统的传递函数模型为 14647 3)(2342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线试分别绘制。 方法一:用step( )函数绘制系统阶跃响应曲线。 程序如下: num=[0 0 1 3 7]; den=[1 4 6 4 1]; t=0::10; step(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-step Response of G(s)=s^2+3s+7/(s^4+4s^3+6s^2+4s+1)') Unit-step Response of G(s)=s 2+3s+7/(s 4+4s 3+6s 2+4s+1) t/s (sec) c (t ) 方法二:用impulse( )函数绘制系统阶跃响应曲线。 程序如下: num=[0 0 0 1 3 7 ]; den=[1 4 6 4 1 0]; t=0::10; impulse(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-impulse Response of G(s)/s=s^2+3s+7/(s^5+4s^4+6s^3+4s^2+s)') 第三章 控制系统的时域分析法 一、知识点总结 1.掌握典型输入信号(单位脉冲、单位阶跃、单位速度、单位加速度、正弦信号)的拉氏变换表达式。 2.掌握系统动态响应的概念,能够从系统的响应中分离出稳态响应分量和瞬态响应分量;掌握系统动态响应的性能评价指标的概念及计算方法(对于典型二阶系统可以直接应用公式求解,非典型二阶系统则应按定义求解)。 解释:若将系统的响应表达成拉普拉氏变换结果(即S 域表达式),将响应表达式进行部分分式展开,与系统输入信号极点相同的分式对应稳态响应;与传递函数极点相同的分式对应系统的瞬态响应。将稳态响应和瞬态响应分式分别进行拉氏逆变换即获得各自的时域表达式。 性能指标:延迟时间、上升时间、峰值时间、调节时间、超调量 3.掌握一阶系统的传递函数形式,在典型输入信号下的时域响应及其响应特征;掌握典型二阶系统的传递函数形式,掌握欠阻尼系统的阶跃响应时域表达及其性能指标的计算公式和计算方法;了解高阶系统的性能分析方法,熟悉主导极点的概念,定性了解高阶系统非主导极点和零点对系统性能的影响。 tr tp ts td 4.熟悉两种改善二阶系统性能的方法和结构形式(比例微分和测速反馈),了解两种方法改善系统性能的特点。 5.掌握系统稳定性分析方法:劳斯判据的判断系统稳定性的判据及劳斯判据表特殊情况的构建方法(首列元素出现0,首列出现无穷大,某一行全为0);掌握应用劳斯判据解决系统稳定裕度问题的方法。了解赫尔维茨稳定性判据。 6.掌握稳态误差的概念和计算方法;掌握根据系统型别和静态误差系数计算典型输入下的稳态误差的方法(可直接应用公式);了解消除稳态误差和干扰误差的方法;了解动态误差系数法。 二、相关知识点例题 例1. 已知某系统的方块图如下图1所示,若要求系统的性能指标为: δδ%=2222%,tt pp=1111,试确定K和τ的值,并计算系统单位阶跃输入下的特征响应量:tt,tt。 图1 解:系统闭环传递函数为:Φ(s)=CC(ss)RR(ss)=KK ss2+(1+KKKK)ss+KK 因此,ωnn=√KK,ζζ=1+KKKK2√KK, δ%=e?ππππ?1?ππ2?ζζ=0.46, t pp=ππωωdd=1ss?ωdd=ωnn?1?ζζ2=3.14 ?ωnn=3.54 K=ωnn2=12.53,τ=2ζζωnn?1KK=0.18 t ss=3ζζωωnn=1.84ss 实验报告 实验名称 线性系统时域响应分析 一、 实验目的 1.熟练掌握step( )函数和impulse( )函数的使用方法,研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2.通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 3.熟练掌握系统的稳定性的判断方法。 二、 实验内容 1.观察函数step( )和impulse( )的调用格式,假设系统的传递函数模型为 1 4647 3)(2 342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线?试分别绘制。 2.对典型二阶系统 2 22 2)(n n n s s s G ωζωω++= 1)分别绘出)/(2s rad n =ω,ζ分别取0,0.25,0.5,1.0和2.0时的单位阶跃响应曲线,分析参数ζ对系统的影响,并计算ζ=0.25时的时域性能指标 ss s p r p e t t t ,,,,σ。 2)绘制出当ζ=0.25, n ω分别取1,2,4,6时单位阶跃响应曲线,分析参数n ω对系统的影响。 3.系统的特征方程式为010532234=++++s s s s ,试用两种判稳方式判别该系统的稳定性。 4.单位负反馈系统的开环模型为 ) 256)(4)(2()(2 ++++= s s s s K s G 试用劳斯稳定判据判断系统的稳定性,并求出使得闭环系统稳定的K 值范围。 三、 实验结果及分析 1.观察函数step( )和impulse( )的调用格式,假设系统的传递函数模型为 1 4647 3)(2342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线?试分别绘制。 方法一: num=[1 3 7]; den=[1 4 6 4 1]; step(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-step Respinse of G(s)=(s^2+3s+7)/(s^4+4s^3+6s^2+4s+1)') 方法二: num=[1 3 7]; den=[1 4 6 4 1 0]; impulse(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-impulse Respinse of G(s)/s=(s^2+3s+7)/(s^5+4s^4+6s^3+4s^2+s)') 实验报告 实验名称:实验1:控制系统的时域分析 课程名称:自控控制原理 专业:电气工程及其自动化 班级:130037 学生姓名:施苏伟 班级学号:13003723 指导教师:杨杨 实验日期:2015 年10 月16日 一、实验目的 1.观察控制系统的时域响应; 2.记录单位阶跃响应曲线; 3.掌握时间响应分析的一般方法; 4.初步了解控制系统的调节过程。 二.实验步骤: 1.将‘实验一代码’这个文件夹拷贝到桌面上; 2.开机进入Matlab6.1 运行界面(其他版本亦可); 3.通过下面方法将当前路径设置为‘实验一代码’这个文件夹所在的路径 4.Matlab 指令窗>>后面输入指令:con_sys; 进入本次实验主界面。 5.分别双击上图中的三个按键,依次完成实验内容。 6.本次实验的相关Matlab 函数: 传递函数G=tf([num],[den])可输入一传递函数,其中num、den 分别表示分子、分母按降幂排列的系数。 三、仿真结果: (一)观察一阶系统G=1/(T+s)的时域响应: T=5s T=8s T=13s 结果分析:一阶系统 G=1/(T+s)的,通过观察曲线发现,随着时间常数T的增大,同种响应要达到相同响应的时间增大,说明T越大,响应越慢。 (二)二阶系统的时域性能分析 (1) 结果分析:自然频率和阻尼比的适当时,通过调节相应的时间,阶跃响应可以得到稳定值。 (2)数据一:自然频率=5.96rad/sec 阻尼比=0.701 数据二:自然频率=8.2964rad/sec 阻尼比=0.701 结果分析:要达到既定范围,自然频率增大阻尼比要随之增大 (3) 实验一线性控制系统时域分析 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时的调节时间。讨论阻尼比变化对系统性能的影响。 课程名称:控制理论指导老师:成绩: 实验名称:控制系统的时域分析实验类型:冋组学生姓名: 、实验目的和要求 1用计算机辅助分析的办法,掌握系统的时域分析方法。 2. 熟悉SimUlink仿真环境。 二、实验内容和原理 (一)实验原理 系统仿真实质上就是对系统模型的求解,对控制系统来说,一般模型可转化成某个微分方程或差分方程表示,因此在仿真过程中,一般以某种数值算法从初态出发,逐步计算系统的响应,最后绘制出系统的响应曲线,进而可分析系统的性能。控制系统最常用的时域分析方法是,当输入信号为单位阶跃和单位冲激函数时,求出系统的输出响应,分别称为单位阶跃响应和单位冲激响应。在MATLAB中,提供了求取连 续系统的单位阶跃响应函数step,单位冲激响应函数impulse,零输入响应函数initial等等。 (二)实验内容 二阶系统,其状态方程模型为 U X I y = [1.9691 6.4493] +[0] U X2 1?画出系统的单位阶跃响应曲线; 2. 画出系统的冲激响应曲线; 3. 当系统的初始状态为x0=[1,0]时,画出系统的零输入响应; 4. 当系统的初始状态为零时,画出系统斜坡输入响应; (三)实验要求 1. 编制MATLAB程序,画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应; 2. 在SimUIink仿真环境中,组成系统的仿真框图,观察单位阶跃响应曲线并记录之。 三、主要仪器设备 计算机一台以及matlab软件,SimUIink仿真环境 四、操作方法与实验步骤 1、程序解决方案: 在MATLAB 中建立文件shiyu.m ,其程序如下: %时域响应函数 fun ction G1 = shiyu( A,B,C,D) 第三章 线性系统的时域分析与校正习题及答案 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 +=+=+ =Φ 实验二 线性系统时域响应分析 一、实验目的 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 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 时,求扰动对两种系统的温度的影响。 第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章 线性系统的时域分析 本章讨论线性系统的运动分析。主要介绍连续系统与离散系统的状态空间模型的求解、状态转移矩阵的性质和计算以及连续系统状态方程的离散化。本章最后介绍基于Matlab的状态空间模型求解与控制系统的运动仿真问题的程序设计与仿真计算。 建立了系统的数学描述之后,接下来要对系统作定量和定性分析。定量分析主要研究系统对给定输入信号的响应问题,也就是对描述系统的状态方程和输出方程的求解问题。定性分析主要研究系统的结构性质,如能控性、能观性、稳定性等。本章先讨论用状态空间模型描述的线性系统的定量分析问题,即状态空间模型的求解问题。根据常微分方程理论求解一个一阶定系数线性微分方程组是很容易的,可是求解一个一阶变系数线性微分方程组却非易事。状态转移矩阵的引入,使得定常系统和时变系统的求解公式具有一个统一的形式。为此,本章将重点讨论状态转移矩阵的定义、性质和计算方法,并在此基础上导出状态方程的求解公式。本章讨论的另一个中心问题是连续系统状态方程的离散化,即建立连续系统的离散系统状态方程。随着计算机在控制系统分析、设计和实时控制中的广泛应用,这个问题显得越来越重要。在离散系统状态方程建立的基础上,本章也将讨论相应的状态方程求解问题,并将导出在形式上与连续系统状态方程的解一致的离散系统状态方程的解。 3.1 线性定常连续系统状态方程的解 在讨论一般线性定常连续系统状态方程的解之前,我们先讨论线性定常齐次状态方程的解,以便引入矩阵指数函数和状态转移矩阵的概念。所谓齐次状态方程,就是指状态方程中不考虑输入项的作用,满足方程解的齐次性的一类状态方程。研究齐次状态方程的解,就是研究系统本身在无外力作用下的自由运动。 3.2 状态转移矩阵及其计算 在状态方程求解过程中,关键是状态转移矩阵Φ(t)的计算。对于线性定常连续系统,该问题又归结到矩阵指数函数e At的计算。上一节已经介绍了基于拉氏反变换技术的矩阵指数函数e At的计算方法,下面讲述计算矩阵指数函数的其他3种常用方法。 3.2.1级数求和法线性系统的时域分析法(第七讲)
matlab实验二线性系统时域响应分析
线性系统的时域分析方法
大作业1(机电控制系统时域频域分析)
线性系统时域分析
自动控制原理_线性系统时域响应分析
第三章控制系统的时域分析法知识点
MATLAB线性系统时域响应分析实验
控制系统的时域分析
自动控制原理实验报告《线性控制系统时域分析》
控制系统的时域分析实验报告
线性系统的时域分析与校正习题及答案
自动控制原理》实验2(线性系统时域响应分析
第3章--线性系统的时域分析--练习与解答
线性系统的时域分析习题答案
第3章线性系统的时域分析习题答案
第3章线性系统的时域分析习题答案
线性系统的时域分析