极限环 - 版本历史

褚天舒：/* 范德珀尔振子 */

2026-08-31T12:51:52Z

范德珀尔振子

褚天舒：/* 范德珀尔振子 */

2026-08-11T14:56:32Z

范德珀尔振子

夏沛雯：/* 不变性 */

2026-06-21T19:15:07Z

不变性

2026年6月21日 (日) 19:14 夏沛雯

2026-06-21T19:14:05Z

显示更改

Jake：/* 应用 */

2026-06-17T05:58:19Z

应用

Jake：/* 待解决的问题 */

2026-06-17T05:57:23Z

待解决的问题

Jake：/* 近似与数值方法 */

2026-06-17T05:55:05Z

近似与数值方法

Jake：/* 分岔与极限环的产生和消亡 */

2026-06-17T05:53:17Z

分岔与极限环的产生和消亡

Jake：/* 极限环失稳后的动力学行为 */

2026-06-17T05:48:39Z

极限环失稳后的动力学行为

Jake：/* 极限环周围的相空间结构：不变流形 */

2026-06-17T05:46:18Z

极限环周围的相空间结构：不变流形

@@ 第84行： / 第84行： @@
 [[Image:VanDerPolPhaseSpace.png|none|'''范德波尔振荡器 Van der Pol oscillator''' 的稳定极限环|thumb|300x300像素]]
-[[范德珀尔振子]]是研究极限环的经典[[非线性系统]]，由荷兰物理学家[[Balthasar van der Pol]]（通常译为范德波尔）于1920年在研究真空电子管振荡电路时提出<ref name="vanderPol">van der Pol, B. (1926). On relaxation-oscillations. ''The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science'', 2(11), 978–992.</ref>：
+[[范德珀尔振子]]是研究极限环的经典[[非线性系统]]，由荷兰物理学家[[Balthasar van der Pol]]（通常译为范德珀尔）于1920年在研究真空电子管振荡电路时提出<ref name="vanderPol">van der Pol, B. (1926). On relaxation-oscillations. ''The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science'', 2(11), 978–992.</ref>：
 :<math display="block">\ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu > 0</math>

@@ 第84行： / 第84行： @@
 [[Image:VanDerPolPhaseSpace.png|none|'''范德波尔振荡器 Van der Pol oscillator''' 的稳定极限环|thumb|300x300像素]]
-[[范德珀尔振子]]是研究极限环的经典[[非线性系统]]，由荷兰物理学家[[Balthasar van der Pol]]在研究电子管振荡电路时提出<ref name="vanderPol">van der Pol, B. (1926). On relaxation-oscillations. ''The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science'', 2(11), 978–992.</ref>：
+[[范德珀尔振子]]是研究极限环的经典[[非线性系统]]，由荷兰物理学家[[Balthasar van der Pol]]（通常译为范德波尔）于1920年在研究真空电子管振荡电路时提出<ref name="vanderPol">van der Pol, B. (1926). On relaxation-oscillations. ''The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science'', 2(11), 978–992.</ref>：
 :<math display="block">\ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, \quad \mu > 0</math>
-原点处线性化矩阵的特征值实部为正，故原点是[[不稳定焦点]]。另一方面，通过估计[[能量函数]]可以证明轨道被限制在有界环形区域内。由[[庞加莱–本迪克松定理]]知该区域内存在极限环，结合唯一性论证可知恰好存在一条，所有非平衡轨道最终趋向该极限环。参数 <math>\mu</math> 较小时振荡近似正弦波；较大时出现[[松弛振荡]]（relaxation oscillation）<ref name="Strogatz" />。
+原点处线性化矩阵的特征值实部为正，故原点是[[不稳定焦点]]。另一方面，通过估计[[能量函数]]可以证明轨道被限制在有界环形区域内。由[[庞加莱–本迪克松定理]]知该区域内存在极限环，结合唯一性论证可知恰好存在一条，所有非平衡轨道最终趋向该极限环。参数 <math>\mu</math>为0时，方程退化为简谐振荡器，系统做等幅正弦振荡，没有能量耗散或注入，能量守恒； 较小时（当[math]\displaystyle{ 0<\mu<<1 }[/math]时）振荡近似正弦波；较大时（当[math]\displaystyle{ \mu>>1 }[/math]时）出现[[松弛振荡]]（relaxation oscillation）<ref name="Strogatz" />。
 === 极限环的排除方法 ===

