常微分方程的数值解法

问题描述

对于一阶的方程：

$$\left\{ \begin{align} &y'= f(x, y) (a \leqslant x \leqslant b) \\ &y(a) = \eta \end{align} \right.$$
所谓的数值解法，就是寻求解 $y(x)$在一系列离散的点上 $a = x_0 < x_1 < x_2 < ... < x_n = b$上的近似值 $y_0, y_1, y_2, ..., y_n$。
相邻两个节点之间的间距 $h_i = x_{i+1} - x_i$是步长，一般是常数，这时节点为 $x_i = x_0 + ih (i = 0, 1, 2, ..., n)$。

欧拉方法

普通的欧拉方法使用欧拉公式计算：
$y_{i+1} = y_i + hf(x_i, y_i)$，此时$y_0 = y(x_0) = y(a) = \eta$

还有一种是改进欧拉公式：它是把欧拉公式与梯形公式联合使用。

$$\left\{ \begin{align} &y_p = y_i + hf(x_i, y_i) \\ &y_c = y_i + hf(x_{i+1}, y_p) \\ &y_{i+1} = \frac{1}{2}(y_p + y_c) \end{align} \right.$$

龙格-库塔方法

这个公式根据精度，有二阶及更高阶的，实际生活中常用四阶龙格-库塔公式。

⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪y