第四章线性控制系统的稳定性

1、2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系1 第四章第四章 线性控制系统的稳定性线性控制系统的稳定性 4.1 4.1 线性系统稳定性的基本概念线性系统稳定性的基本概念 4.2 4.2 传递函数表示的系统稳定性判定传递函数表示的系统稳定性判定 4.3 4.3 状态空间表示的系统稳定性判定状态空间表示的系统稳定性判定 4.4 4.4 本章小结本章小结2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系2稳定是控制系统能够正常运行的首要条件 对系统进行各类品质指标的分析必须在系统稳定的 前提下进行。自动控制理论的基本任务(之一) 分析系统的稳

2、定性问题 提出保证系统稳定的措施 一、稳定性分析的重要性一、稳定性分析的重要性 4.14.1线性系统稳定性的基本概念线性系统稳定性的基本概念 2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系3二、线性系统稳定性分析的理论框架二、线性系统稳定性分析的理论框架 第一第一方法方法第二第二方法方法稳定性分析稳定性分析18921892年俄国数学年俄国数学家李雅普诺夫家李雅普诺夫SISOSISO的代数的代数分析方法分析方法解析解析方法方法RouthRouth判据判据HouwitzHouwitz判据判据根据根据SISOSISO闭环特闭环特征方程的系数判征方程的系数判定定系统的系统

3、的稳定性稳定性根据状态方程根据状态方程A A阵阵判定系统的稳定性判定系统的稳定性2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系4A.Lyapunov(1857-1918)A.Lyapunov(1857-1918)，俄国数学家（俄国数学家（ChebyshevChebyshev 的学生，的学生，MarkovMarkov的同学），的同学），在他的博士论文中，在他的博士论文中，LyapunovLyapunov系统地研究了由系统地研究了由微分方程描述的一般运动微分方程描述的一般运动的稳定性问题，建立了著的稳定性问题，建立了著名的名的LaypunovLaypunov方法，他的

4、方法，他的工作为现代控制及非线性工作为现代控制及非线性控制奠定科基础。控制奠定科基础。三、线性系统稳定性分析的划时代人物三、线性系统稳定性分析的划时代人物 2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系54.2 4.2 传递函数表示的系统稳定性判定传递函数表示的系统稳定性判定本小节是本章的重点，主要介绍以下内容：本小节是本章的重点，主要介绍以下内容：4.2.1 SISO4.2.1 SISO线性定常系统的稳定性问题线性定常系统的稳定性问题4.2.2 Routh4.2.2 Routh稳定判据稳定判据4.2.3 Routh4.2.3 Routh判据的两种特殊情况判据的两

5、种特殊情况4.2.4 Routh4.2.4 Routh判据的推广判据的推广4.2.5 Routh4.2.5 Routh判据的应用判据的应用2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系62022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系72022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系8一、稳定性基本概念一、稳定性基本概念 1 1、稳定性稳定性 任何系统在扰动的作用下都会偏离原平衡状态，产生初任何系统在扰动的作用下都会偏离原平衡状态，产生初始偏差。始偏差。 所谓稳定性，是指系统在扰动消失后，由初始偏差状态恢所谓

6、稳定性，是指系统在扰动消失后，由初始偏差状态恢复到原平衡状态的性能。复到原平衡状态的性能。 4.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性2 2、平衡状态平衡状态系统所受的作用力达到平衡，系统所受的作用力达到平衡，使系统处于稳定（不运动）的状态。使系统处于稳定（不运动）的状态。称为平衡状态。称为平衡状态。ab2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系94.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系104.2.1 SIS

7、O4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系114.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系124.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性临界稳定临界稳定2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系134.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性李雅普诺夫（渐进）稳定性定义：若线性系统在

