第3章控制系统的时域分析ppt课件.ppt

《第3章控制系统的时域分析ppt课件.ppt》由会员分享，可在线阅读，更多相关《第3章控制系统的时域分析ppt课件.ppt（108页珍藏版）》请在三一办公上搜索。

1、第3章 控制系统的时域分析,内 容 提 要,控制系统在典型输入信号作用下的动态过程的品质及稳态性能直接表征了系统的优劣。系统的稳定性是系统正常工作的首要条件，系统的稳定性完全由系统自身的结构和参数决定，而与系统的输入无关；系统的稳态误差是系统的稳态性能指标，它标志着系统的控制精度；系统的时域响应可定性或定量分析系统的动态性能。介绍了如何用MATLAB和SIMULINK进行瞬态响应分析。,知 识 要 点,系统稳定的充分必要条件，Routh判据，误差与稳态误差的定义，静态误差系数及系统的型号，线性定常一阶、二阶系统的时域响应及动态性能的计算，高阶系统的主导极点，偶极子及高阶系统的降阶。,目 录,3

2、.1 线性定常系统的时域响应 3.2 控制系统时域响应的性能指标 3.3 线性定常系统的稳定性 3.4 系统的稳态误差 3.5 一阶系统的时域响应 3.6 二阶系统的时域响应 3.7 高阶系统的瞬态响应 3.8 用MATLAB和SIMULINK进行瞬态响应分析,3.1 线性定常系统的时域响应,对于一单输入单输出n阶线性定常系统，可用一n阶常系数线性微分方程来描述。,系统在输入信号r(t)作用下，输出c(t)随时间变化的规律，即式(3-1)微分方程的解，就是系统的时域响应。由线性微分方程理论知，方程式的解由两部分组成，即 c(t)=c1(t)+c2(t)(3-2)c1(t)对应齐次微分方程的通解

3、 c2(t)非齐次微分方程的一个特解,从系统时域响应的两部分看，稳态分量(特解)是系统在时间t时系统的输出，衡量其好坏是稳态性能指标：稳态误差。系统响应的暂态分量是指从t=0开始到进入稳态之前的这一段过程，采用动态性能指标(瞬态响应指标)，如稳定性、快速性、平稳性等来衡量。,返回,3.2 控制系统时域响应的性能指标,3.2.1 稳态性能指标 采用稳态误差ess来衡量，其定义为：当时间t趋于无穷时，系统输出响应的期望值与实际值之差。即,1上升时间tr：从零时刻首次到达稳态值的时间,即阶跃响应曲线从t=0开始第一次上升到稳态值所需要的时间。,3.2.2 动态性能指标,2.峰值时间tp：从零时刻到达

4、峰值的时间，即 阶跃响应曲线从t=0开始上升到第一个峰值所需要的时间.3.最大超调量Mp：阶跃响应曲线的最大峰值与稳态值的差与稳态值之比，即,4.调整时间ts：阶跃响应曲线进入允许的误差带(一般取稳态值附近5%或2%作为误差带)并不再超出该误差带的最小时间，称为调整时间(或过渡过程时间)。5.振荡次数：在调整时间ts内响应曲线振荡的次数。,返回,3.3 线性定常系统的稳定性,3.3.1 稳定性的概念 若控制系统在足够小的初始偏差的作用下，其过渡过程随时间的推移，逐渐衰减并趋于零，即具有恢复原平衡状态的能力，则称这个系统稳定。否则，称这个系统不稳定。,3.3.2 线性定常系统稳定的充分必要条件,

5、设n阶线性定常系统的微分方程为 对式(3-7)作拉氏变换，得,在式(3-8)中取R(s)=0，得到在初始状态影响下系统的时间响应(即零输入响应)为若pi为系统特征方程D(s)=0的根且当pi各不相同时，有若系统所有特征根pi的实部均为负值，即Repi0 则零输入响应(暂态响应)最终将衰减到零，即 这样的系统就是稳定的。,反之，若特征根中有一个或多个根具有正实部时，则暂态响应将随时间的推移而发散，即这样的系统就是不稳定的。综上所述，系统稳定的充分必要条件是系统特征根的实部均小于零，或系统的特征根均在根平面的左半平面。,3.3.3 劳斯判据,设n阶系统的特征方程为 D(s)=a0sn+a1sn-1

