集智百科 - 复杂系统|人工智能|复杂科学|复杂网络|自组织 - 用户贡献 [zh-cn]

非线性中心

2026-07-12T07:42:28Z

王晔：/* 定理（保守系统的非线性中心） */

#重定向 [[非线性中心]]{{DISPLAYTITLE:非线性中心}}
'''非线性中心'''（'''Nonlinear Center'''）是[[非线性动力学]]中的基本结构之一。当系统[[奇点]]（平衡点）的局部邻域内布满环绕该点的嵌套闭合相轨迹时，该奇点即为非线性中心<ref>Strogatz, S. H. (2015). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering. CRC press.</ref>。由于一般的中心结构缺乏结构稳定性，极易在微弱的非线性扰动下退化为向内或向外螺旋的[[焦点]]，因此严格的非线性中心通常仅存在于受到特定物理法则约束的系统中，最典型的是具备能量守恒的[[保守系统]]或具有[[时间反演对称性]]的[[可逆系统]]。探讨非线性中心在扰动下如何[[分岔]]出[[极限环]]的研究（即'''[[庞加莱-李雅普诺夫中心问题]]'''），是揭示系统自激振荡机制与探索'''[[希尔伯特第十六问题]]'''的核心途径<ref>Ilyashenko, Y. (2002). Centennial history of Hilbert's 16th problem. Bulletin of the American Mathematical Society, 39(3), 301-354.</ref>。

=='''背景'''==
==='''线性中心'''===
在动力系统理论中，[[中心]]（Center）这一概念最初起源于对纯[[线性系统]]相空间轨线的定性分类。1881年，法国数学家[[庞加莱]]（Henri Poincaré）在创立微分方程定性理论时，首次根据系统矩阵的特征值性质，将二维平衡点划分为[[中心]]、[[焦点]]、[[结点]]与[[鞍点]]<ref>Poincaré, H. (1881). Mémoire sur les courbes définies par une équation différentielle. ''Journal de mathématiques pures et appliquées'', 375-422.</ref>，参见[[二维线性动力系统]]。按照该定义，当[[雅可比矩阵]]拥有一对纯虚数特征值 <math>\pm i\omega </math> 时，对应纯[[线性系统]]的平衡点即表现为[[线性中心]]。在相空间中，其轨线呈现为一族完美的同心圆或椭圆闭轨（见图2a），代表着一种无耗散的理想周期振荡状态<ref>Strogatz, S. H. (2018). ''Nonlinear dynamics and chaos''. CRC press.</ref>。

==='''线性化方法'''===
对于一般[[非线性动力系统]]，<math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>，研究[[不动点]]<math>\mathbf{x}^* </math>附近局部动力学行为的标准手段是[[线性化方法]]。该方法通过在不动点处进行[[泰勒级数展开]]并忽略高阶项，将原系统近似简化为[[线性系统]]：

<math>\dot{\mathbf{x}} \approx J(\mathbf{x}^*) (\mathbf{x} - \mathbf{x}^*)</math>

其中<math> J(\mathbf{x}^*) </math>为[[雅可比矩阵]]。该方法利用线性系统的[[特征值]]性质，定性判定原[[非线性系统]]在不动点邻域内的[[拓扑结构]]与[[稳定性]]。

==='''非线性系统的线性化失效'''===
根据[[Hartman-Grobman定理]]<ref>Perko, L. (2013). ''Differential Equations and Dynamical Systems''. Springer.</ref>，局部线性化系统与原非线性系统之间的[[拓扑同胚]]关系，仅在平衡点为[[双曲平衡点]]（特征值实部均非零）时成立。

而当系统处于非双曲的拓扑临界状态时（例如雅可比矩阵存在实部严格为零的特征值，即[[线性中心]]），被忽略的高阶非线性项将对系统的真实动力学行为起决定性作用，此时局部线性化方法失效，无法通过线性近似推断原系统的拓扑性质。

==='''中心的非线性扰动'''===
线性化失效情形下，理想的闭轨结构极度脆弱。高阶非线性项的介入通常会打破原有的局部守恒性，使闭轨演变为渐近吸引或排斥的[[焦点]]。仅在满足特定对称性或守恒律的条件下，系统才能抵御高阶项扰动，形成非线性中心。

为直观展示高阶非线性扰动对临界拓扑系统的决定性影响，可考察以下二维非线性动力系统：
:<math>\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>

[[文件:非线性中心扰动.png|缩略图|右|图1：非线性中心扰动|480px]]

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

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

:<math>\left\{ \begin{array}{l} x = r\cos\theta \\ y = r\sin\theta \end{array} \right.</math>

将极坐标代入原方程并化简，系统可严格转换为如下形式：
:<math>\left\{ \begin{array}{l} \dot{r} = ar^3 \\[6pt] \dot{\theta} = 1 \end{array} \right.</math>
基于极坐标方程，高阶非线性参数 <math>a</math> 对系统相图结构的决定性影响表现如下（见图1）：

* '''稳定焦点'''：当 <math> a < 0 </math> 时，对于任意 <math> r > 0</math>，径向演化 <math> \dot{r} < 0</math>。轨线随时间向内螺旋收缩，原有闭轨被高阶项破坏，原点退化为渐近稳定的[[焦点]]。
*
* '''非线性中心'''：当且仅当 <math> a = 0 </math> 时，系统无径向耗散 <math>（\dot{r} \equiv 0）</math>，角速度恒定。轨线在[[相空间]]中形成一族完美的同心圆闭轨，原点维持中心结构。
*
* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。
*
该经典实例表明，当动力系统处于纯虚数特征值的临界状态时，其相图拓扑极度脆弱，局部的真实稳定性由极微小的高阶非线性扰动（本例中的 <math>a r^3 </math> 项）所主导。

[[文件:非线性系统线性化失效典型情形.png|缩略图|图2：非线性系统线性化失效典型情形|480px]]

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

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

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

=='''非线性动力系统的非线性中心'''==
在前述的线性化失效与临界扰动背景下，要使系统在非线性范畴内维持循环运动，其不动点必须满足更严格的条件。
===定义与实例===
====定义====
在非线性系统 <math>\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x})</math> 的平衡点 <math>\mathbf{x}_0</math> 处，若线性化后的雅可比矩阵具有一对纯虚数特征值 <math>\pm i\omega</math>（即线性部分为“中心”），且加入高阶非线性项后，系统在该平衡点邻域内仍然存在一簇闭合轨道（周期解），则称该平衡点 <math>\mathbf{x}_0</math> 为'''非线性中心'''<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York.</ref>。

====实例====
[[文件:Lotka-Volterra 相图.png|缩略图|右|图3：Lotka-Volterra 模型的非线性中心|260px]]
为了直观展示这一解析定义的物理与生物学意义，我们可以引入经典的 '''[[Lotka-Volterra 捕食者-猎物模型]]'''。该模型用于描述生态系统中猎物 <math>x(t)</math> 与捕食者 <math>y(t)</math> 的相互作用，其非量纲化方程为<ref>Murray, J. D. (2002). Mathematical biology: I. An introduction (Vol. 17). Springer Science & Business Media.</ref>：

