常微分方程的稳定点分类：节点、鞍点、焦点与极限环

要研究一个二维常微分方程组的定性行为，首先必须找到并分类它的稳定点（也称为平衡点或不动点）。这就像给系统的“DNA”做一次扫描，能直接告诉你系统在长时间后会停留在原地、周期运动，还是彻底跑偏。

1. 建立你的方程组与找到稳定点

写下一个标准的二维自治方程组，它描述了系统状态 (x, y) 随时间 t 的变化：

$$\frac{dx}{dt} = P(x,y)$$
$$\frac{dy}{dt} = Q(x,y)$$

其中，P 和 Q 是关于 x 和 y 的函数。

寻找系统的稳定点。解以下方程组：

$$P(x,y) = 0$$
$$Q(x,y) = 0$$

得到的解 (x*, y*) 就是稳定点。在这个点上，dx/dt 和 dy/dt 同时为零，系统状态不再变化。

2. 线性化：从非线性到可解的近似

大多数方程组是非线性的，直接分析困难。我们需要在稳定点附近用线性方程来近似它，这叫线性化。

计算函数 P 和 Q 在稳定点 (x*, y*) 处的雅可比矩阵 J：

$$J = \begin{bmatrix} \frac{\partial P}{\partial x} & \frac{\partial P}{\partial y} \\ \frac{\partial Q}{\partial x} & \frac{\partial Q}{\partial y} \end{bmatrix}_{(x^*, y^*)}$$

这个矩阵的每一个元素都是原方程在稳定点处的偏导数。它决定了系统在稳定点附近被扰动后的“弹回”或“逃逸”趋势。

3. 关键步骤：求解特征值

求解雅可比矩阵 J 的特征值 λ。解特征方程：

$$\det(J - \lambda I) = 0$$

其中 I 是单位矩阵。展开后通常得到一个关于 λ 的一元二次方程：

$$\lambda^2 - \text{tr}(J)\lambda + \det(J) = 0$$

其中 tr(J) 是矩阵 J 的迹（对角线元素之和），det(J) 是矩阵 J 的行列式。

利用求根公式，解出两个特征值 λ₁ 和 λ₂：

$$\lambda_{1,2} = \frac{\text{tr}(J) \pm \sqrt{\text{tr}(J)^2 - 4\det(J)}}{2}$$

这两个特征值的符号和类型，完全决定了稳定点的分类。

4. 根据特征值进行分类

现在，根据 λ₁ 和 λ₂ 的情况，对照以下规则进行分类：

计算判别式 Δ = [tr(J)]² - 4det(J)，以及 det(J) 的值。

1. 节点：特征值 λ₁ 和 λ₂ 是同号实数。

   • 若 λ₁ < 0 且 λ₂ < 0：稳定节点。附近的解曲线会沿着两条特定方向单调地趋近于该点。
   • 若 λ₁ > 0 且 λ₂ > 0：不稳定节点。附近的解曲线会沿着两条特定方向单调地远离该点。

2. 鞍点：特征值 λ₁ 和 λ₂ 是异号实数（即 det(J) < 0）。

   • 这是永远不稳定的点。解曲线会沿着一个特征方向被吸引，同时沿着另一个特征方向被排斥，形成马鞍状的路径。任何微小的扰动最终都会使系统远离。

3. 焦点：特征值 λ₁ 和 λ₂ 是共轭复数（即 Δ < 0）。

   • 若复数的实部为负（即 tr(J) < 0）：稳定焦点。附近的解曲线会旋转并螺旋式地趋近于该点。
   • 若复数的实部为正（即 tr(J) > 0）：不稳定焦点。附近的解曲线会旋转并螺旋式地远离该点。

4. 中心：特征值 λ₁ 和 λ₂ 是纯虚数（即 Δ < 0 且 tr(J) = 0）。

   • 这是一个临界情况。附近的解曲线是围绕该点的闭合曲线，系统呈现周期运动。但这种中心在非线性系统中是“脆弱”的，微小的参数改变可能将其变为焦点或极限环。

5. 进阶概念：极限环与全局行为

注意，上面的分类是局部的，仅描述稳定点附近（小范围）的行为。系统的全局行为可能更复杂。

极限环是一种孤立的周期解，在相平面上表现为一个闭合的环。所有从环内或环外出发的解曲线，最终都会趋向于这个环。它是一个稳定的全局结构，不同于中心点（中心点附近有无穷多个周期解，而极限环是唯一的）。

判断一个系统是否存在极限环非常困难，常用的方法有庞加莱-本迪克松定理（适用于二维系统），该定理指出：如果一个解曲线被限制在一个没有稳定点的有界区域内，那么它最终必须趋向于一个周期解（即极限环）。

6. 操作总结：你的分类清单

拿到一个二维系统后，按以下流程操作：

1. 求出所有稳定点 (x*, y*)。
2. 在每个稳定点处，计算雅可比矩阵 J。
3. 解特征方程，得到两个特征值 λ₁ 和 λ₂。
4. 根据 λ₁ 和 λ₂ 是实数还是复数、是同号还是异号、实部的正负，参照第4节的规则，判定该点是节点、鞍点还是焦点。
5. 结合方程组的具体形式和定理（如庞加莱-本迪克松定理），分析系统是否存在极限环等全局结构。

这个过程将二维微分方程复杂的动力学行为，分解为对几个关键数字（迹、行列式、特征值）的计算和判断，是分析非线性系统的基石。