6、+an-1s+an=a0(s-p1)(s-p2)(s-pn)=0 将上式的系数排成下面的行和列，即为劳斯阵列(劳斯表),sn a0 a2 a4 a6 sn-1 a1 a3 a5 a7 sn-2 b1 b2 b3 b4 sn-3 c1 c2 c3 c4 s2 f1 f2s1 g1s0 h1,其中劳斯判据给出了控制系统稳定的充分条件是，劳斯表中第一列所有元素均大于零。,例3-1 已知三阶系统特征方程为劳斯阵列为 故得出三阶系统稳定的充要条件为各系数大于零，且a1a2a0a3,例 3-2已知系统特征方程 方程无缺项，且系数大于零。列劳斯表：劳斯表中第一列元素大于零，系统是稳定的，即所有特征根均s平面

7、的左半平面。,例 3-3 系统特征方程为各项系数均大于零。列劳斯表：劳斯表中第一列各元素符号不完全一致，系统不稳定。第一列元素符号改变两次，因此系统有两个右半平面的根。,例 3-4 系统特征方程 它有一个系数为负的，有劳斯判据的系统不稳定。但究竟有几个右根，仍需列劳斯表：劳斯表中第一列元素符号改变两次，系统有2个右半平面的根,有两种特殊情况需要说明：1.劳斯表中某一行的第一个元素为零，而该行其它元素并不为零，则在计算下一行第一个元素时，该元素必将趋于无穷大，以至劳斯表的计算无法进行。2.劳斯表中某一行的元素全为零。则表示在s平面内存在一些大小相等符号相反的实根或共轭虚根，系统是不稳定的。,例

8、3-6系统特征方程 列劳斯表 劳斯表中第一列元素符号没有改变，系统没有右半平面的根，但由P(s)=0求得 即系统有一对共轭虚根，系统处于临界稳定，从工程角度来看，临界稳定属于不稳定系统。,例3-7 系统的特征方程为列劳斯表：劳斯表中第一列元素符号改变一次，系统不稳定，且有一个右半平面的根，由P(s)=0得,3.3.4 赫尔维茨判据,设系统的特征方程式为以特征方程式的各项系数组成如下行列式,赫尔维茨判据指出，系统稳定的充分必要条件是在a00的情况下，上述行列式的各阶主子式i均大于零，即,例3-8 系统的特征方程为列出行列式 由赫尔维茨判据，该系统稳定的充分必要条件是：,或写成系统稳定的充分必要条

9、件为a00 a10 a20 a30a1a2-a0a30,例3-9 二阶系统的特征方程为列出行列式 由Hurwitz判据，系统稳定的充分必要条件为a00 a10 a1a20即二阶系统稳定的充分必要条件是特征方程式的所有系数均大于零。,3.3.5 系统参数对稳定性的影响,应用代数判据不仅可以判断系统的稳定性，还可以用来分析系统参数对系统稳定性的影响。例3-10 系统结构图如图3-2所示，试确定系统稳定时K的取值范围解 系统的闭环传递函数 其特征方程式为,列劳斯表 按劳斯判据，要使系统稳定，应有K0，且30-K0，故K的取值范围为0K30,例 3-11 系统结构图如图3-3所示，试分析参数K1，K2

10、，K3和T对系统稳定性的影响。解 系统的闭环传递函数特征方程为,由于特征方程缺项，由劳斯判据知，不论K1，K2，K3和T取何值系统总是不稳定的，称为结构不稳定系统。欲使系统稳定，必须改变系统的结构。如在原系统的前向通道中引入一比例微分环节，如图3-4所示。变结构后系统的闭环传递函数为特征方程为,列劳斯阵列：系统稳定的充分必要条件为 即对于结构不稳定系统，改变系统结构后，只要适当选配参数就可使系统稳定。,3.3.6 相对稳定性和稳定裕量,相对稳定性即系统的特征根在s平面的左半平面且与虚轴有一定的距离，称之为稳定裕量。为了能应用上述的代数判据，通常将s平面的虚轴左移一个距离，得新的复平面s1，即令

11、s1=s+或s=s1-得到以s1为变量的新特征方程式D(s1)=0，再利用代数判据判别新特征方程式的稳定性，若新特征方程式的所有根均在s1平面的左半平面，则说明原系统不但稳定，而且所有特征根均位于-的左侧，称为系统的稳定裕量,例 3-12 检验特征方程式 是否有根在s右半平面，以及有几个根在s=-1垂线的右边。解 列劳斯表：由劳斯判据知，系统稳定，所有特征根均在s的左半平面。,令s=s1-1代入D(s)得s1的特征方程式列劳斯表：劳斯表中第一列元素符号改变一次，表示系统有一个根在s1右半平面，也就是有一个根在s=-1垂线的右边(虚轴的左边)，系统的稳定裕量不到1。,返回,3.4 系统的稳态误差

