12常微分沉方程数值解法.pptVIP

变步长的龙格—库塔方法 公　式 线性多步法 线性多步公式的导出 线性多步公式的导出(续1) 线性多步公式的导出(续2) 线性多步公式的导出(续3) 二步Euler方法 中心差商 梯形公式 梯形公式(续) 梯形公式(见上页)，实际上是Euler方法和隐式Euler方法的算术平均。 梯形公式的精度为二阶。 例：用梯形公式求下列初值问题的解在 改进的Euler方法 改进的Euler方法为Euler方法和梯形公式的结合，也称作预估---校正法。 改进的Euler方法(续1) 嵌套形式 改进的Euler方法(续2) 局部截断误差 称一种数值方法是p阶的，如果其局部截断误差为 。 Euler方法和隐式Euler方法的精度是一阶的。 二步Euler方法的精度是二阶的。 龙格-库塔方法 改进的Euler方法也可写成 二阶龙格-库塔方法 二阶龙格-库塔方法(续1) 要使二阶方法的局部截断误差为 ，四个系数值应满足下列关系式： 二阶龙格-库塔方法(续2) 特例1： 二阶龙格-库塔方法(续3) 特例2： 三阶龙格－库塔方法 四阶龙格－库塔方法 例题分析 两点说明