二阶常微分方程求解器

二阶常微分方程求解器

二阶常微分方程求解器可处理形式为 a·y″ + b·y′ + c·y = g(x) 且具有实常系数的线性常微分方程。它能自动推导特征方程，对阻尼状态（过阻尼、临界阻尼、欠阻尼、无阻尼或不稳定）进行分类，并同时提供符号形式的闭式解和高精度数值解。交互式输出将 y(x) 和 y′(x) 的双曲线时域图与 (y, y′) 的相平面轨迹配对——这一视图能让你一目了然地识别状态：欠阻尼为向内螺旋，过阻尼为向内节点，无阻尼为闭合回路，不稳定为向外螺旋。

什么是具有常系数的二阶线性 ODE？

具有实常系数的二阶线性常微分方程是如下形式的方程：

其中 a ≠ 0，b，c 是实常数，g(x) 是强迫项。给定两个初始条件 y(x₀) = y₀ 和 y′(x₀) = y′₀，这将变成一个在 x₀ 邻域内具有唯一解的初值问题——这可以通过应用于等效一阶系统的 Picard-Lindelöf 定理得出。

如果 g(x) = 0，则方程是齐次的。否则它是非齐次的，完整解可以分解为：

其中 y_h 是相关齐次方程的通解（包含两个任意常数），而 y_p 是完整方程的任一特解。应用两个初始条件可以固定这两个任意常数。

特征方程

在齐次方程中猜测 y = e^(r·x) 可得到特征方程（或辅助方程）：

这是一个二次方程，其判别式 Δ = b² − 4ac 控制了整个定性行为：

三种根情况与阻尼状态

待定系数法（非齐次情况）

当 g(x) 采用以下简单形式之一时，待定系数法通过假设具有相同形式且带有未知系数的试验项，并代入 ODE 求解这些系数，从而提供一个特解：

• 常数 g(x) = k。试验项：y_p = K。若 c = 0 则乘以 x；若 b = 0 且 c = 0，则再乘以 x。
• n 次多项式。试验项：n 次一般多项式。如果常数项或线性项产生共振，则乘以 x 或 x²。
• 指数 g(x) = A·e^(k·x)。试验项：y_p = K·e^(k·x)。如果 k 与特征根重合，则乘以 x（单根）或 x²（二重根）——这就是共振。
• 正弦/余弦 g(x) = A·cos(ω·x) + B·sin(ω·x)。试验项：y_p = K₁·cos(ω·x) + K₂·sin(ω·x)。如果 iω 是根（纯频率共振），则乘以 x。
• 乘积与求和遵循线性叠加原理和乘法法则。

解读相平面

等效的一阶系统为 u = y，v = y′，其中 u′ = v 且 v′ = (g(x) − b·v − c·u)/a。在相平面上绘制 v 对 u 随 x 变化的参数图即可得到相平面轨迹。对于齐次自主系统（g 中不含 x），轨道由其起始点 (y₀, y′₀) 唯一确定，并能一目了然地揭示状态：

• 欠阻尼：轨迹向原点螺旋式缩减。
• 过阻尼：轨迹沿不变直线（慢特征向量）向原点靠近。
• 临界阻尼：退化节点，轨迹切于单一特征向量。
• 无阻尼：环绕原点的闭合椭圆——永久振荡。
• 不稳定：轨迹呈螺旋状或直线向无穷远处发散。

计算示例：受迫阻尼谐振子

考虑方程 y″ + 2·y′ + 5·y = 10，初始条件为 y(0) = 0, y′(0) = 0 —— 这是一个受迫的欠阻尼系统。

1. 特征方程： r² + 2r + 5 = 0 → Δ = 4 − 20 = −16 → r = −1 ± 2i。
2. 齐次解： y_h = e^(−x)·(C₁·cos 2x + C₂·sin 2x)。
3. 特解： 对于常数强迫项 g = 10：尝试 y_p = K，则 5K = 10，得到 y_p = 2。
4. 应用初始条件： y(0) = 0 → C₁ + 2 = 0 → C₁ = −2。y′(0) = 0 → −C₁ + 2C₂ = 0 → C₂ = −1。
5. 最终答案： y(x) = 2 − e^(−x)·(2·cos 2x + sin 2x) —— 随衰减包络振荡，极限 y → 2。

如何使用此计算器