12、,3.4.1 误差及稳态误差的定义 系统的误差e(t)一般定义为被控量的希望值与实际值之差。即 e(t)=被控量的希望值被控量的实际值,对于图3-5所示的反馈控制系统，常用的误差定义有两种 1.输入端定义2.输出端定义当图3-5中反馈为单位反馈时，即H(s)=1时，上述两种定义可统一为,误差响应e(t)与系统输出响应c(t)一样，也包含暂态分量和稳态分量两部分，对于一个稳定系统，暂态分量随着时间的推移逐渐消失，而我们主要关心的是控制系统平稳以后的误差，即系统误差响应的稳态分量稳态误差记为ess。定义稳态误差为稳定系统误差响应e(t)的终值。当时间t趋于无穷时，e(t)的极限存在，则稳态误差为,

13、3.4.2 稳态误差分析,根据误差和稳态误差的定义，系统误差e(t)的象函数 定义 为系统对输入信号的误差传递函数。,由拉普拉斯变换的终值定理计算稳态误差，则代入E(s)表达式得 从上式得出两点结论：1.稳态误差与系统输入信号r(t)的形式有关；2.稳态误差与系统的结构及参数有关。,3.4.3 稳态误差的计算,对于线性系统，响应具有叠加性，不同输入信号作用于系统产生的误差等于每一个输入信号单独作用时产生的误差的叠加。对于图3-7所示系统，控制信号r(t)和扰动信号n(t)同时作用于系统。,1.控制信号r(t)单独作用下，误差 稳态误差essr为,(3-22),2.扰动信号单独作用下，误差 稳态

14、误差,定义 为系统对扰动的误差传递函数。控制系统在给定信号r(t)和扰动信号n(t)同时作用下的稳态误差ess为,例3-13 系统结构图如图所示，当输入r(t)=4t时，求系统的稳态误差ess。解 系统只有在稳定的条件下计算稳态误差才有意义，所以应先判别系统的稳定性。系统的特征方程为列劳斯表,由劳斯判据知，系统稳定条件为系统的误差函数为,由终值定理求得稳态误差计算表明，稳定误差的大小与系统的放大倍数K有关。即K越大，稳定误差ess越小。要减小稳态误差则应增大倍数K，耳闻定性分析却得出，使系统稳定的K只应小于5/4，表明系统的稳态精度和稳态性对放大倍数的要求常是矛盾的。,3.4.4 应用静态误差

15、系数计算给定信号作用下的稳态误差,1.系统的类型 系统的开环传递函数G(s)H(s)可表示为 系统常按开环传递函数中所含有的积分环节个数来分类。把=0，1，2，的系统，分别称为0型，型，型，系统。,2.静态位置误差系数Kp 当系统的输入为单位阶跃信号r(t)=1(t)时，由(3-22)式，有其中，定义为系统静态位置误差系数。,对于0型系统对于型或高于型以上系统,3.静态速度误差系数Kv当系统的输入为单位斜坡信号时r(t)=t1(t)，即则 由式(3-22)有 其中，定义为系统静态速度误差系数。,对于0型系统 对于型系统,对于型或型以上系统,4.静态加速度误差系数Ka,当系统输入为单位加速度信号

16、时，即 则系统稳态误差为其中，定义为系统静态加速度误差系数。,对于0型系统，Ka=0，ess=；对于型系统，Ka=0，ess=；对于型系统，Ka=K，；对于型或型以上系统，Ka=，ess=0。,例3-14 系统结构如图3-9所示，求当输入信号r(t)=2t+t2时，系统的稳态误差ess。首先判别系统的稳定性。由开环传递函数知，闭环特征方程为根据劳斯判据知闭环系统稳定。,3.4.5 干扰信号作用下的稳态误差与系统结构的关系,扰动信号n(t)作用下的系统结构图如图所示 扰动信号n(t)作用下的误差函数为,稳态误差 若，则上式可近似为由上可得，干扰信号作用下产生的稳态误差essn除了与干扰信号的形式

