计算方法Chapter03-常微分方程的数值解法精选.pptVIP

* 转化例子 将 转化为微分方程组。 解：令 ，有： 初始条件： * * 作业 课本 P124： 3 5 10 * * hhhhh * * 二阶龙格-库塔方法 以 k1 和 k2 的加权平均来近似取代 k* 为分析局部截断误差，令 yn = y(xn)，由泰勒公式得： * 补充：二元泰勒展开式 * 用二元泰勒公式展开 将 k1, k2 代入 中可得： * 二阶龙格-库塔方法（续） 2 阶精度 * 四个未知变量，只有三个方程，有无穷多组解 每组解的构成的龙格-库塔方法均为二阶 二阶龙格-库塔方法即为改进的欧拉方法 变形的欧拉法 中 点 方 法 * 三阶龙格-库塔方法 三阶龙格-库塔方法是用三个值 k1, k2, k3 的加权平均来近似取代 k* 要使三阶龙格-库塔方法具有三阶精度，必须使其局部截断误差为 O(h4) 将 k1, k2, k3 代入 yn+1 的表达式中，在 (xn, yn) 处用二元泰勒公式展开，与 y(xn+1) 在 xn 处的泰勒展开式比较 * 三阶龙格-库塔方法（续） 类似二阶龙格-库