计 算 方 法第7章.ppt

第7章 常微分方程的数值解法 7.1 引 言 7.2 欧 拉 方 法 7.2.1 欧拉公式 定义 如果一种方法的局部截断误差为O(hp+1)，则称此种方法具有P阶精度。 图7.1 欧拉折线法 7.2.2 改进欧拉法及局部截断误差 定义 如果一种方法的局部截断误差为O(hp+1)，则称此种方法具有p阶精度。 7.3 龙格-库塔方法 7.3.1 龙格-库塔方法的基本思想 7.3.2 二阶龙格-库塔公式 综上所述，构造二阶龙格库塔公式主要由以下几步产生： 1)在区间[xi，xi+1]上取二点，预报相应点的斜率值； 2)对此两斜率值加权平均作为平均斜率值的近似值； 3)将yi+1与y(xi+1)都在xi处作泰勒展开，为使公式达到二阶精度，比较相应系数，建立有关参数所应满足的方程组； 4)解此方程组得一族二阶龙格-库塔公式。 7.3.3 高阶龙格-库塔公式 7.3.4 步长的自适应问题 7.4 线性多步法 7.4.1 阿当姆斯内插公式及误差 7.4.2 阿当姆斯外推公式及误差 7.5 一阶方程组与高阶方程 7.5.1 一阶方程组