更改

添加157字节 、 2026年4月1日 (星期三)
第109行: 第109行:  
参数 <math>\mu</math> 的大小决定了振荡的形态:当 <math>\mu</math> 较小时,极限环近似为圆,振荡波形接近正弦波;当 <math>\mu</math> 较大时,极限环在相平面中显著拉伸变形,轨道呈现出缓慢漂移与快速跃迁交替的特征,这种现象称为'''松弛振荡'''(relaxation oscillation)。<ref name="Strogatz" />
 
参数 <math>\mu</math> 的大小决定了振荡的形态:当 <math>\mu</math> 较小时,极限环近似为圆,振荡波形接近正弦波;当 <math>\mu</math> 较大时,极限环在相平面中显著拉伸变形,轨道呈现出缓慢漂移与快速跃迁交替的特征,这种现象称为'''松弛振荡'''(relaxation oscillation)。<ref name="Strogatz" />
   −
== 极限环的排除方法==
+
== 极限环的排除方法 ==
   −
=== 梯度系统(Gradient systems) ===
+
判断一个系统是否存在极限环,在一般情形下是困难的。然而,对于某些特殊结构的系统,可以通过分析系统的整体性质来直接排除极限环的存在。以下三种方法各自利用了不同的结构特征,适用范围依次扩大。
   −
梯度系统提供了一类可以直接排除极限环存在性的动力系统。对于二维梯度系统,其动力学方程可以写为 
+
=== 梯度系统 ===
    +
最直接的一类是'''梯度系统'''(gradient system)。若二维系统可以写成
 
: <math>x' = -\nabla V(x)</math>
 
: <math>x' = -\nabla V(x)</math>
 +
其中 <math>V : \mathbb{R}^2 \to \mathbb{R}</math> 为光滑势函数(potential function),则沿任意轨道计算 <math>V</math> 的时间导数,得
 +
: <math>\frac{d}{dt}V(x(t)) = \nabla V(x)\cdot x' = -|\nabla V(x)|^2 \leq 0.</math>
 +
这说明 <math>V</math> 沿轨道单调不增,且在非平衡点处严格递减——系统始终沿势函数下降最快的方向运动。而周期轨道要求系统在一个周期后回到初始状态,<math>V</math> 的值也必须恢复原值,这与上述严格单调性矛盾。因此,梯度系统中不可能存在闭轨道或极限环。<ref name="Strogatz" />
   −
其中 <math>V : \mathbb{R}^2 \to \mathbb{R}</math> 为光滑标量函数,称为势函数(potential function)。在这种系统中,系统的运动方向始终沿着势函数下降最快的方向。
+
=== 李雅普诺夫函数 ===
   −
沿着系统轨道计算函数 <math>V(x(t))</math> 的时间导数可得 
+
梯度系统的条件较为严格,要求系统具有特定的势函数结构。更一般地,若能构造一个光滑函数 <math>V(x)</math>,使其沿轨道的时间导数
 +
: <math>\dot{V}(x) = \nabla V(x)\cdot f(x)</math>
 +
在某一区域内处处保持严格同号(恒正或恒负),则 <math>V</math> 在该区域内沿任意轨道单调变化。由同样的理由,该区域内不可能存在周期轨道。这样的函数称为该区域上的'''李雅普诺夫函数'''(Lyapunov function)。<ref name="Perko" />
   −
: <math>\frac{d}{dt}V(x(t)) = \nabla V(x)\cdot x' = -|\nabla V(x)|^2 \le 0</math>
+
李雅普诺夫函数方法比梯度系统更灵活,但其局限在于需要预先猜测或构造出合适的 <math>V</math>,而这在一般情形下并无系统性方法。
   −
因此,<math>V</math> 在轨道上是单调不增的函数,并且除了在平衡点处外严格下降。这意味着系统在演化过程中始终向势函数较小的方向运动。
+
=== Dulac 准则 ===
−
 
  −
由于周期轨道要求系统状态在一个周期后返回初始点,从而使所有系统量恢复原值,而 <math>V</math> 在轨道上严格单调变化,因此在梯度系统中不可能存在闭轨道或极限环。因此,若一个二维系统可以写成梯度系统的形式,则可以直接排除极限环的存在。<ref name="Strogatz" />
  −
 
  −
=== 李雅普诺夫函数(Lyapunov functions) ===
  −
 
  −
李雅普诺夫函数方法也是排除极限环存在的一种常用工具。若存在一个光滑函数 <math>V(x)</math>,使得沿系统轨道有
  −
 
  −
: <math>\frac{d}{dt}V(x(t)) = \nabla V(x)\cdot f(x)</math>
  −
 
  −
在某一区域内始终保持严格正值或严格负值,则函数 <math>V</math> 在该区域内沿轨道单调变化。
  −
 
  −
由于周期轨道要求系统在一个周期后返回初始状态,因此任何在轨道上单调变化的函数都无法在周期结束时恢复原值。因此,如果存在这样的李雅普诺夫函数,则该区域内不可能存在闭轨道或极限环。
  −
 
  −
因此,通过构造适当的李雅普诺夫函数,可以证明某一区域内不存在极限环,并进一步分析系统的长期行为。<ref name="Perko" />
  −
 
  −
=== Dulac 准则(Dulac criterion) ===
  −
 
  −
Dulac 准则是排除二维系统中闭轨道的重要工具,也称为 Bendixson–Dulac 定理。考虑平面动力系统 
  −
 
  −
: <math>x' = P(x,y), \qquad y' = Q(x,y)</math>
  −
 
  −
其中 <math>P</math> 和 <math>Q</math> 为连续可微函数。若存在一个连续可微函数 <math>B(x,y)</math>,称为 Dulac 函数,使得 
      +
对于无法直接找到李雅普诺夫函数的情形,'''Dulac 准则'''(Dulac criterion,又称 Bendixson–Dulac 定理)提供了另一种排除方法。考虑平面系统
 +
: <math>x' = P(x,y), \qquad y' = Q(x,y),</math>
 +
其中 <math>P</math>、<math>Q</math> 连续可微。若存在连续可微函数 <math>B(x,y)</math>(称为 Dulac 函数),使得
 
: <math>\frac{\partial (BP)}{\partial x} + \frac{\partial (BQ)}{\partial y}</math>
 
: <math>\frac{\partial (BP)}{\partial x} + \frac{\partial (BQ)}{\partial y}</math>
 +
在某个单连通区域内处处同号且不恒为零,则该区域内不存在闭轨道,从而不存在极限环。<ref name="Perko" />
   −
在某个单连通区域内始终保持同号且不恒为零,则该区域内不存在闭轨道,因此也不存在极限环。
+
其思想与李雅普诺夫方法一脉相承:上式本质上是加权向量场 <math>(BP, BQ)</math> 的散度。若散度处处同号,则向量场在该区域内整体表现为纯收缩或纯扩张,轨道无法首尾相接构成闭合回路。当 <math>B \equiv 1</math> 时,该准则退化为经典的 '''Bendixson 判据'''。
   −
该准则的思想来源于向量场散度的性质:如果适当加权后的向量场在区域内始终表现为整体收缩或扩张,则轨道无法形成封闭循环。因此 Dulac 准则常被用于证明某些平面动力系统不存在极限环。<ref name="Perko" />
+
三种方法在适用范围上逐步推广:梯度系统是最特殊的情形,李雅普诺夫函数方法将其推广到更一般的单调量,而 Dulac 准则则通过引入权函数 <math>B</math> 进一步扩大了适用范围。
    
==极限环存在性理论==
 
==极限环存在性理论==
152

个编辑