17、有关外，还与干扰作用点之前(干扰点与误差点之间)的传递函数的结构及参数有关，但与干扰作用点之后的传递函数无关。,3.4.6 改善系统稳态精度的途径,从上面稳态误差分析可知，采用以下途径来改善系统的稳态精度：1.提高系统的型号或增大系统的开环增益，可以保证系统对给定信号的跟踪能力。但同时带来系统稳定性变差，甚至导致系统不稳定。2.增大误差信号与扰动作用点之间前向通道的开环增益或积分环节的个数，可以降低扰动信号引起的稳态误差。但同样也有稳定性问题。3.采用复合控制，即将反馈控制与扰动信号的前馈或与给定信号的顺馈相结合。,返回,3.5一阶系统的时域响应,3.5.1 数学模型 能够用一阶微分方程描述的

18、系统为一阶系统，其传递函数为其中T一阶系统的时间常数,3.5.2 单位阶跃响应,当r(t)=1(t)时，一阶系统的输出c(t)称为单位阶跃响应，记作h(t)。,图3-12 一阶系统的单位阶跃响应,3.5.3 性能指标,1.调整时间ts 经过时间3T4T，响应曲线已达稳态值的95%98%，可以认为其调整过程已完成，故一般取ts=(34)T。2.稳态误差ess 系统的实际输出h(t)在时间t趋于无穷大时，接近于输入值，即3.超调量Mp 一阶系统的单位阶跃响应为非周期响应，故系统无振荡、无超调，Mp=0。,3.5.4 一阶系统的单位脉冲响应,当输入信号r(t)=(t)时，系统的输出称为单位脉冲响应，

19、记为g(t)。当r(t)=(t)，即R(s)=1时，有,返回,3.6 二阶系统的时域响应,3.6.1 二阶系统的数学模型 典型二阶系统的结构图如图3-14所示，其闭环传递函数为或,图3-14 典型二阶系统结构图,其中 系统的阻尼比 n系统的无阻尼自然振荡角频率 系统振荡周期系统的特征方程为特征根为,3.6.2 二阶系统的单位阶跃响应,1.当1时，系统有两个不相等的负实根，称为过阻尼状态。两个不相等的负实根为单位阶跃响应,2.当01时，系统有一对实部为负的共轭复根，称为欠阻尼状态。在欠阻尼状态下，系统的两个闭环极点为一对共轭复极点，即其中，称为阻尼振荡角频率。,单位阶跃响应 也可写成 式中,图3

20、-16 欠阻尼状态下系统单位阶跃响应,上升时间tr其中(2)峰值时间tp,(3)最大超调量Mp(4)调整时间ts 误差带范围为 5%误差带范围为 2%(5)振荡次数N,3.当阻尼比=1时，系统的特征根为两相等的负实根，称为临界阻尼状态。,此时系统在单位阶跃函数作用下，系统的超调量Mp=0，调节时间(对应误差带为5%),图3-18 临界阻尼系统阶跃响应,4.当阻尼比=0时，系统特征根为一对纯虚根，称为无阻尼状态。系统特征根 单位阶跃响应为,3.6.3 二阶系统的单位脉冲响应,1.脉冲响应及脉冲响应函数 当系统输入信号为单位脉冲函数(t)时，系统的响应为单位脉冲响应，记为g(t)。,2.脉冲响应与

21、阶跃响应的关系 系统的单位阶跃响应是该系统单位脉冲响应的积分，或系统的单位脉冲响应是该系统单位阶跃响应的导数。即或,3.二阶系统的单位脉冲响应 当1时当01时当=1时 当=0时,例3-15 原控制系统如图3-23(a)所示，引入速度反馈后的控制系统如图3-23(b)所示，已知在图3-23(b)中，系统单位阶跃响应的超调量Mp%=16.4%，峰值时间tp=1.14s，试确定参数K和Kt，并计算系统在(a)和(b)的单位阶跃响应h(t)。,图3-23 例3-15图,解 对于系统(b)，其闭环传递函数为与典型二阶系统相比较，有（3-55）而已知Mp=16.4%tp=1.14s根据 求得,求得 将 代

22、入(3-55)得 其单位阶跃响应为,对于系统(a)，其闭环传递函数为与典型二阶系统比较有 系统的最大超调量 峰值时间其单位阶跃响应为,返回,3.7 高阶系统的瞬态响应,3.7.1 高阶系统的瞬态响应 n阶系统的闭环传递函数为,当输入为单位阶跃函数r(t)=1(t)，即时，则假设所有闭环零点和极点互不相等且均为实数,当极点中还包含共轭复极点时 进行拉普拉斯反变换可得系统的单位阶跃响应,3.7.2 高阶系统的降阶,主导极点 在整个响应过程中起着主要的决定性作用的闭环极点，我们称它为主导极点。工程上往往只用主导极点估算系统的动态特性。即将系统近似地看成是一阶或二阶系统。,2.偶极子 将一对靠得很近的

