自动控制原理 时域分析法

第三章时域分析法控制系统的数学模型，是分析、研究和设计控制系统的基础，经典控制论中三种分析（时域，根轨迹，频域）、研究和设计控制系统的方法，都是建立在这个基础上。1北京航空航天大学3－1时域分析基础

一、时域分析法的特点它根据系统微分方程，通过拉氏变换，直接求出系统的时间响应。依据响应的表达式及时间响应曲线来分析系统控制性能，并找出系统结构、参数与这些性能之间的关系。这是一种直接方法，而且比较准确，可以提供系统时间响应的全部信息。2北京航空航天大学二、典型初始状态，典型外作用1.典型初始状态通常规定控制系统的初始状态为零状态。即在外作用加于系统之前，被控量及其各阶导数相对于平衡工作点的增量为零，系统处于相对平衡状态。3北京航空航天大学2.典型外作用单位阶跃函数1(t)tf(t)0îíì<³==0t00t1)t(1)t(f其拉氏变换为：s1dte1)s(F)]t(f[L0st===ò¥-其数学表达式为：4北京航空航天大学t

单位斜坡函数0t0t0t)t(1t)t(f<³îíì=.=其拉氏变换为：20sts1dtet)s(F)]t(f[L===ò¥-f(t)0其数学表达式为：5北京航空航天大学单位脉冲函数000)()(=¹îíì¥==ttttfd其数学表达式为：其拉氏变换为：1)()]([==sFtfLò+¥¥-=1)(dttd定义：图中1代表了脉冲强度。单位脉冲作用在现实中是不存在的，它是某些物理现象经数学抽象化的结果。6北京航空航天大学正弦函数其拉氏变换为：220sin)()]([ωsωdteωtsFtfLst+===ò¥-000sin)(<³îíì=ttωttf其数学表达式为：f(t)7北京航空航天大学三、典型时间响应初状态为零的系统，在典型输入作用下输出量的动态过程，称为典型时间响应。8北京航空航天大学1.单位阶跃响应定义：系统在单位阶跃输入[r(t)=1(t)]作用下的响应，常用h(t)表示。若系统的闭环传函为，则h(t)的拉氏变换为故9北京航空航天大学2.单位斜坡响应定义：系统在单位斜坡输入[r(t)=t·1(t)]作用下的响应，常用表示。故则有10北京航空航天大学3.单位脉冲响应定义：系统在单位脉冲输入[r(t)=δ(t)]作用下的响应，常用k(t)表示。注：关于正弦响应，将在第五章里讨论故则有11北京航空航天大学4.三种响应之间的关系由式(3-1-3)可(3-1-1),(3-1-2)写为相应的时域表达式为12北京航空航天大学四、阶跃响应的性能指标t)(th)(pth1ptst误差带013北京航空航天大学1、峰值时间tp：指h(t)曲线中超过其稳态值而达到第一个峰值所需的时间。2、超调量％：指h(t)中对稳态值的最大超出量与稳态值之比。3、调节时间ts：指响应曲线中，h(t)进入稳态值附近5%h()或2%h()误差带，而不再超出的最小时间。4、稳态误差ess：指响应的稳态值与期望值之差。14北京航空航天大学注意事项：15北京航空航天大学3－2一、二阶系统分析与计算定义：由一阶微分方程描述的系统称为一阶系统。一、一阶系统的数学模型及单位阶跃响应16北京航空航天大学一阶系统数学模型微分方程：动态结构图：传递函数：17北京航空航天大学一阶系统单位阶跃响应输入：输出：18北京航空航天大学单位阶跃响应曲线初始斜率：19北京航空航天大学性能指标1.平稳性：2.快速性ts：3.准确性ess：非周期、无振荡，＝020北京航空航天大学举例说明（一阶系统）一阶系统如图所示，试求：当KH＝0.1时，求系统单位阶跃响应的调节时间ts，放大倍数K，稳态误差ess；如果要求ts＝0.1秒，试问系统的反馈系数KH应调整为何值？讨论KH的大小对系统性能的影响及KH与ess的关系。看懂例3-1回答上述各题21二、二阶系统的数学模型及单位阶跃响应定义：由二阶微分方程描述的系统称为二阶系统。22北京航空航天大学二阶系统数学模型二阶系统的微分方程一般式为：23北京航空航天大学二阶系统的反馈结构图24北京航空航天大学二阶系统的传递函数开环传递函数：闭环传递函数：25北京航空航天大学二阶系统的特征方程为解方程求得特征根：当输入为阶跃信号时，则微分方程解的形式为：式中为由r(t)和初始条件确定的待定的系数。

s1,s2完全取决于，n两个参数。26北京航空航天大学此时s1,s2为一对共轭复根，且位于复平面的左半部。特征根分析——（欠阻尼）27北京航空航天大学特征根分析——（临界阻尼）此时s1,s2为一对相等的负实根。

s1=s2=-n28北京航空航天大学特征根分析——（过阻尼）此时s1,s2为两个负实根，且位于复平面的负实轴上。29北京航空航天大学特征根分析——（零阻尼）此时s1,s2为一对纯虚根，位于虚轴上。S1,2=jn30北京航空航天大学特征根分析—（负阻尼）此时s1,s2为一对实部为正的共轭复根，位于复平面的右半部。31北京航空航天大学特征根分析——（负阻尼）此时s1,s2为两个正实根，且位于复平面的正实轴上。32北京航空航天大学二阶系统单位阶跃响应1.过阻尼二阶系统的单位阶跃响应取C(s)拉氏反变换得：33北京航空航天大学过阻尼系统分析衰减项的幂指数的绝对值一个大，一个小。绝对值大的离虚轴远，衰减速度快，绝对值小的离虚轴近，衰减速度慢；衰减项前的系数一个大，一个小；二阶过阻尼系统的动态响应呈非周期性，没有振荡和超调，但又不同于一阶系统；离虚轴近的极点所决定的分量对响应产生的影响大，离虚轴远的极点所决定的分量对响应产生的影响小，有时甚至可以忽略不计。34北京航空航天大学过阻尼系统单位阶跃响应35北京航空航天大学与一阶系统阶跃响应的比较tc(t)0二阶过阻尼系统一阶系统响应136北京航空航天大学二阶过阻尼系统阶跃响应指标分析对于过阻尼二阶系统的响应指标，只着重讨论，它反映了系统响应过渡过程的长短，是系统响应快速性的一个方面，但确定的表达式是很困难的，一般根据（3－1－4）取相对量及经计算机计算后制成曲线或表格。37北京航空航天大学2.欠阻尼二阶系统的单位阶跃响应38北京航空航天大学二阶欠阻尼系统的输出拉氏反变换得：39北京航空航天大学二阶欠阻尼系统输出分析二阶欠阻尼系统的单位阶跃响应由稳态分量和暂态分量组成。稳态分量值等于1，暂态分量为衰减过程，振荡频率为ωd。40北京航空航天大学下图为二阶系统单位阶跃响应的通用曲线。41北京航空航天大学下图根据上图来分析系统的结构参数、对阶跃响应的影响平稳性（％）结论：越大，ωd越小，幅值也越小，响应的振荡倾向越弱，超调越小，平稳性越好。反之，越小，ωd越大，振荡越严重，平稳性越差。42北京航空航天大学当＝0时，为零阻尼响应，具有频率为的不衰减（等幅）振荡。阻尼比和超调量的关系曲线如下图所示43北京航空航天大学结论：对于二阶欠阻尼系统而言，大，

小，系统响应的平稳性好。在一定的情况下，越大，振荡频率

也越高，响应平稳性也越差。44北京航空航天大学快速性从图中看出，对于5％误差带，当时，调节时间最短，即快速性最好。同时，其超调量<5％，平稳性也较好，故称为最佳阻尼比。总结：越大，调节时间越短；当一定时，越大，快速性越好。45北京航空航天大学稳态精度从上式可看出，瞬态分量随时间t的增长衰减到零，而稳态分量等于1，因此，上述欠阻尼二阶系统的单位阶跃响应稳态误差为零。46北京航空航天大学欠阻尼二阶系统

单位阶跃响应性能指标1.上升时间：令，则所以：47北京航空航天大学根据极值定理有：该项不可能为零2.峰值时间48北京航空航天大学取n=1得:49北京航空航天大学3.超调量将峰值时间代入下式得:所以:50北京航空航天大学4.调节时间写出调节时间的表达式相当困难。在分析设设计系统十，经常采用下列近似公式。当阻尼比时51北京航空航天大学三、二阶系统举例2设位置随动系统，其结构图如图所示，当给定输入为单位阶跃时，试计算放大器增益KA＝200，1500，13.5时，输出位置响应特性的性能指标：峰值时间tp，调节时间ts和超调量，并分析比较之。52北京航空航天大学例题解析(1)输入：单位阶跃系统的闭环传递函数：53北京航空航天大学例题解析(2)当KA＝200时系统的闭环传递函数：与标准的二阶系统传递函数对照得：54北京航空航天大学例题解析(3)当KA＝1500时系统的闭环传递函数：与标准的二阶系统传递函数对照得：55北京航空航天大学例题解析(4)当KA＝13.5时系统的闭环传递函数：与标准的二阶系统传递函数对照得：无56北京航空航天大学系统在单位阶跃作用下的响应曲线57北京航空航天大学四改善二阶系统响应的措施1.误差信号的比例－微分控制58北京航空航天大学系统开环传函为：闭环传函为：等效阻尼比：59北京航空航天大学可见，引入了比例－微分控制，使系统的等效阻尼比加大了，从而抑制了振荡，使超调减弱，可以改善系统的平稳性。微分作用之所以能改善动态性能，因为它产生一种早期控制（和称为超前控制），能在实际超调量出来之前，就产生一个修正作用。60北京航空航天大学前面图的相应的等效结构由此知道：61北京航空航天大学和及的大致形状如下一方面，增加项，增大了等效阻尼比，使曲线比较平稳。另一方面，它又使加上了它的微分信号，加速了c(t)的响应速度，但同时削弱了等效阻尼比的平稳作用。62北京航空航天大学总结：引入误差信号的比例－微分控制，能否真正改善二阶系统的响应特性，还需要适当选择微分时间常数。若大一些，使具有过阻尼的形式，而闭环零点的微分作用，将在保证响应特性平稳的情况下，显著地提高系统的快速性。63北京航空航天大学2.输出量的速度反馈控制将输出量的速度信号c(t)采用负反馈形式，反馈导输入端并与误差信号e(t)比较，构成一个内回路，称为速度反馈控制。如下图示。64北京航空航天大学闭环传函为：等效阻尼比：等效阻尼比增大了，振荡倾向和超调量减小，改善了系统的平稳性。65北京航空航天大学3.比例－微分控制和速度反馈控制比较从实现角度看，比例－微分控制的线路结构比较简单，成本低；而速度反馈控制部件则较昂贵。从抗干扰来看，前者抗干扰能力较后者差。从控制性能看，两者均能改善系统的平稳性，在相同的阻尼比和自然频率下，采用速度反馈不足之处是其会使系统的开环增益下降，但又能使内回路中被包围部件的非线性特性、参数漂移等不利影响大大削弱。66北京航空航天大学五高阶系统的时域分析定义：用高阶微分方程描述的系统称高阶系统。由于求高阶系统的时间响应很是困难，所以通常总是将多数高阶系统化为一、二阶系统加以分析。通常对于高阶系统来说，离虚轴最近的一个或两个闭环极点在时间响应中起主导作用，而其他离虚轴较远的极点，他们在时间响应中相应的分量衰减较快，只起次要作用，可以忽略。67北京航空航天大学这时，高阶系统的时域分析就转化为相应的一、二系统。这就是所谓的主导极点的概念，将在第四章中详细介绍。一、二阶系统的极点分布如下：68北京航空航天大学3－3系统稳定性分析本节主要内容：线性定常系统稳定的概念系统稳定的条件和稳定性的判定方法。69北京航空航天大学一、系统稳定的概念是指系统当扰动作用消失后，由初始偏差状态恢复到原平衡状态的性能。若系统能恢复平衡状态，就称该系统是稳定的，若系统在扰动作用消失后不能恢复平衡状态，且偏差越来越大，则称系统是不稳定的。70北京航空航天大学二、稳定性的数学条件设系统的线形化增量方程为：71北京航空航天大学对上式行拉氏变换得：其中：D(s)为系统闭环特征式，也称输出端算子式；M(s)称为输入端算子式。R(s)为输入，C(s)为输出，M0(s)为总的初始条件，与系统的初始状态有关的多项式。或简写为：72北京航空航天大学则有：假定：将C(s)等式右的两项分别展开成部分分式，可得73北京航空航天大学再进行拉氏反变换，得该部分为稳态分量，即微分方程的特解，取决于输入作用。74北京航空航天大学该为瞬态分量，即微分方程的通解，运动规律取决于，由系统的结构参数确定。75北京航空航天大学系统去掉扰动后的恢复能力，应由瞬态分量决定。此时，系统的输入为零。故：稳定性定义可转化为：式中，Ai，Ci均为常值，因此，系统的稳定性仅取决于特征根si的性质。76北京航空航天大学特征根的性质对系统稳定性的影响当si为实根时，即si＝i：77北京航空航天大学78北京航空航天大学特征根与系统稳定性的关系(2)当si为共轭复根时，即si，i+1＝i±

jωi79北京航空航天大学共轭复根情况下系统的稳定性80北京航空航天大学结论：系统稳定的充分必要条件是：系统的特征方程的所有根都具有负实部，或者说都位于S平面的虚轴之左。注：拉氏变换性质中的终值定理的适用条件：S

E(S)在S平面的右半平面解析，就是上面稳定条件的另一种表示，即特征方程的所有根Si位于S平面的虚轴之左。81北京航空航天大学三、稳定性判据判据之一：赫尔维茨（Hurwitz)稳定判据系统稳定的充分必要条件是：特征方程的赫尔维茨行列式Dk（k＝1,2,3,….n）全部为正。82北京航空航天大学赫尔维茨判据系统特征方程的一般形式为：各阶赫尔维茨行列式为：（一般规定）83北京航空航天大学举例：系统的特征方程为：试用赫尔维茨判据判断系统的稳定性。84北京航空航天大学解：第一步：由特征方程得到各项系数２１３５第二步：计算各阶古尔维茨行列式结论：系统不稳定。１０85北京航空航天大学三、稳定性判据判据之二：林纳德－奇帕特（Lienard-Chipard)判据系统稳定的充分必要条件为：１．系统特征方程的各项系数大于零，即２．奇数阶或偶数阶的赫尔维茨行列式大于零。即或必要条件86北京航空航天大学举例：单位负反馈系统的开环传递函数为：试求开环增益Ｋ的稳定域。87北京航空航天大学解：第一步：求系统的闭环特征方程第二步：列出特征方程的各项系数第三步：系统稳定的充分必要条件是：88北京航空航天大学解得：Ｋ＜１４开环增益Ｋ的稳定域为：由此例可见，K越大，系统的稳定性越差。上述判据不仅可以判断系统的稳定性，而且还可根据稳定性的要求确定系统参数的允许范围（即稳定域）。89北京航空航天大学三、稳定性判据判据之三：劳思(Routh)判据系统稳定的充分必要条件是：劳思表中第一列所有元素的计算值均大于零。90北京航空航天大学若系统的特征方程为：则劳思表中各项系数如下图：91北京航空航天大学关于劳思判据的几点说明如果第一列中出现一个小于零的值，系统就不稳定；如果第一列中有等于零的值，说明系统处于临界稳定状态；第一列中数据符号改变的次数等于系统特征方程正实部根的数目，即系统中不稳定根的个数。92北京航空航天大学例1设系统特征方程如下：试用劳斯判据判断该系统的稳定性，并确定正实部根的数目。93北京航空航天大学解：将特征方程系数列成劳斯表结论：系统不稳定；系统特征方程有两个正实部的根。94北京航空航天大学劳思表判据的特殊情况在劳思表的中间某一行中，第一列项为零。在劳思表的中间某一行中，所有各个元素均为零。在这两种情况下，都要进行一些数学处理，原则是不影响劳思判据的结果。95北京航空航天大学例2设系统的特征方程为：试用劳思判据确定正实部根的个数。96北京航空航天大学解：将特征方程系数列成劳斯表由表可见，第二行中的第一列项为零，所以第三行的第一列项出现无穷大。为避免这种情况，可用因子(s+a)乘以原特征式，其中a可为任意正数，这里取a=1。97北京航空航天大学于是得到新的特征方程为：将特征方程系数列成劳斯表：结论：第一列有两次符号变化，故方程有两个正实部根。98北京航空航天大学例3设系统的特征方程为：试用劳思判据确定正实部根的个数。99北京航空航天大学解：将特征方程系数列成劳斯表劳思表中出现全零行，表明特征方程中存在一些大小相等，但位置相反的根。这时，可用全零行上一行的系数构造一个辅助方程，对其求导，用所得方程的系数代替全零行，继续下去直到得到全部劳思表。100北京航空航天大学用行的系数构造系列辅助方程求导得：用上述方程的系数代替原表中全零行，然后按正常规则计算下去，得到101北京航空航天大学102北京航空航天大学表中的第一列各系数中，只有符号的变化，所以该特征方程只有一个正实部根。求解辅助方程，可知产生全零行的根为。再可求出特征方程的其它两个根为。103北京航空航天大学四.结构不稳定及改进措施某些系统，仅仅靠调整参数仍无法稳定，称结构不稳定系统。如下图液位可能控制系统。104北京航空航天大学消除结构不稳定的措施有两种改变积分性质引入比例－微分控制，补上特征方程中的缺项。该系统的闭环特征方程为：系数缺项，显然不满足系统稳定的必要条件，且无论怎么调整系统参数，都不能使系统稳定。105北京航空航天大学1.改变积分性质用反馈包围积分环节，或者包围电动机的传递函数，破坏其积分性质。106北京航空航天大学2.引入比例－微分控制在原系统的前向通路中引入比例－微分控制。107北京航空航天大学其闭环特征方程为：由稳定的充分必要条件：引入比例－微分控制后，补上了特征方程中s的一次项系数。只要适当匹配参数，满足上述条件，系统就可以稳定。108北京航空航天大学3－4稳态误差分析计算一.误差与稳态误差系统的误差e(t)常定义为：e(t)=期望值－实际值误差:(1).e(t)=r(t)-c(t)(2).e(t)=r(t)-b(t)109北京航空航天大学稳态误差定义：稳定系统误差的终值称为稳态系统。当时间t趋于无穷时，e(t)极限存在，则稳态误差为二.稳态误差的计算若e(t)的拉普拉斯变换为E(s)，且110北京航空航天大学在计算系统误差的终值(稳态误差)时，遇到的误差的象函数一般是s的有理分式函数，这时当且仅当的极点均在左半面，就可保证存在，式就成立。注：sE(s)的极点均在左半面的条件中蕴涵了闭环系统稳定的条件。111北京航空航天大学对上述系统，若定义e(t)=r(t)-b(t),则E(s)=R(s)-B(s)112北京航空航天大学称之为系统对输入信号的误差传递函数。称为系统对干扰的误差传递函数。113北京航空航天大学例：系统结构如下图。当输入信号r(t)=1(t),干扰n(t)=1(t)时，求系统的总的稳态误差解：①判别稳定性。由于是一阶系统，所以只要参数大于零，系统就稳定。②求E(s)。114北京航空航天大学根据结构图可以求出：依题意：R(s)=N(s)=1/s,则③应用终值定理得稳态误差115北京航空航天大学三.输入信号作用下的稳态误差与系统结构参数的关系当系统只有输入r(t)作用时，系统的开环传递函数为：RECB116北京航空航天大学将G(s)H(s)写成典型环节串联形式：当sE(s)的极点全部在s平面的座半平面是，可用终值定理求得：上式表明：系统的稳态误差除与输入有关外，只与系统的开环增益K和积分环节的个数有关。117北京航空航天大学1.阶跃输入下的稳态误差要消除阶跃作用下的稳态误差，开环传递函数中至少要有一个积分环节。但是，积分环节多会导致系统不稳定。118北京航空航天大学2.斜坡输入下的稳态误差要消除斜坡作用下的稳态误差，开环传递函数中至少要有两个积分环节。119北京航空航天大学3.等加速输入下的稳态误差要消除等加速作用下的稳态误差，开环传递函数中至少要有三个积分环节。但是，积分环节多会导致系统不稳定。120北京航空航天大学由以上分析可见，要消除系统在幂函数输入信号作用下的稳态误差，则要求增加积分环节的数目，要减小系统的稳态误差，则要求提高开环增益。系统型别是针对系统的开环传递函数中积分环节的个数而言的。=０的系统称为０型系统；＝１的系统称为Ⅰ型系统；＝２的系统称为Ⅱ型系统；121北京航空