常微分方程组求解器

常微分方程组求解器

常微分方程组求解器是一个针对耦合线性及非线性系统的一站式微分方程工具箱。输入 2×2 或 3×3 系数矩阵，该工具即可执行完整的特征值 / 特征向量分析，以 LaTeX 格式生成闭式通解和特解，将原点平衡点分类为鞍点、节点、螺旋点或中心点，并绘制带有动画轨迹的交互式相图。对于非线性平面系统，您可以输入任意右侧项 \(f(x,y)\) 和 \(g(x,y)\)，该工具将生成高精度的 RK4 相图。

什么是常微分方程组？

常微分方程组通过导数将几个关于单一变量（通常是时间 \(t\)）的未知函数耦合在一起。其最简形式为：

当 \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) 且 \(A\) 为常数矩阵时，该系统被称为线性自治系统。此时理论非常完美：整个长期行为完全由 \(A\) 的特征值决定。

线性系统的特征值求解步骤

对于 \(\mathbf{x}' = A\mathbf{x}\)，标准方法是：

1. 计算特征多项式 \(\det(\lambda I - A) = 0\)。
2. 解出特征值 \(\lambda_1, \lambda_2, \dots\)。
3. 对每个特征值，通过求解 \((A - \lambda I) v = 0\) 找到特征向量 \(v\)。
4. 组建通解作为线性组合：\(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\)。
5. 将初始条件 \(\mathbf{x}(0)\) 代入通解以确定常数 \(c_i\)。

2×2 系统的三种情况

迹-行列式平面

对于一个迹为 \(T = a_{11} + a_{22}\)、行列式为 \(D = a_{11} a_{22} - a_{12} a_{21}\) 的 2×2 矩阵，所有分类都可以包含在一个图中：

这就是为什么结果面板显著显示 \(T\)、\(D\) 和 \(\Delta = T^2 - 4D\)——这三个数字足以命名平衡点。

非线性系统与相图

大多数现实世界的 ODEs 是非线性的，没有闭式解。本工具通过四阶龙格-库塔 (RK4) 方法对数值进行积分处理。该方法具有 \(O(h^5)\) 的局部截断误差，是处理平滑向量场的默认首选方案。

相图包含：

• 在 13×13 网格上采样的向量场，显示各点的流动方向。
• 来自您初始条件的轨迹，以红色绘制，并带有显示时间方向的橙色动态游标。
• 来自环形起始点的几条种子流线，提供动态特性的全局视图。
• 对于 2×2 线性系统，显示特征向量轴（虚线青色）——这些是解呈指数级滑动的变分方向。

如何使用此求解器

1. 选择模式 — 通过表单顶部的选项卡选择 2×2 线性、3×3 线性或非线性 2D。
2. 填写系数或方程。点击任何“快速示例”以预填典型系统（稳定节点、中心点、鞍点、单摆、Van der Pol 振荡器等）。
3. 输入初始条件 \((x_0, y_0)\) 和时间跨度 \(T\)。振荡器的典型 \(T\) 值为 6–20，快速衰减的稳定系统为 3–6。
4. 点击求解。将出现包含分类、特征值、特征向量、闭式解（线性模式）、动画相图和时间序列图的完整结果页面。
5. 使用相图下方的按钮重播轨迹，如果您想再次观察游标穿过轨迹的过程。

计算示例 — 阻尼谐振子

阻尼振子 \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\) 可以通过令 \(y = \dot{x}\) 重写为 2D 系统：

对于 \(\omega = 1\) 且 \(\zeta = 0.2\)（欠阻尼情况），矩阵为 \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\)。迹 \(T = -0.4\)，行列式 \(D = 1\)，判别式 \(\Delta = 0.16 - 4 = -3.84 < 0\)，因此我们得到一个稳定螺旋点，特征值为 \(-0.2 \pm 0.9798\,i\)。轨迹向原点螺旋汇聚，时间序列显示指数衰减的正弦曲线。

应用领域

