19、数值积分方法：稳定性、自适应与高精度求解

数值积分方法：稳定性、自适应与高精度求解

1. 初始值问题与基尔霍夫方程

在处理一些电路问题时，会遇到初始值问题。例如，对于包含电感 (L = 1.2 \times 10^{-3} H) 和电容 (C = 5 \times 10^{-3} F) 的双回路电路，基尔霍夫方程如下：
[
\begin{cases}
L\frac{di_1}{dt} + Ri_1 + \frac{q_1 - q_2}{C} = E \
L\frac{di_2}{dt} + Ri_2 + \frac{q_2 - q_1}{C} + \frac{q_2}{C} = 0
\end{cases}
]
同时，还有电流 - 电荷关系：
[
\begin{cases}
\frac{dq_1}{dt} = i_1 \
\frac{dq_2}{dt} = i_2
\end{cases}
]

2. 数值积分方法的稳定性

2.1 稳定性的定义

数值积分方法的稳定性是指局部误差的影响不会灾难性地累积，即全局误差保持有界。若方法不稳定，全局误差会呈指数增长，最终导致数值溢出。需要注意的是，稳定性与精度无关，一个不准确的方法也可能非常稳定。稳定性由三个因素决定：微分方程、求解方法和增量 (h) 的值。除非微分方程是线性的，否则很难预先确定稳定性。

2.2 欧拉方法的稳定性

考虑线性问题 (y’ = -\lambda y)，(y(0) = \beta)（其中 (\lambda) 是正常数），其解析解为 (y(x) = \beta e^