更改

=='''背景'''==

=='''背景'''==

==='''线性中心'''===

==='''线性中心'''===

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

在动力系统理论中，[[中心]]（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> 时，对应纯线性系统的平衡点即表现为[[线性中心]]。在相空间中，其轨线呈现为一族完美的同心圆或椭圆闭轨，代表着一种无耗散的理想周期振荡状态<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>\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>。对于特征值为纯虚数的临界情形，线性化方法失效。此时，平衡点附近的真实拓扑结构取决于[[高阶非线性项]]。

==='''中心的非线性扰动'''===

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

 

=== 非线性系统的线性化失效 ===

根据[[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>

:<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>

* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。

* '''不稳定焦点'''：当 <math>a > 0 </math> 时，径向演化 <math> \dot{r} > 0</math>。轨线随时间向外螺旋发散，原点退化为不稳定焦点。

*  

*  

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

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

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

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