控制系统典型的输入信号

控制系统典型的输入信号

1. 阶跃函数

阶跃函数的定义是

0, t0xr(t)

A, t0式中A为常数。A等于1的阶跃函数称为单位阶跃函数，如图所示。它表示为

xr(t)＝l(t)，或xr(t)=u(t)

单位阶跃函数的拉氏变换为

Xr(s)=L[1(t)]=1/s

在t＝0处的阶跃信号，相当于一个不变的信号突然加到系统上；对于恒值系统，相当于给定值突然变化或者突然变化的扰动量；对于随动系统，相当于加一突变的给定位置信号。

2. 斜坡函数

这种函数的定义是

0, t0 xr(t)A t, t0式中A为常数。该函数的拉氏变换是

Xr(s)=L[At]=A/s2

这种函数相当于随动系统中加入一按恒速变化的位置信号，该恒速度为A。当A＝l时，称为单位斜坡函数，如图所示。

1

3. 抛物线函数

如图 所示，这种函数的定义是

0, t0xr(t)

2A t, t0

式中A为常数。这种函数相当于随动系统中加入一按照恒加速变化的位置信号，该恒加速度为A。抛物线函数的拉氏变换是

Xr(s)=L[At2]=2A/s3

当A＝1/2时，称为单位抛物线函数，即Xr(s)=1/s3。

4. 脉冲函数

这种函数的定义是

A, 0 t (0)xr(t)

0, t0, t (0)式中A为常数，ε为趋于零的正数。脉冲函数的拉氏变换是

AXr(s)LlimA

0当A＝1，ε→0时，称为单位脉冲函数δ(t)，如图 所示。单位脉冲函数的面积等于l，即

2

(t)dt1

在t＝t0处的单位脉冲函数用δ(t-t0)来表示，它满足如下条件

幅值为无穷大、持续时间为零的脉冲纯属数学上的假设，但在系统分析中却很有用处。单位脉冲函数δ(t)可认为是在间断点上单位阶跃函数对时间的导数，即

反之，单位脉冲函数δ(t)的积分就是单位阶跃函数。

控制系统的时域性能指标

对控制系统的一般要求归纳为稳、准、快。工程上为了定量评价系统性能好坏，必须给出控制系统的性能指标的准确定义和定量计算方法。

1 动态性能指标

动态性能指标通常有如下几项：

延迟时间td 阶跃响应第一次达到终值h()的50％所需的时间。

上升时间tr 阶跃响应从终值的10％上升到终值的90％所需的时间；对有振荡的系统，也可定义为从0到第一次达到终值所需的时间。

峰值时间tp 阶跃响应越过稳态值h()达到第一个峰值所需的时间。

调节时间ts 阶跃响到达并保持在终值h()5％误差带内所需的最短时间；有时也用终值的2％误差带来定义调节时间。

超调量％ 峰值h(tp)超出终值h()的百分比，即

％h(tp)h()h()100％

在上述动态性能指标中，工程上最常用的是调节时间ts（描述“快”），超调量％（描述“匀”）以及峰值时间tp。

2 稳态性能指标

稳态误差是时间趋于无穷时系统实际输出与理想输出之间的误差，是系统控制精度或抗干扰能力的一种度量。稳态误差有不同定义，通常在典型输入下进行测定或计算。

3

一阶系统的阶跃响应

一. 一阶系统的数学模型

由一阶微分方程描述的系统，称为一阶系统。一些控制元部件及简单系统如RC网络、发电机、空气加热器、液面控制系统等都是一阶系统。

因为单位阶跃函数的拉氏变换为R(s)=1/s，故输出的拉氏变换式为

C(s)(s)R(s)取C(s)的拉氏反变换得

111T Ts1ssTs11tTc(t)1e或写成

c(t)cssctt

式中，css=1，代表稳态分量；ctte1tT代表暂态分量。当时间t趋于无穷，暂态分

量衰减为零。显然，一阶系统的单位阶跃响应曲线是一条由零开始，按指数规律上升并最终趋于1的曲线，如图所示。响应曲线具有非振荡特征，故又称为非周期响应。