@@ 第34行： / 第34行： @@
 === 不变性 ===
-极限环是相空间中的'''[[不变集]]'''：若系统的初始条件恰好落在该闭轨道上，则对应解将永远沿该轨道周期性演化。<ref name="Strogatz"/>
+极限环是相空间中的'''[[不变集]]'''：若系统的初始条件恰好落在该闭轨道上，则对应解将永远沿该轨道周期性演化<ref name="Strogatz"/>。
 这一性质来源于微分方程解的'''[[存在唯一性定理]]'''（[[Picard定理]]）：[[相空间]]中过同一点只能有一条轨道。因此，若某轨道与极限环有一点相交，则由解的唯一性，该轨道必与极限环重合；反之，极限环外的轨道也不能穿入极限环<ref name="LibreTexts">Lebl, J. ''Notes on Diffy Qs: Differential Equations for Engineers'', §8.4: Limit cycles. Available at: https://www.jirka.org/diffyqs/ (accessed June 2026).</ref>。

@@ 第218行： / 第218行： @@
 科学应用中许多自持振荡系统的仿真中，极限环具有重要意义。典型例子包括：
-* 空气动力学中的极限环振荡<ref name="Thomas">Thomas, J. P., Dowell, E. H., & Hall, K. C. (2002). Nonlinear inviscid aerodynamic effects on transonic divergence, flutter, and limit-cycle oscillations. ''AIAA Journal'', 40(4), 638–646.</ref>
+* [[空气动力学]]中的极限环振荡<ref name="Thomas">Thomas, J. P., Dowell, E. H., & Hall, K. C. (2002). Nonlinear inviscid aerodynamic effects on transonic divergence, flutter, and limit-cycle oscillations. ''AIAA Journal'', 40(4), 638–646.</ref>
 * 神经元动作电位的[[霍奇金–赫胥黎模型]]（Hodgkin–Huxley model）
 * 糖酵解过程中的[[塞尔科夫模型]]（Sel'kov model）<ref name="Selkov">Sel'kov, E. E. (1968). Self-oscillations in glycolysis 1. A simple kinetic model. ''European Journal of Biochemistry'', 4(1), 79–86.</ref>
 * 动物基因表达、激素水平和体温的周期性变化，属于[[昼夜节律]]（circadian rhythm）<ref name="Leloup">Leloup, J.-C., Gonze, D., & Goldbeter, A. (1999). Limit cycle models for circadian rhythms based on transcriptional regulation in ''Drosophila'' and ''Neurospora''. ''Journal of Biological Rhythms'', 14(6), 433–448.</ref>
 * 癌细胞在局限微环境中的迁移<ref name="Bruckner">Brückner, D. B., Fink, A., Schreiber, C., Röttgermann, P. J. F., Rädler, J. O., & Broedersz, C. P. (2019). Stochastic nonlinear dynamics of confined cell migration in two-state systems. ''Nature Physics'', 15(6), 595–601.</ref>
-* 非线性电路系统，包括[[范德波尔模型]]（Van der Pol model）<ref name="Ginoux">Ginoux, J.-M., & Letellier, C. (2012). Van der Pol and the history of relaxation oscillations: Toward the emergence of a concept. ''Chaos'', 22(2), 023120.</ref>
+* [[非线性电路系统]]，包括[[范德波尔模型]]（Van der Pol model）<ref name="Ginoux">Ginoux, J.-M., & Letellier, C. (2012). Van der Pol and the history of relaxation oscillations: Toward the emergence of a concept. ''Chaos'', 22(2), 023120.</ref>
 == 参考文献 ==
 <references/>

