更改

===基于迹—行列式平面的分类===

===基于迹—行列式平面的分类===

对于二维线性系统，局部几何结构可以由矩阵<math>\displaystyle{ A }</math>的特征值完全刻画。设矩阵<math>\displaystyle{ A }</math>的两个特征值为<math>\displaystyle{ \lambda_1,\lambda_2 }</math>，其迹（Trace）与行列式（Determinant）分别为：

对于二维线性系统，局部几何结构可以由矩阵<math>\displaystyle{ A }</math>的特征值完全刻画。设矩阵<math>\displaystyle{ A }</math>的两个特征值为<math>\displaystyle{ \lambda_1,\lambda_2 }</math>，其迹（Trace）与行列式（Determinant）分别为：

如右图Fig.1 所示，以行列式 [math]\Delta[/math] 为横轴，迹 [math]T[/math] 为纵轴，可以将系统所有的定性行为在二维平面上进行直观的几何划分

* '''鞍点（Saddle Point）：''' 当 [math]\Delta < 0[/math]（即左半平面）时，特征值为一正一负的实数。 此时鞍点在两个特征方向上表现出不同的动力学行为。负特征值对应<nowiki>'''</nowiki>稳定方向<nowiki>'''</nowiki>，该方向上的轨迹随时间逐渐趋近平衡点；正特征值对应<nowiki>'''</nowiki>不稳定方向<nowiki>'''</nowiki>，轨迹则随时间远离平衡 点。由于只有严格位于稳定方向上的轨迹才能收敛到平衡点，而附近的微小扰动通常会使轨迹转向不稳定方向，因此鞍点整体上是不稳定的。例如，系统 [math]\dot{x}=x,\ \dot{y}=-y[/math] 的两个特征值分别为 [math]1[/math] 和 [math]-1[/math]，其中 [math]y[/math] 轴为稳定方向，[math]x[/math] 轴为不稳定方向，其余轨迹呈双曲线形状，由此形成如图Fig.2所示的典型鞍点相图。[[文件:Saddle phase portrait.svg|无|缩略图|396x396像素|Fig.2 鞍点相图]]

*'''抛物线分界线：'''

*'''抛物线分界线：'''

*<nowiki>曲线 [math]T^2-4\Delta=0[/math] 是迹—行列式平面中区分实特征值与复特征值的临界线。由特征值公式 [math]\lambda_{1,2}=\dfrac{T\pm\sqrt{T^2-4\Delta}}{2}[/math] 可知，[math]T^2-4\Delta[/math] 的符号直接决定特征值的类型，因此该曲线将 [math]\Delta>0[/math] 的区域划分为结点区和焦点区。</nowiki>

*<nowiki>曲线 [math]T^2-4\Delta=0[/math] 是Fig.1 迹—行列式平面 中区分实特征值与复特征值的临界线。由特征值公式 [math]\lambda_{1,2}=\dfrac{T\pm\sqrt{T^2-4\Delta}}{2}[/math] 可知，[math]T^2-4\Delta[/math] 的符号直接决定特征值的类型，因此该曲线将 [math]\Delta>0[/math] 的区域划分为结点区和焦点区。</nowiki>

*在抛物线'''外部'''，即 [math]\Delta>0[/math] 且 [math]T^2-4\Delta>0[/math] 时，系统具有两个同号的实特征值，对应于'''结点（Node）'''。当 [math]T<0[/math] 时，两个特征值均为负，轨迹沿不同特征方向趋近平衡点，形成稳定结点；当 [math]T>0[/math] 时，两个特征值均为正，轨迹从平衡点向外发散，形成不稳定结点。

*在抛物线'''外部'''，即 [math]\Delta>0[/math] 且 [math]T^2-4\Delta>0[/math] 时，系统具有两个同号的实特征值，对应于'''结点（Node）'''。当 [math]T<0[/math] 时，两个特征值均为负，轨迹沿不同特征方向趋近平衡点，形成稳定结点；当 [math]T>0[/math] 时，两个特征值均为正，轨迹从平衡点向外发散，形成不稳定结点。

*在分界曲线 [math]T^2-4\Delta=0[/math] 上，两个特征值合并为重实根 [math]\lambda_1=\lambda_2=T/2[/math]，系统处于结点与焦点之间的临界状态，通常对应星形结点或退化结点，其具体相图还取决于雅可比矩阵是否具有两个线性无关的特征向量。由此，抛物线不仅区分了轨迹是否具有旋转行为，也给出了结点与焦点之间的动力学分界，如图Fig.3 所示。

*在分界曲线 [math]T^2-4\Delta=0[/math] 上，两个特征值合并为重实根 [math]\lambda_1=\lambda_2=T/2[/math]，系统处于结点与焦点之间的临界状态，通常对应星形结点或退化结点，其具体相图还取决于雅可比矩阵是否具有两个线性无关的特征向量。由此，抛物线不仅区分了轨迹是否具有旋转行为，也给出了结点与焦点之间的动力学分界，如图Fig.3 所示。

 

*'''稳定与不稳定的分界：''' 对于右半平面（[math]\Delta > 0[/math]），横轴 [math]T = 0[/math] 是判断稳定性的关键。

*'''稳定与不稳定的分界：'''在Fig.1 迹—行列式平面 的右半平面，即 [math]\Delta>0[/math] 时，横轴 [math]T=0[/math] 构成稳定系统与不稳定系统之间的分界。这是因为两个特征值实部之和等于矩阵的迹 [math]T[/math]；当特征值为共轭复数时，它们的实部均为 [math]T/2[/math]，因此 [math]T[/math] 的符号直接决定平衡点附近的轨迹是收缩还是发散。