一阶系统的单位阶跃响应

二阶系统的阶跃响应

典型二阶系统方框图，其闭环传递函数为：

2Kv/s(Tms1)KvCsn s2Rs1Kv/s(Tms1)Tms2sKvs22nsn式中

Kv--开环增益；

ωn--无阻尼自然频率或固有频率，n ζ--阻尼比，Kv； Tm1。

2nTms2+2ζωns+ω2n=0

二阶系统的闭环特征方程为 其特征根为

2s1,21n

 4

1. 临界阻尼(ζ=1)

其时域响应为

ct1ent(1nt)

上式包含一个衰减指数项。c(t)为一无超调的单调上升曲线，如图3-8b所示。

(a) (b) (c)

ζ≥1时二阶系统的特征根的分布与单位阶跃响应

2. 过阻尼(ζ＞1)

具有两个不同负实根[s1,s2(21)n]的惯性环节单位阶跃响应拉氏变换

式。其时域响应必然包含二个衰减的指数项，其动态过程呈现非周期性，没有超调和振荡。图为其特征根分布图。

3. 欠阻尼（0<ζ＜1）

图3-9 0<ζ＜1时二阶系统特征根的分布 图3-10 欠阻尼时二阶系统的单位阶跃响应

4. 无阻尼(ζ＝0)

Cs其时域响应为

2ns(s)22n

ct1cosnt

在这种情况下，系统的响应为等幅(不衰减)振荡，

5

图ζ＝0时特征根的分布 图ζ＝0时二阶系统的阶跃响应

5. 负阻尼（ζ＜0）

当ζ＜0时，特征根将位于复平面的虚轴之右，其时域响应中的e的指数将是正的时间函数，因而en为发散的，系统是不稳定的。

显然，ζ≤0时的二阶系统都是不稳定的，而在ζ≥1时，系统动态响应的速度又太慢，所以对二阶系统而言，欠阻尼情况是最有实际意义的。下面讨论这种情况下的二阶系统的动态性能指标。

t欠阻尼二阶系统的动态性能指标

1. 上升时间tr

上升时间tr是指瞬态响应第一次到达稳态值所需的时间。 tr 2dn1由此式可见，阻尼比ζ越小，上升时间tr则越小；ζ越大则tr越大。固有频率ωn越大，

tr越小，反之则tr越大。

2. 峰值时间tp及最大超调量Mp

tp 2dn1

(/12)最大超调量 Mpcmaxc()e最大超调百分数 c%

cmaxc()e(/c()12).100%

3. 调整时间ts

ts(5%)13[3ln(12)]， 00.707 n2n1114ts(2%)[4ln(12)]， 00.707 n2n

图3-13 二阶系统单位阶跃响应的一对包络线 图3-14 调节时间和阻尼比的近似关系