1. 在第一行输入系数 a, b, c。a 必须为非零值（否则方程将变为一阶）。
2. 输入强迫项 g(x)，或者对于齐次问题将其保持为 0。该工具可为常数、最高 2 次的多项式和单个指数 A·e^(k·x)（包括共振情况）推导符号特解。
3. 提供初始条件 (x₀, y₀, y′₀)。由于方程是二阶的，必须同时指定 x₀ 处的 y 和 y′。
4. 选择绘图的 x 范围。求解器将使用 RK4 从 x₀ 开始向 x 的两个方向进行积分。
5. 点击“求解并可视化”。 你将获得特征方程及其在复平面上的根、阻尼状态分类、齐次和特解的符号形式、y 和 y′ 的双曲线图以及相平面轨迹。

常见应用

• 机械弹簧-质量-阻尼系统： m·x″ + c·x′ + k·x = F(t)。过阻尼、临界阻尼和欠阻尼对应于不同的阻尼比 ζ = c/(2·√(m·k))。
• RLC 电路： 串联 RLC 电路遵循 L·Q″ + R·Q′ + Q/C = V(t) —— 结构相同，符号不同。
• 单摆（小角度）： θ″ + (g/L)·θ = 0 产生简谐运动；加入空气阻力则产生阻尼振荡。
• 建筑物的地震响应： 以基础加速度作为强迫项的单自由度结构。
• PID 控制的伺服系统： 闭环误差动力学可简化为二阶 ODE，其阻尼比决定了超调量。
• 带惯性的种群模型： 资本积累滞后的经济增长模型，或响应延迟的生态模型。

数值方法 — 二维系统上的经典龙格-库塔 (RK4)

此工具将 a·y″ + b·y′ + c·y = g(x) 简化为一阶系统：

初始条件为 u(x₀) = y₀, v(x₀) = y′₀。然后对向量状态 (u, v) 应用四阶龙格-库塔法。RK4 的局部截断误差为 O(h⁵)，全局误差为 O(h⁴)；每个方向默认的 400 个子步在非刚性问题中可提供大约六位数的精度。

常见问题解答

什么是具有常系数的二阶线性 ODE？

具有常系数的二阶线性 ODE 的形式为 a·y″ + b·y′ + c·y = g(x)，其中 a, b, c 是实常数，g(x) 是强迫（非齐次）项。给定两个初始条件 y(x₀) = y₀ 和 y′(x₀) = y′₀，解是唯一的。齐次情况 g(x) = 0 总是可以通过特征方程 a·r² + b·r + c = 0 获得闭式解；非齐次情况的解为 y(x) = y_h(x) + y_p(x)。

什么是特征方程？

对于 a·y″ + b·y′ + c·y = 0，代入试探解 y = e^(r·x) 可得到 a·r² + b·r + c = 0 —— 即特征方程或辅助方程。它的根决定了齐次解的形式：两个互异实根给出 y_h = C₁·e^(r₁·x) + C₂·e^(r₂·x)；重根 r 给出 y_h = (C₁ + C₂·x)·e^(r·x)；共轭复根 α ± β·i 给出 y_h = e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x))。

欠阻尼、临界阻尼和过阻尼是什么意思？

这些术语源自弹簧-质量-阻尼器模型 m·x″ + c·x′ + k·x = 0。过阻尼（判别式 > 0，两个实根）意味着系统缓慢返回平衡状态且无振荡。临界阻尼（判别式 = 0，重根）是不产生超调的最快返回方式。欠阻尼（判别式 < 0，复根）产生衰减振荡。无阻尼 (b = 0, c/a > 0) 产生永久的正弦振荡。

什么是待定系数法？

对于简单的强迫项 g(x) —— 常数、多项式、指数、正弦、余弦及其乘积 —— 假设特解 y_p 与 g 具有相同的形式且带有未知系数，通过代入 ODE 并匹配各项来确定这些系数。当 g(x) 与特征根产生共振时，试验项必须乘以 x（或对于二重根乘以 x²）。

什么是相平面？

对于简化为二维系统 (y, y′) 的二阶方程，相平面随 x 的推进绘制 y′ 对 y 的图。相平面中的解曲线能一目了然地揭示状态：欠阻尼为衰减螺旋，过阻尼为向内节点，无阻尼简谐运动为闭合椭圆，不稳定振荡为向外螺旋。它是特征方程根图的几何对应物。

此工具使用什么数值方法？

经典的四阶龙格-库塔 (RK4) 方法应用于等效的一阶系统 u = y, v = y′，其中 u′ = v 且 v′ = (g(x) − b·v − c·u)/a。RK4 的局部截断误差为 O(h⁵)，每个方向默认 400 个子步，在选定窗口内为非刚性方程提供大约六位数的精度。

延伸阅读

• 线性微分方程 — 维基百科
• 特征方程 — 维基百科
• 待定系数法 — 维基百科
• 谐振子 — 维基百科
• 相平面 — 维基百科
• 龙格-库塔法 — 维基百科