计算方法 教学课件 作者 江爱民 主编 张有正 副主编 nn05.pptVIP

上一页 下一页 返回首页 第5章 常微分方程初值问题的数值解法 5.1 数值解法的基本思路 5.2 欧拉方法 5.3 龙格-库塔法 5.4 线性多步法 5.1 数值解法的基本思路 所谓初值问题的常微分方程 是指： 为了便于计算，通常在区间[a,b]上作等步长处理，即 思路 从 出发求出 的方法有 1.化导数为差商 2.数值积分法 对于积分 可以用定积分的数值积分的 方法求得，用不同的数值求积公式，即可获得不同的微分方程数值解的不同公式. 若 若 对数值积分若用梯形公式: 3.泰勒（Taylor）展开法 不考虑余项，则有 两边同时求导得 当n 1时 5.2 欧拉方法 5.2.1 基本公式 考虑一阶常微分方程的初值问题 若采用数值积分法，便有 或 和 欧拉公式 后退欧拉公式 隐形 梯形公式 隐形 5.2.2 欧拉公式的几何解释 5.2.3欧拉方法的误差估计 根据泰勒公式，有 1.欧拉公式 或者写成 2.梯形公式 由泰勒公式有 与 相减 欧拉公式的误差是 梯形公式其每步的截断误差为 比欧拉公式要高,但是梯形公式是隐式的不易计算 5.2.4 改进的欧拉公式 ?先用欧拉公式 对 求出一个初步值 预报值 然后代入式 中，求出校正的 ? 可把它改写成下列形式： 即有 例1 取步长h 0.1，求初值问题 在x 0到x 1各节点上的数值解。 5.3 龙格-库塔法 5.3.1 原理 由微分中值定理 称K为在区间 上的平均斜率 只要平均斜率给出一个算法，就可以得到相应的计算 将改进欧拉法推广为： 相当于 5.3.2 龙格—库塔公式 思路 在区间 上再增加一个或多个预报点，即用三个或更多点上的斜率进行加权平均，作为平均斜率 的近似值，可望得到精度更高的计算公式。 用三个点进行预报： 由泰勒展开 又有 比较两式 代入前面的公式 每步的截断误差为 令p 1/2 中点 ，q 1 终点） 令p 1/3，q 2/3 类似的可得 经典龙格—库塔公式（也称为古典形式）：其每一步的截断误差为 还有库塔公式 例 用四阶龙格—库塔公式解以下微分方程的初值问题，并与改进的欧拉公式、准确值比较（取步长0.2）： 解 对本题的四阶龙格—库塔公式具体形式为 上一页 下一页 返回首页