常微分方程数值解法

问题描述

针对一阶方程的求解问题：

所谓数值求解方法，即是在给定区间 a = x_0 < x_1 < x_2 < ... < x_n = b 内，通过确定一系列离散点上的近似值 y_0, y_1, y_2, ..., y_n 来逼近解 y(x)。
在相邻节点之间，其间隔 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

此外，还存在一种优化后的欧拉公式，该方法结合了传统欧拉公式与梯形公式的应用。

龙格-库塔方法

该表达式依据精度等级，存在二阶及以上形式，在实际应用中最为常见的是四阶龙格-库塔方法。

一阶方程组求解方法

一阶方程组解析

针对如下方程组：

我们可以通过向量符号进行描述：
$ \vec y = (u, v)^T, \quad \vec \eta = (\alpha, \beta)^T,\quad\vec f = (\varphi, \psi