• 机械系统 — 耦合弹簧-质量系统、单摆、陀螺仪。
• 电路系统 — RLC 网络、运放滤波器、状态空间控制。
• 种群动力学 — Lotka-Volterra 捕食者-猎物模型、竞争物种、流行病学 (SIR, SIS)。
• 化学动力学 — 反应网络、BZ 振荡反应。
• 神经科学 — FitzHugh-Nagumo 神经元模型、Hodgkin-Huxley 简化模型。
• 控制理论 — 线性化工厂模型、观测器设计、稳定性裕度。

提示与注意事项

• 如果轨迹迅速发散，请减小时间跨度 T — 不稳定系统可能在几个时间单位内就超出视口范围。
• 对于重特征值，求解器会自动通过解 \((A - \lambda I)w = v\) 找到广义特征向量 \(w\)，因此您无需手动计算即可获得 \(tv\) 项。
• 对于非线性系统，向量场箭头也会将非原点平衡点显示为青色点 — 注意相图中幅度为零的区域。
• 对于 3×3 系统，没有相图（3D 在 2D 页面中难以展示），但闭式解和稳定性判定仍然适用。
• 初始条件和时间跨度与分类无关：修改它们只会改变红色轨迹，不会改变特征值判定。

常见问题

什么是常微分方程组？

常微分方程组 (ODEs) 是一组耦合方程，它们将几个未知函数的导数与单个自变量（通常是时间）联系起来。经典形式是 \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \)，其中 \( \mathbf{x} \) 是状态向量，\(F\) 是向量场。线性系统可以简写为 \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \)，其行为几乎完全由系数矩阵 \(A\) 的特征值决定。

如何通过特征值对 2×2 线性系统的平衡点进行分类？

对于 2×2 系统 \( \mathbf{x}' = A\mathbf{x} \)，原点根据 \(A\) 的迹 \(T\) 和行列式 \(D\) 进行分类：\(D < 0\) 为鞍点（不稳定）；\(D > 0\) 且 \(T^2 > 4D\) 为节点（\(T < 0\) 稳定，\(T > 0\) 不稳定）；\(D > 0\) 且 \(T^2 < 4D\) 为螺旋点（\(T < 0\) 稳定，\(T > 0\) 不稳定，\(T = 0\) 为纯中心点）。临界情况 \(T^2 = 4D\) 产生退化节点。

当特征值为复数时，闭式解是什么样的？

如果 \(A\) 具有复共轭特征值 \( \alpha \pm i\beta \)，复特征向量为 \( v = p + iq \)，则实通解为 \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \)。指数项 \(e^{\alpha t}\) 控制幅度（增长、衰减或恒定），而正弦和余弦项处理旋转。

当矩阵有重特征值时会发生什么？

如果矩阵有重特征值 \(\lambda\) 但只有一个线性无关的特征向量 \(v\)，您还需要一个满足 \( (A - \lambda I) w = v \) 的广义特征向量 \(w\)。通解形式为 \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \)。如果特征空间恰好是二维的，则矩阵在该不变子空间上是单位矩阵的标量倍数，解简化为 \( \mathbf{x}(t) = (c_1 v_1 + c_2 v_2) e^{\lambda t} \)。

该工具可以符号化求解非线性系统吗？

非线性模式使用四阶龙格-库塔 (RK4) 积分器对系统进行数值求解并绘制相图。大多数非线性系统没有闭式解，因此这是标准方法。您仍然可以通过线性化来读取平衡点附近的局部行为，这可以通过 2×2 线性模式处理 — 在不动点处计算雅可比矩阵并将其作为 \(A\) 输入即可。

什么是相图？

相图是 2D 系统在 \(x\)–\(y\) 平面上的解的几何图形。每个解描绘出一条称为轨迹的曲线，轨迹集合与向量场箭头一起揭示了定性行为：解是向内螺旋、向外发散、振荡还是趋于平衡点。相图使系统的全局结构一目了然。

延伸阅读

• 微分方程组 — 维基百科
• 相图 — 维基百科
• 特征值和特征向量 — 维基百科
• 龙格-库塔方法 — 维基百科
• Van der Pol 振荡器 — 维基百科