8、初始扰动的影响下，其动态过程随时间的推移逐渐衰减并趋于零或原平衡工作点，则称系统渐进稳定，简称稳定。反之，若初始扰动的影响下，系统的动态过程随时间的推移而发散，则称系统不稳定。在古典控制理论中的稳定均指渐进稳定！2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系144.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性由稳定性定义可知：1）线性系统的稳定性取决于系统自身的固有特征（结构、参数），与系统的输入信号无关。2）若处于平衡状态的线性定常系统在脉冲信号的作用下，系统的响应最终能够回到平衡状态，则该线性定常系统稳定。2022-3-9北京科

9、技大学自动化学院自动化系北京科技大学自动化学院自动化系154.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性11( )( )( ) ( )( )( ) ( ) ( )Y sC sG s U sG sc tL C sL G s 所所以以对于脉冲响应，我们有：对于脉冲响应，我们有：() 1Lt 显然，系统是否稳定取决于显然，系统是否稳定取决于G(sG(s) )极点在极点在S S平面中的位置。平面中的位置。推论推论1 1：如果当时间趋于无穷时，线性定常系统的如果当时间趋于无穷时，线性定常系统的脉冲脉冲响应函数趋于零响应函数趋于零，则该线性定常系统稳定。，则该线性定常系统

