第二章 连续系统的时域分析习题解答
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
试求系统的零输入响应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
!
6
1H (b) i
+ -
图题2-5
< 1 i
1H
(a) 1 6
F - u $
5
i 1H F 1 6
(a) + ,(b) 1 1
1F
1H + u -
i #
1F
+ u -
i
H 1 4
( 5
图题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 )2
4
(a)
解:(a)f u pu pu u i i pu
i f 4804221
11=+???
?=--+=-
.V )(e 5.0)(125
.05.0184)( 81
t t h p p f u p H t
ε-=?+=+==? .
V )()e 4.0e 4.2()( 6
4.214.0 6723
115.01111311)( )c (.
V )()e 2e 2()( 2
21223235.015.01)( )b (62
22
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(24121
121)( )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 ( 2
1 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 ),试求整个系统的冲激响应h (t )。
1F
+ -
(b) +
&
f
1 3 ;
1
+ u -
(c) + f - 【 1
3
图题2-8
+ f
-
$ 1F 1H + u
-(a)
+ f -
)
1
+ -
(b)
1 + f
-
… 1H
(c) 1
2 + u - 图题2-9
#
y (t )
h 2(t ) h 1(t )
f (t )
解:求和号后的冲激响应为)1()(-+t 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 επδ=-=∑∞
=
图题2-11
(a)
(b)
(c)
(d)
(1)
(2)
t
解:;
)()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 -εεεε-=--=-==--=?τετ
.
)()(sin sin
(sin )()
( )6(; )]1()2([sin )]2()1()[1(sin )
2()1()( )5(;
)()sin cos (e 2
1)(j 4)
e e ()e e (j je 2 )()e e ()
1j (j 21)()e e ()1j (j 21)( )
e e (j 21)( )4(0
211 j j j j j j j j 2t 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 ) = f 1(t )* f 2(t ) 。
`
解:
图题2-13
e -t (t ) {
(b)
)
1()e e (e )( : ;
e e )( , 1
e 2)e e (2e e 2e )( , 1
: , ,0)( : )b ()
4(]e 1[)2(]e 1[)4()2(*)()e ( )4()2(*)(e )( )a (1
2
2
22
121212
1
1 2)4()2(0 0 ][][]
[][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 已知)()1e
()()(t t t t t f t
εε-+=-*,求f (t ) 。
解:微分:)e 1()()1e 0()()e 1()()(0
t 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所示,试求系统的零状态响应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-15
图题2-16
(a) 输入 0
4
(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)( 34
221e e )( 2x 2222x 212
121221x 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 )时,系统的零状态响应。
解:(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)(:2222332333221t 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 εεεεδ------===+=
+++=
?+=+=--===*显然
+
1V -
1 1 2H
1
2 F + 2V -
a b S i L
+C -
+ u - 图题2-18
∫
f (t )2
-
-
.
y (t )
(a)
x 1
x 2
∫
f (t )-3 -
(
y (t )
(b)
∫
-2
x 1
x 2
x 3
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 025.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 212
121t 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 +=++,当激励
)(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 εεεεεετετττ----------------------=-=∴+--=---=-=-=+-+=
+++=
?
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
《机电系统控制基础》大作业一 基于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.掌握典型输入信号(单位脉冲、单位阶跃、单位速度、单位加速度、正弦信号)的拉氏变换表达式。 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
第三章线性系统的时域分析方法 教学目的:通过本章学习,熟悉控制系统动态性能指标定义,掌握线性系统稳定的充要条件和劳斯判椐的应用,以及稳态误差计算方法,掌握一阶、 二阶系统的时域分析方法。 教学重点:掌握系统的动态性能指标,能熟练地应用劳斯判椐判断系统稳定性,二阶系统的动态响应特性分析。 教学难点:高阶系统的的动态响应特性分析。 本章知识结构图: 系统结构图闭环传递函数 一阶标准式 二阶标准式 特征方程稳定性、稳定域 代数判据 误差传递函数误差象函数终值定理稳态误差开环传递函数系统型别、开环增益 公式 静态误差系数 第九讲
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.已知一单位负反馈系统的单位阶跃响应曲线如下图所示,求系统的闭环传递函数。 解答: ①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)
实验一线性控制系统时域分析 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时的调节时间。讨论阻尼比变化对系统性能的影响。
第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系统模型。
专业:电气工程及其自动化 学号: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 ●系统的零输入响应可以用经典法求解,在其中 只有自然响应部分; ●系统的零状态响应也可以用经典法求解,但是 用卷积积分法更加方便。借助于计算机数值计算,可以求出任意信号激励下的响应(数值解)。 ●卷积法要求激励信号是一个有始信号,否则无
实验三 连续时间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)
实验三 二阶系统的时域分析 一、实验目的 1、通过考察系统的过渡过程指标,研究二阶系统的特征参数—阻尼比和自然频率对系统特性的影响,以及系统特征根的位置与过渡过程的关系。 2、学习自己设计实验,安排适当的实验参数,达到以上实验目标。 二、实验内容 根据传递函数2 22 2)(n n n s s s G ωζωω++=的单位阶跃响应,求取过渡过程的质量指标。按表1的形式整理实验数据,分析实验结果,完成实验报告。 此时,系统的特征根为j j s n n βαζωζω±=-±-=2 2,11。 1、令ζ=0.5,取三种不同的n ω,观察根在根平面上的位置,求其过渡过程和它的质量指
标,进行比较。说明当ζ相同时,过渡过程的哪些指标是相同的? 00.2 0.4 0.6 0.8 1 1.2 1.4 ωn 改变,ζ=0.5不变 Tim e (sec) A m p l i t u d e
2、固定n ω,取ζ=0、0. 3、 0.5、0.7、1,观察根在根平面上的位置,求其过渡过程和它的质量指标。总结当ζ不同时,质量指标有哪些变化? 24681012141618 00.20.40.60.811.2 1.41.61.82 Time (sec) A m p l i t u d e
通过上面两图形与表格总结可以得出: n ω影响二阶系统过渡过程中的峰值时间,过渡时间(在ζ不变的情况下,峰值时间随n ω增 大而减小,过渡时间随n ω的增大而减小) ζ影响几乎全部过渡过程指标,其中超调量,衰减比仅与ζ有关(超调量随着ζ的增大而 减小,衰减比随着ζ的增大而增大;在n ω不变的情况下,峰值时间随ζ增大而增大,过渡时间随ζ的增大而减小。) n ω,ζ对系统的稳态误差均没有影响,且均为0.
实验七 控制系统频域分析方法 1.实验目的 (1)熟练掌握Nyquist 图和Bode 图的绘制。 (2)熟练掌握利用Nyquist 图和Bode 图分析系统的性能。 2.实验仪器 (1)Matlab6.5应用软件安装版 一套 (3)PC 机 一台 3. 实验原理 依据MA TLAB 的建模指令,利用MATLAB 对系统仿真,分析系统的频率特性。 4. 实验步骤 (1)建立系统的MATLAB 模型,绘制系统Nyquist 图和Bode 图,分析系统稳定性 (2)求系统的幅值穿越频率和相位穿越频率,分析系统的稳定性。 (3)依据系统框图建立系统模型,利用LTI Viewer 分析系统的稳定性。 (4)绘制离散系统开环传递函数的Nyquist 图和Bode 图,绘制系统单位阶跃响应图。 5. 实验报告内容(选做其中三题) 1、绘制下列各单位反馈系统开环传递函数的Bode 图和Nyquist 图,并根据其稳定裕度判断系统的稳定性。(使用subplot 指令) ) 31)(2s 1)(s 1(10)s (G 1k s +++=)( )101)(s 1(s 10)s (G 2k s ++= )( ) 2.01)(s 1.01(s 10)s (G 32k s ++=)( )101)(s 1.01(s 10)s (G 42k s ++= )( 2、设单位反馈系统的开环传递函数为)12s (s K )s (G 2k ++=n n w s w ξ,其中无阻尼固有频率 Wn=90rad/s ,阻尼比ξ=0.2,试确定是系统稳定的K 的范围。 3、设系统如图7-22所示,试用LTI Viewer 分析系统的稳定性,并求出系统的稳定裕度及单位阶跃响应峰值. 4、设闭环离散系统结构如图7-23所示,其中) 1(10s +=s s G )(,1s =)(H ,绘制T=0.01s,1s 时离散系统开环传递函数的Bode 图和Nyquist 图,以及系统的单位阶跃响应曲线..
专业 班号 组别 指导教师 姓名 学号 实验名称 线性系统时域响应分析 一、实验目的 1.熟练掌握step( )函数和impulse( )函数的使用方法,研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2.通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 3.熟练掌握系统的稳定性的判断方法。 二、实验内容 1.观察函数step( )和impulse( )的调用格式,假设系统的传递函数模型为 1 4647 3)(2342++++++=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)') 方法二:用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)') 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 ,,,,σ。 程序如下: num= [0 0 4]; den1=[1 0 4]; den2=[1 1 4]; den3=[1 2 4];
一、实验目的和要求 1.用计算机辅助分析的办法,掌握系统的时域分析方法。 2.熟悉Simulink 仿真环境。 二、实验内容和原理 (一)实验原理 系统仿真实质上就是对系统模型的求解,对控制系统来说,一般模型可转化成某个微分方程或差分方程表示,因此在仿真过程中,一般以某种数值算法从初态出发,逐步计算系统的响应,最后绘制出系统的响应曲线,进而可分析系统的性能。控制系统最常用的时域分析方法是,当输入信号为单位阶跃和单位冲激函数时,求出系统的输出响应,分别称为单位阶跃响应和单位冲激响应。在MA TLAB 中,提供了求取连续系统的单位阶跃响应函数step ,单位冲激响应函数impulse ,零输入响应函数initial 等等。 (二)实验内容 二阶系统,其状态方程模型为 ?1x -0.5572 -0.7814 1x 1 = + u ?2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 1.画出系统的单位阶跃响应曲线; 2.画出系统的冲激响应曲线; 3.当系统的初始状态为x0=[1,0]时,画出系统的零输入响应; 4.当系统的初始状态为零时,画出系统斜坡输入响应; (三)实验要求 1.编制MA TLAB 程序,画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应; 2.在Simulink 仿真环境中,组成系统的仿真框图,观察单位阶跃响应曲线并记录之。 三、主要仪器设备 计算机一台以及matlab 软件,simulink 仿真环境 四、操作方法与实验步骤 1、程序解决方案:
在MATLAB命令窗口中输入下列命令:并返回系统的传递函数 其输出的曲线如下