152
个编辑
更改
跳到导航
跳到搜索
第1行：
第1行： − + + 第11行：
第12行： + − + − 该轨道可以看作某些其他轨道的极限集，即至少存在一条不同于它本身的轨道，在 <math>t \to +\infty</math> 或 <math>t \to -\infty</math> 时趋近于该闭合轨道。<ref name="Perko">Perko, L. (2001). ''Differential Equations and Dynamical Systems''. Springer.</ref><ref name="Devaney">Devaney, R. (1989). ''An Introduction to Chaotic Dynamical Systems''. CRC Press.</ref> 值得注意的是，一维自治系统中不可能存在极限环，因为一维相空间中的轨道只能单调递增、单调递减或停留于平衡点，无法构成闭合曲线。+ 第36行：
第38行： − + 第52行：
第54行： − + 第134行：
第136行： − + 第150行：
第152行： − 李雅普诺夫函数方法适用于任意 <math>n</math> 维自治系统，比梯度系统更为灵活。然而，该方法的主要局限在于需要预先猜测或构造出合适的 <math>V</math>，而这在一般情形下并无系统性方法。+ 第162行：
第164行： − 综上，三种方法在适用范围上各有侧重：梯度系统与李雅普诺夫函数方法均适用于任意 <math>n</math> 维系统，前者要求系统具有势函数结构，后者则更为灵活但缺乏系统性的构造方法；Dulac 准则通过引入权函数 <math>B</math> 在二维框架内进一步扩大了适用范围，但无法推广至高维。+ 第223行：
第225行： − + − 尽管二维的专属工具在高维失效，部分方法仍可直接推广至任意 <math>n</math> 维系统，分别用于排除和确认极限环的存在。+ − + + − + − − : <math>\dot{\boldsymbol{\xi}} = J(\mathbf{x}^*(t))\boldsymbol{\xi},</math> − 其中 <math>J(\mathbf{x}^*(t))</math> 为系统在 <math>\mathbf{x}^*(t)</math> 处的雅可比矩阵，是以 <math>T</math> 为周期的矩阵函数。<ref name="Chicone">Chicone, C. (2006). ''Ordinary Differential Equations with Applications'' (2nd ed.). Springer.</ref>+ − 该线性系统的'''单值矩阵'''（monodromy matrix）<math>M</math> 定义为线性变分方程在一个周期 <math>T</math> 内的基本解矩阵，其特征值称为'''Floquet 乘子'''（Floquet multipliers），又称'''特征乘子'''（characteristic multipliers）。对于 <math>n</math> 维系统，极限环对应的庞加莱映射有 <math>n-1</math> 个特征乘子，其模决定了极限环的稳定性：<ref name="Kuznetsov" /><ref name="Wiggins" />+ − + − + − + − 注意，由于沿极限环方向本身对应一个模为 1 的 Floquet 乘子（对应轨道的切方向），因此实际决定横截稳定性的乘子共有 <math>n-1</math> 个。在二维情形下，<math>n-1=1</math>，退化为单个特征乘子，与[[#稳定性分类]]一节中的讨论一致。<ref name="Wiggins" />+ − + 第266行：
第267行： − + 第276行：
第277行： − +
无编辑摘要
{{#seo:
{{#seo:
|keywords=稳定极限环,相空间,范德波尔振荡器
|keywords=稳定极限环,相空间,范德波尔振荡器
|description=是一个在相空间中的闭合轨迹，它具有当时间趋于无穷大或时间趋于负无穷大时至少有一条其他轨迹螺旋进入的性质。
|description=极限环是非线性动力系统相空间中的孤立周期轨道，用于描述系统长期稳定的周期性行为。
}}
}}
极限环（limit cycle）是非线性动力系统中的重要概念，用于描述系统在长期演化过程中所呈现的稳定周期行为。当系统的轨道不收敛于平衡点，而是逐渐趋近于一条孤立的闭合轨道时，系统的运动最终表现为沿该闭合轨道的周期运动，这样的闭轨道称为极限环。
极限环（limit cycle）是非线性动力系统中的重要概念，用于描述系统在长期演化过程中所呈现的稳定周期行为。当系统的轨道不收敛于平衡点，而是逐渐趋近于一条孤立的闭合轨道时，系统的运动最终表现为沿该闭合轨道的周期运动，这样的闭轨道称为极限环。
本文首先在一般 n 维动力系统的框架下给出极限环的严格数学定义及其基本性质，随后介绍二维系统中关于极限环的主要判定理论，最后概述高维系统中极限环的若干重要结果及相关研究问题。
本文首先在一般 n 维动力系统的框架下给出极限环的严格数学定义及其基本性质，随后介绍二维系统中关于极限环的主要判定理论，最后概述高维系统中极限环的若干重要结果及相关研究问题。
==定义==
==定义==
极限环的标准定义是相空间中的孤立周期轨道。
'''极限环的定义：'''极限环是 <math>n</math> 维（<math>n \geq 2</math>）动力系统相空间中一条孤立的闭合轨道，且在其某个邻域内不存在其他周期轨道。
'''极限环的定义：'''极限环的标准定义是 <math>n</math> 维（<math>n \geq 2</math>）动力系统相空间中一条孤立周期轨道，且在其某个邻域内不存在其他周期轨道。
在许多动力系统中，极限环同时还是邻近轨道的 α 极限集或 ω 极限集，因此刻画了系统的长期周期行为。，即至少存在一条不同于它本身的轨道，在 <math>t \to +\infty</math> 或 <math>t \to -\infty</math> 时趋近于该闭合轨道。<ref name="Perko">Perko, L. (2001). ''Differential Equations and Dynamical Systems''. Springer.</ref><ref name="Devaney">Devaney, R. (1989). ''An Introduction to Chaotic Dynamical Systems''. CRC Press.</ref> 值得注意的是，一维自治系统中不可能存在极限环，因为一维相空间中的轨道只能单调递增、单调递减或停留于平衡点，无法构成闭合曲线。
为使用更形式化的语言描述这一概念，现考虑一个 <math>n</math> 维（<math>n \geq 2</math>）自治动力系统
为使用更形式化的语言描述这一概念，现考虑一个 <math>n</math> 维（<math>n \geq 2</math>）自治动力系统
极限环是相空间中的'''不变集'''：若系统的初始条件恰好落在该闭轨道上，则对应解将永远沿该轨道周期性演化。<ref name="Strogatz"/>
极限环是相空间中的'''不变集'''：若系统的初始条件恰好落在该闭轨道上，则对应解将永远沿该轨道周期性演化。<ref name="Strogatz"/>
这一性质来源于微分方程解的'''存在唯一性定理'''（Picard定理）：相空间中过同一点只能有一条轨道。因此，一旦轨道进入极限环，它就不可能离开；反之，极限环外的轨道也不能穿入极限环。<ref name="LibreTexts">Lebl, J. ''Differential Equations for Engineers'', §8.4: Limit cycles. LibreTexts. https://math.libretexts.org</ref>
这一性质来源于微分方程解的'''存在唯一性定理'''（Picard定理）：相空间中过同一点只能有一条轨道。因此，若某轨道与极限环有一点相交，则由解的唯一性，该轨道必与该极限环重合<ref name="LibreTexts">Lebl, J. ''Differential Equations for Engineers'', §8.4: Limit cycles. LibreTexts. https://math.libretexts.org</ref>
=== 孤立性 ===
=== 孤立性 ===
* '''半稳定极限环'''（semi-stable limit cycle）：若极限环在某些方向上吸引邻近轨道，而在另一些方向上排斥邻近轨道，则称其为半稳定极限环。在二维系统中，这种情形常表现为极限环的一侧轨道趋近、另一侧轨道远离。半稳定极限环通常出现在系统参数的临界值附近，例如当一个稳定极限环与一个不稳定极限环合并并消失时所对应的'''鞍结分岔'''（saddle-node bifurcation of periodic orbits）。<ref name="Kuznetsov">Kuznetsov, Yu. A. (2004). ''Elements of Applied Bifurcation Theory'' (3rd ed.). Springer.</ref>
* '''半稳定极限环'''（semi-stable limit cycle）：若极限环在某些方向上吸引邻近轨道，而在另一些方向上排斥邻近轨道，则称其为半稳定极限环。在二维系统中，这种情形常表现为极限环的一侧轨道趋近、另一侧轨道远离。半稳定极限环通常出现在系统参数的临界值附近，例如当一个稳定极限环与一个不稳定极限环合并并消失时所对应的'''鞍结分岔'''（saddle-node bifurcation of periodic orbits）。<ref name="Kuznetsov">Kuznetsov, Yu. A. (2004). ''Elements of Applied Bifurcation Theory'' (3rd ed.). Springer.</ref>
[[文件:极限环类型.png|thumb|none稳定、不稳定与半稳定极限环示意图（二维情形），引自<ref name="Strogatz" />|600x600像素]]
[[文件:极限环类型.png|thumb|none|稳定、不稳定与半稳定极限环示意图（二维情形），引自<ref name="Strogatz" />|600x600像素]]
上述分类提供了对极限环稳定性的直观描述，但在一般维数的动力系统中，稳定性的严格判定通常需要借助'''庞加莱映射'''（Poincaré map）。其基本思想是：在相空间中选取一个与极限环横截相交的 <math>(n-1)</math> 维超曲面（称为'''庞加莱截面'''），记录轨道每次穿越该截面的位置，从而将连续时间系统的动力学转化为一个离散映射的不动点问题。<ref name="Wikipedia_PM">Poincaré map. ''Wikipedia''. https://en.wikipedia.org/wiki/Poincaré_map</ref>
上述分类提供了对极限环稳定性的直观描述，但在一般维数的动力系统中，稳定性的严格判定通常需要借助'''庞加莱映射'''（Poincaré map）。其基本思想是：在相空间中选取一个与极限环横截相交的 <math>(n-1)</math> 维超曲面（称为'''庞加莱截面'''），记录轨道每次穿越该截面的位置，从而将连续时间系统的动力学转化为一个离散映射的不动点问题。<ref name="Wikipedia_PM">Poincaré map. ''Wikipedia''. https://en.wikipedia.org/wiki/Poincaré_map</ref>
=== 极限环的排除方法 ===
=== 极限环的排除方法 ===
'''极限环排除方法的核心思想是：'''若能找到一个在系统演化过程中单调变化的量，则系统不可能存在周期轨道。以下三种方法各自利用了不同的结构特征，其中梯度系统与李雅普诺夫函数方法适用于任意 <math>n</math> 维系统，而 Dulac 准则则为二维专属方法。
'''极限环排除方法的核心思想是：'''若能找到一个在系统演化过程中单调变化的量，则系统不可能存在周期轨道。以下三种方法各自利用了不同的结构特征，其中梯度系统与李雅普诺夫函数方法在高维中仍然成立，而 Dulac 准则则为二维专属方法。
====梯度系统====
====梯度系统====
在某一区域内处处保持严格同号（恒正或恒负），则 <math>V</math> 在该区域内沿任意轨道单调变化。由同样的理由，该区域内不可能存在周期轨道。这样的函数称为该区域上的'''李雅普诺夫函数'''（Lyapunov function）。<ref name="Perko" />
在某一区域内处处保持严格同号（恒正或恒负），则 <math>V</math> 在该区域内沿任意轨道单调变化。由同样的理由，该区域内不可能存在周期轨道。这样的函数称为该区域上的'''李雅普诺夫函数'''（Lyapunov function）。<ref name="Perko" />
李雅普诺夫函数方法适用于一般自治系统，比梯度系统更为灵活。然而，该方法的主要局限在于需要预先猜测或构造出合适的 <math>V</math>，而这在一般情形下并无系统性方法。
==== Dulac 准则 ====
==== Dulac 准则 ====
上式本质上是加权向量场 <math>(BP, BQ)</math> 的散度。若散度处处同号，则向量场在该区域内整体表现为纯收缩或纯扩张，轨道无法首尾相接构成闭合回路。当 <math>B \equiv 1</math> 时，该准则退化为经典的 '''Bendixson 判据'''。
上式本质上是加权向量场 <math>(BP, BQ)</math> 的散度。若散度处处同号，则向量场在该区域内整体表现为纯收缩或纯扩张，轨道无法首尾相接构成闭合回路。当 <math>B \equiv 1</math> 时，该准则退化为经典的 '''Bendixson 判据'''。
综上，三种方法在适用范围上各有侧重：梯度系统要求系统具有势函数结构，李雅普诺夫函数则更为灵活但缺乏系统性的构造方法，而 Dulac 准则通过引入权函数 <math>B</math> 在二维框架内进一步扩大了适用范围，但无法推广至高维。
=== 极限环存在性理论 ===
=== 极限环存在性理论 ===
* '''Dulac 准则不可推广'''：Dulac 准则依赖于平面向量场的散度，其核心论证在三维及以上系统中没有对应的推广形式。<ref name="Strogatz" />
* '''Dulac 准则不可推广'''：Dulac 准则依赖于平面向量场的散度，其核心论证在三维及以上系统中没有对应的推广形式。<ref name="Strogatz" />
* '''存在性判定困难'''：在高维系统中，目前尚无通用的极限环存在性判定方法。实际分析通常依赖霍普夫分岔理论（见[[#分岔与极限环]]）或数值方法（见[[#近似与数值方法]]）。
* '''存在性判定困难'''：在高维系统中，目前尚无通用的极限环存在性判定方法。实际分析通常依赖霍普夫分岔理论（见[[#分岔与极限环的产生和消亡]]）或数值方法（见[[#极限环的近似定量分析]]）。
=== 高维系统中的通用分析工具 ===
=== 高维系统中的通用分析工具 ===
尽管二维的专属工具在高维失效，部分方法仍可直接推广至高维系统，分别用于排除和确认极限环的存在。
在'''排除'''极限环存在方面，'''梯度系统'''与'''李雅普诺夫函数'''方法（见[[#极限环的排除方法]]）的核心逻辑不依赖于平面拓扑，在任意维度下同样适用：若能找到一个沿轨道严格单调变化的函数，则系统不可能存在周期轨道。<ref name="Strogatz" /><ref name="Perko" />
在'''排除'''极限环存在方面，'''梯度系统'''与'''李雅普诺夫函数'''方法（见[[#极限环的排除方法]]）的核心逻辑不依赖于平面拓扑，在任意维度下同样适用：若能找到一个沿轨道严格单调变化的函数，则系统不可能存在周期轨道。<ref name="Strogatz" /><ref name="Perko" />
在'''确认'''极限环存在方面，'''霍普夫分岔理论'''（见[[#分岔与极限环]]）同样适用于高维系统：当平衡点的特征值跨越虚轴时，可以严格证明附近产生极限环，这一结论在任意维度下均成立。<ref name="Kuznetsov" />
在'''确认'''极限环存在方面，'''霍普夫分岔理论'''（见[[#分岔与极限环的产生和消亡]]）同样适用于高维系统：当平衡点的特征值跨越虚轴时，可以严格证明附近产生极限环，这一结论在任意维度下均成立。<ref name="Kuznetsov" />
=== 高维系统中的稳定性分析：Floquet 理论 ===
=== 高维系统中的稳定性分析：Floquet 理论 ===
Floquet 理论通过分析极限环附近的线性化系统，给出了任意维度下极限环稳定性的严格定量判据。在任意维度的系统中，极限环的稳定性可通过'''Floquet 理论'''（Floquet theory）进行严格分析。设系统存在一条周期为 <math>T</math> 的极限环 <math>\Gamma</math>，将轨道写为 <math>\mathbf{x}(t) = \mathbf{x}^*(t) + \boldsymbol{\xi}(t)</math>，其中 <math>\mathbf{x}^*(t)</math> 为极限环上的周期解，<math>\boldsymbol{\xi}(t)</math> 为小扰动。对系统在 <math>\Gamma</math> 附近线性化，得到线性变分方程
Floquet 理论通过分析极限环附近的线性化系统，给出了任意维度下极限环稳定性的严格定量判据。
设系统存在一条周期为 <math>T</math> 的极限环 <math>\Gamma</math>。对系统在 <math>\Gamma</math> 附近线性化，得到一个以 <math>T</math> 为周期的线性变分方程。该方程在一个完整周期内的基本解矩阵称为'''单值矩阵'''（monodromy matrix），其特征值称为'''Floquet 乘子'''（Floquet multipliers），又称'''特征乘子'''（characteristic multipliers）。<ref name="Chicone">Chicone, C. (2006). ''Ordinary Differential Equations with Applications'' (2nd ed.). Springer.</ref>
对于 一般自治系统，极限环对应 <math>n-1</math> 个决定横截稳定性的 Floquet 乘子（沿极限环切方向本身对应一个模恒为 1 的乘子，不计入横截稳定性分析）。稳定性由这 <math>n-1</math> 个乘子的模决定：<ref name="Kuznetsov" /><ref name="Wiggins" />
* 若所有特征乘子的模均小于 1，则极限环'''渐近稳定'''；
* 若所有乘子的模均小于 1，则极限环'''渐近稳定'''；
* 若存在模大于 1 的特征乘子，则极限环'''不稳定'''；
* 若存在模大于 1 的乘子，则极限环'''不稳定'''；
* 若存在模等于 1 的特征乘子（且其余均小于 1），则需进一步分析，对应极限环稳定性的临界情形。
* 若存在模等于 1 的乘子（且其余均小于 1），则对应稳定性的临界情形，需进一步分析。
在二维情形下，<math>n-1=1</math>，退化为单个特征乘子，与[[#稳定性分类]]一节中的讨论一致。<ref name="Wiggins" />
=== 极限环周围的相空间结构：不变流形 ===
=== 极限环周围的相空间结构：不变流形 ===
不变流形理论描述了极限环周围相空间的几何结构，揭示了高维系统中特有的鞍型周期轨道现象。
不变流形理论描述了极限环周围相空间的几何结构，揭示了高维系统中特有的鞍型周期轨道现象。
Floquet 理论给出了极限环稳定性的定量判据，而'''不变流形理论'''（invariant manifold theory）则进一步描述了极限环周围相空间的几何结构。在 <math>n</math> 维系统中，极限环 <math>\Gamma</math> 周围的相空间可按特征乘子的模分解为三类不变流形：<ref name="Wiggins" /><ref name="GH">Guckenheimer, J., & Holmes, P. (1983). ''Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields''. Springer.</ref>
Floquet 理论给出了极限环稳定性的定量判据，而'''不变流形理论'''（invariant manifold theory）则进一步描述了极限环周围相空间的几何结构。在 一般自治系统中，极限环 <math>\Gamma</math> 周围的相空间可按特征乘子的模分解为三类不变流形：<ref name="Wiggins" /><ref name="GH">Guckenheimer, J., & Holmes, P. (1983). ''Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields''. Springer.</ref>
* '''稳定流形'''（stable manifold）<math>W^s(\Gamma)</math>：由在 <math>t \to +\infty</math> 时趋近于 <math>\Gamma</math> 的轨道构成，对应模小于 1 的特征乘子；
* '''稳定流形'''（stable manifold）<math>W^s(\Gamma)</math>：由在 <math>t \to +\infty</math> 时趋近于 <math>\Gamma</math> 的轨道构成，对应模小于 1 的特征乘子；
在三维及以上系统中，极限环失稳后可触发倍周期分岔、环面振荡乃至混沌等一系列复杂动力学行为，这是二维系统中不可能出现的现象。
在三维及以上系统中，极限环失稳后可触发倍周期分岔、环面振荡乃至混沌等一系列复杂动力学行为，这是二维系统中不可能出现的现象。
在二维系统中，极限环失稳后的行为相对简单：稳定极限环只能通过分岔转变为不稳定极限环或消亡。而在三维及以上系统中，极限环失稳可触发一系列更为复杂的动力学行为。<ref name="Strogatz" /> 各类分岔机制将在[[#分岔与极限环]]一节中详细介绍，此处重点介绍两种高维系统所特有的失稳路径：
在二维系统中，极限环失稳后的行为相对简单：稳定极限环只能通过分岔转变为不稳定极限环或消亡。而在三维及以上系统中，极限环失稳可触发一系列更为复杂的动力学行为。<ref name="Strogatz" /> 各类分岔机制将在[[#分岔与极限环的产生和消亡]]一节中详细介绍，此处重点介绍两种高维系统所特有的失稳路径：
* '''倍周期分岔'''（period-doubling bifurcation）：极限环失稳时，系统产生一个周期为原来两倍的新极限环。随着参数持续变化，倍周期分岔可以反复发生，形成'''倍周期级联'''（period-doubling cascade），最终在有限的参数范围内导致混沌。这一路径由 Feigenbaum 于1978年系统研究，具有普适的定量规律。<ref name="Feigenbaum">Feigenbaum, M. J. (1978). Quantitative universality for a class of nonlinear transformations. ''Journal of Statistical Physics'', 19(1), 25–52.</ref>
* '''倍周期分岔'''（period-doubling bifurcation）：极限环失稳时，系统产生一个周期为原来两倍的新极限环。随着参数持续变化，倍周期分岔可以反复发生，形成'''倍周期级联'''（period-doubling cascade），最终在有限的参数范围内导致混沌。这一路径由 Feigenbaum 于1978年系统研究，具有普适的定量规律。<ref name="Feigenbaum">Feigenbaum, M. J. (1978). Quantitative universality for a class of nonlinear transformations. ''Journal of Statistical Physics'', 19(1), 25–52.</ref>
==分岔与极限环的产生和消亡==
==分岔与极限环的产生和消亡==
在动力系统中，极限环的产生、消亡与结构变化通常与系统参数的变化密切相关。当参数跨越某一临界值时，系统的定性行为发生突变，这一现象称为'''分岔'''（bifurcation）。上述分岔可分为两类：'''局部分岔'''（霍普夫分岔、鞍结分岔、倍周期分岔、环面分岔等）仅涉及极限环邻域内的局部行为；'''全局分岔'''（同宿分岔、异宿分岔等）则涉及相空间的全局结构。本节仅就分岔与极限环的直接关联作简要介绍，以下各类分岔与极限环关系的核心结论均在满足相应非退化条件下成立，各定理的完整条件与证明详见[[分岔理论]]词条及各类分岔的对应词条。除倍周期分岔仅适用于三维及以上连续自治系统外，以下各类分岔均适用于任意 <math>n</math> 维系统。
在动力系统中，极限环的产生、消亡与结构变化通常与系统参数的变化密切相关。当参数跨越某一临界值时，系统的定性行为发生突变，这一现象称为'''分岔'''（bifurcation）。上述分岔可分为两类：'''局部分岔'''（霍普夫分岔、鞍结分岔、倍周期分岔、环面分岔等）仅涉及极限环邻域内的局部行为；'''全局分岔'''（同宿分岔、异宿分岔等）则涉及相空间的全局结构。本节仅就分岔与极限环的直接关联作简要介绍，以下各类分岔与极限环关系的核心结论均在满足相应非退化条件下成立，各定理的完整条件与证明详见[[分岔理论]]词条及各类分岔的对应词条。除倍周期分岔仅适用于三维及以上连续自治系统外，以下各类分岔均适用于任意 <math>n</math> 维自治系统。
与极限环直接相关的主要分岔类型包括：
与极限环直接相关的主要分岔类型包括：