二阶龙格库塔公式.ppt

第一节 常微分方程 第二节 欧拉方法 第三节 龙格—库塔法 在上一节中，我们得到了一些求微分方程近似解的数值方法，这些方法的局部截断误差较大，精度较低，我们希望得到有更高阶精度的方法。 一阶龙格—库塔方法 如果以y(x)在xi处的斜率作为y(x)在 [xi，xi+1]上的平均斜率k*，即 二阶龙格—库塔方法 在[xi，xi+1]上取两点xi，xi+p（0 p≤1）的斜率值k1，k2的线性组合λ1k1+λ2k2作为k*的近似值(λ1、λ2为待定常数)，此公式一般形式可写成 这就是二阶龙格—库塔法公式。 三阶龙格—库塔公式 为了提高精度,考虑在[xi，xi+1]上取三点xi,xi+p,xi+q的斜率值分别为k1,k2,k3,将k1,k2,k3的线性组合作为平均斜率k*的近似值，其中 k* 这就是欧拉法. 则得 其中k1 = f (xi，yi)，k2为[xi，xi+1]内任意一点xi+p = xi+ ph (0 p≤1) 的斜率f (xi+p，y(xi+p))。 由于y (xi+p)并没有给出,所以先应该求y (xi+p)，仿照改进欧拉公式的构造思想，得到 （8-7） 这样构造出的公式为 k1 k2 k* 公式中含有三个参数λ1,λ2和p,如果我们适当选取参数的值，可以使公式的局部截断误差为O(h3)。 对k1和k2作泰勒展开 代入（8-7）得 （*） 又 y (x)在xi处的二阶泰勒展开式为 当x = xi+1时， ，有 （**） （**） 比较（*）与（**）的系数即可发现， 要使公式（7-7）的局部截断误差满足 ，即要求公式具有二阶 精度只要下列条件成立即可。 （8-8） 满足条件（8-8）的一簇公式统称为 二阶龙格—库塔公式。 特别的，当 塔公式就成为改进欧拉公式。 时，龙格-库 　　改进欧拉公式就是以y(x)在xi和xi+1 处的斜率k1和k2的算术平均 值作为y(x)在[xi，xi+1]上的 平均斜率k*来进行计算的。 若取 时，龙格-库塔 公式就称为变形的欧拉公式，其形式为 （8-9） 此处的 就是欧拉方法预报出的中点 处的近似解；而 等于中点的斜率值 则近似 ； 所以公式可以看作用中点斜率近似 代替平均斜率k*，因此，公式（8-9）也 称作中点公式。 xi+q=xi+qh （0p≤q1） 这时的计算公式为 其中 为了预报点xi+q的斜率k3，一种很自然的 想法是用欧拉法预报，即取 但是，这样做效率比较差。因为在区间 [xi，xi+q]内已有两个斜率可以使用，所以把k1， k2的线性组合作为[xi，xi+q]上平均斜率的近似 值，当然比用欧拉法预报y(x i+q)精度要好。由 此，得到y(xi+q)的预报值。 于是可取 从而得到三点计算公式的形式 （8-10） 下面列出其中的一种形式 类似前面的推导，利用泰勒展开适当选择参数λ1,λ2和p,q,r,s的值使上述公式具有三阶精度,即局部截断误差为O(h4)。这类公式统称为三阶龙格-库塔公式。 （8-11） 高阶龙格-库塔法 为了进一步提高精度，在[xi，xi+1]上可取多个点，预报相应点的斜率值，对这些斜率值加权平均作为平均斜率值。利用泰勒展开，比较相应系数，从而确定在尽可能高的精度下有关参数应满足的条件。 从理论上讲，可以构造任意高阶的龙格-库塔公式,但实践证明,高于四阶的龙格-库塔公式,不但计算量大，而且精确度并不一定提高。在实际计算中，四阶龙格-库塔公式是精度及计算量较理想的公式。 经典的四阶龙格-库塔公式为 （8-12） 龙格-库塔法的步长选择 上面所讨论的龙格-库塔方法是定步长的，单从每一步来说，步长 h 越小，局部截断误差就越小；但随着步长的缩小，不但引起计算量的增加，而且也有可能引起舍入误差的严重积累；而步长h太大又不能达到预期的精度要求，所以怎样选择合适的步长h，这在实际计算中是很重要的。 第八章 常微分方程数值解法 工程技术和自然科学中的许多实际问题，在数学上往往可以归结为求解微分线性方程的形式。高数中我们学过的解析方法主要是用在一些简单和特殊的微分方程求解中，而对于大量的一般形式的微分方程，不能用解析方法求出其精确解，而只能用数值方法求近似解。 在具体求解微分方程时，要有某种定解条件，微分方程与定解条件在一起称为定解问题。定解条件有两种：一种是给出积分曲线在初始点的状态（初始条件），相应问题称作初值问题；另一种