@@ 第206行： / 第206行： @@
 极限环的一般理论至今仍存在若干重要的开放问题，尤其在高维系统中，现有理论框架远未完善。
-* '''[[希尔伯特第十六问题]]'''：1900年，大卫·希尔伯特（David Hilbert）提出的第十六问题第二部分询问：次数不超过 <math>n</math> 的平面多项式系统最多可能存在多少条极限环。<ref name="Hilbert">Hilbert, D. (1902). Mathematical problems. ''Bulletin of the American Mathematical Society'', 8(10), 437–479.</ref> 即使对于二次多项式情形（<math>n=2</math>）该问题至今仍未完全解决；对于一般次数 <math>n</math>，极限环个数的有限性也尚未得到证明。<ref name="Perko" />
+* '''[[希尔伯特第十六问题]]'''：1900年，[[大卫·希尔伯特]]（David Hilbert）提出的第十六问题第二部分询问：次数不超过 <math>n</math> 的平面多项式系统最多可能存在多少条极限环。<ref name="Hilbert">Hilbert, D. (1902). Mathematical problems. ''Bulletin of the American Mathematical Society'', 8(10), 437–479.</ref> 即使对于二次多项式情形（<math>n=2</math>）该问题至今仍未完全解决；对于一般次数 <math>n</math>，极限环个数的有限性也尚未得到证明。<ref name="Perko" />
-* '''高维系统中的存在性理论'''：目前高维系统中确认极限环存在的方法（霍普夫分岔、数值延拓等）均为局部或逐案分析，缺乏系统性的全局存在性理论。建立适用于高维系统的一般性存在性判定框架，是动力系统领域的重要开放问题。<ref name="GH" />
+* '''高维系统中的存在性理论'''：目前高维系统中确认极限环存在的方法（[[霍普夫分岔]]、[[数值延拓]]等）均为局部或逐案分析，缺乏系统性的全局存在性理论。建立适用于高维系统的一般性存在性判定框架，是[[动力系统]]领域的重要开放问题。<ref name="GH" />
-* '''高维系统中的排除方法'''：Dulac 准则在高维系统中无法推广，而李雅普诺夫函数方法虽适用于任意维度，但其构造缺乏系统性方法，实际应用中往往需要针对具体系统逐案设计。<ref name="Strogatz" /><ref name="GH" />
+* '''高维系统中的排除方法'''：[[Dulac 准则]]在高维系统中无法推广，而[[李雅普诺夫函数]]方法虽适用于任意维度，但其构造缺乏系统性方法，实际应用中往往需要针对具体系统逐案设计。<ref name="Strogatz" /><ref name="GH" />
 * '''极限环的个数与分布'''：在高维系统中，极限环的个数问题远比二维复杂，系统可以同时存在多个极限环，且随参数变化可能经历复杂的分岔序列。对高维系统中极限环个数与分布的系统性描述，至今仍是开放问题。<ref name="Kuznetsov" />

@@ 第196行： / 第196行： @@
 极限环的精确解析解在一般情形下难以获得，但对于某些特殊结构的系统可通过近似方法得到定量信息。常用方法包括：