:<math>\left\{ \begin{array}{l} \dot{x} = x(1-y) \\[6pt] \dot{y} = y(x-1) \end{array} \right.</math>

该系统在正象限内存在平衡点 <math>(1, 1)</math>，其雅可比矩阵特征值为 <math>\lambda = \pm i</math>，满足线性中心的前提条件。进一步分析可知，该系统存在严格守恒的能量函数 <math>H(x, y) = (x - \ln x) + (y - \ln y)</math>。由于沿系统轨线 <math>\frac{dH}{dt} \equiv 0</math>，且 <math>H(x,y)</math> 在 <math>(1,1)</math> 处达到全局极小值，高阶非线性项无法破坏轨线的闭合性。因此，该平衡点被确认为一个严格的'''非线性中心'''。在相图上，这表现为捕食者与猎物的种群数量不会趋于静态，而是围绕平衡点呈现出持续的、由初始状态决定的周期性闭轨波动（见图3）。

===性质===
除上述解析定义与直观模型外，非线性中心在动力系统中还具备以下几个核心的几何与拓扑特征：

==== 几何连续性与极限环的区别 ====
[[文件:非线性中心与极限环对比图.png|缩略图|右|图4：非线性中心与极限环对比图|400px]]
非线性中心的核心几何特征在于其闭轨的连续分布性。在中心的一个充分小的邻域内，除平衡点本身外，所有的相轨线均是嵌套的闭合周期轨线，构成一个[[连续统]]（Continuum）。

这一特征将其与[[极限环]]（Limit Cycle）区分开来（见图4）：极限环同为[[闭合周期轨线]]，但它是孤立的，即在极限环的充分小邻域内不存在其他闭合轨线，其内侧或外侧的相邻轨线会以螺旋形式渐近地趋向（或远离）该极限环。

==== 拓扑特征与庞加莱指数 ====
在[[向量场拓扑学]]中，奇点的局部结构可通过[[庞加莱指数]]（Poincaré Index）进行定量刻画<ref>Perko, L. (2013). Differential Equations and Dynamical Systems (3rd ed.). Springer Science & Business Media. ISBN 978-1-4613-0003-8.</ref>。对于[[二维非线性系统]]，沿包围非线性中心且不包含其他奇点的任意简单闭曲线，其向量场旋转的圈数恒为 1。因此，非线性中心的[[庞加莱指数]]为 <math>+1</math>。这一拓扑不变量与[[结点]]、[[焦点]]的指数相同，但与指数为 <math>-1</math> 的[[鞍点]]具有明确的拓扑差异。

==== 周期与振幅的非线性依赖 ====
在纯[[线性系统]]中，中心周围所有闭轨的运转[[周期]]恒为常数 <math>T = 2\pi/\omega</math>，与轨线的振幅大小无关。对于一般的非线性中心，由于高阶非线性项的影响，系统沿闭轨运动的周期 <math>T(r)</math> 通常会依赖于轨线的振幅（即极径 <math>r</math>）。在特殊情况下，若非线性系统不仅存在中心，且周围所有闭轨的周期依然保持恒定，则该奇点被称为“'''[[等时中心]]'''”（Isochronous Center）。等时中心的判定是[[常微分方程定性理论]]中的一项具体研究课题<ref>Chavarriga, J., & Sabatini, M. (1999). A survey of isochronous centers. Qualitative Theory of Dynamical Systems, 1(1), 1-70. </ref>。

==== 高维系统的中心流形推广 ====
对于相空间维度 <math>n \ge 3</math> 的非线性系统，中心的概念可通过'''[[中心流形定理]]'''（Center Manifold Theorem）进行高维推广。若系统在[[不动点]]处的[[雅可比矩阵]]存在一对纯虚数特征值，且其余特征值的实部均不为零，则系统在[[不动点]]附近存在一个二维的[[中心流形]]。高维系统的局部渐近动力学可降维至该[[流形]]上进行分析。若降维后的二维子系统呈现非线性中心特征，则高维系统在对应子空间内同样表现为局部的周期振荡行为<ref>Kuznetsov, Yu. A. (2004). Elements of Applied Bifurcation Theory (3rd edition). Springer.</ref>。

=='''非线性中心判据'''==

'''[[经典中心问题]]'''（Classical Center Problem）由 Poincaré 正式提出<ref>POINCARÉ H. Mémoire sur les courbes définies par une équation différentielle (I)[J]. Journal de Mathématiques Pures et Appliquées, 1881, 7(3): 375-422.</ref>：

: 给定一个[[平面解析系统]]，其线性部分为[[中心]]，如何通过系统的系数判定该平衡点究竟是中心还是[[焦点]]？

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

解决经典中心问题最直接的方法是逐阶计算系统的[[Lyapunov量]]<ref>LYAPUNOV A M. Problème général de la stabilité du mouvement[M]. Princeton: Princeton University Press, 1947.</ref>，通过其消失条件判定中心是否存在。

==='''一般情况的判定——Lyapunov量'''===
====定义====
考虑平面系统
:<math>\left\{ \begin{array}{l} \dot{x} = -y + P(x,y) \\[6pt] \dot{y} = x + Q(x,y) \end{array} \right.</math>
其中 <math>P(x,y), Q(x,y)</math> 为高阶非线性项。通过[[形式级数法]]构造
:<math>F(x,y) = x^2 + y^2 + F_3(x,y) + F_4(x,y) + \cdots</math>
使得沿系统轨线
:<math>\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>L_1, L_2, L_3, \ldots</math> 称为'''[[Lyapunov量]]'''（Lyapunov quantities）或'''[[焦点量]]'''（focal values）。

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

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

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

==='''特殊情况1——保守系统中的非线性中心'''===

==== 定义 ====
在经典力学与非线性动力学中，若一个二阶系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>（其中 <math>\mathbf{x} \in \mathbb{R}^2</math>）存在一个光滑的实值函数 <math>E(\mathbf{x})</math>，使得沿任意轨迹 <math>\mathbf{x}(t)</math> 均有 <math>\frac{dE}{dt} = 0</math>，且 <math>E(\mathbf{x})</math> 在任何开集上不为常数，<math>E(\mathbf{x})</math> 称为'''[[守恒量]]'''，称具有该守恒量的系统为'''[[保守系统]]'''<ref>STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017</ref>。

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

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

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

[[文件:保守系统的非线性中心示例图.png|缩略图|图5：双阱系统能量等高线图及其相平面投影|right|360px]]

==== 示例：双阱势系统 ====
如图5，考虑一个单位质量粒子在双阱势 <math>V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4</math> 中的无阻尼运动。运动方程为 <math>\ddot{x} = x - x^3</math>，等价系统为：
:<math>\dot{x} = y, \quad \dot{y} = x - x^3.</math>
系统具有守恒量（总能量）：
:<math>E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4.</math>
不动点为 <math>(0,0)</math> 和 <math>(\pm 1, 0)</math>。在 <math>(0,0)</math> 处，能量曲面的等高线呈现鞍点特征；而在 <math>(\pm 1, 0)</math> 处，<math>E(x, y)</math> 达到局部极小值。根据定理，这两个不动点均为非线性中心，附近存在一族闭轨（周期轨道）。

==='''特殊情况2——可逆系统中的非线性中心'''===
====定义====
如果存在一个相空间到自身的映射 <math>R(\mathbf{x})</math>，满足 <math>R^2(\mathbf{x}) = \mathbf{x}</math>（即 <math>R</math> 是一个[[对合变换]]），使得系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math> 在变换 <math>t \to -t</math>，<math>\mathbf{x} \to R(\mathbf{x})</math> 下保持不变，则称该系统为'''[[可逆系统]]'''。

最常见的二维情形是系统具有关于 <math>x</math> 轴的反射对称性：即变换 <math>t \to -t</math>，<math>y \to -y</math> 保持系统不变。这等价于要求方程 <math>\dot{x} = f(x,y)</math>，<math>\dot{y} = g(x,y)</math> 中的 <math>f</math> 关于 <math>y</math> 是奇函数（<math>f(x,-y) = -f(x,y)</math>），且 <math>g</math> 关于 <math>y</math> 是偶函数（<math>g(x,-y) = g(x,y)</math>）。

==== 定理（可逆系统的非线性中心）====
假设原点 <math>\mathbf{x}^* = \mathbf{0}</math> 是一个连续可微的可逆系统的'''线性中心'''。若系统是可逆的，则原点也是一个'''非线性中心'''，即所有充分接近原点的轨迹都是闭合的。

'''证明思路：''' 考虑一条从正 <math>x</math> 轴出发靠近原点的轨迹。由于线性部分是中心，该轨迹将绕原点旋转并最终与负 <math>x</math> 轴相交。利用可逆性，将轨迹关于 <math>x</math> 轴反射并反转时间方向，将得到一条具有相同端点但时间箭头反向的相异轨迹，这两条轨迹共同构成一个闭合曲线。

'''反证法：''' 若非线性扰动致使闭合轨迹转变为向内螺旋的焦点，则由可逆性推知其反射逆像必然形成向外螺旋的轨迹。此种轨迹行为与微分方程解的唯一性相冲突，因而在可逆系统中，线性中心在非线性扰动下仍保持为闭合轨迹族，即构成'''非线性中心'''。

该定理表明，[[可逆系统]]所具备的[[时间反演对称性|时间反演和空间反射对称性]]为非线性中心提供了独立于[[能量守恒]]的稳定性机制。

==== 示例：无阻尼钟摆 ====
如图6，无阻尼钟摆的非线性本质通过小角度的近似值<math>\sin\theta\approx\theta</math>被忽略。利用相图的分析钟摆旋转到顶部这一大角度的区域，在缺少阻力和外部吸引力的情况下，钟摆由下式决定：

[[文件:钟摆.png|缩略图|右|图6：无阻尼钟摆|right|280px]]

:<math>\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0.</math>
令 <math>\omega = \sqrt{g/L}</math>，<math>\tau = \omega t</math> 进行无量纲化，得：
:<math>\ddot{\theta} + \sin\theta = 0.</math>
写成相平面系统（<math>v = \dot{\theta}</math>）：
:<math>\dot{\theta} = v, \quad \dot{v} = -\sin\theta.</math>
该系统在变换 <math>\tau \to -\tau</math>，<math>v \to -v</math> 下保持不变，因此是可逆的。不动点 <math>(\theta^*, v^*) = (0,0)</math> 处的雅可比矩阵为 <math>\begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix}</math>，特征值为 <math>\pm i</math>，故为线性中心。根据可逆系统的非线性中心定理，<math>(0,0)</math> 也是一个非线性中心，对应钟摆在最低点的微小摆动。此外，系统也具有守恒量 <math>E = \frac{1}{2}v^2 - \cos\theta</math>，且 <math>(0,0)</math> 是 <math>E</math> 的局部极小值点，这从保守系统的角度也印证了该结论。

== '''参见''' ==

* [[非线性系统]]
* [[线性系统]]
* [[不动点]]
* [[极限环]]
* [[分岔理论]]
* [[中心流形定理]]

=='''参考资料'''==
<references />

== '''外部链接''' ==
* [http://www.scholarpedia.org/article/Equilibrium Scholarpedia: Equilibrium ]
* [http://www.scholarpedia.org/article/Dynamical_systems Scholarpedia: Dynamical Systems]
* [https://www.bohrium.com/sciencepedia/feynman/ordinary_differential_equations_undergraduate-the_Hartman-Grobman_theorem Bohrium-sciencepedia: the Hartman-Grobman theorem]

----
本中文词条由[[用户:王晔|王晔]]编辑，[[用户:陶如意|陶如意]]审校，欢迎在讨论页面留言。
本词条内容源自公开资料，遵守 CC3.0协议。
{{DEFAULTSORT:非线性中心}}

非线性中心

2026-07-12T07:39:49Z

王晔：优化了保守和可逆系统的非线性中心定理证明部分的表述

#重定向 [[非线性中心]]{{DISPLAYTITLE:非线性中心}}
'''非线性中心'''（'''Nonlinear Center'''）是[[非线性动力学]]中的基本结构之一。当系统[[奇点]]（平衡点）的局部邻域内布满环绕该点的嵌套闭合相轨迹时，该奇点即为非线性中心<ref>Strogatz, S. H. (2015). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering. CRC press.</ref>。由于一般的中心结构缺乏结构稳定性，极易在微弱的非线性扰动下退化为向内或向外螺旋的[[焦点]]，因此严格的非线性中心通常仅存在于受到特定物理法则约束的系统中，最典型的是具备能量守恒的[[保守系统]]或具有[[时间反演对称性]]的[[可逆系统]]。探讨非线性中心在扰动下如何[[分岔]]出[[极限环]]的研究（即'''[[庞加莱-李雅普诺夫中心问题]]'''），是揭示系统自激振荡机制与探索'''[[希尔伯特第十六问题]]'''的核心途径<ref>Ilyashenko, Y. (2002). Centennial history of Hilbert's 16th problem. Bulletin of the American Mathematical Society, 39(3), 301-354.</ref>。

=='''背景'''==
==='''线性中心'''===
在动力系统理论中，[[中心]]（Center）这一概念最初起源于对纯[[线性系统]]相空间轨线的定性分类。1881年，法国数学家[[庞加莱]]（Henri Poincaré）在创立微分方程定性理论时，首次根据系统矩阵的特征值性质，将二维平衡点划分为[[中心]]、[[焦点]]、[[结点]]与[[鞍点]]<ref>Poincaré, H. (1881). Mémoire sur les courbes définies par une équation différentielle. ''Journal de mathématiques pures et appliquées'', 375-422.</ref>，参见[[二维线性动力系统]]。按照该定义，当[[雅可比矩阵]]拥有一对纯虚数特征值 <math>\pm i\omega </math> 时，对应纯[[线性系统]]的平衡点即表现为[[线性中心]]。在相空间中，其轨线呈现为一族完美的同心圆或椭圆闭轨（见图2a），代表着一种无耗散的理想周期振荡状态<ref>Strogatz, S. H. (2018). ''Nonlinear dynamics and chaos''. CRC press.</ref>。

==='''线性化方法'''===
对于一般[[非线性动力系统]]，<math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>，研究[[不动点]]<math>\mathbf{x}^* </math>附近局部动力学行为的标准手段是[[线性化方法]]。该方法通过在不动点处进行[[泰勒级数展开]]并忽略高阶项，将原系统近似简化为[[线性系统]]：

<math>\dot{\mathbf{x}} \approx J(\mathbf{x}^*) (\mathbf{x} - \mathbf{x}^*)</math>

其中<math> J(\mathbf{x}^*) </math>为[[雅可比矩阵]]。该方法利用线性系统的[[特征值]]性质，定性判定原[[非线性系统]]在不动点邻域内的[[拓扑结构]]与[[稳定性]]。

==='''非线性系统的线性化失效'''===
根据[[Hartman-Grobman定理]]<ref>Perko, L. (2013). ''Differential Equations and Dynamical Systems''. Springer.</ref>，局部线性化系统与原非线性系统之间的[[拓扑同胚]]关系，仅在平衡点为[[双曲平衡点]]（特征值实部均非零）时成立。

而当系统处于非双曲的拓扑临界状态时（例如雅可比矩阵存在实部严格为零的特征值，即[[线性中心]]），被忽略的高阶非线性项将对系统的真实动力学行为起决定性作用，此时局部线性化方法失效，无法通过线性近似推断原系统的拓扑性质。

==='''中心的非线性扰动'''===
线性化失效情形下，理想的闭轨结构极度脆弱。高阶非线性项的介入通常会打破原有的局部守恒性，使闭轨演变为渐近吸引或排斥的[[焦点]]。仅在满足特定对称性或守恒律的条件下，系统才能抵御高阶项扰动，形成非线性中心。

为直观展示高阶非线性扰动对临界拓扑系统的决定性影响，可考察以下二维非线性动力系统：
:<math>\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>

[[文件:非线性中心扰动.png|缩略图|右|图1：非线性中心扰动|480px]]

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

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

:<math>\left\{ \begin{array}{l} x = r\cos\theta \\ y = r\sin\theta \end{array} \right.</math>

将极坐标代入原方程并化简，系统可严格转换为如下形式：
:<math>\left\{ \begin{array}{l} \dot{r} = ar^3 \\[6pt] \dot{\theta} = 1 \end{array} \right.</math>
基于极坐标方程，高阶非线性参数 <math>a</math> 对系统相图结构的决定性影响表现如下（见图1）：

* '''稳定焦点'''：当 <math> a < 0 </math> 时，对于任意 <math> r > 0</math>，径向演化 <math> \dot{r} < 0</math>。轨线随时间向内螺旋收缩，原有闭轨被高阶项破坏，原点退化为渐近稳定的[[焦点]]。
*
* '''非线性中心'''：当且仅当 <math> a = 0 </math> 时，系统无径向耗散 <math>（\dot{r} \equiv 0）</math>，角速度恒定。轨线在[[相空间]]中形成一族完美的同心圆闭轨，原点维持中心结构。
*
* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。
*
该经典实例表明，当动力系统处于纯虚数特征值的临界状态时，其相图拓扑极度脆弱，局部的真实稳定性由极微小的高阶非线性扰动（本例中的 <math>a r^3 </math> 项）所主导。

[[文件:非线性系统线性化失效典型情形.png|缩略图|图2：非线性系统线性化失效典型情形|480px]]

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

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

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

=='''非线性动力系统的非线性中心'''==
在前述的线性化失效与临界扰动背景下，要使系统在非线性范畴内维持循环运动，其不动点必须满足更严格的条件。
===定义与实例===
====定义====
在非线性系统 <math>\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x})</math> 的平衡点 <math>\mathbf{x}_0</math> 处，若线性化后的雅可比矩阵具有一对纯虚数特征值 <math>\pm i\omega</math>（即线性部分为“中心”），且加入高阶非线性项后，系统在该平衡点邻域内仍然存在一簇闭合轨道（周期解），则称该平衡点 <math>\mathbf{x}_0</math> 为'''非线性中心'''<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York.</ref>。

====实例====
[[文件:Lotka-Volterra 相图.png|缩略图|右|图3：Lotka-Volterra 模型的非线性中心|260px]]
为了直观展示这一解析定义的物理与生物学意义，我们可以引入经典的 '''[[Lotka-Volterra 捕食者-猎物模型]]'''。该模型用于描述生态系统中猎物 <math>x(t)</math> 与捕食者 <math>y(t)</math> 的相互作用，其非量纲化方程为<ref>Murray, J. D. (2002). Mathematical biology: I. An introduction (Vol. 17). Springer Science & Business Media.</ref>：

:<math>\left\{ \begin{array}{l} \dot{x} = x(1-y) \\[6pt] \dot{y} = y(x-1) \end{array} \right.</math>

该系统在正象限内存在平衡点 <math>(1, 1)</math>，其雅可比矩阵特征值为 <math>\lambda = \pm i</math>，满足线性中心的前提条件。进一步分析可知，该系统存在严格守恒的能量函数 <math>H(x, y) = (x - \ln x) + (y - \ln y)</math>。由于沿系统轨线 <math>\frac{dH}{dt} \equiv 0</math>，且 <math>H(x,y)</math> 在 <math>(1,1)</math> 处达到全局极小值，高阶非线性项无法破坏轨线的闭合性。因此，该平衡点被确认为一个严格的'''非线性中心'''。在相图上，这表现为捕食者与猎物的种群数量不会趋于静态，而是围绕平衡点呈现出持续的、由初始状态决定的周期性闭轨波动（见图3）。

===性质===
除上述解析定义与直观模型外，非线性中心在动力系统中还具备以下几个核心的几何与拓扑特征：

==== 几何连续性与极限环的区别 ====
[[文件:非线性中心与极限环对比图.png|缩略图|右|图4：非线性中心与极限环对比图|400px]]
非线性中心的核心几何特征在于其闭轨的连续分布性。在中心的一个充分小的邻域内，除平衡点本身外，所有的相轨线均是嵌套的闭合周期轨线，构成一个[[连续统]]（Continuum）。

这一特征将其与[[极限环]]（Limit Cycle）区分开来（见图4）：极限环同为[[闭合周期轨线]]，但它是孤立的，即在极限环的充分小邻域内不存在其他闭合轨线，其内侧或外侧的相邻轨线会以螺旋形式渐近地趋向（或远离）该极限环。

==== 拓扑特征与庞加莱指数 ====
在[[向量场拓扑学]]中，奇点的局部结构可通过[[庞加莱指数]]（Poincaré Index）进行定量刻画<ref>Perko, L. (2013). Differential Equations and Dynamical Systems (3rd ed.). Springer Science & Business Media. ISBN 978-1-4613-0003-8.</ref>。对于[[二维非线性系统]]，沿包围非线性中心且不包含其他奇点的任意简单闭曲线，其向量场旋转的圈数恒为 1。因此，非线性中心的[[庞加莱指数]]为 <math>+1</math>。这一拓扑不变量与[[结点]]、[[焦点]]的指数相同，但与指数为 <math>-1</math> 的[[鞍点]]具有明确的拓扑差异。

==== 周期与振幅的非线性依赖 ====
在纯[[线性系统]]中，中心周围所有闭轨的运转[[周期]]恒为常数 <math>T = 2\pi/\omega</math>，与轨线的振幅大小无关。对于一般的非线性中心，由于高阶非线性项的影响，系统沿闭轨运动的周期 <math>T(r)</math> 通常会依赖于轨线的振幅（即极径 <math>r</math>）。在特殊情况下，若非线性系统不仅存在中心，且周围所有闭轨的周期依然保持恒定，则该奇点被称为“'''[[等时中心]]'''”（Isochronous Center）。等时中心的判定是[[常微分方程定性理论]]中的一项具体研究课题<ref>Chavarriga, J., & Sabatini, M. (1999). A survey of isochronous centers. Qualitative Theory of Dynamical Systems, 1(1), 1-70. </ref>。

==== 高维系统的中心流形推广 ====
对于相空间维度 <math>n \ge 3</math> 的非线性系统，中心的概念可通过'''[[中心流形定理]]'''（Center Manifold Theorem）进行高维推广。若系统在[[不动点]]处的[[雅可比矩阵]]存在一对纯虚数特征值，且其余特征值的实部均不为零，则系统在[[不动点]]附近存在一个二维的[[中心流形]]。高维系统的局部渐近动力学可降维至该[[流形]]上进行分析。若降维后的二维子系统呈现非线性中心特征，则高维系统在对应子空间内同样表现为局部的周期振荡行为<ref>Kuznetsov, Yu. A. (2004). Elements of Applied Bifurcation Theory (3rd edition). Springer.</ref>。

=='''非线性中心判据'''==

'''[[经典中心问题]]'''（Classical Center Problem）由 Poincaré 正式提出<ref>POINCARÉ H. Mémoire sur les courbes définies par une équation différentielle (I)[J]. Journal de Mathématiques Pures et Appliquées, 1881, 7(3): 375-422.</ref>：

: 给定一个[[平面解析系统]]，其线性部分为[[中心]]，如何通过系统的系数判定该平衡点究竟是中心还是[[焦点]]？

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

解决经典中心问题最直接的方法是逐阶计算系统的[[Lyapunov量]]<ref>LYAPUNOV A M. Problème général de la stabilité du mouvement[M]. Princeton: Princeton University Press, 1947.</ref>，通过其消失条件判定中心是否存在。

==='''一般情况的判定——Lyapunov量'''===
====定义====
考虑平面系统
:<math>\left\{ \begin{array}{l} \dot{x} = -y + P(x,y) \\[6pt] \dot{y} = x + Q(x,y) \end{array} \right.</math>
其中 <math>P(x,y), Q(x,y)</math> 为高阶非线性项。通过[[形式级数法]]构造
:<math>F(x,y) = x^2 + y^2 + F_3(x,y) + F_4(x,y) + \cdots</math>
使得沿系统轨线
:<math>\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>L_1, L_2, L_3, \ldots</math> 称为'''[[Lyapunov量]]'''（Lyapunov quantities）或'''[[焦点量]]'''（focal values）。

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

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

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

==='''特殊情况1——保守系统中的非线性中心'''===

==== 定义 ====
在经典力学与非线性动力学中，若一个二阶系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>（其中 <math>\mathbf{x} \in \mathbb{R}^2</math>）存在一个光滑的实值函数 <math>E(\mathbf{x})</math>，使得沿任意轨迹 <math>\mathbf{x}(t)</math> 均有 <math>\frac{dE}{dt} = 0</math>，且 <math>E(\mathbf{x})</math> 在任何开集上不为常数，<math>E(\mathbf{x})</math> 称为'''[[守恒量]]'''，称具有该守恒量的系统为'''[[保守系统]]'''<ref>STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017</ref>。

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

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

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

[[文件:保守系统的非线性中心示例图.png|缩略图|图5：双阱系统能量等高线图及其相平面投影|right|360px]]

==== 示例：双阱势系统 ====
如图5，考虑一个单位质量粒子在双阱势 <math>V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4</math> 中的无阻尼运动。运动方程为 <math>\ddot{x} = x - x^3</math>，等价系统为：
:<math>\dot{x} = y, \quad \dot{y} = x - x^3.</math>
系统具有守恒量（总能量）：
:<math>E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4.</math>
不动点为 <math>(0,0)</math> 和 <math>(\pm 1, 0)</math>。在 <math>(0,0)</math> 处，能量曲面的等高线呈现鞍点特征；而在 <math>(\pm 1, 0)</math> 处，<math>E(x, y)</math> 达到局部极小值。根据定理，这两个不动点均为非线性中心，附近存在一族闭轨（周期轨道）。

==='''特殊情况2——可逆系统中的非线性中心'''===
====定义====
如果存在一个相空间到自身的映射 <math>R(\mathbf{x})</math>，满足 <math>R^2(\mathbf{x}) = \mathbf{x}</math>（即 <math>R</math> 是一个[[对合变换]]），使得系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math> 在变换 <math>t \to -t</math>，<math>\mathbf{x} \to R(\mathbf{x})</math> 下保持不变，则称该系统为'''[[可逆系统]]'''。

最常见的二维情形是系统具有关于 <math>x</math> 轴的反射对称性：即变换 <math>t \to -t</math>，<math>y \to -y</math> 保持系统不变。这等价于要求方程 <math>\dot{x} = f(x,y)</math>，<math>\dot{y} = g(x,y)</math> 中的 <math>f</math> 关于 <math>y</math> 是奇函数（<math>f(x,-y) = -f(x,y)</math>），且 <math>g</math> 关于 <math>y</math> 是偶函数（<math>g(x,-y) = g(x,y)</math>）。

==== 定理（可逆系统的非线性中心）====
假设原点 <math>\mathbf{x}^* = \mathbf{0}</math> 是一个连续可微的可逆系统的'''线性中心'''。若系统是可逆的，则原点也是一个'''非线性中心'''，即所有充分接近原点的轨迹都是闭合的。

'''证明思路：''' 考虑一条从正 <math>x</math> 轴出发靠近原点的轨迹。由于线性部分是中心，该轨迹将绕原点旋转并最终与负 <math>x</math> 轴相交。利用可逆性，将轨迹关于 <math>x</math> 轴反射并反转时间方向，将得到一条具有相同端点但时间箭头反向的相异轨迹，这两条轨迹共同构成一个闭合曲线。

'''反证法：''' 若非线性扰动致使闭合轨迹转变为向内螺旋的焦点，则由可逆性推知其反射逆像必然形成向外螺旋的轨迹。此种轨迹行为与微分方程解的唯一性相冲突，因而在可逆系统中，线性中心在非线性扰动下仍保持为闭合轨迹族，即构成'''非线性中心'''。

该定理表明，[[可逆系统]]所具备的[[时间反演对称性|时间反演和空间反射对称性]]为非线性中心提供了独立于[[能量守恒]]的稳定性机制。

==== 示例：无阻尼钟摆 ====
如图6，无阻尼钟摆的非线性本质通过小角度的近似值<math>\sin\theta\approx\theta</math>被忽略。利用相图的分析钟摆旋转到顶部这一大角度的区域，在缺少阻力和外部吸引力的情况下，钟摆由下式决定：

[[文件:钟摆.png|缩略图|右|图6：无阻尼钟摆|right|280px]]

:<math>\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0.</math>
令 <math>\omega = \sqrt{g/L}</math>，<math>\tau = \omega t</math> 进行无量纲化，得：
:<math>\ddot{\theta} + \sin\theta = 0.</math>
写成相平面系统（<math>v = \dot{\theta}</math>）：
:<math>\dot{\theta} = v, \quad \dot{v} = -\sin\theta.</math>
该系统在变换 <math>\tau \to -\tau</math>，<math>v \to -v</math> 下保持不变，因此是可逆的。不动点 <math>(\theta^*, v^*) = (0,0)</math> 处的雅可比矩阵为 <math>\begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix}</math>，特征值为 <math>\pm i</math>，故为线性中心。根据可逆系统的非线性中心定理，<math>(0,0)</math> 也是一个非线性中心，对应钟摆在最低点的微小摆动。此外，系统也具有守恒量 <math>E = \frac{1}{2}v^2 - \cos\theta</math>，且 <math>(0,0)</math> 是 <math>E</math> 的局部极小值点，这从保守系统的角度也印证了该结论。

== '''参见''' ==

* [[非线性系统]]
* [[线性系统]]
* [[不动点]]
* [[极限环]]
* [[分岔理论]]
* [[中心流形定理]]

=='''参考资料'''==
<references />

== '''外部链接''' ==
* [http://www.scholarpedia.org/article/Equilibrium Scholarpedia: Equilibrium ]
* [http://www.scholarpedia.org/article/Dynamical_systems Scholarpedia: Dynamical Systems]
* [https://www.bohrium.com/sciencepedia/feynman/ordinary_differential_equations_undergraduate-the_Hartman-Grobman_theorem Bohrium-sciencepedia: the Hartman-Grobman theorem]

----
本中文词条由[[用户:王晔|王晔]]编辑，[[用户:陶如意|陶如意]]审校，欢迎在讨论页面留言。
本词条内容源自公开资料，遵守 CC3.0协议。
{{DEFAULTSORT:非线性中心}}

2026-06-08T12:01:46Z

王晔：王晔上传文件:非线性中心扰动.png的新版本

非线性中心

2026-06-08T11:56:02Z

王晔：/* 几何连续性与极限环的区别 */

#重定向 [[非线性中心]]{{DISPLAYTITLE:非线性中心}}
'''非线性中心'''（'''Nonlinear Center'''）是[[非线性动力学]]中的基本结构之一。当系统[[奇点]]（平衡点）的局部邻域内布满环绕该点的嵌套闭合相轨迹时，该奇点即为非线性中心<ref>Strogatz, S. H. (2015). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering. CRC press.</ref>。由于一般的中心结构缺乏结构稳定性，极易在微弱的非线性扰动下退化为向内或向外螺旋的[[焦点]]<ref>Andronov, A. A., & Pontryagin, L. (1937). Systèmes grossiers. In Dokl. Akad. Nauk. SSSR (Vol. 14, No. 5, pp. 247-250).</ref>，因此严格的非线性中心通常仅存在于受到特定物理法则约束的系统中，最典型的是具备能量守恒的[[保守系统]]或具有[[时间反演对称性]]的[[可逆系统]]<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear oscillations, dynamical systems, and bifurcations of vector fields. Springer Science & Business Media.</ref>。探讨非线性中心在扰动下如何[[分岔]]出[[极限环]]的研究（即'''[[庞加莱-李雅普诺夫中心问题]]'''），是揭示系统自激振荡机制与探索'''[[希尔伯特第十六问题]]'''的核心途径<ref>Ilyashenko, Y. (2002). Centennial history of Hilbert's 16th problem. Bulletin of the American Mathematical Society, 39(3), 301-354.</ref>。

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

==='''非线性系统的线性化失效'''===
对于一般的非[[线性系统]] <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>，[[不动点]] <math>\mathbf{x}^*</math> 附近的局部动力学行为通常可通过线性化系统来近似判定。根据'''[[哈特曼-格罗伯曼定理]]'''（'''Hartman-Grobman Theorem'''），仅当平衡点[[特征值]]实部非零（[[双曲平衡点]]）时，即它在任何方向上都有明确的吸引或排斥倾向，非线性系统的局部拓扑才与其线性化系统[[同胚]]<ref>Hartman,P.(1960)."On the local linearization of differential equations." Proc.AMS.</ref>。对于特征值为纯虚数的临界情形，线性化方法失效。此时，平衡点附近的真实拓扑结构取决于[[高阶非线性项]]。

==='''中心的非线性扰动'''===
为直观展示高阶非线性项如何破坏线性中心，考虑经典的二维非线性[[动力系统]]：
:<math>\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>

[[文件:非线性中心扰动.png|缩略图|右|图1：非线性中心扰动|480px]]

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

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

:<math>\left\{ \begin{array}{l} x = r\cos\theta \\ y = r\sin\theta \end{array} \right.</math>

将极坐标代入原方程并化简，系统可严格转换为如下形式：
:<math>\left\{ \begin{array}{l} \dot{r} = ar^3 \\[6pt] \dot{\theta} = 1 \end{array} \right.</math>
基于极坐标方程，高阶非线性参数 <math>a</math> 对系统相图结构的决定性影响表现如下（见图1）：

* '''稳定焦点'''：当 <math> a < 0 </math> 时，对于任意 <math> r > 0</math>，径向演化 <math> \dot{r} < 0</math>。轨线随时间向内螺旋收缩，原有闭轨被高阶项破坏，原点退化为渐近稳定的[[焦点]]。
*
* '''非线性中心'''：当且仅当 <math> a = 0 </math> 时，系统无径向耗散 <math>（\dot{r} \equiv 0）</math>，角速度恒定。轨线在[[相空间]]中形成一族完美的同心圆闭轨，原点维持中心结构。
*
* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。
*
该经典实例表明，当动力系统处于纯虚数特征值的临界状态时，其相图拓扑极度脆弱，局部的真实稳定性完全由极微小的高阶非线性扰动（本例中的 <math>a r^3 </math> 项）所主导。

[[文件:非线性系统线性化失效典型情形.png|缩略图|图2：非线性系统线性化失效典型情形]]

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

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

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

=='''非线性动力系统的非线性中心'''==
在前述的线性化失效与临界扰动背景下，要使系统在非线性范畴内维持循环运动，其不动点必须满足更严格的条件。在非线性系统 <math>\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x})</math> 的平衡点 <math>\mathbf{x}_0</math> 处，若线性化后的雅可比矩阵具有一对纯虚数特征值 <math>\pm i\omega</math>（即线性部分为“中心”），且加入高阶非线性项后，系统在该平衡点邻域内仍然存在一簇闭合轨道（周期解），则称该平衡点 <math>\mathbf{x}_0</math> 为'''非线性中心'''<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York.</ref>。

为了直观展示这一解析定义的物理与生物学意义，我们可以引入经典的 '''Lotka-Volterra 捕食者-猎物模型'''。该模型用于描述生态系统中猎物 <math>x(t)</math> 与捕食者 <math>y(t)</math> 的相互作用，其非量纲化方程为<ref>Murray, J. D. (2002). Mathematical biology: I. An introduction (Vol. 17). Springer Science & Business Media.</ref>：

[[文件:Lotka-Volterra 相图.png|缩略图|右|图3：Lotka-Volterra 模型的非线性中心。相图显示猎物（x）与捕食者（y）数量围绕平衡点 (1,1) 呈周期性循环运动。|320px]]

:<math>\left\{ \begin{array}{l} \dot{x} = x(1-y) \\[6pt] \dot{y} = y(x-1) \end{array} \right.</math>

该系统在正象限内存在平衡点 <math>(1, 1)</math>，其雅可比矩阵特征值为 <math>\lambda = \pm i</math>，满足线性中心的前提条件。进一步分析可知，该系统存在严格守恒的能量函数 <math>H(x, y) = (x - \ln x) + (y - \ln y)</math>。由于沿系统轨线 <math>\frac{dH}{dt} \equiv 0</math>，且 <math>H(x,y)</math> 在 <math>(1,1)</math> 处达到全局极小值，高阶非线性项无法破坏轨线的闭合性。因此，该平衡点被确认为一个严格的'''非线性中心'''。在相图上，这表现为捕食者与猎物的种群数量不会趋于静态，而是围绕平衡点呈现出持续的、由初始状态决定的周期性闭轨波动。

除上述解析定义与直观模型外，非线性中心在动力系统中还具备以下几个核心的几何与拓扑特征：

=== 几何连续性与极限环的区别 ===
非线性中心的核心几何特征在于其闭轨的连续分布性。在中心的一个充分小的邻域内，除平衡点本身外，所有的相轨线均是嵌套的闭合周期轨线，构成一个连续统（Continuum）。这一特征将其与[[极限环]]（Limit Cycle）区分开来：极限环同为闭合周期轨线，但它是孤立的，即在极限环的充分小邻域内不存在其他闭合轨线。中心结构通常代表系统处于无耗散的局部保守状态，而极限环多出现于耗散系统中，是能量输入与耗散达到动态平衡的结果。
[[文件:非线性中心与极限环对比图.png|缩略图|右|非线性中心与极限环对比图]]

=== 拓扑特征与庞加莱指数 ===
在向量场拓扑学中，奇点的局部结构可通过庞加莱指数（Poincaré Index）进行定量刻画。对于二维非线性系统，沿包围非线性中心且不包含其他奇点的任意简单闭曲线，其向量场旋转的圈数恒为 1。因此，非线性中心的庞加莱指数为 <math>+1</math>。这一拓扑不变量与[[结点]]、[[焦点]]的指数相同，但与指数为 <math>-1</math> 的[[鞍点]]具有明确的拓扑差异。

=== 周期与振幅的非线性依赖 ===
在纯线性系统中，中心周围所有闭轨的运转周期恒为常数 <math>T = 2\pi/\omega</math>，与轨线的振幅大小无关。对于一般的非线性中心，由于高阶非线性项的影响，系统沿闭轨运动的周期 <math>T(r)</math> 通常会依赖于轨线的振幅（即极径 <math>r</math>）。在特殊情况下，若非线性系统不仅存在中心，且周围所有闭轨的周期依然保持恒定，则该奇点被称为“'''等时中心'''”（Isochronous Center）。等时中心的判定是常微分方程定性理论中的一项具体研究课题<ref>Chicone, C. (1999). Ordinary Differential Equations with Applications. Springer Science & Business Media.</ref>。

=== 高维系统的中心流形推广 ===
对于相空间维度 <math>n \ge 3</math> 的非线性系统，中心的概念可通过'''[[中心流形定理]]'''（Center Manifold Theorem）进行高维推广。若系统在不动点处的雅可比矩阵存在一对纯虚数特征值，且其余特征值的实部均不为零，则系统在不动点附近存在一个二维的中心流形。高维系统的局部渐近动力学可降维至该流形上进行分析。若降维后的二维子系统呈现非线性中心特征，则高维系统在对应子空间内同样表现为局部的周期振荡行为<ref>Kuznetsov, Yu. A. (2004). Elements of Applied Bifurcation Theory (3rd edition). Springer.</ref>。

=='''非线性中心判据'''==

'''经典中心问题'''（Classical Center Problem）由 Poincaré 正式提出<ref>POINCARÉ H. Mémoire sur les courbes définies par une équation différentielle (I)[J]. Journal de Mathématiques Pures et Appliquées, 1881, 7(3): 375-422.</ref>：

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

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

解决经典中心问题最直接的方法是逐阶计算系统的 Lyapunov 量<ref>LYAPUNOV A M. Problème général de la stabilité du mouvement[M]. Princeton: Princeton University Press, 1947.</ref>，通过其消失条件判定中心是否存在。

===一般情况的判定——Lyapunov量===
====定义====
考虑平面系统
:<math>\left\{ \begin{array}{l} \dot{x} = -y + P(x,y) \\[6pt] \dot{y} = x + Q(x,y) \end{array} \right.</math>
其中 <math>P(x,y), Q(x,y)</math> 为高阶非线性项。通过形式级数法构造
:<math>F(x,y) = x^2 + y^2 + F_3(x,y) + F_4(x,y) + \cdots</math>
使得沿系统轨线
:<math>\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>L_1, L_2, L_3, \ldots</math> 称为'''Lyapunov量'''（Lyapunov quantities）或'''焦点量'''（focal values）。

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

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

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

===特殊情况1——保守系统中的非线性中心===

==== 定义 ====
在经典力学与非线性动力学中，若一个二阶系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>（其中 <math>\mathbf{x} \in \mathbb{R}^2</math>）存在一个光滑的实值函数 <math>E(\mathbf{x})</math>，使得沿任意轨迹 <math>\mathbf{x}(t)</math> 均有 <math>\frac{dE}{dt} = 0</math>，且 <math>E(\mathbf{x})</math> 在任何开集上不为常数，<math>E(\mathbf{x})</math> 称为'''守恒量'''，称具有该守恒量的系统为'''保守系统'''<ref>STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017</ref>。

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

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

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

[[文件:保守系统的非线性中心示例图.png|缩略图|图3：双阱系统能量等高线图及其相平面投影|right|360px]]

==== 示例：双阱势系统 ====
考虑一个单位质量粒子在双阱势 <math>V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4</math> 中的无阻尼运动。运动方程为 <math>\ddot{x} = x - x^3</math>，等价系统为：
:<math>\dot{x} = y, \quad \dot{y} = x - x^3.</math>
系统具有守恒量（总能量）：
:<math>E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4.</math>
不动点为 <math>(0,0)</math> 和 <math>(\pm 1, 0)</math>。在 <math>(0,0)</math> 处，能量曲面的等高线呈现鞍点特征；而在 <math>(\pm 1, 0)</math> 处，<math>E(x, y)</math> 达到局部极小值。根据定理，这两个不动点均为非线性中心，附近存在一族闭轨（周期轨道）。

===特殊情况2——可逆系统中的非线性中心===
====定义====
如果存在一个相空间到自身的映射 <math>R(\mathbf{x})</math>，满足 <math>R^2(\mathbf{x}) = \mathbf{x}</math>（即 <math>R</math> 是一个对合变换），使得系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math> 在变换 <math>t \to -t</math>，<math>\mathbf{x} \to R(\mathbf{x})</math> 下保持不变，则称该系统为'''可逆系统'''。

最常见的二维情形是系统具有关于 <math>x</math> 轴的反射对称性：即变换 <math>t \to -t</math>，<math>y \to -y</math> 保持系统不变。这等价于要求方程 <math>\dot{x} = f(x,y)</math>，<math>\dot{y} = g(x,y)</math> 中的 <math>f</math> 关于 <math>y</math> 是奇函数（<math>f(x,-y) = -f(x,y)</math>），且 <math>g</math> 关于 <math>y</math> 是偶函数（<math>g(x,-y) = g(x,y)</math>）。

==== 定理（可逆系统的非线性中心）====
假设原点 <math>\mathbf{x}^* = \mathbf{0}</math> 是一个连续可微可逆系统的'''线性中心'''。若系统是可逆的，则原点也是一个'''非线性中心'''，即所有充分接近原点的轨迹都是闭合的。

'''证明思路：''' 考虑一条从正 <math>x</math> 轴出发靠近原点的轨迹。由于线性部分是中心，该轨迹将绕原点旋转并最终与负 <math>x</math> 轴相交。利用可逆性，将轨迹关于 <math>x</math> 轴反射并反转时间方向，将得到一条具有相同端点但箭头相反的轨迹。这两条轨迹共同构成一个闭合曲线。因此原点被一族闭轨包围。

'''判定逻辑：'''
该定理表明，可逆系统所具备的时间反演对称性为非线性中心提供了独立于能量守恒的稳定性机制。即便系统不满足保守条件，只要其在时间反演与空间反射的复合变换下保持不变，线性中心便具有抵御非线性扰动的能力。其几何解释为：若非线性扰动致使闭合轨迹转变为向内螺旋的焦点，则由可逆性推知其反射逆像必然形成向外螺旋的轨迹。此种轨迹行为与微分方程解的唯一性相冲突，因而在可逆系统中，线性中心在非线性扰动下仍保持为闭合轨迹族，即构成'''非线性中心'''。

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

[[文件:钟摆.png|缩略图|右|图4：无阻尼钟摆|right|280px]]

:<math>\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0.</math>
令 <math>\omega = \sqrt{g/L}</math>，<math>\tau = \omega t</math> 进行无量纲化，得：
:<math>\ddot{\theta} + \sin\theta = 0.</math>
写成相平面系统（<math>v = \dot{\theta}</math>）：
:<math>\dot{\theta} = v, \quad \dot{v} = -\sin\theta.</math>
该系统在变换 <math>\tau \to -\tau</math>，<math>v \to -v</math> 下保持不变，因此是可逆的。不动点 <math>(\theta^*, v^*) = (0,0)</math> 处的雅可比矩阵为 <math>\begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix}</math>，特征值为 <math>\pm i</math>，故为线性中心。根据可逆系统的非线性中心定理，<math>(0,0)</math> 也是一个非线性中心，对应钟摆在最低点的微小摆动。此外，系统也具有守恒量 <math>E = \frac{1}{2}v^2 - \cos\theta</math>，且 <math>(0,0)</math> 是 <math>E</math> 的局部极小值点，这从保守系统的角度也印证了该结论。

== 参见 ==

* [[非线性系统]]
* [[线性系统]]
* [[不动点]]
* [[极限环]]
* [[分岔理论]]
* [[中心流形定理]]

=='''参考资料'''==
<references />

== 外部链接 ==
* [http://www.scholarpedia.org/article/Equilibrium Scholarpedia: Equilibrium ]
* [http://www.scholarpedia.org/article/Dynamical_systems Scholarpedia: Dynamical Systems]
* [https://www.bohrium.com/sciencepedia/feynman/ordinary_differential_equations_undergraduate-the_Hartman-Grobman_theorem Bohrium-sciencepedia: the Hartman-Grobman theorem]

----
本中文词条由[[用户:王晔|王晔]]编辑，[[用户:陶如意|陶如意]]审校，欢迎在讨论页面留言。
本词条内容源自公开资料，遵守 CC3.0协议。
{{DEFAULTSORT:非线性中心}}

非线性中心

2026-06-08T11:54:57Z

王晔：

#重定向 [[非线性中心]]{{DISPLAYTITLE:非线性中心}}
'''非线性中心'''（'''Nonlinear Center'''）是[[非线性动力学]]中的基本结构之一。当系统[[奇点]]（平衡点）的局部邻域内布满环绕该点的嵌套闭合相轨迹时，该奇点即为非线性中心<ref>Strogatz, S. H. (2015). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering. CRC press.</ref>。由于一般的中心结构缺乏结构稳定性，极易在微弱的非线性扰动下退化为向内或向外螺旋的[[焦点]]<ref>Andronov, A. A., & Pontryagin, L. (1937). Systèmes grossiers. In Dokl. Akad. Nauk. SSSR (Vol. 14, No. 5, pp. 247-250).</ref>，因此严格的非线性中心通常仅存在于受到特定物理法则约束的系统中，最典型的是具备能量守恒的[[保守系统]]或具有[[时间反演对称性]]的[[可逆系统]]<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear oscillations, dynamical systems, and bifurcations of vector fields. Springer Science & Business Media.</ref>。探讨非线性中心在扰动下如何[[分岔]]出[[极限环]]的研究（即'''[[庞加莱-李雅普诺夫中心问题]]'''），是揭示系统自激振荡机制与探索'''[[希尔伯特第十六问题]]'''的核心途径<ref>Ilyashenko, Y. (2002). Centennial history of Hilbert's 16th problem. Bulletin of the American Mathematical Society, 39(3), 301-354.</ref>。

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

==='''非线性系统的线性化失效'''===
对于一般的非[[线性系统]] <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>，[[不动点]] <math>\mathbf{x}^*</math> 附近的局部动力学行为通常可通过线性化系统来近似判定。根据'''[[哈特曼-格罗伯曼定理]]'''（'''Hartman-Grobman Theorem'''），仅当平衡点[[特征值]]实部非零（[[双曲平衡点]]）时，即它在任何方向上都有明确的吸引或排斥倾向，非线性系统的局部拓扑才与其线性化系统[[同胚]]<ref>Hartman,P.(1960)."On the local linearization of differential equations." Proc.AMS.</ref>。对于特征值为纯虚数的临界情形，线性化方法失效。此时，平衡点附近的真实拓扑结构取决于[[高阶非线性项]]。

==='''中心的非线性扰动'''===
为直观展示高阶非线性项如何破坏线性中心，考虑经典的二维非线性[[动力系统]]：
:<math>\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>

[[文件:非线性中心扰动.png|缩略图|右|图1：非线性中心扰动|480px]]

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

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

:<math>\left\{ \begin{array}{l} x = r\cos\theta \\ y = r\sin\theta \end{array} \right.</math>

将极坐标代入原方程并化简，系统可严格转换为如下形式：
:<math>\left\{ \begin{array}{l} \dot{r} = ar^3 \\[6pt] \dot{\theta} = 1 \end{array} \right.</math>
基于极坐标方程，高阶非线性参数 <math>a</math> 对系统相图结构的决定性影响表现如下（见图1）：

* '''稳定焦点'''：当 <math> a < 0 </math> 时，对于任意 <math> r > 0</math>，径向演化 <math> \dot{r} < 0</math>。轨线随时间向内螺旋收缩，原有闭轨被高阶项破坏，原点退化为渐近稳定的[[焦点]]。
*
* '''非线性中心'''：当且仅当 <math> a = 0 </math> 时，系统无径向耗散 <math>（\dot{r} \equiv 0）</math>，角速度恒定。轨线在[[相空间]]中形成一族完美的同心圆闭轨，原点维持中心结构。
*
* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。
*
该经典实例表明，当动力系统处于纯虚数特征值的临界状态时，其相图拓扑极度脆弱，局部的真实稳定性完全由极微小的高阶非线性扰动（本例中的 <math>a r^3 </math> 项）所主导。

[[文件:非线性系统线性化失效典型情形.png|缩略图|图2：非线性系统线性化失效典型情形]]

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

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

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

=='''非线性动力系统的非线性中心'''==
在前述的线性化失效与临界扰动背景下，要使系统在非线性范畴内维持循环运动，其不动点必须满足更严格的条件。在非线性系统 <math>\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x})</math> 的平衡点 <math>\mathbf{x}_0</math> 处，若线性化后的雅可比矩阵具有一对纯虚数特征值 <math>\pm i\omega</math>（即线性部分为“中心”），且加入高阶非线性项后，系统在该平衡点邻域内仍然存在一簇闭合轨道（周期解），则称该平衡点 <math>\mathbf{x}_0</math> 为'''非线性中心'''<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York.</ref>。

为了直观展示这一解析定义的物理与生物学意义，我们可以引入经典的 '''Lotka-Volterra 捕食者-猎物模型'''。该模型用于描述生态系统中猎物 <math>x(t)</math> 与捕食者 <math>y(t)</math> 的相互作用，其非量纲化方程为<ref>Murray, J. D. (2002). Mathematical biology: I. An introduction (Vol. 17). Springer Science & Business Media.</ref>：

[[文件:Lotka-Volterra 相图.png|缩略图|右|图3：Lotka-Volterra 模型的非线性中心。相图显示猎物（x）与捕食者（y）数量围绕平衡点 (1,1) 呈周期性循环运动。|320px]]

:<math>\left\{ \begin{array}{l} \dot{x} = x(1-y) \\[6pt] \dot{y} = y(x-1) \end{array} \right.</math>

该系统在正象限内存在平衡点 <math>(1, 1)</math>，其雅可比矩阵特征值为 <math>\lambda = \pm i</math>，满足线性中心的前提条件。进一步分析可知，该系统存在严格守恒的能量函数 <math>H(x, y) = (x - \ln x) + (y - \ln y)</math>。由于沿系统轨线 <math>\frac{dH}{dt} \equiv 0</math>，且 <math>H(x,y)</math> 在 <math>(1,1)</math> 处达到全局极小值，高阶非线性项无法破坏轨线的闭合性。因此，该平衡点被确认为一个严格的'''非线性中心'''。在相图上，这表现为捕食者与猎物的种群数量不会趋于静态，而是围绕平衡点呈现出持续的、由初始状态决定的周期性闭轨波动。

除上述解析定义与直观模型外，非线性中心在动力系统中还具备以下几个核心的几何与拓扑特征：

=== 几何连续性与极限环的区别 ===
非线性中心的核心几何特征在于其闭轨的连续分布性。在中心的一个充分小的邻域内，除平衡点本身外，所有的相轨线均是嵌套的闭合周期轨线，构成一个连续统（Continuum）。这一特征将其与[[极限环]]（Limit Cycle）区分开来：极限环同为闭合周期轨线，但它是孤立的，即在极限环的充分小邻域内不存在其他闭合轨线。中心结构通常代表系统处于无耗散的局部保守状态，而极限环多出现于耗散系统中，是能量输入与耗散达到动态平衡的结果。
[[文件:非线性中心与极限环对比图|缩略图|右|非线性中心与极限环对比图]]

=== 拓扑特征与庞加莱指数 ===
在向量场拓扑学中，奇点的局部结构可通过庞加莱指数（Poincaré Index）进行定量刻画。对于二维非线性系统，沿包围非线性中心且不包含其他奇点的任意简单闭曲线，其向量场旋转的圈数恒为 1。因此，非线性中心的庞加莱指数为 <math>+1</math>。这一拓扑不变量与[[结点]]、[[焦点]]的指数相同，但与指数为 <math>-1</math> 的[[鞍点]]具有明确的拓扑差异。

=== 周期与振幅的非线性依赖 ===
在纯线性系统中，中心周围所有闭轨的运转周期恒为常数 <math>T = 2\pi/\omega</math>，与轨线的振幅大小无关。对于一般的非线性中心，由于高阶非线性项的影响，系统沿闭轨运动的周期 <math>T(r)</math> 通常会依赖于轨线的振幅（即极径 <math>r</math>）。在特殊情况下，若非线性系统不仅存在中心，且周围所有闭轨的周期依然保持恒定，则该奇点被称为“'''等时中心'''”（Isochronous Center）。等时中心的判定是常微分方程定性理论中的一项具体研究课题<ref>Chicone, C. (1999). Ordinary Differential Equations with Applications. Springer Science & Business Media.</ref>。

=== 高维系统的中心流形推广 ===
对于相空间维度 <math>n \ge 3</math> 的非线性系统，中心的概念可通过'''[[中心流形定理]]'''（Center Manifold Theorem）进行高维推广。若系统在不动点处的雅可比矩阵存在一对纯虚数特征值，且其余特征值的实部均不为零，则系统在不动点附近存在一个二维的中心流形。高维系统的局部渐近动力学可降维至该流形上进行分析。若降维后的二维子系统呈现非线性中心特征，则高维系统在对应子空间内同样表现为局部的周期振荡行为<ref>Kuznetsov, Yu. A. (2004). Elements of Applied Bifurcation Theory (3rd edition). Springer.</ref>。

=='''非线性中心判据'''==

'''经典中心问题'''（Classical Center Problem）由 Poincaré 正式提出<ref>POINCARÉ H. Mémoire sur les courbes définies par une équation différentielle (I)[J]. Journal de Mathématiques Pures et Appliquées, 1881, 7(3): 375-422.</ref>：

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

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

解决经典中心问题最直接的方法是逐阶计算系统的 Lyapunov 量<ref>LYAPUNOV A M. Problème général de la stabilité du mouvement[M]. Princeton: Princeton University Press, 1947.</ref>，通过其消失条件判定中心是否存在。

===一般情况的判定——Lyapunov量===
====定义====
考虑平面系统
:<math>\left\{ \begin{array}{l} \dot{x} = -y + P(x,y) \\[6pt] \dot{y} = x + Q(x,y) \end{array} \right.</math>
其中 <math>P(x,y), Q(x,y)</math> 为高阶非线性项。通过形式级数法构造
:<math>F(x,y) = x^2 + y^2 + F_3(x,y) + F_4(x,y) + \cdots</math>
使得沿系统轨线
:<math>\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>L_1, L_2, L_3, \ldots</math> 称为'''Lyapunov量'''（Lyapunov quantities）或'''焦点量'''（focal values）。

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

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

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

===特殊情况1——保守系统中的非线性中心===

==== 定义 ====
在经典力学与非线性动力学中，若一个二阶系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>（其中 <math>\mathbf{x} \in \mathbb{R}^2</math>）存在一个光滑的实值函数 <math>E(\mathbf{x})</math>，使得沿任意轨迹 <math>\mathbf{x}(t)</math> 均有 <math>\frac{dE}{dt} = 0</math>，且 <math>E(\mathbf{x})</math> 在任何开集上不为常数，<math>E(\mathbf{x})</math> 称为'''守恒量'''，称具有该守恒量的系统为'''保守系统'''<ref>STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017</ref>。

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

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

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

[[文件:保守系统的非线性中心示例图.png|缩略图|图3：双阱系统能量等高线图及其相平面投影|right|360px]]

==== 示例：双阱势系统 ====
考虑一个单位质量粒子在双阱势 <math>V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4</math> 中的无阻尼运动。运动方程为 <math>\ddot{x} = x - x^3</math>，等价系统为：
:<math>\dot{x} = y, \quad \dot{y} = x - x^3.</math>
系统具有守恒量（总能量）：
:<math>E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4.</math>
不动点为 <math>(0,0)</math> 和 <math>(\pm 1, 0)</math>。在 <math>(0,0)</math> 处，能量曲面的等高线呈现鞍点特征；而在 <math>(\pm 1, 0)</math> 处，<math>E(x, y)</math> 达到局部极小值。根据定理，这两个不动点均为非线性中心，附近存在一族闭轨（周期轨道）。

===特殊情况2——可逆系统中的非线性中心===
====定义====
如果存在一个相空间到自身的映射 <math>R(\mathbf{x})</math>，满足 <math>R^2(\mathbf{x}) = \mathbf{x}</math>（即 <math>R</math> 是一个对合变换），使得系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math> 在变换 <math>t \to -t</math>，<math>\mathbf{x} \to R(\mathbf{x})</math> 下保持不变，则称该系统为'''可逆系统'''。

最常见的二维情形是系统具有关于 <math>x</math> 轴的反射对称性：即变换 <math>t \to -t</math>，<math>y \to -y</math> 保持系统不变。这等价于要求方程 <math>\dot{x} = f(x,y)</math>，<math>\dot{y} = g(x,y)</math> 中的 <math>f</math> 关于 <math>y</math> 是奇函数（<math>f(x,-y) = -f(x,y)</math>），且 <math>g</math> 关于 <math>y</math> 是偶函数（<math>g(x,-y) = g(x,y)</math>）。

==== 定理（可逆系统的非线性中心）====
假设原点 <math>\mathbf{x}^* = \mathbf{0}</math> 是一个连续可微可逆系统的'''线性中心'''。若系统是可逆的，则原点也是一个'''非线性中心'''，即所有充分接近原点的轨迹都是闭合的。

'''证明思路：''' 考虑一条从正 <math>x</math> 轴出发靠近原点的轨迹。由于线性部分是中心，该轨迹将绕原点旋转并最终与负 <math>x</math> 轴相交。利用可逆性，将轨迹关于 <math>x</math> 轴反射并反转时间方向，将得到一条具有相同端点但箭头相反的轨迹。这两条轨迹共同构成一个闭合曲线。因此原点被一族闭轨包围。

'''判定逻辑：'''
该定理表明，可逆系统所具备的时间反演对称性为非线性中心提供了独立于能量守恒的稳定性机制。即便系统不满足保守条件，只要其在时间反演与空间反射的复合变换下保持不变，线性中心便具有抵御非线性扰动的能力。其几何解释为：若非线性扰动致使闭合轨迹转变为向内螺旋的焦点，则由可逆性推知其反射逆像必然形成向外螺旋的轨迹。此种轨迹行为与微分方程解的唯一性相冲突，因而在可逆系统中，线性中心在非线性扰动下仍保持为闭合轨迹族，即构成'''非线性中心'''。

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

[[文件:钟摆.png|缩略图|右|图4：无阻尼钟摆|right|280px]]

:<math>\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0.</math>
令 <math>\omega = \sqrt{g/L}</math>，<math>\tau = \omega t</math> 进行无量纲化，得：
:<math>\ddot{\theta} + \sin\theta = 0.</math>
写成相平面系统（<math>v = \dot{\theta}</math>）：
:<math>\dot{\theta} = v, \quad \dot{v} = -\sin\theta.</math>
该系统在变换 <math>\tau \to -\tau</math>，<math>v \to -v</math> 下保持不变，因此是可逆的。不动点 <math>(\theta^*, v^*) = (0,0)</math> 处的雅可比矩阵为 <math>\begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix}</math>，特征值为 <math>\pm i</math>，故为线性中心。根据可逆系统的非线性中心定理，<math>(0,0)</math> 也是一个非线性中心，对应钟摆在最低点的微小摆动。此外，系统也具有守恒量 <math>E = \frac{1}{2}v^2 - \cos\theta</math>，且 <math>(0,0)</math> 是 <math>E</math> 的局部极小值点，这从保守系统的角度也印证了该结论。

== 参见 ==

* [[非线性系统]]
* [[线性系统]]
* [[不动点]]
* [[极限环]]
* [[分岔理论]]
* [[中心流形定理]]

=='''参考资料'''==
<references />

== 外部链接 ==
* [http://www.scholarpedia.org/article/Equilibrium Scholarpedia: Equilibrium ]
* [http://www.scholarpedia.org/article/Dynamical_systems Scholarpedia: Dynamical Systems]
* [https://www.bohrium.com/sciencepedia/feynman/ordinary_differential_equations_undergraduate-the_Hartman-Grobman_theorem Bohrium-sciencepedia: the Hartman-Grobman theorem]

----
本中文词条由[[用户:王晔|王晔]]编辑，[[用户:陶如意|陶如意]]审校，欢迎在讨论页面留言。
本词条内容源自公开资料，遵守 CC3.0协议。
{{DEFAULTSORT:非线性中心}}

文件:非线性中心扰动.png

2026-06-08T11:51:13Z

王晔：

文件:非线性系统线性化失效典型情形.png

2026-06-08T11:50:47Z

王晔：王晔上传文件:非线性系统线性化失效典型情形.png的新版本

文件:非线性中心与极限环对比图.png

2026-06-08T11:47:20Z

王晔：

2026-05-24T15:58:32Z

王晔：/* 非线性系统的线性化失效 */

#重定向 [[非线性中心]]{{DISPLAYTITLE:非线性中心}}
'''非线性中心'''（'''Nonlinear Center'''）是[[非线性动力学]]理论中的基本拓扑结构之一。当系统奇点（平衡点）的局部邻域内布满环绕该点的嵌套闭合相轨迹时，该奇点即为非线性中心<ref>Strogatz, S. H. (2015). Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering. CRC press.</ref>。由于一般的中心结构缺乏结构稳定性，极易在微弱的非线性扰动下退化为向内或向外螺旋的[[焦点]]<ref>Andronov, A. A., & Pontryagin, L. (1937). Systèmes grossiers. In Dokl. Akad. Nauk. SSSR (Vol. 14, No. 5, pp. 247-250).</ref>，因此严格的非线性中心通常仅存在于受到特定物理法则约束的系统中，最典型的是具备能量守恒的[[保守系统]]或具有时间反演对称性的[[可逆系统]]<ref>Guckenheimer, J., & Holmes, P. (1983). Nonlinear oscillations, dynamical systems, and bifurcations of vector fields. Springer Science & Business Media.</ref>。对非线性中心及其在扰动下分岔出[[极限环]]的研究（即'''庞加莱-李雅普诺夫中心问题'''），是揭示系统自激振荡机制与探索'''希尔伯特第十六问题'''的核心途径<ref>Ilyashenko, Y. (2002). Centennial history of Hilbert's 16th problem. Bulletin of the American Mathematical Society, 39(3), 301-354.</ref>。

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

==='''非线性系统的线性化失效'''===
对于一般的非线性系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>，[[不动点]] <math>\mathbf{x}^*</math> 附近的局部动力学行为通常可通过线性化系统来近似判定。根据'''[[哈特曼-格罗伯曼定理]]'''（'''Hartman-Grobman Theorem'''），仅当平衡点特征值实部非零（双曲平衡点）时，即它在任何方向上都有明确的吸引或排斥倾向，非线性系统的局部拓扑才与其线性化系统同胚<ref>Hartman,P.(1960)."On the local linearization of differential equations." Proc.AMS.</ref>。对于特征值为纯虚数的临界情形，线性化方法失效。此时，平衡点附近的真实拓扑结构取决于高阶非线性项。

==='''中心的非线性扰动'''===
为直观展示高阶非线性项如何破坏线性中心，考虑经典的二维非线性[[动力系统]]：
:<math>\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>\lambda = \pm i</math>。
仅从线性化分析来看，原点表现为一个标准的[[线性中心]]。

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

[[文件:线性中心扰动.png|缩略图|右|图1：非线性中心扰动|480px]]

:<math>\left\{ \begin{array}{l} x = r\cos\theta \\ y = r\sin\theta \end{array} \right.</math>

将极坐标代入原方程并化简，系统可严格转换为如下形式：
:<math>\left\{ \begin{array}{l} \dot{r} = ar^3 \\[6pt] \dot{\theta} = 1 \end{array} \right.</math>
基于极坐标方程，高阶非线性参数 <math>a</math> 对系统相图结构的决定性影响表现如下：

* '''稳定焦点'''：当 <math> a < 0 </math> 时，对于任意 <math> r > 0</math>，径向演化 <math> \dot{r} < 0</math>。轨线随时间向内螺旋收缩，原有闭轨被高阶项破坏，原点退化为渐近稳定的[[焦点]]。
*
* '''非线性中心'''：当且仅当 <math> a = 0 </math> 时，系统无径向耗散 <math>（\dot{r} \equiv 0）</math>，角速度恒定。轨线在[[相空间]]中形成一族完美的同心圆闭轨，原点维持中心结构。
*
* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。
*
该经典实例表明，当动力系统处于纯虚数特征值的临界状态时，其相图拓扑极度脆弱，局部的真实稳定性完全由极微小的高阶非线性扰动（本例中的 <math>a r^3 </math> 项）所主导。

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

[[文件:非线性系统线性化失效典型情形.png|缩略图|图2：非线性系统线性化失效典型情形]]

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

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

=='''非线性中心'''==
在非线性系统 <math>\dot{x}=f(x)</math> 的平衡点 <math>x_0</math> 处，线性化后的雅可比矩阵具有一对纯虚数特征值 <math>\pm i\omega</math> (即线性部分为“中心”)，且加入非线性项后，系统在该平衡点领域内仍然存在一簇闭合轨道（周期解），则称该平衡点 <math>x_0</math> 为非线性中心<ref>Guckenheimer,J.& Holmes,P.(1983).Nonlinear Oscillations,Dynamical Systems,and Bifurcations of Vector Fields.Springer-Verlag,New York.</ref>。

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

=='''非线性中心判据'''==

'''经典中心问题'''（Classical Center Problem）由 Poincaré 正式提出<ref>POINCARÉ H. Mémoire sur les courbes définies par une équation différentielle (I)[J]. Journal de Mathématiques Pures et Appliquées, 1881, 7(3): 375-422.</ref>：

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

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

解决经典中心问题最直接的方法是逐阶计算系统的 Lyapunov 量<ref>LYAPUNOV A M. Problème général de la stabilité du mouvement[M]. Princeton: Princeton University Press, 1947.</ref>，通过其消失条件判定中心是否存在。

===一般情况的判定——Lyapunov量===
====定义====
考虑平面系统
:<math>\left\{ \begin{array}{l} \dot{x} = -y + P(x,y) \\[6pt] \dot{y} = x + Q(x,y) \end{array} \right.</math>
其中 <math>P(x,y), Q(x,y)</math> 为高阶非线性项。通过形式级数法构造
:<math>F(x,y) = x^2 + y^2 + F_3(x,y) + F_4(x,y) + \cdots</math>
使得沿系统轨线
:<math>\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>L_1, L_2, L_3, \ldots</math> 称为'''Lyapunov量'''（Lyapunov quantities）或'''焦点量'''（focal values）。

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

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

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

===特殊情况1——保守系统中的非线性中心===

==== 定义 ====
在经典力学与非线性动力学中，若一个二阶系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math>（其中 <math>\mathbf{x} \in \mathbb{R}^2</math>）存在一个光滑的实值函数 <math>E(\mathbf{x})</math>，使得沿任意轨迹 <math>\mathbf{x}(t)</math> 均有 <math>\frac{dE}{dt} = 0</math>，且 <math>E(\mathbf{x})</math> 在任何开集上不为常数，<math>E(\mathbf{x})</math> 称为'''守恒量'''，称具有该守恒量的系统为'''保守系统'''<ref>STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆, 译. 2版. 北京: 机械工业出版社, 2017</ref>。

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

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

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

[[文件:保守系统的非线性中心示例图.png|缩略图|图3：双阱系统能量等高线图及其相平面投影|right|360px]]

==== 示例：双阱势系统 ====
考虑一个单位质量粒子在双阱势 <math>V(x) = -\frac{1}{2}x^2 + \frac{1}{4}x^4</math> 中的无阻尼运动。运动方程为 <math>\ddot{x} = x - x^3</math>，等价系统为：
:<math>\dot{x} = y, \quad \dot{y} = x - x^3.</math>
系统具有守恒量（总能量）：
:<math>E(x, y) = \frac{1}{2}y^2 - \frac{1}{2}x^2 + \frac{1}{4}x^4.</math>
不动点为 <math>(0,0)</math> 和 <math>(\pm 1, 0)</math>。在 <math>(0,0)</math> 处，能量曲面的等高线呈现鞍点特征；而在 <math>(\pm 1, 0)</math> 处，<math>E(x, y)</math> 达到局部极小值。根据定理，这两个不动点均为非线性中心，附近存在一族闭轨（周期轨道）。

===特殊情况2——可逆系统中的非线性中心===
====定义====
如果存在一个相空间到自身的映射 <math>R(\mathbf{x})</math>，满足 <math>R^2(\mathbf{x}) = \mathbf{x}</math>（即 <math>R</math> 是一个对合变换），使得系统 <math>\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})</math> 在变换 <math>t \to -t</math>，<math>\mathbf{x} \to R(\mathbf{x})</math> 下保持不变，则称该系统为'''可逆系统'''。

最常见的二维情形是系统具有关于 <math>x</math> 轴的反射对称性：即变换 <math>t \to -t</math>，<math>y \to -y</math> 保持系统不变。这等价于要求方程 <math>\dot{x} = f(x,y)</math>，<math>\dot{y} = g(x,y)</math> 中的 <math>f</math> 关于 <math>y</math> 是奇函数（<math>f(x,-y) = -f(x,y)</math>），且 <math>g</math> 关于 <math>y</math> 是偶函数（<math>g(x,-y) = g(x,y)</math>）。

==== 定理（可逆系统的非线性中心）====
假设原点 <math>\mathbf{x}^* = \mathbf{0}</math> 是一个连续可微可逆系统的'''线性中心'''。若系统是可逆的，则原点也是一个'''非线性中心'''，即所有充分接近原点的轨迹都是闭合的。

'''证明思路：''' 考虑一条从正 <math>x</math> 轴出发靠近原点的轨迹。由于线性部分是中心，该轨迹将绕原点旋转并最终与负 <math>x</math> 轴相交。利用可逆性，将轨迹关于 <math>x</math> 轴反射并反转时间方向，将得到一条具有相同端点但箭头相反的轨迹。这两条轨迹共同构成一个闭合曲线。因此原点被一族闭轨包围。

'''判定逻辑：'''
该定理表明，可逆系统所具备的时间反演对称性为非线性中心提供了独立于能量守恒的稳定性机制。即便系统不满足保守条件，只要其在时间反演与空间反射的复合变换下保持不变，线性中心便具有抵御非线性扰动的能力。其几何解释为：若非线性扰动致使闭合轨迹转变为向内螺旋的焦点，则由可逆性推知其反射逆像必然形成向外螺旋的轨迹。此种轨迹行为与微分方程解的唯一性相冲突，因而在可逆系统中，线性中心在非线性扰动下仍保持为闭合轨迹族，即构成'''非线性中心'''。

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

[[文件:钟摆.png|缩略图|右|图4：无阻尼钟摆|right|280px]]

:<math>\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0.</math>
令 <math>\omega = \sqrt{g/L}</math>，<math>\tau = \omega t</math> 进行无量纲化，得：
:<math>\ddot{\theta} + \sin\theta = 0.</math>
写成相平面系统（<math>v = \dot{\theta}</math>）：
:<math>\dot{\theta} = v, \quad \dot{v} = -\sin\theta.</math>
该系统在变换 <math>\tau \to -\tau</math>，<math>v \to -v</math> 下保持不变，因此是可逆的。不动点 <math>(\theta^*, v^*) = (0,0)</math> 处的雅可比矩阵为 <math>\begin{pmatrix}0 & 1 \\ -1 & 0\end{pmatrix}</math>，特征值为 <math>\pm i</math>，故为线性中心。根据可逆系统的非线性中心定理，<math>(0,0)</math> 也是一个非线性中心，对应钟摆在最低点的微小摆动。此外，系统也具有守恒量 <math>E = \frac{1}{2}v^2 - \cos\theta</math>，且 <math>(0,0)</math> 是 <math>E</math> 的局部极小值点，这从保守系统的角度也印证了该结论。

== 参见 ==

* [[非线性系统]]
* [[复杂系统]]
* [[极限环]]
* [[分岔理论]]
* [[吸引子]]
* [[自组织]]
* [[中心流形定理]]

=='''参考资料'''==
<references />

== 外部链接 ==
* [http://www.scholarpedia.org/article/Equilibrium Scholarpedia: Equilibrium ]
* [http://www.scholarpedia.org/article/Dynamical_systems Scholarpedia: Dynamical Systems]
{{DEFAULTSORT:非线性中心}}