非线性中心

在非线性动力学中，非线性中心（Nonlinear Center）是指相平面上某一孤立不动点附近的所有轨迹均为闭合曲线（周期轨道）的动力学现象。在线性系统中，若雅可比矩阵的特征值为纯虚数，则该不动点为线性中心；然而，对于非线性系统，特征值为纯虚数并不能保证该点一定是中心——微小的非线性项就可能将中心转变为稳定或不稳定的焦点（Focus）。

为了判定非线性系统的中心稳定性，需要借助系统的附加结构性质。两类重要的系统——保守系统(Conservative system)与可逆系统（Reversible system）——提供了判定非线性中心的有力工具。这两种方法分别从能量守恒的几何性质和时间反演的对称性出发，确立了非线性中心存在的充分条件。

保守系统中的非线性中心

定义

在经典力学与非线性动力学中，若一个二阶系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math]（其中 [math]\displaystyle{ \mathbf{x} \in \mathbb{R}^2 }[/math]）存在一个光滑的实值函数 [math]\displaystyle{ E(\mathbf{x}) }[/math]，使得沿任意轨迹 [math]\displaystyle{ \mathbf{x}(t) }[/math] 均有 [math]\displaystyle{ \frac{dE}{dt} = 0 }[/math]，且 [math]\displaystyle{ E(\mathbf{x}) }[/math] 在任何开集上不为常数，则称该系统为保守系统，[math]\displaystyle{ E(\mathbf{x}) }[/math] 称为守恒量或第一积分。

定理（保守系统的非线性中心）

考虑系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math]，其中 [math]\displaystyle{ \mathbf{x} = (x, y) \in \mathbb{R}^2 }[/math]，[math]\displaystyle{ \mathbf{f} }[/math] 是连续可微的。假设存在一个守恒量 [math]\displaystyle{ E(\mathbf{x}) }[/math]，且 [math]\displaystyle{ \mathbf{x}^* }[/math] 是一个孤立不动点（即在其某个邻域内没有其他不动点）。若 [math]\displaystyle{ \mathbf{x}^* }[/math] 是 [math]\displaystyle{ E(\mathbf{x}) }[/math] 的一个局部极小值（或局部极大值），则所有充分接近 [math]\displaystyle{ \mathbf{x}^* }[/math] 的轨迹都是闭合的，即 [math]\displaystyle{ \mathbf{x}^* }[/math] 是一个非线性中心。

证明思路： 由于能量沿轨迹守恒，每条轨迹都位于 [math]\displaystyle{ E(\mathbf{x}) }[/math] 的某一等高线上。在局部极值点附近，等高线是闭合的。又因假设 [math]\displaystyle{ \mathbf{x}^* }[/math] 是孤立不动点，故在该点的小邻域内，等高线上不存在其他不动点，因此轨迹必定沿着闭合等高线运动，形成周期轨道。

图1：双阱系统能量等高线图及其相平面投影

示例：双阱势系统

考虑一个单位质量粒子在双阱势 [math]\displaystyle{ V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4 }[/math] 中的无阻尼运动。运动方程为 [math]\displaystyle{ \ddot{x} = x - x^3 }[/math]，等价系统为：

[math]\displaystyle{ \dot{x} = y, \quad \dot{y} = x - x^3. }[/math]

系统具有守恒量（总能量）：

[math]\displaystyle{ E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4. }[/math]

不动点为 [math]\displaystyle{ (0,0) }[/math] 和 [math]\displaystyle{ (\pm 1, 0) }[/math]。在 [math]\displaystyle{ (0,0) }[/math] 处，能量曲面的等高线呈现鞍点特征；而在 [math]\displaystyle{ (\pm 1, 0) }[/math] 处，[math]\displaystyle{ E(x, y) }[/math] 达到局部极小值。根据定理，这两个不动点均为非线性中心，附近存在一族闭轨（周期轨道）。

可逆系统中的非线性中心

定义

