第三章线性系统的时域分析法

第三章线性系统的时域分析法 3.1 典型输入信号和时域指标 3.2 一阶系统的时域响应 3.3 二阶系统的时域响应 3.4 高阶系统分析 3.5 稳定性分析 3.6 稳态误差分析

3.1典型输入信号和时域指标 分析和设计控制系统的首要工作是确定系统的数模，一旦获得系统的数学模型，就可以采用几种不同的方法去分析系统的性能。线性系统：时域分析法，根轨迹法，频率法非线性系统：描述函数法，相平面法采样系统： z变换法状态空间法多输入多输出系统：

对线性系统，时域分析法的要点是： （1）建立数模（微分方程式，传递函数）（2）选择合适的输入函数（典型信号）。取决于系统常见工作状态，同时，在所有的可能的输入信号中，选取最不利的信号作为系统的典型输入信号。（3）求出系统输出随时间变化的关系 C(s) = G(s)R(s) c(t) = L−1[C(s)] （4）根据时间响应确定系统的性能，包括稳定性快速性和准确性等方面指标，看这些指标是否符合生产工艺的要求。

1(t) f(t) (t) t t 0 0 t 0 目前，常用的典型外作用有以下几种： 1 （1）单位阶跃函数其数学表达式为（2）单位斜坡函数其数学表达式为（3）单位脉冲函数其数学表达式为

f(t) f(t) t t 0 0 （4）单位匀加速函数其数学表达式为（5）正弦函数其数学表达式为 f(t) = Asinωt

任何一个实际控制系统的时间响应，都由过渡过程和稳态过程两部分组成:任何一个实际控制系统的时间响应，都由过渡过程和稳态过程两部分组成: （1）过渡过程：系统从刚加入输入信号后，到系统输出量达到稳态值前的响应过程，称为过渡过程或动态过程。在这一期间，由于系统具有惯性、摩擦以及其它一些原因，输出量不可能完全复现输入量的变化。根据结构和参数选择情况，过渡过程表现为衰减、发散或等幅振荡形式，如图所示。显然，一个可以运行的控制系统，其过渡过程必须是衰减的（稳定的）。（2）稳态过程：时间 t 趋于无穷大时的响应过程，稳态过程表征输出量最复现输入量的程度，用稳态性能描述。

c(t) c(t) c(t) t t t 0 0 0

0.05 或 0.02 p 1 c(t) 0.5 td tr tp ts t 0 用tr,tp, p,ts四个性能指标来衡量瞬态响应的好坏。

R r(t) c(t) C R(s) 1 Ts C(s) + ﹣ 3.2 一阶系统的时域分析凡是可用一阶微分方程描述的系统，称为一阶系统。 T＝RC，时间常数。其典型结构图及传递函数为：