**当 [math]T < 0[/math]（下半平面），系统是'''渐近稳定'''的（稳定结点或稳定焦点），能量逐渐耗散，轨迹最终汇聚于原点。

**当当 [math]T<0[/math] 时，两个特征值的实部均小于零，系统为'''渐近稳定系统'''，对应稳定结点或稳定焦点。轨迹随时间逐渐趋近平衡点：在稳定结点中，轨迹沿特征方向直接汇聚；在稳定焦点中，轨迹则以螺旋方式向内收缩。

**当 [math]T > 0[/math]（上半平面），系统是'''不稳定'''的（不稳定结点或不稳定焦点），轨迹向外发散。

**当 [math]T = 0[/math]（正半横轴），特征值为纯虚数，系统能量守恒，形成围绕原点的闭合周期轨道，即'''中心（Center）'''。下文提到的“相爱相杀”模型正属于此类。

**当 [math]T=0[/math] 且 [math]\Delta>0[/math] 时，两个特征值为纯虚数 [math]\lambda_{1,2}=\pm i\sqrt{\Delta}[/math]。在线性系统中，轨迹既不趋近平衡点，也不向外发散，而是形成围绕原点的闭合周期轨道，对应中心（Center）。此时系统通常存在一个保持不变的二次型；当该二次型具有明确的物理含义时，可以将其理解为系统能量守恒。下文介绍的“相爱相杀”模型属于这一类周期性动力学情形，其典型相图如 Fig.4 所示。[[文件:Fig4 center phase portrait.svg|无|缩略图|325x325像素|Fig.4 中心相图]]

*'''退化情况：''' 抛物线本身对应着特征值相等的退化结点或星点；而纵轴（[math]\Delta = 0[/math]）对应着存在零特征值的非孤立不动点（即存在一条由平衡点组成的直线）。

*'''退化情况：''' 在Fig.1 迹—行列式平面 中，抛物线 [math]T^2-4\Delta=0[/math] 对应两个特征值相等的临界情形，此时 [math]\lambda_1=\lambda_2=T/2[/math]。若雅可比矩阵具有两个线性无关的特征向量，则平衡点表现为'''星形结点（Star Node）'''，轨迹沿不同方向直接趋近或远离平衡点；若矩阵只有一个独立特征向量，则形成'''退化结点'''或非正常结点（Improper Node），一般轨迹在接近平衡点时逐渐与唯一的特征方向相切。两类结点的稳定性仍由 [math]T[/math] 的符号决定：当 [math]T<0[/math] 时为稳定结点，当 [math]T>0[/math] 时为不稳定结点。

*<nowiki>纵轴 [math]\Delta=0[/math] 对应至少存在一个零特征值。对于二维齐次线性系统，若雅可比矩阵的秩为 1，其零空间是一条直线，直线上的每一点均满足 [math]\dot{\boldsymbol{x}}=0[/math]，因而构成一条由平衡点组成的</nowiki>'''非孤立平衡点'''集合；其余轨迹则沿非零特征值对应的方向趋近或远离这条平衡点直线。若雅可比矩阵为零矩阵，则整个相平面上的点都是平衡点。需要注意的是，对于一般非线性系统，零特征值只表明平衡点是非双曲的，不能仅根据线性化结果断定一定存在平衡点直线，还需要结合高阶项或中心流形进一步判断。'''星形结点'''、'''退化结点'''和'''非孤立平衡点'''的典型相图如 Fig.5 所示。

==经典实例：二维谐振子系统<ref name=":4" />==

==经典实例：二维谐振子系统<ref name=":4" />==

从该演化方程可以看出，要完全确定系统未来的状态，仅仅知道位置<math>\displaystyle{ x }</math>这一物理量是不够的，还必须掌握系统当下的速度<math>\displaystyle{ v }</math>。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维'''[[相平面]]（Phase plane），即<math>\displaystyle{ (x,v) }</math>空间。'''

从该演化方程可以看出，要完全确定系统未来的状态，仅仅知道位置<math>\displaystyle{ x }</math>这一物理量是不够的，还必须掌握系统当下的速度<math>\displaystyle{ v }</math>。因此，描述该系统的自然状态空间不再是一维的直线，而是由位置与速度构成的二维'''[[相平面]]（Phase plane），即<math>\displaystyle{ (x,v) }</math>空间。'''

==跨学科案例：罗密欧与朱丽叶的恋爱动力学<ref name=":4" />==

==跨学科案例：罗密欧与朱丽叶的恋爱动力学<ref name=":4" /> ==

 

 

 

二维线性系统不仅可以描述机械振动、电路演化等物理过程，也可以作为抽象模型刻画个体之间的相互作用。为了说明线性系统分类在不同学科中的普适性，非线性动力学研究者 Steven Strogatz 引入了一个经典的“恋爱动力学”（Love Dynamics）模型，用以描述两个人情感状态之间的线性反馈关系。

二维线性系统不仅可以描述机械振动、电路演化等物理过程，也可以作为抽象模型刻画个体之间的相互作用。为了说明线性系统分类在不同学科中的普适性，非线性动力学研究者 Steven Strogatz 引入了一个经典的“恋爱动力学”（Love Dynamics）模型，用以描述两个人情感状态之间的线性反馈关系。