-* '''[[松弛振荡]]'''（relaxation oscillation）：适用于具有快慢时间尺度分离特征的系统（典型情形为 Van der Pol 振子大参数极限 <math>\mu \gg 1</math>）。轨道在相平面上缓慢漂移与快速跃迁交替，极限环周期可估算为 <math>T \approx \mu(3 - 2\ln 2)</math>。<ref name="Strogatz" />
+* '''[[松弛振荡]]'''（relaxation oscillation）：适用于具有快慢时间尺度分离特征的系统（典型情形为[[Van der Pol 振子]]大参数极限 <math>\mu \gg 1</math>）。轨道在相平面上缓慢漂移与快速跃迁交替，极限环周期可估算为 <math>T \approx \mu(3 - 2\ln 2)</math>。<ref name="Strogatz" />
-* '''[[平均法]]'''（averaging method）：适用于弱非线性振荡（<math>\varepsilon \ll 1</math>）。通过对快速振荡取平均，将原系统化为振幅慢变方程，其不动点对应极限环振幅，稳定性由线性化特征值判定。对于 Van der Pol 振子，一阶近似给出振幅 <math>a = 2</math>。<ref name="Perko" /><ref name="GH" />
+* '''[[平均法]]'''（averaging method）：适用于弱非线性振荡（<math>\varepsilon \ll 1</math>）。通过对快速振荡取平均，将原系统化为振幅慢变方程，其不动点对应极限环振幅，稳定性由线性化特征值判定。对于 [[Van der Pol 振子]]，一阶近似给出振幅 <math>a = 2</math>。<ref name="Perko" /><ref name="GH" />
 * '''[[数值延拓方法]]'''（numerical continuation）：从已知极限环出发，沿参数方向系统性追踪极限环分支，同步检测分岔点。可处理任意强非线性系统，不依赖小参数假设。常用软件包括 AUTO、MATCONT 与 PyDSTool。<ref name="Kuznetsov" /><ref name="GH" />

@@ 第176行： / 第176行： @@
 == 分岔与极限环的产生和消亡 ==
-在动力系统中，极限环的产生、消亡与结构变化通常与系统参数的变化密切相关。当参数跨越某一临界值时，系统的定性行为发生突变，这一现象称为'''[[分岔]]'''（bifurcation）。分岔可分为两类：'''[[局部分岔]]'''（[[霍普夫分岔]]、[[鞍结分岔]]、[[倍周期分岔]]、[[环面分岔]]等）仅涉及极限环邻域内的局部行为；'''[[全局分岔]]'''（[[同宿分岔]]、[[异宿分岔]]等）则涉及相空间的全局结构。本节仅就分岔与极限环的直接关联作简要介绍，以下各类分岔的核心结论均在满足相应非退化条件下成立，各定理的完整条件与证明详见[[分岔理论]]词条及各类分岔的对应词条。除倍周期分岔仅适用于三维及以上连续自治系统外，以下各类分岔均适用于任意 <math>n</math> 维系统。
+在[[动力系统]]中，极限环的产生、消亡与结构变化通常与系统参数的变化密切相关。当参数跨越某一临界值时，系统的定性行为发生突变，这一现象称为'''[[分岔]]'''（bifurcation）。分岔可分为两类：'''[[局部分岔]]'''（[[霍普夫分岔]]、[[鞍结分岔]]、[[倍周期分岔]]、[[环面分岔]]等）仅涉及极限环邻域内的局部行为；'''[[全局分岔]]'''（[[同宿分岔]]、[[异宿分岔]]等）则涉及相空间的全局结构。本节仅就分岔与极限环的直接关联作简要介绍，以下各类分岔的核心结论均在满足相应非退化条件下成立，各定理的完整条件与证明详见[[分岔理论]]词条及各类分岔的对应词条。除[[倍周期分岔]]仅适用于三维及以上连续自治系统外，以下各类分岔均适用于任意 <math>n</math> 维系统。
-[[文件:超临界霍普夫分岔示意.png|thumb|none|400px|超临界霍普夫分岔示意图。横轴为参数 <math>\mu</math>，纵轴为振幅 <math>A</math>。当 <math>\mu < 0</math> 时平衡点稳定；当 <math>\mu > 0</math> 时平衡点失稳，系统产生稳定极限环。]]
+[[文件:超临界霍普夫分岔示意.png|thumb|none|400px|[[超临界霍普夫分岔]]示意图。横轴为参数 <math>\mu</math>，纵轴为振幅 <math>A</math>。当 <math>\mu < 0</math> 时平衡点稳定；当 <math>\mu > 0</math> 时平衡点失稳，系统产生稳定极限环。]]
 与极限环直接相关的主要分岔类型包括：