c(t) t 3.2.1 单位阶跃响应当输入信号r(t)=1(t)时，系统的响应c(t)称作其单位阶跃响应。 1.0 响应曲线在[0，）的时间区间中始终不会超过其稳态值，把这样的响应称为非周期响应。 0.982 0.95 0.865 0.632 ０　 T2T 3T 4T

c(t) 1.0 0.982 0.95 0.865 0.632 t ０　 T2T 3T 4T 一阶系统响应具备两个重要的特点： ① 可以用时间常数T去度量系统输出量的数值。 ② 响应曲线的初始斜率等于1/T。一阶系统的瞬态响应指标调整时间ts 定义：︱c(ts) 1 ︱=  ( 取5%或2%)

c(t) t 0 3.2.2 单位斜坡响应[ r(t) = t ] T r(t)= t c(t) = t ﹣T + Te﹣t/T T 稳态响应是一个与输入斜坡函数斜率相同但在时间上迟后了一个时间常数T的斜坡函数。

c(t) 表明过渡过程结束后，其稳态输出与单位斜坡输入之间，在位置上仍有误差，一般叫做跟踪误差。比较阶跃响应曲线和斜坡响应曲线： 1.0 T r(t)= t t c(t) 0 T t 0 在阶跃响应中，输出量与输入量之间的位置误差随时间而减小，最终趋于0，而在初始状态下，位置误差最大，响应曲线的斜率也最大；在斜坡响应中，输出量与输入量之间的位置误差随时间而增大，最终趋于常值T，在初始状态下，位置误差和响应曲线的斜率均等于0。其原因在稳态误差的计算中说明。

h(t) 1/T 0.368/T 0.135/T 0.05/T t 0 T 2T 3T 3.2.3　单位脉冲响应 [R(s)=1] 它恰是系统的闭环传函，这时输出称为脉冲响应函数，以h(t)标志。求系统闭环传函提供了实验方法，以单位脉冲输入信号作用于系统，测定出系统的单位脉冲响应，可以得到闭环传函。对应

线性定常系统的重要性质 1.当系统输入信号为原来输入信号的导数时，这时系统的输出则为原来输出的导数。 2. 在零初始条件下，当系统输入信号为原来输入信号时间的积分时，系统的输出则为原来输出对时间的积分，积分常数由零初始条件决定。

C(s) n2 R(s) s(s+2n) + ﹣ 3.3 二阶系统的时域分析 3.3.1二阶系统单位阶跃响应 1. 二阶系统的数学模型标准化二阶系统的结构图为：　闭环传递函数为二阶系统有两个结构参数(阻尼比)和n(无阻尼振荡频率)。二阶系统的性能分析和描述，都是用这两个参数表示的。

L R c(t) r(t) C 例如RLC电路微分方程式为：对于不同的二阶系统，阻尼比和无阻尼振荡频率的含义是不同的。

2　二阶系统的闭环极点 j  0 二阶系统的闭环特征方程，即 s 2 + 2n s + n2 = 0 其两个特征根为：上述二阶系统的特征根表达式中，随着阻尼比 的不同取值，特征根有不同类型的值，或者说在s平面上有不同的分布规律。分述如下： (1)  >1 时，特征根为一对不等值的负实根，位于s 平面的负实轴上，使得系统的响应表现为过阻尼的。 s1 s2

j j  0  0 (2) =1时，特征根为一对等值的负实根，位于s 平面的负实轴上，使得系统的响应表现为临界阻尼的。 (3)0 < < 1时，特征根为一对具有负实部的共轭复根，位于s平面的左半平面上，使得系统的响应表现为欠阻尼的。 jd s1 n  n s1= s2 = n s2

j j   0 0 (4)  = 0 时，特征根为一对幅值相等的虚根，位于s平面的虚轴上，使得系统的响应表现为无阻尼的等幅振荡过程。 (5)  < 0 时，特征根位于s平面的右半平面，使得系统的响应表现为幅值随时间增加而发散。 jn

j j s1 s2 jd s1 n jn j j    n 0 0 s1=s2   0 0 s2 阻尼比取不同值时，二阶系统根的分布

3. 单位阶跃响应 由式，其输出的拉氏变换为式中s1，s2是系统的两个闭环特征根。对上式两端取拉氏反变换，可以求出系统的单位阶跃响应表达式。阻尼比在不同的范围内取值时，二阶系统的特征根在s 平面上的位置不同，二阶系统的时间响应对应有不同的运动规律。下面分别加以讨论。

j jd s1 n    n 0 s2 （1）欠阻尼情况

c(t) c(t) t t 0 0 1 欠阻尼二阶系统的单位阶跃响应由两部分组成：稳态分量为1，表明系统在1(t)作用下不存在稳态位置误差；瞬态响应是阻尼正弦项，其振荡频率为阻尼振荡频率，而其幅值则按指数曲线衰减，两者均由参数和n决定。无阻尼情况 c(t) = 1  cosnt ( t > 0)

c(t) t 0 1 （2）临界阻尼情况 s1,2= n 此时响应是稳态值为1 的非周期上升过程，其变化率 t = 0，变化率为0； t > 0变化率为正，c(t) 单调上升； t →∞，变化率趋于0。整个过程不出现振荡。

（3）过阻尼情况 响应特性包含两个单调衰减的指数项，且它们的代数和不会超过1，因而响应是非振荡的。（不同于一阶系统）

横坐标nt，曲线只是的函数。  =0,0.1,0.2,0.4,0.6,0.8,1,2

c(t) t 0 3.3.2 欠阻尼二阶系统的动态性能指标 0.05 或 0.02 p 1 0.5 td tr tp ts 用tr,tp, p,ts四个性能指标来衡量瞬态响应的好坏。

(1) 上升时间tr：从零上升至第一次到达稳态值所需的时间，是系统响应速度的一种度量。tr 越小，响应越快。 (2) 峰值时间tp：响应超过稳态值，到达第一个峰值所需的时间。

(3) 超调量p：响应曲线偏离阶跃曲线最大值，用百分比表示。

p只是 的函数，其大小与自然频率n无关。  = 0.2 p = 52.7%  = 0.4 p = 25.4%  = 0.6 p = 9.5%  = 0.707 p = 4.3% p (4) 调节时间ts：响应曲线衰减到与稳态值之差不超过5%所需要的时间。  c(t)  c()     c() ( t  ts )

工程上，当0.1 < < 0.9时，通常用下列二式近似计算调节时间。 △ = 5% △ = 2%

必须 (1) ωn一定，使tr tp  使 ts (一定范围 ) 必须必须 (2) 一定，使tr tp tsωn  (3) p只由 决定必须总结：各性能指标之间是有矛盾的。

K C(s) R(s) s(Ts+1) + ﹣ 例3-1　单位负反馈随动系统如图所示 (1) 确定系统特征参数与实际参数的关系。 (2) 若K = 16(rad/s)、T = 0.25(s)，试计算系统的动态性能指标。解: (1) 系统的闭环传递函数为与典型二阶系统比较可得： K/T= n2 1/T = 2n

(2) K = 16，T = 0.25时 ( =0.05 )

1.3 c(t) 1 0.1 t 0 例3-2 已知单位负反馈系统的单位阶跃响应曲线如图所示，试求系统的开环传递函数。解：由系统的单位阶跃响应曲线，直接求出超调量和峰值时间。 p = 30% tp = 0.1 求解上述二式，得到 = 0.357，n= 33.6 (rad/s)。于是二阶系统的开环传递函数为

3.3.5　二阶系统性能的改善 1. 误差的比例－微分控制具有误差比例－微分控制的二阶系统如图所示 + + Td s C(s) n2 R(s) s(s+2n) + ﹣ 系统的开环传递函数为闭环传递函数为式中d为系统的有效阻尼比。

1 t 0 上式表明，比例－微分控制的二阶系统不改变系统的自然频率，但是可以增大系统的有效阻尼比以抑制振荡。此时，相当于为系统增加了一个闭环零点。若令Z=1/Td，上式可以表示为比例－微分控制的二阶系统有时称为有零点的二阶系统。与没有零点的二阶系统相比，超调量会增大一些。 c1(t) 有零点的二阶系统。 c(t) 没有零点的二阶系统。 c1(t) c(t)

(1) (2)

C(s) R(s) n2 s(s+2n) Kf s + + ﹣ ﹣ 2.输出量的速度反馈控制系统的闭环传递函数为：式中　　　　　　为系统的有效阻尼比。显然，输出量的速度反馈控制也可以在不改变系统的自然频率基础上，增大系统的有效阻尼比，减小超调量。

与比例微分控制不同的是，输出量的速度反馈控制没有附加零点的影响，两者对系统动态性能的改善程度是不同的。与比例微分控制不同的是，输出量的速度反馈控制没有附加零点的影响，两者对系统动态性能的改善程度是不同的。 3.两种控制方案的比较都为系统提供了一个参数选择的自由度，兼顾了系统响应的快速性和平稳性。但是，二者改善系统性能的机理及其应用场合是不同的。简述如下：（1）微分控制的附加阻尼作用产生于系统输入端误差信号的变化率，而速度反馈控制的附加阻尼作用来源于系统输出量的变化率。微分控制为系统提供了一个实零点，可以缩短系统的初始响应时间，但在相同阻尼程度下，将比速度反馈控制产生更大的阶跃响应超调量。

(2) 比例控制位于系统的输入端，微分作用对输入噪声有明显的放大作用。当输入端噪声严重时，不宜选用比例－微分控制。同时，由于微分器的输入信号是低能量的误差信号，要求比例－微分控制具有足够的放大作用，为了不明显恶化信噪比，需选用高质量的前置放大器。输出速度反馈控制，是从高能量的输出端向低能量的输入端传递信号，无需增设放大器，并对输入端噪声有滤波作用，适合于任何输出可测的控制场合。

C(s) R(s) G(s) H(s) + ﹣ 3.4　高阶系统的时域分析 3.4.1　高阶系统的阶跃响应控制系统的基本结构如图所示。 G(s)，H(s)一般是复变量s 的多项式之比，故上式可记为其闭环传递函数为

根据能量的有限性，分子多项式的阶次m不高于分母多项式的阶次n。对上式进行因式分解，可以表示为根据能量的有限性，分子多项式的阶次m不高于分母多项式的阶次n。对上式进行因式分解，可以表示为式中0 < k <1。即系统有q 个实极点和r 对共轭复数极点。取拉氏变换，并设全部初始条件为零，得到系统单位阶跃响应的时间表达式：

式中　　　　　；k=arccosk；Ak、Bk是与C(s)在对应闭环极点上的留数有关的常数。上式表明，如果系统的所有闭环极点都具有负实部，系统时间响应的各暂态分量都将随时间的增长而趋近于零，这时称高阶系统是稳定的。 3.4.2　闭环主导极点 1）高阶系统瞬态响应各分量的衰减快慢由 pi ， kn决定，也即闭环极点负实部的绝对值越大，相应的分量衰减越快。 2）各分量所对应的系数由系统的零极点分布决定。当某一极点越靠近零点，而远离其他极点和原点，则相应系数越小，该瞬态分量的影响就越小；

当某一极点远离零点，越靠近其他极点和原点，则相应系数越大，该瞬态分量的影响就越大；当某一极点远离零点，越靠近其他极点和原点，则相应系数越大，该瞬态分量的影响就越大；一个零点和一个极点距离非常近，把这一对零极点称为偶极子。 3）系统的零极点共同决定了系统瞬态响应曲线的形状。根据上述，把系数很小的分量，远离虚轴衰减很快的分量常常忽略，高阶系统就可用低阶系统来近似估计。 4）对系统瞬态响应起主导作用的极点，称为主导极点。应用闭环主导极点的概念，可以把一些高阶系统近似为一阶或二阶系统，以实现对高阶系统动态性能的近似评估。

一般情况，高阶系统具有振荡性，所以主导极点常常是一对共轭复数极点。找到了一对共轭复数极点，高阶系统的动态性能就可以应用二阶系统的性能指标来近似估计。

c(t) t 0 例3-3 已知闭环传递函数为试求阶跃响应。解： c(t) = 1 1.1et + 0.11e10t ≈1 1.1et 主导极点是 s = 1，这时系统传递函数近似为 1

j  1.25 1 10 0 例3-4 已知闭环传递函数为试求阶跃响应。解： c(t) = 1 0.22et 0.78e10t

第三章 线性系统的时域分析法