实验4-Laypunov方程求解

《现代控制理论》实验四院系：学生姓名：学号：

一：原理1.李雅普诺夫稳定性概念忽略输入后，非线性时变系统的状态方程如下（1）式中，为维状态向量；为时间变量；为维函数，其展开式为假定方程的解为，和分别为初始状态向量和初始时刻，。平衡状态如果对于所有，满足（2）的状态称为平衡状态（又称为平衡点）。平衡状态的各分量不再随时间变化。若已知状态方程，令所求得的解，便是平衡状态。对于线性定常系统，其平衡状态满足，如果非奇异，系统只有惟一的零解，即存在一个位于状态空间原点的平衡状态。至于非线性系统，的解可能有多个，由系统状态方程决定。控制系统李雅普诺夫意义下的稳定性是关于平衡状态的稳定性，反映了系统在平衡状态附近的动态行为。鉴于实际线性系统只有一个平衡状态，平衡状态的稳定性能够表征整个系统的稳定性。对于具有多个平衡状态的非线性系统来说，由于各平衡状态的稳定性一般并不相同，故需逐个加以考虑，还需结合具体初始条件下的系统运动轨迹来考虑。本节主要研究平衡状态位于状态空间原点（即零状态）的稳定性问题，因为任何非零状态均可以通过坐标变换平移到坐标原点，而坐标变换又不会改变系统的稳定性将上式展开，的每一元素都是的线性组合，因而可写成矩阵多项式故可以显式表出与的关系当式（8-74）成立时，对于任意，均有，系统渐近稳定。只要有一个特征值的实部大于零，对于,便无限增长，系统不稳定。如果只有一个（或一对，且均不能是重根）特征值的实部等于零，其余特征值实部均小于零，便含有常数项或三角函数项，则系统是李雅普诺夫意义下稳定的。4.李雅普诺夫稳定性直接判别法李雅普诺夫第二法（直接法）是利用李雅普诺夫函数直接对平衡状态稳定性进行判断，无需求出系统状态方程的解，它对各种控制系统均适用。根据物理学原理，若系统贮存的能量（含动能与位能）随时间推移而衰减，系统迟早会到达平衡状态。实际系统的能量函数表达式相当难找，因此李雅普诺夫引入了广义能量函数，称之为李雅普诺夫函数。它与及有关，是一个标量函数，记以；若不显含，则记以。考虑到能量总大于零，故为正定函数，能量衰减特性用表示。遗憾的是至今仍未形成构造李雅普诺夫函数的通用方法，需要凭经验与技巧。实践表明，对于大多数系统，可先尝试用二次型函数作为李雅普诺夫函数。1.标量函数定号性 （1）正定性标量函数在域中对所有非零状态有且，称在域内正定。如是正定的。（2）负定性标量函数在域中对所有非零有且，称在域内负定。如是负定的。如果是负定的，-则一定是正定的。（3）负（正）半定性，且在域内某些状态处有，而其它状态处均有（），则称在域内负（正）半定。设为负半定，则为正半定。如为正半定。（4）不定性在域内可正可负，则称不定。如是不定的。关于正定性的提法是：标量函数在域中，对于及所有非零状态有,且，则称在域内正定。的其它定号性提法类同。二次型函数是一类重要的标量函数，记（6）其中，为对称矩阵，有。显然满足，其定号性由赛尔维斯特准则判定。当的各顺序主子行列式均大于零时，即（7）为正定矩阵，则正定。当的各顺序主子行列式负、正相间时，即（8）为负定矩阵，则负定。若主子行列式含有等于零的情况，则为正半定或负半定。不属以上所有情况的不定。下面不打算对李雅普诺夫第二法中诸稳定性定理在数学上作严格证明，而只着重于物理概念的阐述和应用。2.李雅普诺夫第二法诸稳定性定理设系统状态方程为，其平衡状态满足,不失一般性，把状态空间原点作为平衡状态，并设系统在原点邻域存在对的连续的一阶偏导数。定理1若①正定，②负定；则原点是渐近稳定的。负定表示能量随时间连续单调地衰减，故与渐近稳定性定义叙述一致。定理2若①正定；②负半定，且在非零状态不恒为零；则原点是渐近稳定的。负半定表示在非零状态存在，但在从初态出发的轨迹上，不存在的情况，于是系统将继续运行至原点。状态轨迹仅是经历能量不变的状态，而不会维持在该状态。定理3若①正定；②负半定，且在非零状态恒为零；则原点是李雅普诺夫意义下稳定的。沿状态轨迹能维持，表示系统能维持等能量水平运行，使系统维持在非零状态而不运行至原点。定理4若①正定；②正定；则原点是不稳定的。正定表示能量函数随时间增大，故状态轨迹在原点邻域发散。参考定理2可推论：正定，当正半定，且在非零状态不恒为零时，则原点不稳定。应注意到，李雅普诺夫函数[正定的]的选取是不惟一的，但只要找到一个满足定理所述条件，便可对原点的稳定性作出判断，并不因选取的不同而有所影响。不过至今尚无构造李雅普诺夫函数的通用方法，这是应用李雅普诺夫稳定性理论的主要障碍。如果选取不当，会导致不定的结果，这时便作不出确定的判断，需要重新选取。以上定理按照连续单调衰减的要求来确定系统稳定性，并未考虑实际稳定系统可能存在衰减振荡的情况，因此其条件是偏于保守的，故借稳定性定理判稳定者必稳定，李雅普诺夫第二法诸稳定性定理所述条件都是充分条件。具体分析时，先构造一个李雅普诺夫函数，通常选二次型函数，求其导数，再将状态方程代入，最后根据的定号性判别稳定性。至于如何判断在非零状态下是否有恒为零的情况，可按如下方法进行：令，将状态方程代入，若能导出非零解，表示对，的条件是成立的；若导出的是全零解，表示只有原点满足的条件。二：程序：A=[01;-1-1];Q=eye(size(A,1));p=lyap