根据以上分析，二阶振荡系统特征参数ζ和ωn与瞬态性能指标(δ

6

4. 振荡次数μ

在调整时问ts之内，输出c(t)波动的次数称为振荡次数μ，显然

式中 tfts tf2d2n12，称为阻尼振荡的周期时间。

s1 222TS2TS1这一系统的单位阶跃响应瞬态特性指标为： 最大超调百分数

%e上升时间

(/12)100%4.3%

tr调整时间

n124.7T

ts2%8.43T(用近似式求得为8T) ts5%4.14T(用近似式求得为6T)

有一位置随动系统其中Kk＝4。求该系统的(１)固有频率；(2)阻尼比；(3)超调量和调整时间；(4)如果要求实现工程最佳参数ζ＝l／2，开环放大系数kk值应是多少？

【解】系统的闭环传递函数为

sKk Kk4

s2sKk与二阶系统标准形式的传递函数

2n s22s2nsn对比得：(1) 固有频率

nKk42

12n0.25

(2) 阻尼比 由2n1得 (/ (3) 超调 %e12)n100%47%

(4) 调整时间ts5%当要求3n6s

122,Kkn0.5

12时，由2n1 得 n可见该系统要满足工程最佳参数的要求，须降低开环放大系数Kk的值。但是，降低Kk值将增大系统的误差。

7

劳斯稳定判据

将系统的特征方程式写成如下标准式

a0sna1sn1a2sn2an1san0 将各系数组成如下排列的劳斯表

snsn1sn2sn3s2s1soa0a1b1c1e1f1g1a2a3b2c2e2a4a5b3c3a6a7b4c4

表中的有关系数为

b1b2b3a1a2a0a3

a1a1a4a0a5

a1a1a6a0a7

a1



系数bi的计算，一直进行到其余的b值全部等于零为止。

c1c2c3b1a3a1b2

b1b1a5a1b3

b1b1a7a1b4

b1

这一计算过程，一直进行到 n行为止。为了简化数值运算，可以用一个正整数去除或乘某一行的各项，这时并不改变稳定性的结论。

(l) 第一列所有系数均不为零的情况 第一列所有系数均不为零时，劳斯判据指出，特征方程式的实部为正实数根的数目等于劳斯表中第一列的系数符号改变的次数。方程式的根全部在复平面的左半平面的充分必要条件是，方程式的各项系数全部为正值，并且劳斯表的第一列都具有正号。

8

例如， 三阶系统的特征方程式为

a0s3a1s2a2sa30

列出劳斯表为

s3s2s1s0则系统稳定的充分必要条件是

a0a1a1a2a0a3a1a3a2a3

a00，a10，a20，a30，(a1a2a0a3)0

系统的特征方程为

s2ss3s4s50

试用劳斯判据判断系统的稳定性。

解 计算劳斯表中各元素的数值，并排列成下表

32s5s4s3s2129325143550

130s1s0由上表可以看出，第一列各数值的符号改变了两次，由+2变成-1，又由-1改变成+9。因此该系统有两个正实部的根，系统是不稳定的。

(2) 某行第一列的系数等于零而其余项中某些项不等于零的情况 在计算劳斯表中的各元素的数值时，如果某行的第一列的数值等于零，而其余的项中某些项不等于零，那么可以用一有限小的数值ε来代替为零的那一项，然后按照通常方法计算阵列中其余各项。如果零(ε)上面的系数符号与零(ε)下面的系数符号相反，则表明这里有一个符号变化。

例如，对于下列特征方程式

s42s3s22s10

劳斯表为

s4s3s2s1s0111220(0)1

221现在观察第一列中的各项数值。当ε趋近于零时，22的值是一很大的负值，因此可

以认为第一列中的各项数值的符号改变了两次。由此得出结论，该系统特征方程式有两个根具有正实部，系统是不稳定的。

如果零(ε)上面的系数符号与零(ε)下面的系数符号不变，则表示系统有纯虚根。例如，对下列特征方程式

9

s32s2s20

劳斯表为

s3s2s1s011222

可以看出，第一列各项中ε的上面和下面的系数符号不变，故有一对虚根。将特征方程式分解，有

(s21)(s2)0 解得根为

p1,2j1， p32

(3) 某行所有各项系数均为零的情况 如果劳斯表中某一行的各项均为零，或只有等于零的一项， 这表示在 s 平面内存在一些大小相等但符号相反的特征根。 在这种情况下，可利用全零行的上一行各系数构造一个辅助方程，式中s均为偶次。将辅助方程对s求导，用所得的导数方程系数代替全零行，然后继续计算下去。至于这些大小相等，符号相反的根，可以通过解辅助方程得到。

系统特征方程式为

s62s58s412s320s216s160

试用劳斯判据判断系统的稳定性。 解 劳斯表中的s～s各项为

s663182016

s5 s4s33212160 168000304由上表可以看出，s行的各项全部为零。为了求出s-s各项，将s行的各项组成辅助方程为

A(s)s46s28

将辅助方程A(s)对s求导数得

dA(s)4s312s ds用上式中的各项系数作为s行的各项系数，并计算以下各行的各项系数，得劳斯表为

3s6ss 4ss3s2 s1182016212160

168412384 38s0从上表的第一列可以看出，各项符号没有改变，因此可以确定在右半平面没有特征方程

10

式的根。另外，由于s行的各项皆为零，这表示有共轭虚根。这些根可由辅助方程求出。本例中的辅助方程式是

3s46s280

由之求得特征方程式的大小相等符号相反的虚根为

p1,2j2， p3,4j2， p5,61j2

稳态误差及其计算

误差本身是时间t的函数，在时域中以et表示。稳定系统误差的终值称为稳态误差ess，即为误差信号的稳态分量，则稳态误差为

ess系统的误差传递函数

lime(t)limsE(s)

ts0Es1 Rs1G(s)H(s)Rs故 Es 1G(s)H(s)

将系统误差的拉氏变换E(s)代入(3-38)，得稳态误差的计算公式为 esslims0sR(s)

1G(s)H(s)控制系统的型别

控制系统的一般开环传递函数可以写成

G(s)H(s)Kk(Tis1)sN(Tjs1)j1i1nNm

式中Kk为开环放大系数或称为开环传递系数；Ti、Tj为时间常数；N表示开环传递函数中串联的积分环节个数。这是一个很重要的结构参数。根据N的数值，可将系统分为几种不同类型。N＝0的系统称为0型系统；N＝1的系统称为I型系统；N=2的系统称为II型系统。当N＞2时，要使系统稳定是很困难的。因此，一般采用的是0型、I型和II型系统。

典型输入下系统的稳态误差

对于不同输入函数，下面分析系统的稳态误差。

1. 单位阶跃输入下的稳态误差

1)下的系统稳态误差，由式(3-40)得 ss11 esslim