10、稳定。系统是稳定的。系统是稳定的。lim ( )0tc t 系统仍能回到原有的平衡状态系统仍能回到原有的平衡状态简证：简证：2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系164.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性若系统闭环传递函数的所有极点全部 位于S左半平面，则系统稳定。则脉冲响应为：则脉冲响应为： 令系统的闭环传递函数含有令系统的闭环传递函数含有q q个实数极点和个实数极点和r r对复数对复数 极点：极点：)2()()()()(22111+P+P+PnknkkrkjqjimissPsZsKssGwwxf+rkrkknk

11、tkknktkqjtpjteCteBeAtgnkknkkj112211cos1sin)(xwxwwxwx显然只有当系统闭环传递函数的所有极点全部位于显然只有当系统闭环传递函数的所有极点全部位于S S左左半平面时，半平面时，g(t)|g(t)|t t 0 0 成立，即系统才稳定。成立，即系统才稳定。 2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系17j0j0j0j0j04.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系184.2.1 SISO4.2.1 SISO

12、线性定常系统的稳定性线性定常系统的稳定性如果当时间趋于无穷时，线性定常系如果当时间趋于无穷时，线性定常系统的统的阶跃响应函数趋于某一个常数阶跃响应函数趋于某一个常数，则该线性，则该线性定常系统稳定。定常系统稳定。 这个推论的证明请同学们自行完成。这个推论的证明请同学们自行完成。 2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系19二、二、SISOSISO系统阶跃响应的稳定问题系统阶跃响应的稳定问题实根情况：实根情况：2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系20虚根情况：虚根情况：2022-3-9北京科技大学自动化学院自动化系北京科

13、技大学自动化学院自动化系214.2.1 SISO4.2.1 SISO线性定常系统的稳定性线性定常系统的稳定性临界稳定：当系统的极点有在虚轴上时，系统的输出将会出现等幅振荡的状态，称之为临界稳定状态。稳定裕度的概念：S S平面平面2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系22三、三、SISOSISO线性定常系统的稳定性分析方法：线性定常系统的稳定性分析方法： 求脉冲响应求脉冲响应 求阶跃响应求阶跃响应 求系统的闭环特征根求系统的闭环特征根不易求不易求其它简单的判定方法其它简单的判定方法? ?2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动

14、化系234.2.2 Routh4.2.2 Routh稳定判据稳定判据(Rouths(Rouths stability criterion) stability criterion) Routh Routh表表将闭环特征方程将闭环特征方程的各项系数，按的各项系数，按右面的格式排成右面的格式排成RouthRouth表。表。102113212321343212753116420fSeeSdddScccSabbbSaaaaSaaaaSnnnn 000122110 + + + + + + aaSaSaSaSannnnn系统闭环特征方程系统闭环特征方程130211aaaaab150412aaaaab132

15、111bab acb153121bab acb2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系24系统渐进稳定的必要条件是系统渐进稳定的必要条件是特征方程的系数特征方程的系数均大于零均大于零。如果劳斯表中如果劳斯表中第一列的系数第一列的系数均为均为正值正值，则其特征方程式，则其特征方程式的根都在的根都在S S的左半平面，相应的系统是稳定的。的左半平面，相应的系统是稳定的。如果劳斯表中如果劳斯表中第一列系数的符号有变化第一列系数的符号有变化，则符号的变化，则符号的变化次数等于该特征方程式的根在次数等于该特征方程式的根在S S的右半平面上的个数，相的右半平面上的个数，相

16、应的系统为不稳定。应的系统为不稳定。劳劳斯斯稳稳定定判判据据表中表中这样可求得这样可求得n+1n+1行系数行系数 121211141713131512121311170613150412130211,eeddefbbaabcbbaabcbbaabcaaaaabaaaaabaaaaab 2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系250103 . 25175 .41423+SSS例例4.2-14.2-1试用劳斯判据判别系统的稳定性试用劳斯判据判别系统的稳定性。列劳斯表列劳斯表401423103 . 25 .380103 . 25 .4105171SSSS由于该表第

17、一列系数的符号变化了两次，所以该方程中有由于该表第一列系数的符号变化了两次，所以该方程中有二个根在二个根在S S的右半平面，因而系统是不稳定的。的右半平面，因而系统是不稳定的。已知某一调速系统的闭环特征方程式为：已知某一调速系统的闭环特征方程式为：2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系26 4.2.3 Routh4.2.3 Routh判据的两种特殊情况判据的两种特殊情况劳斯表某一行中的第一项元素等于劳斯表某一行中的第一项元素等于0 0，而该行的其余各项，而该行的其余各项不等于不等于0 0或没有其余项。或没有其余项。 以一个很小的正数以一个很小的正数 来代替

18、为来代替为0 0的这项，据此算出其的这项，据此算出其余的各项，完成劳斯表的排列。余的各项，完成劳斯表的排列。解决的办法解决的办法若劳斯表第一列中系数的符号有变化，其变化的次数若劳斯表第一列中系数的符号有变化，其变化的次数就等于该方程在就等于该方程在S S右半平面上根的数目，相应的右半平面上根的数目，相应的 如果第一列如果第一列 上面的系数与下面的系数符号相同，则上面的系数与下面的系数符号相同，则表示该方程中有一对共轭虚根存在，相应的系统为表示该方程中有一对共轭虚根存在，相应的系统为结结论论2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系27已知系统的闭环特征方程式为

19、已知系统的闭环特征方程式为02223+SSS试判别相应系统的稳定性。试判别相应系统的稳定性。2)(022110123SSSS由于表中第一列由于表中第一列 上面元素的符号与其下面元素的符号相同，上面元素的符号与其下面元素的符号相同，所以该闭环特征方程中有一对共轭虚根存在，相应的系统为所以该闭环特征方程中有一对共轭虚根存在，相应的系统为临界稳定系统临界稳定系统( (这在工业上属于不稳定的系统这在工业上属于不稳定的系统) )。例4.2-2列劳斯表列劳斯表2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系28劳斯表某一行元素全为劳斯表某一行元素全为0 0。这表示相应方程中含有

20、一些大。这表示相应方程中含有一些大小相等符号相反的实根或共轭虚根小相等符号相反的实根或共轭虚根( (关于原点对称的根关于原点对称的根) )。 利用系数全为利用系数全为0 0行的上一行系数构造一个辅助多项行的上一行系数构造一个辅助多项式，并以这个辅助多项式导数的系数来代替表中系数为式，并以这个辅助多项式导数的系数来代替表中系数为全全0 0的行。从而完成劳斯表的排列。的行。从而完成劳斯表的排列。解决解决办法办法关于原点对称的根可以通过求解这个辅助方程式得到，而关于原点对称的根可以通过求解这个辅助方程式得到，而且其根的数目总是偶数的。且其根的数目总是偶数的。若劳斯表第一列中系数的符号有变化，其变化的

21、次数就等若劳斯表第一列中系数的符号有变化，其变化的次数就等于该方程在于该方程在S S右半平面上根的数目，相应的右半平面上根的数目，相应的如果第一列如果第一列上的元素没有符号变化，则表示该方程中有共上的元素没有符号变化，则表示该方程中有共轭纯虚根存在，相应的系统为轭纯虚根存在，相应的系统为结结论论4.2.3 Routh4.2.3 Routh判据的两种特殊情况判据的两种特殊情况2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系2916038166248000161220161221620810123456SSSSSSS由于第一列的系数均由于第一列的系数均为正值，表明该方程

22、为正值，表明该方程在在S S右半平面上没有右半平面上没有特征根。特征根。该系统处于临界该系统处于临界稳定状态稳定状态。 已知系统的闭环特征方程式为已知系统的闭环特征方程式为试判别相应系统的稳定性。试判别相应系统的稳定性。例4.2-3列劳斯表列劳斯表ssdssdF248)(3+2,2jj令令F(sF(s)=0)=0，求得：，求得：sssF16122)(24+=2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系304.2.4 Routh4.2.4 Routh判据的推广判据的推广实际系统希望实际系统希望S S左半平面上的根距离左半平面上的根距离虚轴有一定的距离。这种系统在系

23、统虚轴有一定的距离。这种系统在系统参数发生一定变化时仍能保持稳定。参数发生一定变化时仍能保持稳定。 此法可以估计一个稳定系统的所有闭环特征根中此法可以估计一个稳定系统的所有闭环特征根中最靠近虚轴的根离虚轴有多远，从而了解系统稳定最靠近虚轴的根离虚轴有多远，从而了解系统稳定的的“程度程度”稳定裕度。稳定裕度。 令令s=ss=s1 1-a-a，代入原系统的闭环特征方程中，得到，代入原系统的闭环特征方程中，得到以以s s1 1 为变量的特征方程式，然后用劳斯判据去判别该为变量的特征方程式，然后用劳斯判据去判别该方程中是否有根位于垂线方程中是否有根位于垂线s s1 1=-a=-a右侧。右侧。1sa02

24、022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系31用劳斯判据检验下列特征方程用劳斯判据检验下列特征方程例例4.2-44.2-4列劳斯表列劳斯表041310223+SSS是否有根在是否有根在S S的右半平面上，并检验有几个根在的右半平面上，并检验有几个根在 的右方。的右方。1S42 .121081304101320123SSSS第一列全为正，所有的根均位于左半平面，系统稳定。第一列全为正，所有的根均位于左半平面，系统稳定。s s1 1s s-1 -10 0j j2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系32令令S=Z-1S=Z-1代入

25、特征方程：代入特征方程：04) 1( 3) 1(10) 1(223+ZZZ014223+ZZZ式中有负号，显然有根在式中有负号，显然有根在1S的右方的右方。列劳斯表列劳斯表12114120123SSSS第一列的系数符号变第一列的系数符号变化了一次，表示原方化了一次，表示原方程有一个根在垂直直程有一个根在垂直直线线 的右方。的右方。1S041310223+SSS2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系334.2.5 Routh4.2.5 Routh判据的应用判据的应用例例4.2-54.2-51 1 系统参数稳定范围的确定系统参数稳定范围的确定已知某调速系统的特

26、征方程式为已知某调速系统的特征方程式为 0)1 (16705175 .4123+KSSS求该系统稳定的求该系统稳定的K K值范围。值范围。)1 (167005 .41)1 (16705175 .410)1 (16705 .41051710123KSKSKSS+由劳斯判据可知，若系统稳定，则劳斯表中第一列的系数由劳斯判据可知，若系统稳定，则劳斯表中第一列的系数必须全为正值：必须全为正值：+0)1 (16700)1 ( 2 .40517KK9 .111K列劳斯表列劳斯表2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系34第一列均为正值，第一列均为正值，S S全部位全部位

27、于左半平面，故系统稳定。于左半平面，故系统稳定。已知一单位反馈控制系统如下图所示，试回答：已知一单位反馈控制系统如下图所示，试回答： )(sR)(sCsKt）（10)5(20+sss)(sGc时，闭环系统是否稳定？时，闭环系统是否稳定？ ssKsGpc) 1()(+时，闭环系统的稳定条件是什么？时，闭环系统的稳定条件是什么？1)(sGc例4.2-61)(sGc20152075020155010123SSSS时，闭环系统时，闭环系统的的特征方程为特征方程为: :020501523=+SSS2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系35闭环特征方程为：闭环特征方程

28、为： 020205015234+ppKSKSSSssKsGpc) 1()(+开环传递函数：开环传递函数： )10)(5() 1(20)()(2+SSSSKsGsGpc列列劳劳斯斯表表432101502015200750202015750202015 2015(75020)/1520pppppPppsKpsKKsKKKKsKsK2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系36 因此，利用劳斯稳定判据可确定系统一因此，利用劳斯稳定判据可确定系统一个或两个可调参数对系统稳定性的影响。个或两个可调参数对系统稳定性的影响。 欲使系统稳定第一列的系数必须全为正值欲使系统稳定

29、第一列的系数必须全为正值 0pK5 .37020750ppKK020525015152075001520750)151520750(20pppppKKKKK5 .26pK5 .260pK2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系37+)sT)(sT( sk1121+系统闭环特征方程为：系统闭环特征方程为：0221321+kss )TT(sTT02121+kTkTTT01121+kTT为稳定条件为稳定条件例4.2-7系统结构图如下所示，确定系统参数稳定的条件。系统结构图如下所示，确定系统参数稳定的条件。 S3 T1T2 1S2 T1+T2 kS1 0S0 k 0

30、121212(T +T )-kTTT +T2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系38当当K=2K=2时，时，RouthRouth表的第三、表的第三、五行元素全为五行元素全为0 0。系统将有对称于系统将有对称于原点的闭环特征原点的闭环特征根。根。2 2 求特殊情况下系统的闭环特征根求特殊情况下系统的闭环特征根例例4.2-84.2-8已知某系统的闭环特征方程为：已知某系统的闭环特征方程为： 5432SSKS2S +S+1=0+试确定使系统有对称于原点的闭环特征根的试确定使系统有对称于原点的闭环特征根的K K值，并求出此时值，并求出此时系统的所有闭环特征根系统的

31、所有闭环特征根。543210S1K1S121SK-20S31K-2S3S142S2S10+ ，进而得，进而得列劳斯表列劳斯表2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系394.3 4.3 状态空间表示的系统稳定性判定状态空间表示的系统稳定性判定定理定理5.1:5.1: 线性定常系统线性定常系统 xAxbuycxdu+平衡状态平衡状态 渐近稳定的充要条件是矩阵渐近稳定的充要条件是矩阵A A的所有特征值均具有负实部的所有特征值均具有负实部. . 0ex 证明：证明：，由其齐次解，由其齐次解 可知：若可知：若A A的的特征特征0( )( )Atx te x t则当则当

32、 有界，有界， 0(t)0(t)。0 x( )x t值均具有负实部。值均具有负实部。 系统输出稳定：系统输出稳定：如果系统对于有界输入如果系统对于有界输入u u 所引起的输出所引起的输出y y是有是有 界的界的. .则称系统为输出稳定则称系统为输出稳定. .定理定理5.25.2：线性定常系统线性定常系统 输出稳定的充要条件是传输出稳定的充要条件是传 函函 的极点全部位于的极点全部位于s s的左半平面的左半平面. . ( , , )A b c1( )()G Sc SIAb2022-3-9北京科技大学自动化学院自动化系北京科技大学自动化学院自动化系40 设系统的状态空间表达式为设系统的状态空间表达式为: :1011,0011xxuyx + 试分析系统的状态稳定性与输出稳定性试分析系统的状态稳定性与输出稳定性. . 1) 1)有有A A的特征方程的特征方