-* '''[[霍普夫分岔]]'''（Hopf bifurcation）：当平衡点的一对复特征值实部跨越零且满足横截条件与非退化条件时，平衡点失稳并在附近产生极限环。超临界情形产生稳定极限环，亚临界情形产生不稳定极限环，后者通常与一个已存在的稳定极限环共存。霍普夫分岔是高维系统中为数不多的可以严格证明极限环存在的情形之一。<ref name="Kuznetsov" /><ref name="GH" />
+* '''[[霍普夫分岔]]'''（Hopf bifurcation）：当平衡点的一对复特征值实部跨越零且满足横截条件与非退化条件时，平衡点失稳并在附近产生极限环。[[超临界]]情形产生稳定极限环，[[亚临界]]情形产生不稳定极限环，后者通常与一个已存在的稳定极限环共存。霍普夫分岔是高维系统中为数不多的可以严格证明极限环存在的情形之一。<ref name="Kuznetsov" /><ref name="GH" />
-* '''[[极限环鞍结分岔]]'''（saddle-node bifurcation of periodic orbits）：一条稳定极限环与一条不稳定极限环在参数变化下逐渐靠近，在临界参数值处合并为半稳定极限环后同时消亡。当系统存在双稳态时，这一机制可导致'''[[迟滞现象]]'''（hysteresis）。<ref name="Kuznetsov" /><ref name="Strogatz" />
+* '''[[极限环鞍结分岔]]'''（saddle-node bifurcation of periodic orbits）：一条稳定极限环与一条不稳定极限环在参数变化下逐渐靠近，在临界参数值处合并为半稳定极限环后同时消亡。当系统存在[[双稳态]]时，这一机制可导致'''[[迟滞现象]]'''（hysteresis）。<ref name="Kuznetsov" /><ref name="Strogatz" />
-* '''[[同宿与异宿分岔]]'''（homoclinic／heteroclinic bifurcation）：极限环趋近于平衡点的连接轨道时，其周期趋向无穷大并最终消亡。在三维系统中，若同宿轨道趋近于[[鞍焦点]]且满足 [[Shilnikov 条件]]（鞍焦点的不稳定特征值实部大于稳定特征值实部的绝对值），系统附近可出现可数无穷条周期轨道及混沌动力学，称为 '''[[Shilnikov混沌]]'''。<ref name="Shilnikov">Shilnikov, L. P. (1965). A case of the existence of a denumerable set of periodic motions. ''Soviet Mathematics Doklady'', 6, 163–166.</ref><ref name="GH" />
+* '''[[同宿与异宿分岔]]'''（homoclinic／heteroclinic bifurcation）：极限环趋近于平衡点的连接轨道时，其周期趋向无穷大并最终消亡。在三维系统中，若[[同宿轨道]]趋近于[[鞍焦点]]且满足 [[Shilnikov 条件]]（鞍焦点的不稳定特征值实部大于稳定特征值实部的绝对值），系统附近可出现可数无穷条周期轨道及混沌动力学，称为 '''[[Shilnikov混沌]]'''。<ref name="Shilnikov">Shilnikov, L. P. (1965). A case of the existence of a denumerable set of periodic motions. ''Soviet Mathematics Doklady'', 6, 163–166.</ref><ref name="GH" />
 * '''[[倍周期分岔]]'''（period-doubling bifurcation）〔仅适用于三维及以上连续自治系统〕：极限环失稳时产生周期加倍的新极限环，反复发生形成'''[[倍周期级联]]'''（period-doubling cascade），最终导致混沌。[[Feigenbaum]]于1978年证明这一过程具有普适的定量规律，收敛速率由 '''[[Feigenbaum 常数]]''' <math>\delta \approx 4.669</math> 描述。<ref name="Feigenbaum">Feigenbaum, M. J. (1978). Quantitative universality for a class of nonlinear transformations. ''Journal of Statistical Physics'', 19(1), 25–52.</ref><ref name="Kuznetsov" />

