非线性中心

重定向至：

• 非线性中心

背景

线性中心

在研究动力系统的平衡点时，常利用雅可比矩阵进行局部线性化分析。若特征值是一对纯虚数 [math]\displaystyle{ \pm i\omega }[/math]，对应线性系统的平衡点即为线性中心，其相空间轨线表现为一族同心圆或椭圆闭轨。系统在相空间的演化表现为沿着同一条闭合曲线周而复始地运转。线性中心在系统稳定性判定中属于李雅普诺夫稳定（Lyapunov stable），但不是渐进稳定（Asymptotically stable）。

非线性系统的线性化失效

对于一般的非线性系统 [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}) }[/math]，不动点 [math]\displaystyle{ \mathbf{x}^* }[/math] 附近的局部动力学行为通常可通过线性化系统来近似判定。根据哈特曼-格罗伯曼定理（Hartman-Grobman Theorem），仅当平衡点特征值实部非零（双曲平衡点）时，非线性系统的局部拓扑才与其线性化系统同胚。对于特征值为纯虚数的临界情形，线性化方法失效。此时，平衡点附近的真实拓扑结构取决于高阶非线性项。

中心的非线性扰动

为直观展示高阶非线性项如何破坏线性中心，考虑经典的二维非线性动力系统：

[math]\displaystyle{ \left\{ \begin{array}{l} \dot{x} = -y + ax(x^2+y^2) \\[6pt] \dot{y} = x + ay(x^2+y^2) \end{array} \right. }[/math]

对该系统在坐标原点 (0,0) 处进行雅可比矩阵线性化，其近似线性系统的特征值为纯虚数 [math]\displaystyle{ \lambda = \pm i }[/math]。 仅从线性化分析来看，原点表现为一个标准的线性中心。

但要揭示其非线性局部拓扑结构，需利用极坐标变换

图1：非线性中心扰动

[math]\displaystyle{ \left\{ \begin{array}{l} x = r\cos\theta \\ y = r\sin\theta \end{array} \right. }[/math]

将极坐标代入原方程并化简，系统可严格转换为如下形式：

[math]\displaystyle{ \left\{ \begin{array}{l} \dot{r} = ar^3 \\[6pt] \dot{\theta} = 1 \end{array} \right. }[/math]

基于极坐标方程，高阶非线性参数 [math]\displaystyle{ a }[/math] 对系统相图结构的决定性影响表现如下：

• 稳定焦点：当 [math]\displaystyle{ a \lt 0 }[/math] 时，对于任意 [math]\displaystyle{ r \gt 0 }[/math]，径向演化 [math]\displaystyle{ \dot{r} \lt 0 }[/math]。轨线随时间向内螺旋收缩，原有闭轨被高阶项破坏，原点退化为渐近稳定的焦点。
• 非线性中心：当且仅当 [math]\displaystyle{ a = 0 }[/math] 时，系统无径向耗散 [math]\displaystyle{ （\dot{r} \equiv 0） }[/math]，角速度恒定。轨线在相空间中形成一族完美的同心圆闭轨，原点维持中心结构。
• 不稳定焦点：当 [math]\displaystyle{ a \gt 0 }[/math] 时，径向演化 [math]\displaystyle{ \dot{r} \gt 0 }[/math]。轨线随时间向外螺旋发散，原点退化为不稳定焦点。

该经典实例表明，当动力系统处于纯虚数特征值的临界状态时，其相图拓扑极度脆弱，局部的真实稳定性完全由极微小的高阶非线性扰动（本例中的 [math]\displaystyle{ a r^3 }[/math] 项）所主导。

此外，以下两类典型的非双曲边界情形也极易发生定性改变：

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

退化结点

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

非孤立不动点

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

非线性中心

在非线性系统 [math]\displaystyle{ \dot{x}=f(x) }[/math] 的平衡点 [math]\displaystyle{ x_0 }[/math] 处，线性化后的雅可比矩阵具有一对纯虚数特征值 [math]\displaystyle{ \pm i\omega }[/math] (即线性部分为“中心”)，且加入非线性项后，系统在该平衡点领域内仍然存在一簇闭合轨道（周期解），则称该平衡点 [math]\displaystyle{ x_0 }[/math] 为非线性中心。

与焦点的区分：

• 非线性中心：平衡点邻域内所有充分靠近的轨线均为围绕该点的闭合曲线，既不趋近也不远离平衡点。
• 焦点：轨线呈螺旋状——稳定焦点向内收敛，不稳定焦点向外发散。
• 本质区别：中心要求邻域内所有轨线闭合；焦点不存在闭合轨线。

非线性中心判据

经典中心问题

经典中心问题（Classical Center Problem）由 Poincaré 正式提出：

> 对于线性部分为中心的非线性平面系统，在什么条件下，原点邻域仍存在一族闭合轨道（即保持非线性中心）？

线性化分析无法区分中心与焦点，必须考察高阶非线性项的累积效应。

解决经典中心问题最直接的方法是逐阶计算系统的 Lyapunov 量，通过其消失条件判定中心是否存在。

Lyapunov量与焦点量

定义

考虑平面系统

[math]\displaystyle{ \left\{ \begin{array}{l} \dot{x} = -y + P(x,y) \\[6pt] \dot{y} = x + Q(x,y) \end{array} \right. }[/math]

其中 [math]\displaystyle{ P(x,y), Q(x,y) }[/math] 为高阶非线性项。通过形式级数法构造

[math]\displaystyle{ F(x,y) = x^2 + y^2 + F_3(x,y) + F_4(x,y) + \cdots }[/math]

使得沿系统轨线

[math]\displaystyle{ \frac{dF}{dt} = L_1(x^2+y^2)^2 + L_2(x^2+y^2)^3 + L_3(x^2+y^2)^4 + \cdots }[/math]

系数 [math]\displaystyle{ L_1, L_2, L_3, \ldots }[/math] 称为Lyapunov量（Lyapunov quantities）或焦点量（focal values）。

判定规则

• 若所有 [math]\displaystyle{ L_k = 0 }[/math]（<[math]\displaystyle{ k = 1, 2, 3, \ldots }[/math]），则原点为非线性中心。
• 若第一个非零量为 [math]\displaystyle{ L_k \neq 0 }[/math]，则原点为k阶细焦点（fine focus of order k），其稳定性由 [math]\displaystyle{ L_k }[/math] 符号决定：
  • [math]\displaystyle{ L_k \lt 0 }[/math]：k阶稳定焦点；
  • [math]\displaystyle{ L_k \gt 0 }[/math]：k阶不稳定焦点。

说明

Lyapunov量的计算随阶数指数增长，通常只能计算前几项。对于具体系统，逐阶验证所有 [math]\displaystyle{ L_k=0 }[/math] 往往不现实，因此需要寻找更高效的充分条件。

幸运的是，对于两类重要的动力系统——保守系统和可逆系统——无需逐阶计算 Lyapunov 量，即可直接判定非线性中心的存在。

保守系统中的非线性中心

定义

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

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

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

示例：双阱势系统

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

图4：无阻尼钟摆

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