如果存在一个相空间到自身的映射 [math]\displaystyle{ R(\mathbf{x}) }[/math]，满足 [math]\displaystyle{ R^2(\mathbf{x}) = \mathbf{x} }[/math]（即 [math]\displaystyle{ R }[/math] 是一个对合变换），使得系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math] 在变换 [math]\displaystyle{ t \to -t }[/math]，[math]\displaystyle{ \mathbf{x} \to R(\mathbf{x}) }[/math] 下保持不变，则称该系统为可逆系统。 最常见的二维情形是系统具有关于 [math]\displaystyle{ x }[/math] 轴的反射对称性：即变换 [math]\displaystyle{ t \to -t }[/math]，[math]\displaystyle{ y \to -y }[/math] 保持系统不变。这等价于要求方程 [math]\displaystyle{ \dot{x} = f(x,y) }[/math]，[math]\displaystyle{ \dot{y} = g(x,y) }[/math] 中的 [math]\displaystyle{ f }[/math] 关于 [math]\displaystyle{ y }[/math] 是奇函数（[math]\displaystyle{ f(x,-y) = -f(x,y) }[/math]），且 [math]\displaystyle{ g }[/math] 关于 [math]\displaystyle{ y }[/math] 是偶函数（[math]\displaystyle{ g(x,-y) = g(x,y) }[/math]）。

定理（可逆系统的非线性中心）

假设原点 [math]\displaystyle{ \mathbf{x}^* = \mathbf{0} }[/math] 是一个连续可微可逆系统的线性中心。若系统是可逆的，则原点也是一个非线性中心，即所有充分接近原点的轨迹都是闭合的。 证明思路： 考虑一条从正 [math]\displaystyle{ x }[/math] 轴出发靠近原点的轨迹。由于线性部分是中心，该轨迹将绕原点旋转并最终与负 [math]\displaystyle{ x }[/math] 轴相交。利用可逆性，将轨迹关于 [math]\displaystyle{ x }[/math] 轴反射并反转时间方向，将得到一条具有相同端点但箭头相反的轨迹。这两条轨迹共同构成一个闭合曲线。因此原点被一族闭轨包围。

示例：无阻尼钟摆

无阻尼钟摆的非线性本质通过小角度的近似值[math]\displaystyle{ sin𝜃 ≈ 𝜃 }[/math]被忽略。利用相图的分析钟摆旋转到顶部这一大角度的区域，在缺少阻力和外部吸引力的情况下，钟摆由下式决定：

[math]\displaystyle{ \frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0. }[/math]

令 [math]\displaystyle{ \omega = \sqrt{g/L} }[/math]，[math]\displaystyle{ \tau = \omega t }[/math] 进行无量纲化，得：

[math]\displaystyle{ \ddot{\theta} + \sin\theta = 0. }[/math]

写成相平面系统（[math]\displaystyle{ v = \dot{\theta} }[/math]）：

[math]\displaystyle{ \dot{\theta} = v, \quad \dot{v} = -\sin\theta. }[/math]

该系统在变换 [math]\displaystyle{ \tau \to -\tau }[/math]，[math]\displaystyle{ v \to -v }[/math] 下保持不变，因此是可逆的。不动点 [math]\displaystyle{ (\theta^*, v^*) = (0,0) }[/math] 处的雅可比矩阵为 [math]\displaystyle{ \begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix} }[/math]，特征值为 [math]\displaystyle{ \pm i }[/math]，故为线性中心。根据可逆系统的非线性中心定理，[math]\displaystyle{ (0,0) }[/math] 也是一个非线性中心，对应钟摆在最低点的微小摆动。此外，系统也具有守恒量 [math]\displaystyle{ E = \frac{1}{2}v^2 - \cos\theta }[/math]，且 [math]\displaystyle{ (0,0) }[/math] 是 [math]\displaystyle{ E }[/math] 的局部极小值点，这从保守系统的角度也印证了该结论。

图2：无阻尼钟摆