第二章 连续系统的时域分析习题解答
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
试求系统的零输入响应y x (t )(t 0)。 解:,e e
)( ,3 ,1 )1(32121t t
A A t y p p --+=-=-=
. 0 , )e 12(1)(1
21444200 ,
)e ()( , 2 ,0 )3(. 0 , 2sin e 5.0)(905.00cos 240)sin (cos 21cos 0 ,
)2cos(e )( , 2j 2 ,0 )2(;0 , e 5.1e 5.3)(5
.15.3312 232132323123213 ,2123213233232132213 ,213212121 t t t y A A A A A A A A A A t A A t y p p t t t y A A A A A A A A A A A A t A A t y p p t t y A A A A A A t t t t t t ------+-=????
??-=-==???
???+-=-=+=++=-===????
???
-===???
???=+-=+=++=±-==-=????-==????--=+= 2-5 已知图题2-5各电路零输入响应分别为:
.
0 ,V sin e
6cos e
2)( (b); 0 ,V e 4e 6)( )a (3343x x t t t t u t t u t
t
t t ----+=-=
求u (0-)、i (0-)。
解:;V 246)0()0( )a (x =-==+-u u
.
0)66(1.0)0()0( V;202)0()0( )b (A
3
512)1618(61)0()0( x x x x =+-===+===++-==+-+-+-i i u u i i 2-6 图题2-6所示各电路:
(a) 已知i (0-) = 0,u (0-) = 5V ,求u x (t ); (b) 已知u (0-) = 4V ,i (0-) = 0,求i x (t ); (c) 已知i (0-) = 0,u (0-) = 3V ,求u x (t ) .
解:0650650)( )(2
=++?=++?=p p p
p p a Z
(b)
图题2-5
(a) 1 6
F
(a)
(b)
(c)
图题2-6
.
0 ,V e 10e 15)( ,10 ,1532050)
0()0(' ,V 5)0(e e )(3 , 2 32212
121322121---------=-==????--=+=?==
=+=?-=-=? t t u A A A A A A C i u u A A t u p p t t x x x t t x
.
0 ,V e e 4 ,4 ,1 ,15035)0(' ,V 3)0( ,e e ,
4,1,045045: )c (. 0 ,A sin e 4)( ,2 ,4 sin cos 4cos 04401)0(' , 0)0( )
cos(e )(11 , 11 0
22011110)( )b (421421212212
1212
121212-------------==-=-=+?-==+=-=-==++?=++=-==???
?--==?=+?-==+=?--=+-=?=++?=+++?= t u A A u u A A u p p p p p p t t t i A A A A A A A A i i A t A t i j p j p p p p p p t t x x x t t x t x x x t x 同理/πY 2-7 已知三个连续系统的传输算子H (p )分别为:
. )
2(1
3 )3( ; )84()12( )2( ; )3)(1(42
)1(2
2+++++-+++p p p p p p p p p p 试求各系统的单位冲激响应h (t )。 解:; )()e e ()(3
111)( )1(3t t h p p p H t t ε--+=?+++=
.
)()e 4
1e 2541()(241
)2(5.241)( )3(;
)()2sin e 875.02cos e 8181()( 2)2(2
875.0)2(8181)( 5.1,81)84(1)21()81(8481)( )2(222222
2222t t t h p p p p H t t t t h p p p p H B A p p p p B p A p p B Ap p p H t t t t εε-----+=?+-++=-+-=?++?-++-=?-==?++--+-=++++-= 2-8 求图题2-8所示各电路中关于u (t )的冲激响应h (t )
解:(a)f u pu pu u i i pu
i f 480
422111=+???
?=--+=- e 5.0)(125
.05.0184)( 81
t h p p f u p H t
ε-=?+=+==? (a)
.
V )()e 4.0e 4.2()( 64.214.0 6723115.01111311)( )c (.
V )()e 2e 2()( 2
21223235.015.01)( )b (6222
t t h p p p p p p
p p p H t t h p p p p p p
p p H t t t t εε-----=?+++-=++=+++?+=-=?+-+=++=
++=
2-9 求图题2-9所示各电路关于u (t )的冲激响应h (t )与阶跃响应g (t )。
解
:
2
cos 21)(0 2cos 42)(21)()( ),
(2
sin 42)(21)( 2
1
)
21(24121121)( )a ( t _0 22
2
t t t
t d h t g t t t t h
p p p p p p p p H ετεττεδ=??????+==-=?+-=+=++=-?
),(e 41)(21)(2
141211212111)( )b ( 21 t t t h p p p p p p H t εδ-+=?++=++=++=. )()e e ()(0
]e e [)()( )e e 2()(1
12211122)( )c (()e 211()(0e 21)(21)()( 2 2_
0 2 21 2
1 _0 t t
][t t t d h t g t h p p p
p
p p H t t t t d h t g t
t t t t εεττεεεετττττ---
------
--=+-==-=?+-+=+-+=
-=-==??2-10 如图题2-10所示系统,已知两个子系统的冲激响应分别为h 1(t )(t 1),h 2(t )(t ),
解:求和号后的冲激响应为)1()(-+t t δδ,于是整个系统的冲激响应为:
(b)
u
(c) 图题2-8
u
(a)
(b)
(c)
图题2-9
图题2-10
y (t )
f (t )
)1()()(-+=t t t h εε
2-11 各信号波形如题图2-11所示,试计算下列卷积,并画出其波形。
. )(')( )3( ; )()( )2( ; )()( )1(41 31 21t f t f t f t f t f t f ***
解:
.
)3()3(2
1)1()1(23)
1()1(23)3()3(21 )3()3(21)1(
)1()1()1(21 )
1()1(2
1)1
()1()3()3(21 )
1()1()(')( )3();
6()6
(
2
1)5()5(21)4()4(21)3()3( )2()2(2
1)1()1(21)(21 )
4()3()2()()( )2();
4()4(2
1)2()2( )
()2()2()4()4(2
1 )
2()2()()( )1(2()2(21)()2()2(21)(11411113111211-----+++-++=-----+++---+++-++=--+=--+--------+------=-+---=--+---+++-++=-++=--+-++=t t t t t t t t t t t t t t t t t t t t t f t f t f t f t t t t t t t t t t t t t t t f t f t f t f t f t t t t t t t t t t t f t f t f t f t t t t t t t f εεεεεεεεεεεεεεεεεεεεεεεεε***
2-12 求下列各组信号的卷积积分。
;
)2()1()( , )]1()([sin )( )5(; )(sin )( , )(e )( )4( ; )(e )( , )(e )( )3(; )(e )( , )()( )2( ; )1()( , )()( )1(2121 2212121++-=--=======-==----t t t f t t t t f t t t f t t f t t f t t f t t f t t f t t f t t f t t t t δδεεπεεεεεεεε . )( sin )( , )()( )6(20
1t t T t f nT t t f n επδ=-=∑∞
=
解:;
)()e e ()()e e (1
21)( )3(;
)()e 1()()d e ()( )2(;
)1()1()( )1(22 t
t t t y t t t y t t t y t --t t --t -t -εεεε-=--=-==--=?τετ
图题2-11
(a)
(b)
(c)
(d)
(1)
(2)
t
.)(
)
(sin
sin
(
sin
)
(
)
(
)6(
;)]
1
(
)2
(
[
sin
)]
2
(
)1
(
)[
1
(
sin
)2
(
)1
(
)(
)5(
;)(
)
sin
cos
(e
2
1
)(
j4
)
e
e
(
)
e
e(j
je
2
)(
)
e
e(
)1
j(j2
1
)(
)
e
e(
)1
j
(j2
1
)(
)
e
e(
j2
1
)(
)4(
2
1
1
j
j
j
j
j
j
j
j
2
t
t
T
t
T
T
nT
t
f
t
y
t
t
t
t
t
t
t
f
t
f
t
y
t
t
t
t
t
t
t
y
t
f
n
n
-t
t-
t
t-
t
-t
t-
-t
t
-t
t-
t
ε
π
ε
π
π
ε
ε
π
ε
ε
π
ε
ε
ε
ε
=
=
-
=
+
-
+
+
-
-
-
-
=
+
+
-
=
+
-
=
-
+
+
-
=
-
-
-
-
-
-
=
?
-
=
∑
∑∞
∞
=
=
2-13求图示各组波形的卷积积分y(t) = f1(t)* f2(t) 。
解:
)
1(
)
e
e(
e
)(
:
;
e
e
)(
,
1
e2
)
e
e(2
e
e2
e
)(
,
1
:
,
,0
)
(
:
)b(
)4
(
]
e
1[
)2
(
]
e
1[
)4
(
)2
(
*
)(
)
e
(
)4
(
)2
(
*
)
(
e
)(
)a(
1
2
2
2
2
1
2
1
2
1
2
1
1
2
)4
(
)2
(
]
[
]
[
]
[
]
[
t
t
y
d
t
y
t
d
d
t
y
t
f
t
t
t
t
t
t
t
d
t
y
t
t
t
t
t
t
t
t
t
t
-
-
+
=
=
=
>
-
=
-
+
=
+
=
<
∴
≠
-∞
-
-
-
-
-
=
-
-
-
-
=
-
-
-
=
+
∞
-
+
+
+
+
+
∞
-
-
-
-
-
-
-
?
?
?
?
-
-
ε
τ
τ
τ
ε
ε
δ
δ
ε
δ
δ
τ
τ
ε
τ
τ
τ
τ
τ
或
时
时
可借助于图解分段计算
不可用微分积分性质
注意Θ
2-14已知)(
)1
e
(
)(
)(t
t
t
t
t
f tε
ε-
+
=-
*,求f(t) 。
解:微分:)
(
)
e
1(
)
(
)1
e
0(
)
(
)
e
1(
)
(
)
(0t
t
t
t
f t
tε
δ
ε
ε-
--
=
-
+
+
-
=
*
再微分:)
(
)
(
e
)
e
1(
)
(
e
)
(
)
(0t
f
t
t
t
t
f t
t=
=
-
+
=-
-ε
ε
δ
*.
2-15某LTI系统的激励f(t)和冲激响应h(t)如图题2-15所
图题2-13
(a)
(b)
示,试求系统的零状态响应y f (t ),并画出波形。 解:][}{)2()(*)]2()([2
1)( 0 ----=
?
-t t d t y t
f δδττετετ )4()]4()2([)2(4
1)2()]2()([41)2()()2()]2()([41222][}{--------+--=---+--=t t t t t t t t t t t t t t εεεεεεδδεεε*
2-16 图题2-16表示一个LTI 系统的输入-输出关系。试求出该系统的冲激响应。
解:)2(2)(2)()2(2)(2)(--=?--=t t t h t f t f t y δδ
2-17 已知某系统的微分方程为)(3)(2)(3)(t f t f t y t y +'=+'+'',0-初始条件
2)0( , 1)0(='=--y y ,试求:
(1) 系统的零输入响应y x (t );
(2) 激励f (t )(t )时,系统的零状态响应y f (t )和全响应y (t ); (3) 激励f (t ) e 3t (t )时,系统的零状态响应y f (t )和全响应y (t )。 解:(1) 算子方程为:)()3()()2)(1(t f p t y p p +=++
)
()e 2
5e 223()()()( )
()e 2
1e 223()()()( )()e e 2()(2
112233
)( )2(; 0 ,e 3e 4)( 3
4
221e e )( 2x 222
2x 21
2121221x t t y t y t y t t t h t y t t h p p p p p p H t t y A A A A A A A A t y t t t t t t f f t t t t εεεε------------+=+=+-==-=?+-+=+++=
-=????-==????--=+=?+=∴* )
()e
4e 5()()()( )()e e ()(e )()( )3(2x 23t t y t y t y t t t h t y t
t
t t t f f εεε------=+=-==*
2-18 图题2-18所示的系统,求当激励f (t ) e t (t )时,系统的零状态响应。
图题2-15
图题2-16
(a) 输入
解:(a) 令f (t ) = (t ),则y (t ) = h (t ),
)
(e )32()](e [)](e 3)(2[)()(e 3)(2)()1
32()()1112(
)(2)( ,)(:21221t t t t t t y t t t p t p p p
t h x x t h x t f px x t t t t f εεεδεδδδ-----=-=-=+-=+-+=?-=-==*显然 (b) 令f (t ) = (t ),则y (t ) = h (t ),
)
()e e ()](e [)(e )()(e )(21)(2
31
)()1()( ,23)(:2222
332333221t t t t y t t p t f p p p t h x p x x t h x px t f x p px x t t t t t f εεεεδ------===+=+++=
?+=+=--===*显然
2-19 图题2-19所示电路,t < 0时S 在位置a 且电路已达稳态;t = 0时将S 从a 板到b, 求t > 0时的零输入响应u x (t )、零状态响应u f (t )和全响应u (t )。
解:i ) 先求零状态响应u f (t ):
)
()e e 1(2)()](1)][e 2e 5.0()([ )(),
()e 2e 5.0()()(2
2
5.05.012225.05.0211121)(25.0 0
25.025.0}{t t d t t u t t t h p p p p p p
p p H t t t
t t f εεττετδεδττ--------=-+++-=++-=++++-=+-+-+=+-+=
?
ii ) 求零输入响应u x (t ):t t
A A t u 225.01x e e
)(--+=
s
V 5.2)0(')0('1)0('2)0(' ,2
1)0('A
1)0( ,V 1)0( ,V 2111)0(x x x x x x x x /-=+?=?-=-=?-=-==+?=++++++++C L C L C L u i u u i i u u
).()e e
1(1)()()(V 1)0( ;0 ,V e e
)(1125.05.2225.0x 25.0x 212121t t u t u t u u t t u A A A A A A t
t t t
f ε-------+=+==+=????==????--=-+=∴ 2-20 已知某系统的微分方程为)(3)(')(2)(' 3)(" t f t f t y t y t y +=++,当激励
图题2-19
C 图题2-18
f (t )y (t )
(a)
f (t )y (t )
(b)
)(t f =)(e 4t t ε-时,系统的全响应)()e 6
1e 27e 314()(42t t y t t t ε-----=; 试求零输入
响应y x (t )与零状态响应y f (t )、自由响应与强迫响应、暂态响应与稳态响应。 解:
.
, )();
()e 2
7e 314(: );(e 61:)( )()e 3e 4()()()()( )()e 3
221e 61( )
()]e 1(e 2
1)e 1(e 32[)(]e 2e 2[e )(),()e e 2()( ,2
112233)(242x 24223 0 )(2)(422}{不含稳态响应全为暂态自由响应强迫响应零状态响应零状态响应t y t t t t y t y t y t e t t d t y t t h p p p p p p H t t t t t t t t t t t t t
t t t t f f εεεεεετετττ----------------------=-=∴+--=---=-=-=+-+=
+++=
?
第三章 线性系统的时域分析法 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.有一位置随动系统,其结构图如下图所示,其中K = 4。求该系统的:1)自然 k 振荡角频率;2)系统的阻尼比;3)超调量和调节时间;4)如果要求 <0.707 , 值。 应怎样改变系统参数 K k 2.已知受控对象的开环传递函数为
(1)单位反馈时,计算单位脉冲响应的输出。 (2)试采用速度反馈方法,使得系统的阻尼比ζ=05.,确定速度反馈系数τ的值,并计算性能改善后的动态性能。 解 (1)单位反馈时,闭环传递函数为 其单位脉冲响应为 响应曲线为等幅振荡的,所以该系统仅作单位反馈,不能实现调节作用。 (2)增加速度反馈如图所示。 闭环传递函数为 ζωτ=,所以 阻尼比ζ=05.,则有2 n τ=?= 20.50.95 此时,系统阶跃响应的超调量为 调节时间为 3.已知速度反馈控制系统如图所示,要求系统的超调量为20%,峰值时间为1秒,试计算相应的前向增益K与速度反馈系数K 的值。如果保持K值不变,Kf为零时,计算超调量增大值。
解上述系统的闭环传递函数为 比较二阶系统的标准式有 给定的性能指标为 上述指标与系统特征参数ζ和ωn的关系为: 解得 所以: 当K=125.,Kf=0时,也就是没有速度反馈时,闭环传递函数成为: 阻尼比:
超调量增大为: 4.对下图所示系统,试求K为何值时,阻尼比ζ=0.5。并求此时系统单位阶跃响应的最大超调量和调整时间。 解:系统开环传函为: 系统闭环传函为: 最大超调量: 调整时间
5. 系统结构如图,欲使超调量бp =0. 2, 过渡过程时间t s =1秒(Δ=0.02), 试确定K 和τ的值。 答案: ()2222(2)2n n n K s s K s K s ωτζωωΦ==+++++ 0.456ζ= 8.77 n ω= 277n K ω== 0.078τ= 6. 题图所示机械系统,当受到 F =40N 力的作用时,位移量xt ()的阶跃响应如图所示,试确定机械系统的参数m ,k, f 的值。 解: 图示机械系统的传递函数为 由图所示稳态值()1c ∞=,由终值定理 得到 K=40N/m 由超调量: 峰值时间:
第二章连续系统的时域分析 求响应:经典法:已知f(t)、x{0} 全响应y(t)= y f(t)+y x(t) 卷积积分法:先求n(t),已知f(t) y f(t)=h(t) f(t) 主要内容: 一经典法求LTI系统的响应: 齐次解自由响应瞬态零输入 特解强迫响应稳态(阶跃、周期)零状态二冲击响应与阶跃响应:(定义、求解方法仍为经典法)三卷积积分:(定义、图示法求卷积) 四卷积积分的性质:
§2.1 LTI 系统的响应(经典法) 一 常系数线性微分方程的经典解 n 阶:y )(n (t)+ a n-1y )1(-n (t)+…+ a 1y )1((t)+ a 0y(t) = b m f )(m (t)+ b m-1 f )1(-m (t)+……+ b 1 f )1((t)+ b 0f(t) 全解:y(t)=齐次解y h (t)+ 特解y p (t) 1 齐次解:y h (t)=∑=n i t e i C i 1 λ(形式取决于特征根) 特征方程: λ)(n (t)+ a n-1λ)1(-n (t)+… + a 1 λ(t)+ a 0=0 特征根:决定齐次解的函数形式,表2-1 如为2个单实根λ1、λ2, y h (t )=e C t 11 λ +e C t 22 λ 如为2重根(λ+1)2=0,λ= - 1,y h (t)=C 1te -t +C 0e -t 系数C i :求得全解后,由初始条件确定 2 特解: 函数形式:由激励的函数形式决定,与特征根有关系,表2-2 如:f(t)为常数 )(t ε, y p (t)=P 0 f(t)=t 2, y p (t)= P 2t 2+ P 1t+ P 0 f(t)=e -t ,λ= - 2,不等 y p (t)=P e -t f(t)= e -t ,λ= - 1,相等 y p (t)=P 1te -t +P 0e -t 系数P i :由原微分方程求出 3 全解:y(t)= y h (t)+ y p (t)=∑=n i t e i C i 1 λ+ y p (t) 此时利用y(0),y ‘(0),求出系数C i
第三章线性系统的时域分析方法 教学目的:通过本章学习,熟悉控制系统动态性能指标定义,掌握线性系统稳定的充要条件和劳斯判椐的应用,以及稳态误差计算方法,掌握一阶、 二阶系统的时域分析方法。 教学重点:掌握系统的动态性能指标,能熟练地应用劳斯判椐判断系统稳定性,二阶系统的动态响应特性分析。 教学难点:高阶系统的的动态响应特性分析。 本章知识结构图: 系统结构图闭环传递函数 一阶标准式 二阶标准式 特征方程稳定性、稳定域 代数判据 误差传递函数误差象函数终值定理稳态误差开环传递函数系统型别、开环增益 公式 静态误差系数 第九讲
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.掌握系统动态响应的概念,能够从系统的响应中分离出稳态响应分量和瞬态响应分量;掌握系统动态响应的性能评价指标的概念及计算方法(对于典型二阶系统可以直接应用公式求解,非典型二阶系统则应按定义求解)。 解释:若将系统的响应表达成拉普拉氏变换结果(即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.已知一单位负反馈系统的单位阶跃响应曲线如下图所示,求系统的闭环传递函数。 解答: ①max ()100100()X X %%e %X δ-∞=?=?∞ 由 2.1820.090.6082e ξ-==?= ②0.8 4.946m n t ω==?= ③2222224.4648.9222 6.01424.46 6.01424.46 n B n n W K s s s s s s ωωω=?=?=++++++ 2.已知系统如下图所示,求系统的单位阶跃响应,并判断系统的稳定性。 解答: ()() ()210 1101061010.511B s s W s s s s s +==+++++ 3.16n ω==, 260.95n ξωξ=?
( )()1sin n t c X t ξωωθ-= ,arctg θ= ()31 3.2sin 0.98718.19t e t -=-+? (5分) 系统根为 1,2632P j -= =-±,在左半平面,所以系统稳定。 3.一阶系统的结构如下图所示。试求该系统单位阶跃响应的调节时间t s ;如果要求t s (5%)≤ 0.1(秒),试问系统的反馈系数应取何值? (1)首先由系统结构图写出闭环传递函数 得 T =0.1(s ) 因此得调节时间 t s =3T =0.3(s),(取5%误差带) (2)求满足t s (5%) ≤0.1(s )的反馈系数值。 假设反馈系数K t (K t >0) ,那么同样可由结构图写出闭环传递函数 由闭环传递函数可得 T = 0.01/K t 100()10()100()0.1110.1c B r X s s W s X s s s ===++?1001/()1000.0111t B t t K s W s K s s 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 (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:控制系统的时域分析 课程名称:自控控制原理 专业:电气工程及其自动化 班级: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)
第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用计算机辅助分析的办法,掌握系统的时域分析方法。 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)
实验三 连续时间LTI 系统的时域分析 一、实验目的 1、学会使用符号法求解连续系统的零输入响应和零状态响应 2、学会使用数值法求解连续系统的零状态响应 3、学会求解连续系统的冲激响应和阶跃响应 二、实验原理及实例分析 1、连续时间系统零输入响应和零状态响应的符号求解 连续时间系统可以使用常系数微分方程来描述,其完全响应由零输入响应和零状态响应组成。MATLAB 符号工具箱提供了dsolve 函数,可以实现对常系数微分方程的符号求解,其调用格式为: dsolve(‘eq1,eq2…’,’cond1,cond2,…’,’v’) 其中参数eq 表示各个微分方程,它与MATLAB 符号表达式的输入基本相同,微分和导数的输入是使用Dy ,D2y ,D3y 来表示y 的一价导数,二阶导数,三阶导数;参数cond 表示初始条件或者起始条件;参数v 表示自变量,默认是变量t 。通过使用dsolve 函数可以求出系统微分方程的零输入响应和零状态响应,进而求出完全响应。 [实例1]试用Matlab 命令求齐次微分方程0)()(2)(='+''+'''t y t y t y 的零输入响应,已知起始条件为2)0(,1)0(,1)0(=''='=---y y y 。
3、连续时间系统冲激响应和阶跃响应的求解 在连续时间LTI系统中,冲激响应和阶跃响应是系统特性的描述。在MATLAB中,对于冲激响应和阶跃响应的数值求解,可以使用控制工具箱中提供的函数impulse和step来求解。 ) , ( ) , ( t sys step y t sys impulse y = = 其中t表示系统响应的时间抽样点向量,sys表示LTI系统模型。
实验二 线性系统时域响应分析 一、实验目的 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控制系统时域响应的基本概念,典型输入信号及意义; 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章 线性系统的时域分析 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 应满足什么关系?
专业:电气工程及其自动化 学号:07050443 05 姓名: 实验一 二阶系统时域分析 一、 实验目的 1. 研究二阶系统的两个重要参数ξ、n ω与系统结构之间的关系。 2. 观察系统在阶跃输入作用下的响应,运用基本理论,分析系统过度过程特点及各种参数对其学习过程的影响,从而找出改善系统动态性能的方法,并在实验中加以验证。 3. 学习二阶系统阶跃响应的测试方法。 4. 掌握开环传递函数与闭环传递函数之间的对应关系,以及ξ、n ω与传递函数系数之间的关系。 二、 实验内容 选择适当的元器件建立单位负反馈二阶系统。 开环传递函数由积分环节和惯性环节构成:()() 1S T S T K S G 21+= 令T T T 21==。 1. 设1T = 改变K 值,使阻尼比ξ,分别为0、0.5、0.7、1、1.5;观察并记录在单位阶跃信号作用下,不同阻尼比时,系统输出响应曲线,并测量系统的超调量σ%、上升时间r t 、峰值时间p t 、调节时间s t 。 (1)当阻尼比ξ无限大时: (2)当阻尼比ξ=0.5时:
(3)当阻尼比ξ=0.7时: (4)当阻尼比ξ=1时: (5)当阻尼比ξ=1.5时:
2. 设定K 值 使ξ=0.707,改变时间常数T ,观察并记录在单位阶跃信号作用下,系统输出曲线,并测量系统的超调量σ%、上升时间r t 、峰值时间p t 、调节时间s t 。并与(1)的结果加以比较。 (1) 当T=0.1时: (2) 当T=1时:
(3) 当T=1.5时: 3. 改变时间常数 使1T 不等于2T ,观察并记录输出波形的变化情况。 (1) 当1T 1=,2T 2=时: (2) 当2T 1=,1T 2=时:
第二章 连续时间系统的时域分析 §2-1 引 言 线性连续时间系统的时域分析,就是一个建立和求解线性微分方程的过程。 一、建立数学模型 主要应用《电路分析》课程中建立在KCL 和KVL 基础上的各种方法。 线性时不变系统的微分方程的一般形式可以为: )()(...)()()()(...)()(0111101111t e b t e dt d b t e dt d b t e dt d b t r a t r dt d a t r dt d a t r dt d m m m m m m n n n n n ++++=++++------ 二、求解(时域解) 1、时域法 将响应分为通解和特解两部分: 1) 通解:通过方程左边部分对应的特征方程所得 到的特征频率,解得的系统的自然响应(或自由响应); 2) 特解:由激励项得到系统的受迫响应;
3)代入初始条件,确定通解和特解中的待定系数。 经典解法在激励信号形式简单时求解比较简单,但是激励信号形式比较复杂时求解就不容易,这时候很难确定特解的形式。 2、卷积法(或近代时域法,算子法) 这种方法将响应分为两个部分,分别求解: 1)零输入响应:系统在没有输入激励的情况下,仅仅由系统的初始状态引起的响应 r )(t ; zi 2)零状态响应: 状态为零(没有初始储能)的条件下,仅仅由输入信号引起的响应 r )(t 。 zs ●系统的零输入响应可以用经典法求解,在其中 只有自然响应部分; ●系统的零状态响应也可以用经典法求解,但是 用卷积积分法更加方便。借助于计算机数值计算,可以求出任意信号激励下的响应(数值解)。 ●卷积法要求激励信号是一个有始信号,否则无
第3章 线性系统的时域分析 本章讨论线性系统的运动分析。主要介绍连续系统与离散系统的状态空间模型的求解、状态转移矩阵的性质和计算以及连续系统状态方程的离散化。本章最后介绍基于Matlab的状态空间模型求解与控制系统的运动仿真问题的程序设计与仿真计算。 建立了系统的数学描述之后,接下来要对系统作定量和定性分析。定量分析主要研究系统对给定输入信号的响应问题,也就是对描述系统的状态方程和输出方程的求解问题。定性分析主要研究系统的结构性质,如能控性、能观性、稳定性等。本章先讨论用状态空间模型描述的线性系统的定量分析问题,即状态空间模型的求解问题。根据常微分方程理论求解一个一阶定系数线性微分方程组是很容易的,可是求解一个一阶变系数线性微分方程组却非易事。状态转移矩阵的引入,使得定常系统和时变系统的求解公式具有一个统一的形式。为此,本章将重点讨论状态转移矩阵的定义、性质和计算方法,并在此基础上导出状态方程的求解公式。本章讨论的另一个中心问题是连续系统状态方程的离散化,即建立连续系统的离散系统状态方程。随着计算机在控制系统分析、设计和实时控制中的广泛应用,这个问题显得越来越重要。在离散系统状态方程建立的基础上,本章也将讨论相应的状态方程求解问题,并将导出在形式上与连续系统状态方程的解一致的离散系统状态方程的解。 3.1 线性定常连续系统状态方程的解 在讨论一般线性定常连续系统状态方程的解之前,我们先讨论线性定常齐次状态方程的解,以便引入矩阵指数函数和状态转移矩阵的概念。所谓齐次状态方程,就是指状态方程中不考虑输入项的作用,满足方程解的齐次性的一类状态方程。研究齐次状态方程的解,就是研究系统本身在无外力作用下的自由运动。 3.2 状态转移矩阵及其计算 在状态方程求解过程中,关键是状态转移矩阵Φ(t)的计算。对于线性定常连续系统,该问题又归结到矩阵指数函数e At的计算。上一节已经介绍了基于拉氏反变换技术的矩阵指数函数e At的计算方法,下面讲述计算矩阵指数函数的其他3种常用方法。 3.2.1级数求和法
实验三 连续时间LTI 系统的时域分析 一、实验目的 1.学会用MA TLAB 求解连续系统的零状态响应; 2. 学会用MATLAB 求解冲激响应及阶跃响应; 3.学会用MA TLAB 实现连续信号卷积的方法; 二、实验原理 1.连续时间系统零状态响应的数值计算 我们知道,LTI 连续系统可用如下所示的线性常系数微分方程来描述, () ()0 ()()N M i j i j i j a y t b f t ===∑∑ 在MA TLAB 中,控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数lsim 。其调用格式 y=lsim(sys,f,t) 式中,t 表示计算系统响应的抽样点向量,f 是系统输入信号向量,sys 是LTI 系统模型,用来表示微分方程,差分方程或状态方程。其调用格式 sys=tf(b,a) 式中,b 和a 分别是微分方程的右端和左端系数向量。例如,对于以下方程: ''''''''''''32103210()()()()()()()()a y t a y t a y t a y t b f t b f t b f t b f t +++=+++ 可用32103210[,,,];[,,,];a a a a a b b b b b == (,)sys tf b a = 获得其LTI 模型。 注意,如果微分方程的左端或右端表达式中有缺项,则其向量a 或b 中的对应元素应为零,不能省略不写,否则出错。 例3-1 已知某LTI 系统的微分方程为 y’’(t)+ 2y’(t)+100y(t)=f(t) 其中,' (0)(0)0,()10sin(2)y y f t t π===,求系统的输出y(t). 解:显然,这是一个求系统零状态响应的问题。其MATLAB 计算程序如下: ts=0;te=5;dt=0.01; sys=tf([1],[1,2,100]); t=ts:dt:te; f=10*sin(2*pi*t); y=lsim(sys,f,t); plot(t,y); xlabel('Time(sec)'); ylabel('y(t)'); 2.连续时间系统冲激响应和阶跃响应的求解 在MATLAB 中,对于连续LTI 系统的冲激响应和阶跃响应,可分别用控制系统工具箱提供的函数impluse 和step 来求解。其调用格式为 y=impluse(sys,t)