二阶系统的阶跃响应doc.docx
- 文档编号:3886012
- 上传时间:2022-11-26
- 格式:DOCX
- 页数:30
- 大小:26.65KB
二阶系统的阶跃响应doc.docx
《二阶系统的阶跃响应doc.docx》由会员分享,可在线阅读,更多相关《二阶系统的阶跃响应doc.docx(30页珍藏版)》请在冰豆网上搜索。
二阶系统的阶跃响应doc
第3章辅导
控制系统典型的输入信号
1.阶跃函数
阶跃函数的定义是
0,t0
xr(t)
A,t0
式中A为常数。
A等于1的阶跃函数称为单位阶跃函数,如图所示。
它表示为
xr(t)=l(t),或xr(t)=u(t)
单位阶跃函数的拉氏变换为
Xr(s)=L[1(t)]=1/s
在t=0处的阶跃信号,相当于一个不变的信号突然加到系统上;对于恒值系统,相当
于给定值突然变化或者突然变化的扰动量;对于随动系统,相当于加一突变的给定位置信号。
2.斜坡函数
这种函数的定义是
xr(t)
0,t0
At,t0
式中A为常数。
该函数的拉氏变换是
2
Xr(s)=L[At]=A/s
这种函数相当于随动系统中加入一按恒速变化的位置信号,该恒速度为A。
当A=l时,
称为单位斜坡函数,如图所示。
1
3.抛物线函数
如图所示,这种函数的定义是
0,t0
xr(t)
At
2
t0
式中A为常数。
这种函数相当于随动系统中加入一按照恒加速变化的位置信号,该恒
加速度为A。
抛物线函数的拉氏变换是
23
Xr(s)=L[At]=2A/s
当A=1/2时,称为单位抛物线函数,即Xr(s)=1/s
3。
4.脉冲函数
这种函数的定义是
A
0t(0)
xr(t)
0,t0,t(0)
式中A为常数,ε为趋于零的正数。
脉冲函数的拉氏变换是
Xrs)Llim
(
0
A
A
当A=1,ε→0时,称为单位脉冲函数δ(t),如图所示。
单位脉冲函数的面积等于l,
即
2
(t)dt1
在t=t0处的单位脉冲函数用δ(t-t0)来表示,它满足如下条件
幅值为无穷大、持续时间为零的脉冲纯属数学上的假设,但在系统分析中却很有用处。
单位脉冲函数δ(t)可认为是在间断点上单位阶跃函数对时间的导数,即
反之,单位脉冲函数δ(t)的积分就是单位阶跃函数。
控制系统的时域性能指标
对控制系统的一般要求归纳为稳、准、快。
工程上为了定量评价系统性能好坏,必须给
出控制系统的性能指标的准确定义和定量计算方法。
1动态性能指标
动态性能指标通常有如下几项:
延迟时间
t阶跃响应第一次达到终值h()的50%所需的时间。
d
上升时间
t阶跃响应从终值的10%上升到终值的90%所需的时间;对有振荡的系统,
r
也可定义为从0到第一次达到终值所需的时间。
峰值时间tp阶跃响应越过稳态值h()达到第一个峰值所需的时间。
调节时间
t阶跃响到达并保持在终值h()5%误差带内所需的最短时间;有时也用
s
终值的2%误差带来定义调节时间。
超调量%峰值h(t)超出终值h()的百分比,即
p
h(tp)h()
%100
h()
%
在上述动态性能指标中,工程上最常用的是调节时间ts(描述“快”),超调量%(描
述“匀”)以及峰值时间tp。
2稳态性能指标
稳态误差是时间趋于无穷时系统实际输出与理想输出之间的误差,是系统控制精度或抗
干扰能力的一种度量。
稳态误差有不同定义,通常在典型输入下进行测定或计算。
3
一阶系统的阶跃响应
一.一阶系统的数学模型
由一阶微分方程描述的系统,称为一阶系统。
一些控制元部件及简单系统如RC网络、
发电机、空气加热器、液面控制系统等都是一阶系统。
因为单位阶跃函数的拉氏变换为R(s)=1/s,故输出的拉氏变换式为
C(s)(s)R(s)
1
Ts
1
1
s
1
s
T
Ts
1
取C(s)的拉氏反变换得
c(t)1e
1
T
t
或写成
c(t)
cssc
tt
式中,css=1,代表稳态分量;
ctte
1
T
t
代表暂态分量。
当时间t趋于无穷,暂态分
量衰减为零。
显然,一阶系统的单位阶跃响应曲线是一条由零开始,按指数规律上升并最终
趋于1的曲线,如图所示。
响应曲线具有非振荡特征,故又称为非周期响应。
一阶系统的单位阶跃响应
二阶系统的阶跃响应
典型二阶系统方框图,其闭环传递函数为:
s
CsK/s(Ts1)K
vmv
2
Rs1K/s(Ts1)Ts
sK
v
mm
v
2
s
2
n
2
2nsn
式中
Kv--开环增益;
ωn--无阻尼自然频率或固有频率,
n
K
v
T
m
;
ζ--阻尼比,
2
1
nT
m
。
二阶系统的闭环特征方程为
s2+2ζωns+ω2n=0
其特征根为
4
2
s1
1,2
n
1.临界阻尼(ζ=1)
其时域响应为
t
ct1e(1nt)
n
上式包含一个衰减指数项。
c(t)为一无超调的单调上升曲线,如图3-8b所示。
(a)(b)(c)
ζ≥1时二阶系统的特征根的分布与单位阶跃响应
2.过阻尼(ζ>1)
2具有两个不同负实根[1,s
(1)n]
s的惯性环节单位阶跃响应拉氏变换
2
式。
其时域响应必然包含二个衰减的指数项,其动态过程呈现非周期性,没有超调和振荡。
图为其特征根分布图。
3.欠阻尼(0<ζ<1)
图3-90<ζ<1时二阶系统特征根的分布图3-10欠阻尼时二阶系统的单位阶跃响应
4.无阻尼(ζ=0)
2
Cs
2
s(s
n
2
)n
其时域响应为
ct1cos
t
n
在这种情况下,系统的响应为等幅(不衰减)振荡,
5
图ζ=0时特征根的分布图ζ=0时二阶系统的阶跃响应
5.负阻尼(ζ<0)
当ζ<0时,特征根将位于复平面的虚轴之右,其时域响应中的e的指数将是正的时间
函数,因而
t
e为发散的,系统是不稳定的。
n
显然,ζ≤0时的二阶系统都是不稳定的,而在ζ≥1时,系统动态响应的速度又太慢,
所以对二阶系统而言,欠阻尼情况是最有实际意义的。
下面讨论这种情况下的二阶系统的动
态性能指标。
欠阻尼二阶系统的动态性能指标
1.上升时间tr
上升时间tr是指瞬态响应第一次到达稳态值所需的时间。
t
r
d
n
1
2
由此式可见,阻尼比ζ越小,上升时间tr则越小;ζ越大则tr越大。
固有频率ωn越大,
tr越小,反之则tr越大。
2.峰值时间tp及最大超调量Mp
t
p
2
1
d
n
2
(/1)
最大超调量
Mp
c
max
c()e
cc()
2
max(/1)
最大超调百分数.100%
%ec
c()
3.调整时间ts
t
s
(5%)
1
13
2
[3
ln(1
)]
2
nn
,
0
0.707
t
s
(2%)
1
14
2
[
4
ln(1
)]
2
nn
,
0
0.707
图3-13二阶系统单位阶跃响应的一对包络线图3-14调节时间和阻尼比的近似关系
6
根据以上分析,二阶振荡系统特征参数ζ和ωn与瞬态性能指标(δ
4.振荡次数μ
在调整时问ts之内,输出c(t)波动的次数称为振荡次数μ,显然
t
s
t
f
式中
22
t,称为阻尼振荡的周期时间。
f
2
1
d
n
s
2T
1
2STS
2
2
1
这一系统的单位阶跃响应瞬态特性指标为:
最大超调百分数
2
(/1)
e100%4.3%%
上升时间
t
r4.7T
2
1
n
调整时间
ts2%8.43T(用近似式求得为8T)
ts5%4.14T(用近似式求得为6T)
有一位置随动系统其中Kk=4。
求该系统的(1)固有频率;
(2)阻尼比;(3)超调量和调
整时间;(4)如果要求实现工程最佳参数ζ=l/2,开环放大系数
k值应是多少?
k
【解】系统的闭环传递函数为
s
K
k
24
K
k
ssK
k
与二阶系统标准形式的传递函数
s
2
s
2
n
2
2nsn
对比得:
(1)固有频率nKk42
1
(2)阻尼比由2n1得0.25
2
n
2n
(/1)
(3)超调%e100%47%
3
(4)调整时间ts5%6s
n
当要求
1
2
12
时,由2n1得n,K0.5
kn
2
可见该系统要满足工程最佳参数的要求,须降低开环放大系数
K的值。
但是,降低Kk
k
值将增大系统的误差。
7
劳斯稳定判据
将系统的特征方程式写成如下标准式
nasasasa
n1n2
as
012n1n
0
将各系数组成如下排列的劳斯表
n
s
n
s
1
a
0
a
1
a
2
a
3
a
4
a
5
a
6
a
7
n
s
2
b
1
b
2
b
3
b
4
n
s
3
c
1
c
2
c
3
c
4
2
s
e
1
e
2
1
s
f
1
o
s
g
1
表中的有关系数为
aa
12
b
1a
1
aa
03
aa
14
b
2a
1
aa
05
aaaa
1607
b
3a
1
系数
b的计算,一直进行到其余的b值全部等于零为止。
i
ba
13
c
1b
1
ab
12
ba
15
c
2b
1
ab
13
ba
17
c
3b
1
ab
14
这一计算过程,一直进行到n行为止。
为了简化数值运算,可以用一个正整数去除或
乘某一行的各项,这时并不改变稳定性的结论。
(l)第一列所有系数均不为零的情况第一列所有系数均不为零时,劳斯判据指出,
特征方程式的实部为正实数根的数目等于劳斯表中第一列的系数符号改变的次数。
方程式的
根全部在复平面的左半平面的充分必要条件是,方程式的各项系数全部为正值,并且劳斯表
8
的第一列都具有正号。
例如,三阶系统的特征方程式为
32
a0sasasa
123
0
列出劳斯表为
3
s
a
0
a
2
2
s
a
1
a
3
1
s
aa
12
a
1
aa
03
0
s
a
3
则系统稳定的充分必要条件是
a0,a10,a20,a30,(a1a2a0a3)0
0
系统的特征方程为
5sss2s
43
s234
50
试用劳斯判据判断系统的稳定性。
解计算劳斯表中各元素的数值,并排列成下表
5
s
114
4
s
235
3
s
130
2
s
950
1
s
32
0
s
5
由上表可以看出,第一列各数值的符号改变了两次,由+2变成-1,又由-1改变成+9。
因此该系统有两个正实部的根,系统是不稳定的。
(2)某行第一列的系数等于零而其余项中某些项不等于零的情况在计算劳斯表中的
各元素的数值时,如果某行的第一列的数值等于零,而其余的项中某些项不等于零,那么可
以用一有限小的数值ε来代替为零的那一项,然后按照通常方法计算阵列中其余各项。
如果
零(ε)上面的系数符号与零(ε)下面的系数符号相反,则表明这里有一个符号变化。
例如,对于下列特征方程式
4sss
32
s22
10
劳斯表为
4
s
111
3
s
220
2
s
(0)1
1
s
2
2
0
s
1
现在观察第一列中的各项数值。
当ε趋近于零时,
2
2的值是一很大的负值,因此可
以认为第一列中的各项数值的符号改变了两次。
由此得出结论,该系统特征方程式有两个根
具有正实部,系统是不稳定的。
如果零(ε)上面的系数符号与零(ε)下面的系数符号不变,则表示系统有纯虚根。
例如,
9
对下列特征方程式
3s2s
s2
20
劳斯表为
3
s
11
2
1
s
s
22
0
s
2
可以看出,第一列各项中ε的上面和下面的系数符号不变,故有一对虚根。
将特征方程
式分解,有
(
2s
s1)(
2)
0
解得根为
p1j1,p32
2
(3)某行所有各项系数均为零的情况如果劳斯表中某一行的各项均为零,或只有等于
零的一项,这表示在s平面内存在一些大小相等但符号相反的特征根。
在这种情况下,
可利用全零行的上一行各系数构造一个辅助方程,式中s均为偶次。
将辅助方程对s求导,
用所得的导数方程系数代替全零行,然后继续计算下去。
至于这些大小相等,符号相反的根,
可以通过解辅助方程得到。
系统特征方程式为
6sssss
5432
s28122016
160
试用劳斯判据判断系统的稳定性。
解劳斯表中的
6
s~
3
s各项为
6
s
182016
5
s
212160
4
s
168
3
s
000
由上表可以看出,
3
s行的各项全部为零。
为了求出
3
s-
0
s各项,将
4
s行的各项组成辅
助方程为
A(s)
4s2
s6
8
将辅助方程A(s)对s求导数得
dA(s)3
4s
ds
12s
用上式中的各项系数作为
3
s行的各项系数,并计算以下各行的各项系数,得劳斯表为
6
s
182016
s
s
212160
4
s
168
3
s
412
2
s
1
s
0
s
3
4
3
8
8
10
从上表的第一列可以看出,各项符号没有改变,因此可以确定在右半平面没有特征方程
式的根。
另外,由于
3
s行的各项皆为零,这表示有共轭虚根。
这些根可由辅助方程求出。
本例中的辅助方程式是
4s
2
s6
80
由之求得特征方程式的大小相等符号相反的虚根为
p1j2,p3,4j2,p5,61j2
2
稳态误差及其计算
误差本身是时间t的函数,在时域中以et表示。
稳定系统误差的终值称为稳态误差
e,
ss
即为误差信号的稳态分量,则稳态误差为
e
ss
lim
t
et)
(lim
s0
sE(s)
系统的误差传递函数
Es1
Rs1G(s)H(s)
故
Es
Rs
1G(s)H(s)
将系统误差的拉氏变换E(s)代入(3-38),得稳态误差的计算公式为
e
ss
sR(s)
lim
01G(s)H(s)
s
控制系统的型别
控制系统的一般开环传递函数可以写成
m
K
k
(Ts
i
1)
G(s)H(s)
n
i
1
N
N
s
(Ts1)
j
j1
式中
K为开环放大系数或称为开环传递系数;Ti、Tj为时间常数;N表示开环传递函
k
数中串联的积分环节个数。
这是一个很重要的结构参数。
根据N的数值,可将系统分为几
种不同类型。
N=0的系统称为0型系统;N=1的系统称为I型系统;N=2的系统称为II
型系统。
当N>2时,要使系统稳定是很困难的。
因此,一般采用的是0型、I型和II型系
统。
典型输入下系统的稳态误差
对于不同输入函数,下面分析系统的稳态误差。
1.单位阶跃输入下的稳态误差
单位阶跃输入(
Rs
1
s
)下的系统稳态误差,由式(3-40)得
e
ss
s
lim
01G(s)H(s)
s
1
s
1
1
G(s)
H
(s)
11
定义
k
p
lim
G(s)H(
s0
s)
k称为位置误差系数,则
p
1
e
ssK
1
p
0型系统的稳态误差为
1
e
ssK
1
k
I型或高于I型的系统的位置稳态误差为
e
ss
0
2.单位斜坡输入下的稳态误差
单位斜坡输入(
1
Rs)的系统稳态误差
2
s
e
ss
s11
lim
2
s0GsHsssGsHs
0GsHsssGsHs
1()()lim()()
s0
定义
K
lim
sG(s)H
s0
(s)
K称为速度误差系数。
则
ess
1
K
对于0型系统
m
K
k
(Ts
i
1)
k
i1
n
lim
s0s0Ts
(
j
j1
1)
0
所以ess
对于I型系统
m
K
k
(Ts
i
1)
k
i1
lim
s
n
s
0
(1)
sTs
j
j1
K
k
所以
11
e
sskK
k
对于II型或更高型系统
12
m
K
k
(Ts
i
1)
k
i1
n
lims
s0s2(Ts1)
j
j1
所以e0
ss
,0型系统不能跟踪斜坡输入;单位反馈的Ⅰ型系统能跟踪斜坡输入,但总有一
定误差,
3.单位抛物线输入下的稳态误差
定义
K
a
2
limsG(s)H(s)
s0
K称为加速度误差系数
a
则
1
e
ssK
a
m
K
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 系统 阶跃 响应 doc