第8行: |
第8行: |
| {{short description|Part of mathematics that addresses the stability of solutions}} | | {{short description|Part of mathematics that addresses the stability of solutions}} |
| | | |
− | 在数学上,<font color="#ff8000">稳定性理论 Stability theory</font>被用于研究<font color="#ff8000">微分方程Differential equation</font>解的稳定性和<font color="#ff8000">动力系统 Dynamical system</font>在初始条件的微小扰动下轨迹的稳定性问题。例如,<font color="#ff8000">热传导方程 Heat equation</font>是一个稳定的偏微分方程,因为<font color="#ff8000">极大值原理 Maximum principle</font>的存在,初始数据的微小扰动会导致温度随之产生微小的变化。在偏微分方程中,人们可以使用 <math>Lp</math> 范数或 <math>sup</math> 范数来度量函数之间的距离,而在微分几何中,人们可以使用Gromov–Hausdorff距离来度量空间之间的距离<ref> Duplij S . Gromov–Hausdorff Distance[M]. 2003.</ref>。 | + | 在数学上,'''稳定性理论 Stability theory'''被用于研究'''微分方程Differential equation'''解的稳定性和'''动力系统 Dynamical system'''在初始条件的微小扰动下轨迹的稳定性问题。例如,'''热传导方程 Heat equation'''是一个稳定的偏微分方程,因为'''极大值原理 Maximum principle'''的存在,初始数据的微小扰动会导致温度随之产生微小的变化。在偏微分方程中,人们可以使用 <math>Lp</math> 范数或 <math>sup</math> 范数来度量函数之间的距离,而在微分几何中,人们可以使用Gromov–Hausdorff距离来度量空间之间的距离<ref> Duplij S . Gromov–Hausdorff Distance[M]. 2003.</ref>。 |
| | | |
| | | |
− | 在动力系统中,如果一条<font color="#ff8000">轨道 Orbit</font>上任意点的前向轨道都处于一个足够小的邻域内,或者这条轨道整体停留在一个邻域(一般是较小的邻域,也有可能是较大的邻域)内,则称该轨道的状态为<font color="#ff8000">李雅普诺夫稳定 Lyapunov stable</font>。有各种标准来证明轨道的稳定性或不稳定性。在适当的条件下,这个问题可以简化为一个涉及矩阵<font color="#ff8000">特征值 Eigenvalue</font>的问题,关于这类矩阵特征值的问题已被大量研究并且该领域已经比较成熟。一种更一般的方法涉及<font color="#ff8000">李雅普诺夫函数 Lyapunov function</font>。在实践中,很多<font color="#ff8000">稳定性判据 Stability criterion</font>都可以使用,我们可以使用其中的任何一个作为判断系统稳定性的准则。 | + | 在动力系统中,如果一条'''轨道 Orbit'''上任意点的前向轨道都处于一个足够小的邻域内,或者这条轨道整体停留在一个邻域(一般是较小的邻域,也有可能是较大的邻域)内,则称该轨道的状态为'''李雅普诺夫稳定 Lyapunov stable'''。有各种标准来证明轨道的稳定性或不稳定性。在适当的条件下,这个问题可以简化为一个涉及矩阵'''特征值 Eigenvalue'''的问题,关于这类矩阵特征值的问题已被大量研究并且该领域已经比较成熟。一种更一般的方法涉及'''李雅普诺夫函数 Lyapunov function'''。在实践中,很多'''稳定性判据 Stability criterion'''都可以使用,我们可以使用其中的任何一个作为判断系统稳定性的准则。 |
− | [[File:Stability_Diagram.png|thumb|550px|稳定性图将<font color="#ff8000">庞加莱映射 Poincaré map</font> 根据其特征划分为稳定或不稳定区间。如图可见,图中下半部分区域中系统的稳定性增加。<ref>[http://www.egwald.ca/linearalgebra/lineardifferentialequationsstabilityanalysis.php Egwald Mathematics - Linear Algebra: Systems of Linear Differential Equations: Linear Stability Analysis] Accessed 10 October 2019.</ref>|链接=Special:FilePath/Stability_Diagram.png]] | + | [[File:Stability_Diagram.png|thumb|550px|稳定性图将'''庞加莱映射 Poincaré map''' 根据其特征划分为稳定或不稳定区间。如图可见,图中下半部分区域中系统的稳定性增加。<ref>[http://www.egwald.ca/linearalgebra/lineardifferentialequationsstabilityanalysis.php Egwald Mathematics - Linear Algebra: Systems of Linear Differential Equations: Linear Stability Analysis] Accessed 10 October 2019.</ref>|链接=Special:FilePath/Stability_Diagram.png]] |
| | | |
| | | |
− | 常微分方程在经历了长期的求精确解的努力后逐渐停滞,庞加莱在分析的基础上引入几何方法,开创了常微分方程定性理论,同时在分析中引入几何方法,搭建起分析与几何之间的沟通桥梁,带来了微分方程研究的新突破。李雅普诺夫则在庞加莱定性分析的基础上 ,转而进入了新的稳定性研究。如今 ,李雅普诺夫稳定性理论被普遍认为是微分方程定性理论的基本成就之一。不仅有精确的定义 ,更有严格的分析证明 ,将微分方程及稳定性理论的研究推向了新的高度。庞加莱被公认是19世纪后四分之一和二十世纪初的领袖数学家,是对于数学和它的应用具有全面知识的最后一个人,他在数学方面的杰出工作对20世纪和当今的数学造成极其深远的影响。<font color="#ff8000">庞加莱映射 Poincaré map</font>是由相空间中轨道运动定义的一种映射,是当轨道反复穿越同一截面时,反映后继点对先行点依赖关系的映射<ref>Perko. Differential equations and dynamical systems[M]. Springer, 2001.</ref>。一个连续非线性动力系统的求解是非常困难的,庞加莱给出了相图分析法。在相图中虽然不能定量地知道物理量随时间的变化,但可以定性地得到轨线的形态类型及其拓扑结构,从而了解动力系统运动的全局图像。为了更清楚了解高维相空间运动的形态,在连续运动的轨线上用一个截面(称庞加莱截面)将其横截,轨线在截面上穿过的情况就可以简捷地判断运动的形态。对于庞加莱映射是稳定的还是不稳定的判断则取决于其特征,如图所示,在相空间区间中向下的方向上稳定性增加。 | + | 常微分方程在经历了长期的求精确解的努力后逐渐停滞,庞加莱在分析的基础上引入几何方法,开创了常微分方程定性理论,同时在分析中引入几何方法,搭建起分析与几何之间的沟通桥梁,带来了微分方程研究的新突破。李雅普诺夫则在庞加莱定性分析的基础上 ,转而进入了新的稳定性研究。如今 ,李雅普诺夫稳定性理论被普遍认为是微分方程定性理论的基本成就之一。不仅有精确的定义 ,更有严格的分析证明 ,将微分方程及稳定性理论的研究推向了新的高度。庞加莱被公认是19世纪后四分之一和二十世纪初的领袖数学家,是对于数学和它的应用具有全面知识的最后一个人,他在数学方面的杰出工作对20世纪和当今的数学造成极其深远的影响。'''庞加莱映射 Poincaré map'''是由相空间中轨道运动定义的一种映射,是当轨道反复穿越同一截面时,反映后继点对先行点依赖关系的映射<ref>Perko. Differential equations and dynamical systems[M]. Springer, 2001.</ref>。一个连续非线性动力系统的求解是非常困难的,庞加莱给出了相图分析法。在相图中虽然不能定量地知道物理量随时间的变化,但可以定性地得到轨线的形态类型及其拓扑结构,从而了解动力系统运动的全局图像。为了更清楚了解高维相空间运动的形态,在连续运动的轨线上用一个截面(称庞加莱截面)将其横截,轨线在截面上穿过的情况就可以简捷地判断运动的形态。对于庞加莱映射是稳定的还是不稳定的判断则取决于其特征,如图所示,在相空间区间中向下的方向上稳定性增加。 |
| | | |
| ==动力系统概述== | | ==动力系统概述== |
| | | |
− | 微分方程和动力系统定性理论的许多部分关心系统或者方程解的渐近性质及其轨迹,这也意味着系统经过很长时间后会发生什么 <ref>Palis J , Melo W D . Geometric theory of dynamical systems[J]. Springer-Verlag, 1982, 10.1007/978-1-4612-5703-5.</ref>。系统最简单的行为表现为<font color="#ff8000">平衡点 Equilibrium points</font>或不动点,以及<font color="#ff8000">周期轨道 Periodic orbit</font>。如果我们已经很好地理解了一个特定的轨道,那么很自然地就会问下一个问题:初始条件的一个微小变化对于系统来说是否仍会保持类似的行为。稳定性理论解决了以下问题:附近的轨道是否会无限靠近给定的轨道?已知的轨道会收敛到给定的轨道吗?在前一种情况下,轨道被称为是<font color="#ff8000">稳定 Stable</font>的;在后一种情况下,轨道是<font color="#ff8000">渐近稳定 Asymptotically stable </font>的,并且收敛到给定的轨道称为<font color="#ff8000">吸引子 Attractor</font><ref>Zaslavsky G M . The simplest case of a strange attractor[J]. Physics Letters A, 1978, 69(3):145-147.</ref>。 | + | 微分方程和动力系统定性理论的许多部分关心系统或者方程解的渐近性质及其轨迹,这也意味着系统经过很长时间后会发生什么 <ref>Palis J , Melo W D . Geometric theory of dynamical systems[J]. Springer-Verlag, 1982, 10.1007/978-1-4612-5703-5.</ref>。系统最简单的行为表现为'''平衡点 Equilibrium points'''或不动点,以及'''周期轨道 Periodic orbit'''。如果我们已经很好地理解了一个特定的轨道,那么很自然地就会问下一个问题:初始条件的一个微小变化对于系统来说是否仍会保持类似的行为。稳定性理论解决了以下问题:附近的轨道是否会无限靠近给定的轨道?已知的轨道会收敛到给定的轨道吗?在前一种情况下,轨道被称为是'''稳定 Stable'''的;在后一种情况下,轨道是'''渐近稳定 Asymptotically stable '''的,并且收敛到给定的轨道称为'''吸引子 Attractor'''<ref>Zaslavsky G M . The simplest case of a strange attractor[J]. Physics Letters A, 1978, 69(3):145-147.</ref>。 |
| | | |
| 对于一个一阶常微分方程自治系统的平衡解<math>f_e</math>: | | 对于一个一阶常微分方程自治系统的平衡解<math>f_e</math>: |
第29行: |
第29行: |
| 稳定性意味着在微小的扰动下轨迹不会发生太大的变化。相反的情况中,微小的扰动会使得轨迹发生较大变化,即附近的轨道与给定的轨道互相排斥,这也是一种有趣的现象。一般来说,在某些方向对初始状态的扰动使得轨道渐近地接近给定轨道,而在其他方向的扰动则使得轨道远离给定轨道。也可能存在对初始状态在某些方向的扰动使得轨道行为变得比较复杂(比如既不会收敛也不会完全逃逸),从而稳定性理论不能对于这样的动力学状态给予充分的预测信息。 | | 稳定性意味着在微小的扰动下轨迹不会发生太大的变化。相反的情况中,微小的扰动会使得轨迹发生较大变化,即附近的轨道与给定的轨道互相排斥,这也是一种有趣的现象。一般来说,在某些方向对初始状态的扰动使得轨道渐近地接近给定轨道,而在其他方向的扰动则使得轨道远离给定轨道。也可能存在对初始状态在某些方向的扰动使得轨道行为变得比较复杂(比如既不会收敛也不会完全逃逸),从而稳定性理论不能对于这样的动力学状态给予充分的预测信息。 |
| | | |
− | 稳定性理论的关键思想之一是用轨道附近系统的线性化,来分析轨道在扰动下的定性行为。特别地,在 n 维<font color="#ff8000">相空间 Phase space</font>的光滑动力系统的每个平衡点上,都存在一个 n×n 的矩阵 A,其特征值刻画了邻近点的动力学行为<font color="#ff8000">(Hartman-Grobman 定理 Hartman–Grobman theorem)</font>。更确切地说,如果矩阵所有的特征值都是负实数或实部为负的复数,那么这个平衡点就是一个稳定的吸引子,并且附近的点以指数速率收敛到它,参考<font color="#ff8000">李雅普诺夫稳定性 Lyapunov stability</font>和<font color="#ff8000">指数稳定性 Exponential stability</font>。如果所有的特征值都不是纯虚数(或零) ,那么吸引方向和排斥方向都与矩阵 A 的特征空间有关,其特征值的实部分别为负和正。对于更复杂的轨道上的扰动情形,也有类似的表述。 | + | 稳定性理论的关键思想之一是用轨道附近系统的线性化,来分析轨道在扰动下的定性行为。特别地,在 n 维'''相空间 Phase space'''的光滑动力系统的每个平衡点上,都存在一个 n×n 的矩阵 A,其特征值刻画了邻近点的动力学行为'''(Hartman-Grobman 定理 Hartman–Grobman theorem)'''。更确切地说,如果矩阵所有的特征值都是负实数或实部为负的复数,那么这个平衡点就是一个稳定的吸引子,并且附近的点以指数速率收敛到它,参考'''李雅普诺夫稳定性 Lyapunov stability'''和'''指数稳定性 Exponential stability'''。如果所有的特征值都不是纯虚数(或零) ,那么吸引方向和排斥方向都与矩阵 A 的特征空间有关,其特征值的实部分别为负和正。对于更复杂的轨道上的扰动情形,也有类似的表述。 |
| | | |
| ==不动点稳定性== | | ==不动点稳定性== |
第67行: |
第67行: |
| 如我们所知,线性系统是一类数学模型,指的是由线性运算子组成的系统,也就是说,这类系统首先满足线性的特性<ref>Luenberger D G . Observing the State of a Linear System[J]. IEEE Transactions on Military Electronics, 2007, 8(2):74-80.</ref>。相较于非线性系统,线性系统的特性比较简单。 | | 如我们所知,线性系统是一类数学模型,指的是由线性运算子组成的系统,也就是说,这类系统首先满足线性的特性<ref>Luenberger D G . Observing the State of a Linear System[J]. IEEE Transactions on Military Electronics, 2007, 8(2):74-80.</ref>。相较于非线性系统,线性系统的特性比较简单。 |
| | | |
− | 根据系统矩阵A是否随时间变化,引入<font color="#ff8000">自治系统 autonomous system</font>的概念后,可以把线性系统分为自治的和非自治的,对于线性系统一般也可以称为定常的和时变的,也就是说: | + | 根据系统矩阵A是否随时间变化,引入'''自治系统 autonomous system'''的概念后,可以把线性系统分为自治的和非自治的,对于线性系统一般也可以称为定常的和时变的,也就是说: |
| | | |
| (1)自治的线性系统就是定常线性系统。 | | (1)自治的线性系统就是定常线性系统。 |
第79行: |
第79行: |
| | | |
| | | |
− | 对于一个如下的<font color="#ff8000">自治系统 autonomous system</font> | + | 对于一个如下的'''自治系统 autonomous system''' |
| :<math>x' = Ax,</math> | | :<math>x' = Ax,</math> |
| | | |
第92行: |
第92行: |
| 可以这样描述:最初的原点({{Math|0 ∈ '''R'''<sup>''n''</sup>}} ) 是该动力系统的平衡点。当且仅当对于{{Math|''A''}}的所有特征值 {{Math|''λ''}} 有 {{Math|Re(''λ'') < 0}} 时,这个解是随着{{Math|''t'' → ∞}}是渐近稳定的(未来趋势)。类似地,当且仅当对于 {{Math|''A''}} 的所有特征值 {{Math|''λ''}} 有{{Math|Re(''λ'') > 0}} 时,系统随着{{Math|''t'' → -∞}}是渐近稳定的(负号表示方向指向过去趋势)。如果存在一个{{Math|''A''}}的特征值 {{Math|''λ''}} 使得 {{Math|Re(''λ'') > 0}},则该解在{{Math|''t'' → ∞}}时是不稳定的。 | | 可以这样描述:最初的原点({{Math|0 ∈ '''R'''<sup>''n''</sup>}} ) 是该动力系统的平衡点。当且仅当对于{{Math|''A''}}的所有特征值 {{Math|''λ''}} 有 {{Math|Re(''λ'') < 0}} 时,这个解是随着{{Math|''t'' → ∞}}是渐近稳定的(未来趋势)。类似地,当且仅当对于 {{Math|''A''}} 的所有特征值 {{Math|''λ''}} 有{{Math|Re(''λ'') > 0}} 时,系统随着{{Math|''t'' → -∞}}是渐近稳定的(负号表示方向指向过去趋势)。如果存在一个{{Math|''A''}}的特征值 {{Math|''λ''}} 使得 {{Math|Re(''λ'') > 0}},则该解在{{Math|''t'' → ∞}}时是不稳定的。 |
| | | |
− | 为了判定线性系统原点的稳定性,可以使用劳斯-赫尔维茨稳定性判据<font color="#ff8000">Routh–Hurwitz stability criterion</font>,来将这一结果应用在实践中。矩阵的特征值是其特征多项式的根。如果所有根的实部都是严格负的,那么一个具有实系数的单变量多项式称为赫尔维茨多项式 <font color="#ff8000">Hurwitz polynomial</font> 。劳斯-赫尔维茨定理 <font color="#ff8000">Routh–Hurwitz theorem</font>通过一种避免计算根的算法来描述赫尔维茨多项式的特征。 | + | 为了判定线性系统原点的稳定性,可以使用劳斯-赫尔维茨稳定性判据'''Routh–Hurwitz stability criterion''',来将这一结果应用在实践中。矩阵的特征值是其特征多项式的根。如果所有根的实部都是严格负的,那么一个具有实系数的单变量多项式称为赫尔维茨多项式 '''Hurwitz polynomial''' 。劳斯-赫尔维茨定理 '''Routh–Hurwitz theorem'''通过一种避免计算根的算法来描述赫尔维茨多项式的特征。 |
| | | |
| ===非线性自治系统=== | | ===非线性自治系统=== |
第107行: |
第107行: |
| | | |
| | | |
− | 设{{Math|''J''<sub>''p''</sub>(''v'')}}为向量场 {{Math|''v''}}在点{{Math|''p''}}的{{Math|''n''×''n''}}<font color="#ff8000">雅可比矩阵 Jacobian matrix</font>。如果 {{Math|''J''}} 的所有特征值都具有严格负的实部,则系统的解是渐近稳定的。这个条件可以用劳斯-赫尔维茨判据<font color="#ff8000">Routh–Hurwitz stability criterion</font>来检验。 | + | 设{{Math|''J''<sub>''p''</sub>(''v'')}}为向量场 {{Math|''v''}}在点{{Math|''p''}}的{{Math|''n''×''n''}}'''雅可比矩阵 Jacobian matrix'''。如果 {{Math|''J''}} 的所有特征值都具有严格负的实部,则系统的解是渐近稳定的。这个条件可以用劳斯-赫尔维茨判据'''Routh–Hurwitz stability criterion'''来检验。 |
| | | |
| ==一般动力系统的李雅普诺夫函数== | | ==一般动力系统的李雅普诺夫函数== |