[理学]数值分析课件第9章1-4节.pptVIP

第9章 常微分方程初值问题数值解法 9.1 引 言 9.2 简单的数值方法与基本概念 9.2.2 梯形方法 9.2.3 单步法的局部截断误差与阶 9.2.4 改进的欧拉公式 9.3 龙格-库塔方法 9.3.1 显式龙格-库塔法的一般形式 9.3.2 二阶显式R-K方法 9.3.3 三阶与四阶显式R-K方法 9.3.4 变步长的龙格-库塔方法 9.4 单步法的收敛性与稳定性 9.4.2 绝对稳定性与绝对稳定域 若要使公式阶数提高，就必须使右端积分的数值求积公式 精度提高，必然要增加求积节点. 为此可将（3.3）的右端用求积公式表示为 点数 越多，精度越高， 上式右端相当于增量函数 ， 为得到便于计算的显式方法，可类似于改进欧拉法，将公式表示为 （3.4） 其中 （3.3） （3.5） 这里 均为常数. (3.4)和(3.5)称为 级显式龙格-库塔（Runge-Kutta）法， 简称R-K方法. 当 时，就是欧拉法， 此时方法的阶为 . 当 时，改进欧拉法(3.1)，(3.2)也是其中的一种. （3.4） 下面将证明阶 . 要使公式(3.4),(3.5)具有更高的阶 ,就要增加点数 . 下面就 推导R-K方法. （3.4） （3.5） 对 的R-K方法，计算公式如下 （3.6） 这里 均为待定常数. 希望适当选取这些系数，使公式阶数 尽量高. 根据局部截断误差的定义，（3.6）的局部截断误差为 （3.7） 这里 . 为得到 的阶 ，要将上式各项在 处做泰 勒展开， 由于 是二元函数，故要用到二元泰勒展开， 其中 各项展开式为 （3.8） 将以上结果代入局部截断误差公式则有 要使公式（3.6）具有 阶，必须使 （3.6） 即 （3.9）的解是不惟一的. 令 ，则得 这样得到的公式称为二阶R-K方法， 如取 ，则 这就是改进欧拉法（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） 只要将 按二元函数泰勒展开，使 ， 可得待定参数满足方程 （3.12） 这是8个未知数6个方程的方程组，解也不是惟一的. 所以 这是一簇公式. 满足条件（3.12）的公式（3.11）统称为三阶R-K公式. 一个常见的公式为 此公式称为库塔三阶方法. 继续上述过程，经过较复杂的数学演算，可以导出各 种四阶龙格-库塔公式，下列经典公式是其中常用的一个： 可以证明其截断误差为 . 四阶龙格-库塔方法的每一步需要计算4次函数值 ， （3.13） 例3 解 设取步长 ，从 到 用四阶龙格- 库塔方法求解初值问题 这里 ，经典的四阶龙格-库塔公式为 表9-3列出了计算结果，同时列出了相应的精确解. 比较例3和例2的计算结果，显然龙格-库塔方法的精度高. 但由于这里放大了步长 ，所以表9-3和表9-2 所耗费的计算量几乎相同. 龙格-库塔方法的推导基于泰勒展开方法，因而它要求所求的解