稳定性

*13)正确理解李雅普诺夫稳定性的有关概念; 14)掌握李雅普诺夫第二法,初步掌握寻求系统李雅普诺夫函数判定系统稳定性的方法。

stability

Lyapunov

这个章节不是状态空间的主线内容,相对独立,初学放到最后也不影响的。

注意这个定语,不是“线性系统”了,是范围更大的“控制系统”。

到了现代控制理论的这里,“稳定”有了精确的数学定义和不同的稳定状态。前面经典控制里面我们认为系统最后靠近一个值就是稳定的(粗略的讲)。现在遇到的问题是比如极限环自振,说他不稳吧,他和那种系统跑飞崩溃还不一样,说他稳吧,这还在动。所以有了稳有了不同的稳法。

关于稳定的一些提法,最核心的是:李雅普诺夫稳定性、渐进稳定性、大范围渐进稳定。还有像极限定义一样的对稳定的数学描述。

当然稳定不是为了玩弄概念,我们是为了在拿到一个控系统系统后,研究其最根本的稳不稳问题。稳了再谈后面设计校正优化的事情。

李雅普诺夫建立了基于状态空间描述的稳定性理论,提出了依赖于线性系统微分方程的解来判断稳定性的第一方法(称为间接法〉和利用经验和技巧来构造李雅普诺夫函数借以判断稳定性的第二方法(称为直接法)。

李雅普诺夫提出的这一理论是确定系统稳定性的更一般的理论,不仅适用于单变量,线性﹑定常系统,还适用于多变量、非线性、时变系统。

各种稳定性的定义

这里的稳定针对的系统可以上升到非线性时变的情况。系统的状态方程:

x˙=f(x,t)\boldsymbol{\dot{x}} = \boldsymbol{f}(\boldsymbol{x},t)

首先要想想平衡状态(平衡点)该怎么用数学表达呢?状态不再随时间运动就是平衡了,数学表达就是

xe˙=f(xe,t)=0\boldsymbol{\dot{x_e}} = \boldsymbol{f}(\boldsymbol{x_e},t) = \boldsymbol{0}

线性定常系统,平衡状态有$$ \boldsymbol{Ax_e} = \boldsymbol{0} $$,根据数学原理,如果系统矩阵$$A$$非奇异(满秩,可逆,行列式不为0),那么经有唯一0解。也就是说平衡状态(平衡点)在状态空间的零点。

从一个高中物理的小例子出发,一个滑块在圆弧轨道里运动的问题。

滑块A的轨道光滑,可以想想A一直做来回运动,永不停止。滑块B绝对是不稳的,只要动一下就回不到平衡点了。滑块C的轨道是有摩擦力的,最终能量消耗完了,就停在最低点了,就是稳定点。

分别对应了李雅普诺夫意义下的稳定、不稳定、大范围渐近稳定。

李雅普诺夫稳定性

状态空间轨迹在稳定点的附近运动。

系统做不衰减的震荡也可以认为是稳定的。或者稳定极限环。

数学定义和极限很像,不过状态空间里的抽象距离,用范数取讲这个事情

\Vert \Vert,敲公式太麻烦了。后面有空在说吧。

总之这个数学定义相当严谨的,要求状态初始点在状态空间里距离稳定点不远,时间往大走,状态点都离稳定点不远,大概就这么个事情。

要求初始点这个事情在全局稳定里就理解为啥这里要这个前提了。

对于非线性系统,李雅普诺夫稳定性可以进行一阶近似丢掉高阶项进行分析,因为讲的是领域嘛,可以这么搞。而且很合理。稳定点线性化九的就是领域的事情,李雅普诺夫稳定也是讲的领域的事情,这就正好对上了。

渐近稳定性

状态空间轨迹最终停在平衡点了。经典控制理论里面的稳定全是渐近稳定。

也可以看出渐进稳定,首先得李雅普诺夫稳定。并且要求最后状态点就停在稳定点了,当然也包含了状态初始点在稳定点附近这个前提。

大范围渐进稳定性

整个状态空间,不管状态初始点在哪里,轨迹都能运动到平衡点。

非线性系统稳定和初始值也有关系(具体可以随便拿个相图看看),线性系统倒是只和结构有关。这里稳定性是针对所有系统去谈,所以前面强调初始点要离平衡点不远,如果任意初始点,那就是大范围稳定了。

在前面非线性系统,用二阶相图近似分析的时候,体会一下。

BIBS稳定性

任意有界的初始状态,在任意有界的输入下,状态均有界

线性定常系统,BIBO判断主要依据系统矩阵A,特征值在左半开平面,则BIBS稳定。

BIBO稳定性

任意有界的初始状态,在任意有界的输入下,输出均有界

线性定常系统,BIBO判断主要依据传递函数,极点在左半开平面,则BIBO稳定。

对于线性定常系统,如果渐进稳定,则必然BIBS稳定,BIBO稳定。

但是如果李雅普诺夫稳定,系统不一定BIBS稳定,BIBO稳定。

李雅普诺夫稳定性间接判别法

第一法是用状态方程的解来判断系统稳定性的方法。

对于线性系统x˙=Ax\dot{x} = Ax来讲线性定常系统的特征值判据:系统矩阵A的全部特征值在复平面左半部

这是渐近稳定充要条件:

Re(λi)0    (i=1,2,,n)Re(\lambda_i) \leq 0 \ \ \ \ (i=1,2,\cdots,n)

即系统矩阵A的全部特征值位于复平面左半部。这个经典理论相通。

对于非线性系统来说,首先可以使用线性化的手段。实际上在前面非线性分析里已经有过这个行为了。用二阶系统相图近似分析平衡点稳定性。

因此,对于非线性系统的稳定性分析,先考虑如何进行线性化。一元函数经典泰勒。二元也是,二元泰勒

{x1˙=f1(x1,x2)x2˙=f2(x1,x2)\begin{array}{} \left\{ \begin{align} \dot{x_1} &= f_1(x_1,x_2) \\ \dot{x_2} &= f_2(x_1,x_2) \\ \end{align} \right. \end{array}

先得找到平衡点,毕竟是在平衡点附近做线性化。如果平衡点在xe=(x1,x2)=(0,0)x_e = (x_1,x_2) =(0,0)

那么线性化之后的系统矩阵就是

\begin{array} {} \dot{x_1} = x_2 - x_1 ( x_1^2 + x_2^2 ) \ \dot{x_2} = -x_1-x_2 ( x_1^2 + x_2^2 ) \end{array}

Last updated