23、闭环零、极点称为偶极子。工程上，当某极点和某零点之间的距离比它们的模值小一个数量级，就可认为这对零极点为偶极子。闭环传递函数中，如果零、极点数值上相近，则可将该零点和极点一起消掉，称之为偶极子相消。,3.7.3 零极点对阶跃响应的影响,1.零点对阶跃响应的影响 假设系统中增加一个闭环实零点，即系统中增加了一个串联环节 且闭环零点z位于复平面的左半平面，,上式拉普拉斯反变换 可见，增加一个闭环左实零点以后，系统阶跃响应增加了一项，该项的值与c(t)的变化率成正比，与该零点离虚轴的距离成反比。显然，该零点的增加将使系统响应过程加快，超调量增大，系统对输入作用的反应灵敏了。,反之，如果增加的闭环零点

24、位于复平面的右半平面，即，则这将使系统响应过程变慢，超调量减小，系统对输入作用的反应变滞呆了。,2.极点对阶跃响应的影响,假设系统增加一个闭环左实极点-|p|，系统在单位阶跃信号作用下输出 取拉普拉斯反变换得,可以看出：系统中增加一个闭环左实极点，系统的过渡过程将变慢，超调量将减小，系统的反应变得较为滞呆。对于闭环传递函数存在右极点的情况，系统时域响应是发散的，系统不稳定,返回,3.8 用MATLAB和SIMULINK进行瞬态响应分析,3.8.1 单位脉冲响应 当输入信号为单位脉冲函数(t)时，系统输出为单位脉冲响应，MATLAB中求取脉冲响应的函数为impulse()，其调用格式为 y，x，

25、t=impulse(num，den，t)或 impulse(num，den)式中G(s)=num/den；t为仿真时间；y为时间t的输出响应；x为时间t的状态响应。,例3-16 试求下列系统的单位脉冲响应MATLAB命令为：t=0：0.1：40；num=1；den=1，0.3，1；impulse(num，den，t)；grid；title(Unit-impulse Response of G(s)=1/(s2+0.3s+1)其响应结果如图所示。,例3-17 系统传递函数为求取其单位脉冲响应的MATLAB命令为t=0：0.1：10；num=1；den=1，1，1；y，x，t=impulse(nu

26、m，den，t)plot(t，y)；grid xlabel(t)；ylable(y)；其响应结果如图所示。,3.8.2 单位阶跃响应,当输入为单位阶跃信号时，系统的输出为单位阶跃响应，在MATLAB中可用step()函数实现，其调用格式为y，x，t=step(num，den，t)或step(num，den),例3-18 求系统传递函数为num=1；den=1，0.5，1；t=0：0.1：10；y，x，t=step(num，den，t)；plot(t，y)；grid；xlabel(Time sec t)；ylabel(y)响应曲线如图3-26所示,图3-26 单位阶跃响应,3.8.3 斜坡响应,

27、在MATLABA中没有斜坡响应命令，因此，需要利用阶跃响应命令来求斜坡响应。根据单位斜坡响应输入是单位阶跃输入的积分。当求传递函数为的斜坡响应时，可先用除得，再利用阶跃响应命令即可求得斜坡响应。,例3-19 已知闭环系统传递函数 对单位斜坡输入 则,系统单位斜坡响应的MATLAB命令：num=1；den=1，0.3，1，0；t=0：0.1：10；c=step(num，den，t)；plot(t，c)；grid；xlabel(t sec)；ylabel(Input and Output)其响应结果如图所示。,3.8.4 任意函数作用下系统的响应,用线性仿真函数lsim来实现，其调用格式为 y，x

28、=lsim(num，den，u，t)式中；y(t)为系统输出响应；x(t)为系统状态响应；u为系统输入信号；t为 仿真时间。,例3-20 反馈系统如图3-28(a)所示，系统输入信号为图3-28(b)所示的三角波，求取系统输出响应。,MATLAB实现指令numg=10，20；deng=1，10，0；num，den=cloop(numg，deng，-1)；v1=0：0.1：2；v2=1.9：-0.1：-2；v3=-1.9：0.1：0；t=0：0.1：8；u=v1，v2，v3；y，x=lsim(num，den，u，t)；plot(t，y，t，u)；xlabel(Time sec)；,3.8.5 Simulink中的时域响应举例,例3-21图3-30的Simulink的仿真框图可演示系统对典型信号的时间响应曲线，图中给出阶跃响应曲线。,返回,