计算方法第五章1-3节.ppt

第五章 常微分方程初值问题数值解法 5.1 引 言 5.2 简单的数值方法与基本概念 5.2.2 梯形方法 5.2.3 单步法的局部截断误差与阶 5.2.4 改进的欧拉公式 5.3 Runge-Kutta方法 5.3.1 Runge-Kutta法的一般形式 5.3.2 二阶R-K方法 5.3.3 三阶与四阶R-K方法 5.3.4 变步长的龙格-库塔方法 若要使公式阶数提高，就必须使右端积分的数值求积公式 精度提高，必然要增加求积节点. 为此可将（3.3）的右端用求积公式表示为 点数 越多，精度越高， 上式右端相当于增量函数 ， 为得到便于计算的显式方法，可类似于改进Euler法，将公式表示为 （3.4） 其中 （3.3） （3.5） 这里 均为常数. (3.4)和(3.5)称为 级显式龙格-库塔（Runge-Kutta）法， 简称R-K方法. 当 时，就是Euler法， 此时方法的阶为 . 当 时，改进Euler法(3.1)，(3.2)也是其中的一种. （3.4） 下面将证明阶 . 要使公式(3.4),(3.5)具有更高的阶 ,就要增加点数 . 下面就 推导R-K方法. （3.4） （3.5） 对 的R-K方法，计算公式如下 （3.6） 这里 均为待定常数. 希望适当选取这些系数，使公式阶数 尽量高. 根据局部截断误差的定义，（3.6）的局部截断误差为 （3.7） 这里 . 为得到 的阶 ，要将上式各项在 处做Talor 展开， 由于 是二元函数，故要用到二元Talor展开， 其中 各项展开式为 （3.8） 将以上结果代入局部截断误差公式则有 要使公式（3.6）具有 阶，必须使 （3.6） 即 （3.9）的解是不惟一的. 令 ，则得 这样得到的公式称为二阶R-K方法， 如取 ，则 这就是改进Euler法（3.1）. （3.9） 若取 则 . 称为中点公式， （3.10）也可表示为 得计算公式 （3.10） 相当于数值积分的中矩形公式. 的R-K公式(3.6)的局部误差不可能提高到 . （3.6） 把 多展开一项，从（3.8）的 看到展开式中 的项是不能通过选择参数消掉的. 实际上要使 的项为零，需增加3个方程，要确定4个 参数 ，这是不可能的. 故 的显式R-K方法的阶只能是 ，而不能得 到三阶公式. 要得到三阶R-K方法，必须 . （3.11） 其中 及 均为待定参数. 此时（3.4）,（3.5）的公式表示为 公式（3.11）的局部截断误差为 （3.4） （3.5） 只要将 按二元函数Taylor展开，使 ， 可得待定参数满足方程 （3.12） 这是8个未知数6个方程的方程组，解也不是惟一的. 所以 这是一簇公式. 满足条件（3.12）的公式（3.11）统称为三阶R-K公式. 一个常见的公式为 此公式称为库塔三阶方法. 继续上述过程，经过较复杂的数学演算，可以导出各 种四阶R-K公式，下列经典公式是其中常用的一个： 可以证明其截断误差为 . 四阶Runge-Kutta方法的每一步需要计算4次函数值 ， （3.13） 例3 解 设取步长 ，从 到 用四阶龙格- 库塔方法求解初值问题 这里 ，经典的四阶龙格-库塔公式为 本章要着重考察的一阶方程的初值问题 （1.1） （1.2） 只要函数 适当光滑——譬如关于 满足利普 希茨(Lipschitz）条件 （1.3） 理论上就可以保证初值问题（1.1），（1.2）的解 存在并且惟一. 所谓数值解法，就是寻求解 在一系列离散节点 上的