非线性中心定理

在非线性动力学中，非线性中心（Nonlinear Center）是指相平面上某一孤立不动点附近的所有轨迹均为闭合曲线（周期轨道）的动力学现象。

图1：不动点的类型与稳定性

图2：不动点分类树

非线性系统的线性化失效

对于一般的非线性系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math]，不动点 [math]\displaystyle{ \mathbf{x}^* }[/math] 附近的局部动力学行为通常可通过线性化系统来近似判定。根据哈特曼-格罗伯曼定理（Hartman-Grobman Theorem），若不动点是双曲不动点（即雅可比矩阵的所有特征值实部均不为零），则非线性系统与其线性化系统在不动点附近是局部拓扑等价的，线性化分析能够正确反映系统的稳定性等定性行为。

但当不动点为非双曲不动点时（即特征值为纯虚数），该定理不再适用。此时，高阶非线性项主导了不动点邻域内的动力学行为，微小的扰动即可能导致线性化预测的彻底失效。以下三类典型的非双曲边界情形极易发生定性改变：

图3：非线性系统线性化失效典型情形

中心

当雅可比矩阵的特征值为一对共轭纯虚数（[math]\displaystyle{ \lambda = \pm i\omega }[/math]）时，线性化系统预测不动点为线性中心（Linear Center），其周围环绕一族闭合轨迹。但非线性项的引入可能使特征值获得非零实部，使闭合轨迹族退化为螺旋收敛的稳定焦点（Stable Focus）或螺旋发散的不稳定焦点（Unstable Focus）。因此，该线性中心是结构不稳定（Structurally Unstable）的。

退化结点

当雅可比矩阵具有重特征值且仅对应单个独立特征向量时，线性化系统预测不动点为退化结点（Degenerate Node）。虽然其稳定性（稳定或不稳定）在非线性扰动下通常得以保持，但其局部轨迹的几何形状可能发生改变（例如从结点变为焦点）。

非孤立不动点

当雅可比矩阵存在零特征值时，线性化系统预测存在一整条直线（或平面）的不动点。这类情形称为非孤立不动点（Non-Isolated Fixed Point）。然而，任意微小的非线性扰动通常会使这一连续统塌缩为若干个孤立的不动点，或完全消失。这类情形对应于系统处于分岔的临界状态，其动力学行为由高阶非线性项主导。

综上，对于特征值为纯虚数的中心情形，线性化分析无法给出可靠结论。为判定非线性系统的中心稳定性，需要借助系统的附加结构性质。两类重要的系统——保守系统与可逆系统——提供了判定非线性中心的有力工具。这两种方法分别从能量守恒的几何性质和时间反演的对称性出发，确立了非线性中心存在的充分条件。

保守系统中的非线性中心

定义

在经典力学与非线性动力学中，若一个二阶系统 [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] 是孤立不动点，故在该点的小邻域内，等高线上不存在其他不动点，因此轨迹必定沿着闭合等高线运动，形成周期轨道。

内在联系： 该定理揭示了保守系统结构与非线性中心之间的因果联系。在一般非线性系统中，线性化分析所得的中心极易因非线性项的扰动而失去闭合性，演化为螺旋收敛或发散的焦点。然而，保守系统固有的能量守恒性质排除了此种演变：若轨迹以螺旋方式趋近或远离不动点，则其守恒量将随时间单调变化，这与守恒量沿轨迹恒定的定义相矛盾。由此推知，只要不动点对应于守恒量的局部极值，其邻近轨迹便被约束于闭合的等高线上，从而保证了非线性中心的存在。

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

示例：双阱势系统

考虑一个单位质量粒子在双阱势 [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]被忽略。利用相图的分析钟摆旋转到顶部这一大角度的区域，在缺少阻力和外部吸引力的情况下，钟摆由下式决定：

图5：无阻尼钟摆

[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] 的局部极小值点，这从保守系统的角度也印证了该结论。