@@ 第170行： / 第170行： @@
 === 极限环失稳后的动力学行为 ===
-'''在三维及以上系统中，极限环失稳后可触发倍周期分岔、环面振荡乃至混沌等一系列复杂动力学行为，这是二维系统中不可能出现的现象。'''
+'''在三维及以上系统中，极限环失稳后可触发[[倍周期分岔]]、[[环面振荡]]乃至[[混沌]]等一系列复杂动力学行为，这是二维系统中不可能出现的现象。'''
-在二维系统中，极限环失稳后只能通过分岔转变为不稳定极限环或消亡。而在三维及以上系统中，极限环失稳可触发更为复杂的行为，常见路径包括[[倍周期分岔]]（可通过倍周期级联进入混沌）以及[[环面分岔]]（产生准周期运动，环面进一步破裂后可进入混沌，即[[Ruelle–Takens情景]]）。这些现象共同揭示了高维系统动力学的根本复杂性：极限环不再是系统长期行为的终点，而可能只是通往更复杂动力学的一个过渡阶段。无论是二维还是高维系统，极限环的产生与消亡均与系统参数的变化密切相关，这一现象将在下一节中系统介绍。<ref name="Strogatz" /><ref name="Kuznetsov" />
+在二维系统中，极限环失稳后只能通过[[分岔转变]]为不稳定极限环或消亡。而在三维及以上系统中，[[极限环失稳]]可触发更为复杂的行为，常见路径包括[[倍周期分岔]]（可通过[[倍周期级联]]进入[[混沌]]）以及[[环面分岔]]（产生[[准周期运动]]，环面进一步破裂后可进入[[混沌]]，即[[Ruelle–Takens情景]]）。这些现象共同揭示了[[高维系统动力学]]的根本复杂性：极限环不再是系统长期行为的终点，而可能只是通往更复杂动力学的一个过渡阶段。无论是二维还是高维系统，极限环的产生与消亡均与系统参数的变化密切相关，这一现象将在下一节中系统介绍。<ref name="Strogatz" /><ref name="Kuznetsov" />
 == 分岔与极限环的产生和消亡 ==

@@ 第166行： / 第166行： @@
 这一分解揭示了高维系统中极限环的一个二维所没有的现象：'''[[鞍型周期轨道]]'''（saddle-type periodic orbit）——极限环在某些法向方向上吸引轨道，同时在另一些方向上排斥轨道。这类极限环在实际系统中可充当不同吸引域之间的边界。<ref name="Kuznetsov" /><ref name="GH" />
-[[文件:鞍型周期轨道不变流形图（三维示例）.png|thumb|none|350px|[[鞍型周期轨道]][[不变流形]]图（三维示例）|鞍型周期轨道 Γ（蓝色）的稳定流形 <math>W^s(\Gamma)</math>（绿色）上的轨迹在 <math>t\to+\infty</math> 时趋近 Γ，而不稳定流形 <math>W^u(\Gamma)</math>（红色）上的轨迹在 <math>t\to+\infty</math> 时远离 Γ，两者的 Floquet 乘子满足 <math>|\lambda_s|<1</math>，<math>|\lambda_u|>1</math>。]]
+[[文件:鞍型周期轨道不变流形图（三维示例）.png|thumb|none|350px|[[鞍型周期轨道]][[不变流形]]图（三维示例）|鞍型周期轨道[math]\Gamma[/math]（蓝色）的稳定流形 <math>W^s(\Gamma)</math>（绿色）上的轨迹在 <math>t\to+\infty</math> 时趋近[math]\Gamma[/math]，而不稳定流形 <math>W^u(\Gamma)</math>（红色）上的轨迹在 <math>t\to+\infty</math> 时远离[math]\Gamma[/math]，两者的 [[Floquet 乘子]]满足 <math>|\lambda_s|<1</math>，<math>|\lambda_u|>1</math>。]]
 === 极限环失稳后的动力学行为 ===