为什么龙格-库塔法比欧拉法更精确：高阶截断误差消除

1. 理解问题：什么是数值微分方程求解？

面对一个无法求出精确解的常微分方程（例如 dy/dx = f(x, y)），你需要用数值方法一步步地计算 y 在离散点 x0, x1, x2... 上的近似值。核心目标是控制每一步计算引入的误差，让最终结果尽可能接近真实解。

2. 基础：欧拉法的核心与缺陷

实现最简单的欧拉法，只需一个公式：y_{n+1} = y_n + h * f(x_n, y_n)。这里 h 是步长（每一步 x 增加的大小）。
理解它的缺陷：这个公式相当于用当前点 (x_n, y_n) 处的斜率 f(x_n, y_n) 作为整个步长 h 范围内的平均斜率。如果真实解是弯曲的（曲率大），用起点斜率代表全程，误差就会在每一步累积，并且迅速变大。定义这个误差的阶：欧拉法的局部截断误差与 h 的平方成正比（即 O(h^2)），整体误差与 h 的一次方成正比（即 O(h)）。

3. 核心思想：用多个“试算”来提升斜率估计精度

龙格-库塔法（以经典的4阶为例）的根本改进在于：它不使用单一斜率，而是通过在区间内多个点进行“试算”，得到一组斜率，然后对它们加权平均，从而获得一个在整个步长 h 内更具代表性的平均斜率。

4. 执行经典四阶龙格-库塔法（RK4）的步骤

假设你要从 x_n 推进到 x_{n+1} = x_n + h，已知当前值 y_n。

1. 计算起点斜率 k1：k1 = f(x_n, y_n)。
2. 计算中点斜率 k2（用 k1 试探）：
   • 先用 k1 走半个步长，估算中点处的 y 值：y_temp = y_n + (h/2) * k1。
   • 再计算该中点处的斜率：k2 = f(x_n + h/2, y_temp)。
3. 计算中点斜率 k3（用 k2 试探）：
   • 用 k2 重新估算中点 y 值：y_temp2 = y_n + (h/2) * k2。
   • 再计算该中点处的斜率：k3 = f(x_n + h/2, y_temp2)。
4. 计算终点斜率 k4（用 k3 试探）：
   • 用 k3 走完整个步长，估算终点 y 值：y_temp3 = y_n + h * k3。
   • 再计算该终点处的斜率：k4 = f(x_n + h, y_temp3)。
5. 加权平均这四个斜率，得到最优化的平均斜率 k：
   k = (k1 + 2*k2 + 2*k3 + k4) / 6
6. 更新 y 值：y_{n+1} = y_n + h * k。

观察：这个过程通过两次中点试探（k2, k3）和一次终点试探（k4），综合了步长起点、中点和终点的曲率信息。

5. 定量对比：截断误差如何被消除？

分析误差来源。任何数值方法的解都可以展开为 h 的幂级数。一个精确的方法，其级数展开项应与原方程的解析解展开项完全匹配。

• 欧拉法：它的计算公式只精确匹配了解析解泰勒展开到 h^1 的项。因此，其局部截断误差（即一步之内未被匹配的剩余部分）主要来自 h^2 项，即 O(h^2)。
• RK4：其精心设计的四步计算和加权系数，使得它的计算公式能够精确匹配解析解泰勒展开到 h^4 的所有项。因此，其局部截断误差主要来自 h^5 项，即 O(h^5)。

本质：RK4通过增加计算复杂度（每一步计算4次函数 f），换来了对更高阶泰勒展开项的匹配，从而消除了欧拉法中主要的低阶误差项（h^2, h^3, h^4 项）。这就是“高阶截断误差消除”的含义。

6. 实际影响：精度与步长的关系

牢记：局部截断误差阶数越高，在总计算量相同的情况下，最终能达到的整体精度就越高。

• 使用欧拉法时，为了将误差缩小为原来的 1/10，你需要将步长 h 缩小 10 倍。这意味着你需要计算 10 倍多的步数，计算量线性增加。
• 使用RK4时，为了将误差缩小为原来的 1/10，理论上你只需要将步长 h 缩小 10^{1/4} ≈ 1.78 倍（因为整体误差与 h^4 成正比）。这意味着步数只需增加约 1.78 倍，计算量增加很少。

结论：在达到同等精度的前提下，RK4通常比欧拉法高效得多。或者，在相同的步长下，RK4的精度远高于欧拉法。这就是为什么在工程和科学计算中，当需要较高精度时，RK4是比欧拉法更常用的选择。

7. 选择与权衡

认清两者的核心权衡：

• 欧拉法：每一步计算快（只算1次 f），但精度低，需要非常小的步长，导致总步数多，总误差可能仍然很大。
• RK4：每一步计算慢（算4次 f），但精度高，允许使用大得多的步长，总步数少，在达到高精度要求时总计算成本可能更低。

执行你的选择：对于对精度要求不高或函数 f 计算极其昂贵的情况，欧拉法或低阶方法可能更合适。对于大多数需要可靠精度的通用情况，RK4是平衡精度与计算效率的黄金标准。