“极限环”的版本间的差异
(→性质) |
(→典型实例) |
||
| 第64行: | 第64行: | ||
因此,极限环描述的不是某一特殊初始条件下的解,而是系统在大范围初始条件下共同趋向的'''渐近行为'''。在生物节律、电子振荡器、化学振荡反应等领域,极限环正是系统自发产生稳定周期行为的数学根源。<ref name="Wikipedia_LC" /> | 因此,极限环描述的不是某一特殊初始条件下的解,而是系统在大范围初始条件下共同趋向的'''渐近行为'''。在生物节律、电子振荡器、化学振荡反应等领域,极限环正是系统自发产生稳定周期行为的数学根源。<ref name="Wikipedia_LC" /> | ||
| − | == 典型实例 | + | == 典型实例 == |
| − | |||
| − | + | === 简单构造的极限环 === | |
| + | 极限环可以通过显式构造二维系统的径向与角向动力学来得到。考虑如下系统: | ||
<math display="block"> | <math display="block"> | ||
\begin{cases} | \begin{cases} | ||
| 第76行: | 第76行: | ||
</math> | </math> | ||
| − | 其中 <math>r^2 = x^2 + y^2</math> | + | 其中 <math>r^2 = x^2 + y^2</math>。将系统改写为极坐标形式,得: |
| − | |||
| − | |||
<math display="block"> | <math display="block"> | ||
| 第86行: | 第84行: | ||
\end{cases} | \end{cases} | ||
</math> | </math> | ||
| + | 两个方程可独立分析:角变量 <math>\theta</math> 以常数速度单调增加,轨道因此始终绕原点旋转;径向变量 <math>r</math> 满足一维方程 <math>\dot{r} = r(1 - r^2)</math>,在 <math>r = 0</math> 与 <math>r = 1</math> 处存在平衡点。对 <math>r = 1</math> 的稳定性分析表明:当 <math>0 < r < 1</math> 时 <math>\dot{r} > 0</math>,轨道向外增长;当 <math>r > 1</math> 时 <math>\dot{r} < 0</math>,轨道向内收缩。因此,只要 <math>r \neq 0</math>,径向分量都将单调收敛至 <math>r = 1</math>。 | ||
| − | + | 综合两个方向的行为:轨道在趋向单位圆的过程中持续旋转,在相平面中呈螺旋形,并最终收敛到 <math>r = 1</math> 上的周期运动。该单位圆即为系统唯一的稳定极限环,除原点外的所有轨道均以其为 ω 极限集。<ref name="Strogatz" /> | |
| − | |||
| − | |||
| − | |||
| − | |||
| − | + | === Van der Pol 振子 === | |
| − | ===Van der Pol 振子=== | ||
| − | |||
| − | Van der Pol | + | [[Image:VanDerPolPhaseSpace.png|right|250px|thumb|'''<font color="#ff8000">范德波尔振荡器 Van der Pol oscillator</font>''' 的稳定极限环]] |
| + | Van der Pol 振子是研究极限环的经典非线性系统之一。该方程最初由荷兰物理学家巴尔塔萨·范德波尔(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>其方程为: | ||
<math display="block"> | <math display="block"> | ||
\ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, | \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, | ||
</math> | </math> | ||
| − | + | 其中 <math>\mu > 0</math> 为参数,控制非线性阻尼的强度。将其改写为二维一阶系统: | |
| − | 其中 <math>\mu > 0</math> | ||
| − | |||
| − | |||
| − | |||
<math display="block"> | <math display="block"> | ||
\begin{cases} | \begin{cases} | ||
| 第113行: | 第103行: | ||
\end{cases} | \end{cases} | ||
</math> | </math> | ||
| + | 该系统在原点处存在唯一平衡点。对原点处的线性化矩阵进行分析,其特征值实部为正,因此原点是不稳定焦点,附近轨道向外螺旋发散。<ref name="Strogatz" /> | ||
| − | + | 另一方面,通过估计能量函数可以证明,系统轨道被限制在相平面的某个有界环形区域内,且该区域内不含平衡点。由'''庞加莱-本迪克松定理'''(Poincaré–Bendixson theorem),该区域内必然存在闭轨道;再结合唯一性论证,可知系统恰好存在一条极限环,所有非平衡轨道最终都趋向该闭轨道。<ref name="Strogatz" /><ref name="Perko">Perko, L. (2001). ''Differential Equations and Dynamical Systems'' (3rd ed.). Springer.</ref> | |
| − | |||
| − | |||
| − | |||
| − | |||
| − | + | 参数 <math>\mu</math> 的大小决定了振荡的形态:当 <math>\mu</math> 较小时,极限环近似为圆,振荡波形接近正弦波;当 <math>\mu</math> 较大时,极限环在相平面中显著拉伸变形,轨道呈现出缓慢漂移与快速跃迁交替的特征,这种现象称为'''松弛振荡'''(relaxation oscillation)。<ref name="Strogatz" /> | |
== 极限环的排除方法== | == 极限环的排除方法== | ||
2026年4月1日 (三) 21:50的版本
极限环(Limit Cycle)是非线性动力系统中的重要概念,用于描述系统在长期演化过程中所呈现的周期性行为。当系统的轨道不收敛于平衡点,而是趋于一条闭合曲线时,即形成周期运动。在本文中,主要讨论二维自治系统(二维相空间)中的极限环。在这一情形下,极限环可以理解为相平面中的一条孤立闭轨道,即在该闭轨道的邻域内不存在其他周期轨道。该概念最早由 Henri Poincaré 提出,是微分方程定性理论的重要组成部分,并在分析二维系统的长期行为中具有核心作用。后续内容将基于二维相空间的框架,介绍极限环的定义、性质及其判定方法以及相关应用。
定义
在研究二维相空间的动力系统时,极限环(limit cycle)通常定义为孤立的周期轨道或闭合轨道。
该轨道是某一其他轨道的极限集,即至少存在另一条轨道,其在 [math]\displaystyle{ t \to +\infty }[/math] 或 [math]\displaystyle{ t \to -\infty }[/math] 时以螺旋方式趋近于该轨道。[1][2]
为使用更形式化的语言描述这一概念,现考虑一个二维动力系统
- [math]\displaystyle{ x'(t) = V(x(t)) }[/math]
其中 [math]\displaystyle{ V : \mathbb{R}^2 \to \mathbb{R}^2 }[/math] 为光滑向量场。[3]
该系统的一条轨道(trajectory)是满足上述微分方程的光滑函数 [math]\displaystyle{ x(t) }[/math],其取值于 [math]\displaystyle{ \mathbb{R}^2 }[/math];而该函数在相空间中的像称为其轨道(orbit),即 [math]\displaystyle{ \mathbb{R}^2 }[/math] 的一个子集。[1]
如果一条轨道不是常数解,并且存在 [math]\displaystyle{ T\gt 0 }[/math] 使得 [math]\displaystyle{ x(t+T)=x(t) }[/math] 对所有 [math]\displaystyle{ t \in \mathbb{R} }[/math] 成立,则称该轨道为周期轨道(periodic trajectory),其在相空间中的像称为闭轨道(closed orbit)或循环(cycle)。[3]
在二维系统中,极限环可以看作是一类孤立的周期轨道,即其附近不存在其他周期轨道。
性质
极限环作为二维动力系统中的孤立闭轨道,具有一系列重要的几何与动力学性质。这些性质从几何结构出发,逐步揭示了极限环如何约束其邻域内轨道的行为,并最终决定系统的长期演化趋势。
拓扑性质
极限环是相平面上的一条简单闭曲线。若尔当曲线定理(Jordan curve theorem)指出,任何简单闭曲线都将平面划分为两个互不相交的区域,即内部区域与外部区域,且两者只以该曲线为公共边界。[4]这一性质对极限环邻域内的轨道行为有根本性的约束:内部轨道与外部轨道处于拓扑上完全隔离的两个区域,无法在不穿越极限环本身的情况下相互转化。[5]
此外,极限环内部至少包含系统的一个平衡点(即满足 [math]\displaystyle{ V'(p) = 0 }[/math] 的点 [math]\displaystyle{ p }[/math])。若极限环内部只有一个平衡点,则该点必为汇、源或中心,而不能是鞍点。[4]
不变性
极限环是相空间中的不变集:若系统的初始条件恰好落在该闭轨道上,则对应解将永远沿该轨道周期性演化。[3]
这一性质来源于微分方程解的存在唯一性定理(Picard定理):相空间中过同一点只能有一条轨道。因此,一旦轨道进入极限环,它就不可能离开;反之,极限环外的轨道也不能穿入极限环。[6]
孤立性
极限环与一般闭轨道的本质区别在于其孤立性——在极限环任意足够小的邻域内,不存在其他闭轨道。[3]这一点可与可积哈密顿系统作对比:在哈密顿系统中,闭轨道以连续族的形式出现(轨道之间由能量守恒量连续参数化),任意一条闭轨道的邻域内都存在无穷多条其他闭轨道;而极限环则是耗散系统的产物,其孤立性正反映了系统能量耗散与补充之间的精确平衡。[7]
孤立性的一个直接推论是:极限环邻域内的轨道不再是闭合曲线,而只能呈螺旋状趋近或远离极限环。这就引出了极限环的稳定性问题。
稳定性分类
根据邻近轨道的渐近行为,极限环可分为三类:[3]
- 稳定极限环(stable limit cycle):内外两侧的邻近轨道在 [math]\displaystyle{ t \to +\infty }[/math] 时均螺旋趋近于该极限环,是系统自发产生稳定周期振荡的典型来源。
- 不稳定极限环(unstable limit cycle):内外两侧的邻近轨道在 [math]\displaystyle{ t \to +\infty }[/math] 时均螺旋远离该极限环;若将时间反向([math]\displaystyle{ t \to -\infty }[/math]),则轨道趋近。
- 半稳定极限环(semi-stable limit cycle):一侧轨道趋近,另一侧轨道远离。半稳定极限环通常出现在系统参数的临界值处,对应一个稳定极限环与一个不稳定极限环合并消亡的鞍结分岔(saddle-node bifurcation of periodic orbits)。[8]
上述分类在直觉上清晰,但要对稳定性给出严格的定量判断,需借助庞加莱映射(Poincaré map)。其做法是:在相平面内作一条与极限环横截相交的截线,记录轨道每次穿越截线的位置,从而将连续系统的稳定性问题转化为一个离散映射的不动点稳定性问题。[9]若该映射在极限环对应不动点处的导数(即特征乘子,characteristic multiplier)的模小于1,则极限环渐近稳定;若模大于1,则不稳定。值得注意的是,这一稳定性判断与所选截线的具体位置无关。[10]
作为极限集
以上性质合在一起,揭示了极限环在系统长期行为中的核心作用。对于稳定极限环而言,其吸引域内所有轨道最终都将收敛至该极限环——无论初始条件如何,系统的长期行为都趋向同一个周期振荡。在动力系统的语言中,稳定极限环是其吸引域内轨道的ω极限集(即轨道在 [math]\displaystyle{ t \to +\infty }[/math] 时的极限集合);相应地,不稳定极限环则是对应轨道的α极限集(即 [math]\displaystyle{ t \to -\infty }[/math] 时的极限集合)。[3]
因此,极限环描述的不是某一特殊初始条件下的解,而是系统在大范围初始条件下共同趋向的渐近行为。在生物节律、电子振荡器、化学振荡反应等领域,极限环正是系统自发产生稳定周期行为的数学根源。[4]
典型实例
简单构造的极限环
极限环可以通过显式构造二维系统的径向与角向动力学来得到。考虑如下系统: [math]\displaystyle{ \begin{cases} \dot{x} = x(1 - r^2) - y, \\ \dot{y} = y(1 - r^2) + x, \end{cases} }[/math]
其中 [math]\displaystyle{ r^2 = x^2 + y^2 }[/math]。将系统改写为极坐标形式,得:
[math]\displaystyle{ \begin{cases} \dot{r} = r(1 - r^2), \\ \dot{\theta} = 1. \end{cases} }[/math] 两个方程可独立分析:角变量 [math]\displaystyle{ \theta }[/math] 以常数速度单调增加,轨道因此始终绕原点旋转;径向变量 [math]\displaystyle{ r }[/math] 满足一维方程 [math]\displaystyle{ \dot{r} = r(1 - r^2) }[/math],在 [math]\displaystyle{ r = 0 }[/math] 与 [math]\displaystyle{ r = 1 }[/math] 处存在平衡点。对 [math]\displaystyle{ r = 1 }[/math] 的稳定性分析表明:当 [math]\displaystyle{ 0 \lt r \lt 1 }[/math] 时 [math]\displaystyle{ \dot{r} \gt 0 }[/math],轨道向外增长;当 [math]\displaystyle{ r \gt 1 }[/math] 时 [math]\displaystyle{ \dot{r} \lt 0 }[/math],轨道向内收缩。因此,只要 [math]\displaystyle{ r \neq 0 }[/math],径向分量都将单调收敛至 [math]\displaystyle{ r = 1 }[/math]。
综合两个方向的行为:轨道在趋向单位圆的过程中持续旋转,在相平面中呈螺旋形,并最终收敛到 [math]\displaystyle{ r = 1 }[/math] 上的周期运动。该单位圆即为系统唯一的稳定极限环,除原点外的所有轨道均以其为 ω 极限集。[3]
Van der Pol 振子
Van der Pol 振子是研究极限环的经典非线性系统之一。该方程最初由荷兰物理学家巴尔塔萨·范德波尔(Balthasar van der Pol)在研究电子管振荡电路时提出,[11]其方程为: [math]\displaystyle{ \ddot{x} - \mu (1 - x^2)\dot{x} + x = 0, }[/math] 其中 [math]\displaystyle{ \mu \gt 0 }[/math] 为参数,控制非线性阻尼的强度。将其改写为二维一阶系统: [math]\displaystyle{ \begin{cases} \dot{x} = y, \\ \dot{y} = \mu (1 - x^2)y - x. \end{cases} }[/math] 该系统在原点处存在唯一平衡点。对原点处的线性化矩阵进行分析,其特征值实部为正,因此原点是不稳定焦点,附近轨道向外螺旋发散。[3]
另一方面,通过估计能量函数可以证明,系统轨道被限制在相平面的某个有界环形区域内,且该区域内不含平衡点。由庞加莱-本迪克松定理(Poincaré–Bendixson theorem),该区域内必然存在闭轨道;再结合唯一性论证,可知系统恰好存在一条极限环,所有非平衡轨道最终都趋向该闭轨道。[3][1]
参数 [math]\displaystyle{ \mu }[/math] 的大小决定了振荡的形态:当 [math]\displaystyle{ \mu }[/math] 较小时,极限环近似为圆,振荡波形接近正弦波;当 [math]\displaystyle{ \mu }[/math] 较大时,极限环在相平面中显著拉伸变形,轨道呈现出缓慢漂移与快速跃迁交替的特征,这种现象称为松弛振荡(relaxation oscillation)。[3]
极限环的排除方法
梯度系统(Gradient systems)
梯度系统提供了一类可以直接排除极限环存在性的动力系统。对于二维梯度系统,其动力学方程可以写为
- [math]\displaystyle{ x' = -\nabla V(x) }[/math]
其中 [math]\displaystyle{ V : \mathbb{R}^2 \to \mathbb{R} }[/math] 为光滑标量函数,称为势函数(potential function)。在这种系统中,系统的运动方向始终沿着势函数下降最快的方向。
沿着系统轨道计算函数 [math]\displaystyle{ V(x(t)) }[/math] 的时间导数可得
- [math]\displaystyle{ \frac{d}{dt}V(x(t)) = \nabla V(x)\cdot x' = -|\nabla V(x)|^2 \le 0 }[/math]
因此,[math]\displaystyle{ V }[/math] 在轨道上是单调不增的函数,并且除了在平衡点处外严格下降。这意味着系统在演化过程中始终向势函数较小的方向运动。
由于周期轨道要求系统状态在一个周期后返回初始点,从而使所有系统量恢复原值,而 [math]\displaystyle{ V }[/math] 在轨道上严格单调变化,因此在梯度系统中不可能存在闭轨道或极限环。因此,若一个二维系统可以写成梯度系统的形式,则可以直接排除极限环的存在。[3]
李雅普诺夫函数(Lyapunov functions)
李雅普诺夫函数方法也是排除极限环存在的一种常用工具。若存在一个光滑函数 [math]\displaystyle{ V(x) }[/math],使得沿系统轨道有
- [math]\displaystyle{ \frac{d}{dt}V(x(t)) = \nabla V(x)\cdot f(x) }[/math]
在某一区域内始终保持严格正值或严格负值,则函数 [math]\displaystyle{ V }[/math] 在该区域内沿轨道单调变化。
由于周期轨道要求系统在一个周期后返回初始状态,因此任何在轨道上单调变化的函数都无法在周期结束时恢复原值。因此,如果存在这样的李雅普诺夫函数,则该区域内不可能存在闭轨道或极限环。
因此,通过构造适当的李雅普诺夫函数,可以证明某一区域内不存在极限环,并进一步分析系统的长期行为。[1]
Dulac 准则(Dulac criterion)
Dulac 准则是排除二维系统中闭轨道的重要工具,也称为 Bendixson–Dulac 定理。考虑平面动力系统
- [math]\displaystyle{ x' = P(x,y), \qquad y' = Q(x,y) }[/math]
其中 [math]\displaystyle{ P }[/math] 和 [math]\displaystyle{ Q }[/math] 为连续可微函数。若存在一个连续可微函数 [math]\displaystyle{ B(x,y) }[/math],称为 Dulac 函数,使得
- [math]\displaystyle{ \frac{\partial (BP)}{\partial x} + \frac{\partial (BQ)}{\partial y} }[/math]
在某个单连通区域内始终保持同号且不恒为零,则该区域内不存在闭轨道,因此也不存在极限环。
该准则的思想来源于向量场散度的性质:如果适当加权后的向量场在区域内始终表现为整体收缩或扩张,则轨道无法形成封闭循环。因此 Dulac 准则常被用于证明某些平面动力系统不存在极限环。[1]
极限环存在性理论
庞加莱–本迪克松定理(Poincaré–Bendixson theorem)
李纳定理(Liénard theorem)
每一个闭合轨迹在其内部都包含一个系统的驻点 Stationary Point,即点[math]\displaystyle{ p }[/math]当有[math]\displaystyle{ V(p)=0 }[/math].本迪克森-杜拉克定理 Bendixson–Dulac theorem和庞加莱-本迪克森定理 Poincaré–Bendixson theorem分别预言了二维非线性动力系统极限环的缺失或存在。
极限环的近似定量分析
松弛振荡
弱非线性振荡
摄动与平均法
稳定、不稳定和半稳定极限环
当所有相邻轨迹在时间趋于无穷大时都接近极限环时,称之为稳定或吸引极限环(ω-极限环)。当时间趋近于负无穷时,若所有相邻轨线都逼近它,则它是不稳定的极限环(α-极限环)。如果存在一个相邻轨迹在时间趋于无穷大时螺旋进入极限环,而另一个轨迹在时间趋于负无穷大时螺旋进入极限环,则称其为半稳定极限环。
此外,还存在既不稳定也不不稳定亦不半稳定的情形,例如:邻近轨迹可能从外部接近极限环,而极限环内部由一族其他闭轨道逼近(这些轨道本身不是极限环)。
稳定极限环是吸引子(attractor)的一个例子。它对应于自我维持的振荡(oscillation):闭轨道描述了系统的周期行为,而任何偏离该轨道的微小扰动都会使系统重新回到该极限环上。
待解决的问题
一般来说,寻找极限环是一个非常困难的问题。平面上一个多项式微分方程的极限环个数问题是希尔伯特第十六题(Hilbert's sixteenth problem)第二部分的主要研究对象。
例如,在平面上是否存在一个系统 [math]\displaystyle{ x' = V(x) }[/math],其中 [math]\displaystyle{ V }[/math] 的两个分量都是二次多项式,使得该系统具有超过四个极限环。
应用
科学应用中许多自持振荡系统的仿真中,极限环具有重要意义。一些典型例子包括:
- 空气动力学中的极限环振荡[12]
- 神经元动作电位的霍奇金–赫胥黎模型(Hodgkin–Huxley model)
- 糖酵解过程中的塞尔科夫模型(Sel'kov model)[13]
- 癌细胞在局限微环境中的迁移过程可表现出极限环振荡[16]
- 一些非线性电路系统表现出极限环振荡[17],这启发了范德波尔模型(Van der Pol model)
参考文献
- ↑ 1.0 1.1 1.2 1.3 1.4 Perko, L. (2001). Differential Equations and Dynamical Systems. Springer. 引用错误:无效
<ref>标签;name属性“Perko”使用不同内容定义了多次 - ↑ Devaney, R. (1989). An Introduction to Chaotic Dynamical Systems. CRC Press.
- ↑ 3.00 3.01 3.02 3.03 3.04 3.05 3.06 3.07 3.08 3.09 3.10 3.11 Strogatz, S. H. (1994). Nonlinear Dynamics and Chaos. Westview Press.
- ↑ 4.0 4.1 4.2 Limit cycle. Wikipedia. https://en.wikipedia.org/wiki/Limit_cycle
- ↑ Coddington, E. A., & Levinson, N. (1955). Theory of Ordinary Differential Equations. McGraw-Hill.
- ↑ Lebl, J. Differential Equations for Engineers, §8.4: Limit cycles. LibreTexts. https://math.libretexts.org
- ↑ Jordan, D. W., & Smith, P. (1977). Nonlinear Ordinary Differential Equations. Clarendon Press.
- ↑ Kuznetsov, Yu. A. (2004). Elements of Applied Bifurcation Theory (3rd ed.). Springer.
- ↑ Poincaré map. Wikipedia. https://en.wikipedia.org/wiki/Poincaré_map
- ↑ Wiggins, S. (2003). Introduction to Applied Nonlinear Dynamical Systems and Chaos (2nd ed.). Springer.
- ↑ van der Pol, B. (1926). On relaxation-oscillations. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 2(11), 978–992.
- ↑ Thomas, Jeffrey P.; Dowell, Earl H.; Hall, Kenneth C. (2002). "Nonlinear Inviscid Aerodynamic Effects on Transonic Divergence, Flutter, and Limit-Cycle Oscillations". AIAA Journal, 40(4): 638. DOI: 10.2514/2.1720.
- ↑ Sel'kov, E. E. (1968). "Self-Oscillations in Glycolysis 1. A Simple Kinetic Model". European Journal of Biochemistry, 4(1): 79–86. DOI: 10.1111/j.1432-1033.1968.tb00175.x.
- ↑ Leloup, Jean-Christophe; Gonze, Didier; Goldbeter, Albert (1999). "Limit Cycle Models for Circadian Rhythms Based on Transcriptional Regulation in Drosophila and Neurospora". Journal of Biological Rhythms, 14(6): 433–448. DOI: 10.1177/074873099129000948.
- ↑ Roenneberg, Till et al. (2008). "Modelling Biological Rhythms". Current Biology, 18(17): R826–R835. DOI: 10.1016/j.cub.2008.07.017.
- ↑ Brückner, David B. et al. (2019). "Stochastic nonlinear dynamics of confined cell migration in two-state systems". Nature Physics, 15(6): 595–601. DOI: 10.1038/s41567-019-0445-4.
- ↑ Ginoux, Jean-Marc; Letellier, Christophe (2012). "Van der Pol and the history of relaxation oscillations: Toward the emergence of a concept". Chaos, 22(2): 023120. DOI: 10.1063/1.3670008.