东南大学数值分析第六章微常分方程数值解法.docVIP

第六章 常微分方程数值解法 ——RK4法、AB4法 ******（学号） *****（姓名） 上机题目要求见教材P307,23题。 一、算法原理 题目要求采用RK4法和AB4法求解最简单的常微分方程初值问题 为求解式(1)，采用离散化方法，就是寻求解在区间上的一系列点 上的近似值。记表示相邻两个节点的间距，称为步长。 求微分方程数值解的主要问题： (1) 如何将微分方程离散化,并建立求其数值解的递推公式; (2) 递推公式的局部截断误差、数值数与精确解的误差估计； (3) 递推公式的稳定性与收敛性. Runge-Kutta方法 基本思想：通过在多预报几个点求斜率，并将其加权平均作为的近似值，以此构造更高精度的计算公式。如果每步计算四次函数 的值，完全类似的，可以导出局部截断误差为的四阶Runge-Kutta公式(RK4)： Adams显式公式 Runge-Kutta方法是单步法，计算时，只用到， 而已知信息、等没有被直接利用。可以设想如果充分利用已知信息，，…来计算，那么不但有可能提高精度，而且大大减少了计算量，这就是构造所谓线性多步法的基本思想。我们知道常微分初值问题(1)与积分 是等价的。下面我们讨论最简单的多步法——Admas方法。为了导出求微分方程(1)的数值解法，我们将(3)右边积分的被积函数用插值多项式来近似替代，可以得到4阶精度的显式公式