第三章常微分方程的差分办法(17-18).pptxVIP

;第三章 常微分方程的差分方法;3、3 龙格-库塔方法; 我们称 为区间 上的平均斜率，这样只要对平均斜率 提供一种算法，相应地我们便导出一种计算格式。; 这个处理过程启示我们，如果设法在 内多预测几个点的斜率值，然后将它们加权平均作为平均斜率，则有可能构造出具有更高精度的计算公式。这就是龙格-库塔方法的基本思想。 ;2、 二阶龙格-库塔方法(用两个点的斜率平均作为平均斜率) ;于是我们就得到如下计算格式：; 代 及 入（12）式得;可知欲使上式有二阶精度，只要成立 ; 如果取 ，则上述格式称为变形的欧拉格式，亦称为中点格式。;3、三阶龙格-库塔方法(用三个点的斜率平均作为平均斜率) ;这时有计算格式;令：p=1/2,q=1 ;4、四阶龙格-库塔公式(用四个点的斜率平均作为平均斜率) ;例：运用四阶经典龙格-库塔方法计算 的解在x=0.4处的近似值。取步长h=0.2。 ;解：四阶经典龙格-库塔公式 ;由于取步长h=0.2，所以： ;， 得： ;;; 同积分的数值计算一样，微分方程的数值解法也需要选择步长。同样，我们可以采取步长加倍或折半的办法选择步长，即通过检查步长折半前后的两种计算结果的偏差： 来判断选取的步长是否合适，具体可以分为两种情况来处理。;3、4 亚当姆斯格式; 同样，我们也可导出如下隐式的二阶、三阶及四阶亚当姆斯格式：; 仿照改进的欧拉格式的构造方法，将显式及隐式两种亚当姆斯格式相匹配，可构成下列亚当姆斯预报－校正系统：; 我们可以方便地估计出亚当姆斯预报校正系统的截断误差，从而依据这种估计将该系统 就可改进为如下精度更高的计算方案： 预报 改进 校正 改进 ;本章小节