常微分方程初值问题数值解法 - 副本.ppt

第9章 常微分方程初值问题数值解法 9.2 简单的数值方法与基本概念 9.2.2 梯形方法 9.2.3 单步法的局部截断误差与阶 9.2.4 改进的欧拉公式 9.3 龙格-库塔方法 9.3.2 二阶显式R-K方法 9.3.3 三阶与四阶显式R-K方法 9.3.4 变步长的龙格-库塔方法 将以上结果代入（3.7）则有 要使公式（3.6）具有 阶，必须使 （3.9） 即 （3.9）的解是不唯一的. 令 ，则得 这样得到的公式称为二阶R-K方法，如取 ，则 这就是改进欧拉法（3.1）. 若取 ，则 . 得计算公 式 （3.10） 称为中点公式，相当于数值积分的中矩形公式. （3.10）也可表示为 的R-K公式(3.6)的局部误差不可能提高到 . 把 多展开一项，从（3.8）的 看到展开式中 的项是不能通过选择参数消掉的. 实际上要使 的项为零，需增加3个方程，要确定4个 参数 ，这是不可能的. 故 的显式R-K方法的阶只能是 ，而不能得 到三阶公式. 要得到三阶显式R-K方法，必须 . 此时（3.4）， （3.5）的公式表示为 （3.11） 其中 及 均为待定参数，公式 （3.11）的局部截断误差为 只要将 按二元函数泰勒展开，使 ，可 得待定参数满足方程 （3.12） 这是8个未知数6个方程的方程组，解也不是唯一的. 可以得 到很多公式. 满足条件（3.12）的公式（3.11）统称为三阶R-K公式. 一个常见的公式为 此公式称为库塔三阶方法. 继续上述过程，经过较复杂的数学演算，可以导出各 种四阶龙格-库塔公式，下列经典公式是其中常用的一个： 四阶龙格-库塔方法的每一步需要计算四次函数值 ， 可以证明其截断误差为 . 9.1 引 言 本章要着重考察的一阶方程的初值问题 （1.1） （1.2） 只要函数 适当光滑——譬如关于 满足利普 希茨(Lipschitz）条件 （1.3） 理论上就可以保证初值问题（1.1），（1.2）的解 存在并且唯一. 所谓数值解法，就是寻求解 在一系列离散节点 上的近似值 . 相邻两个节点的间距 称为步长. 如不特别说明，总是假定 为定数， 这时节点为 . 初值问题（1.1），（1.2）的数值解法的基本特点是 采取“步进式”，即求解过程顺着节点排列的次序一步一步 地向前推进. 描述这类算法，只要给出用已知信息 计算 的递推公式. 首先，要对方程（1.1）离散化，建立求数值解的递 推公式. 一类是计算 时只用到前一点的值 ，称为单步 法. 另一类是用到 前面 点的值 ， 称为 步法. 其次，要研究公式的局部截断误差和阶，数值解 与 精确解 的误差估计及收敛性，还有递推公式的计算 稳定性等问题. 9.2.1 欧拉法与后退欧拉法 在 平面上，微分方程（1.1）的解 称作它 的积分曲线. 积分曲线上一点 的切线斜率等于函数 的 值. 如果按函数 在 平面上建立一个方向场，那 么，积分曲线上每一点的切线方向均与方向场在该点的方 向相一致. 基于上述几何解释，从初始点 出发，先依 方向场在该点的方向推进到 上一点 ，然后再从 依方向场的方向推进到 上一点 ，循此前进做出 一条折线 （图9-1). 一般地，设已做出该折线的顶点 ，过 依 方向场的方向再推进到