1G(s)H(s)s1G(s)H(s)s0单位阶跃输入(Rs定义

11

kplimG(s)H(s)

s0kp称为位置误差系数，则

ess0型系统的稳态误差为 ess1

1Kp1

1KkI型或高于I型的系统的位置稳态误差为

ess0

2. 单位斜坡输入下的稳态误差

单位斜坡输入(Rs1)的系统稳态误差 s2s1esslim2ss01G(s)H(s)lims01

sG(s)H(s)定义

KlimsG(s)H(s)

s0K称为速度误差系数。

则 ess对于0型系统

1 Kklims0Kk(Tis1)s0(Tjs1)j1i1nm0

所以 ess 对于I型系统

mklimss0Kk(Tis1)s(Tjs1)j1i1nKk

所以 ess对于II型或更高型系统

11 kKk 12

klimss0Kk(Tis1)s2(Tjs1)j1i1nm

所以 ess0

，0型系统不能跟踪斜坡输入；单位反馈的Ⅰ型系统能跟踪斜坡输入，但总有一

定误差，

3. 单位抛物线输入下的稳态误差

定义

Kalims2G(s)H(s)

s0Ka称为加速度误差系数

则 ess1 Ka对于0型系统 kalims0s2Kk(Tis1)s0(Tjs1)j1i1nm0

所以 ess 对于Ⅰ型系统

mkalims2s0Kk(Tis1)s1(Tjs1)j1i1n0

所以 ess 对于Ⅱ型系统

mkalims2s0Kk(Tis1)s1(Tjs1)j1i1nKk

所以 ess1 kk由此可知，0型和Ⅰ型系统都不能跟踪抛物线输入，Ⅱ型系统能跟踪抛物线输入，但存在稳态误差。

典型输入信号作用下的稳态误差和误差系数

系统型别 静态误差系数 阶跃输入r(t)位置误差eA1(t) A1Kp 斜坡输入r(t)At AKv Kp Kv Ka ss速度误差essAt2加速度输入r(t)2A 加速度误差essKa 13

0 KA 0 0 1KI  K 0 0 II   K 0 。

减小稳态误差的方法

1. 引入给定量顺馈 2. 引入扰动量顺馈

  A K 0 A K

14