第五章 常微分方程初值问题数值解.ppt

5.2.5 Euler两步格式 改进的Euler格式中 下面构造与梯形方法在精度上相匹配的显式方法 预测 校正 在点xn处为 则有 设y(xn-1),y(xn) 的近似值yn-1,yn已知，用它带入上式右端进行计算，并取计算结果yn+1作为y(xn+1)的近似值,则可得 ( 5.7 ) Euler两步格式 取 局部截断误差（设 ） 前面介绍过的数值方法,无论是Euler方法,后退的Euler方法, 梯形方法，还是改进的Euler方法,它们都是单步法,其特点是在计算yn+1时只用到前一步的信息yn 可是公式(5.7)中除了yn外,还用到更前一步的信息yn-1,即调用了前两步的信息,故称其为Euler两步法 ( 5.8 ) 由前面知对于(5.8) ( 5.9 ) 对于(5.9),假定 则 从而有 预测 校正 可得事后估计式 可以期望，利用这样估计出的误差作为计算结果的一种补偿，有可能使精度得到改善 比 更接近准确值 比 更接近准确值 设以pn和cn分别表示第n步的预估值和校正值 预测 改进 校正 改进 设以pn和cn分别表示第n步的预估值和校正值 预测 改进 校正 改进 计算 计算 注：算法开始之前必须先给出y0,y1和p1-c1 §5.3 Runge-Kutta方法 定义5.1：如果一种方法的局部截断误差为 则称该方法具有p阶精度. 5.3.2 Runge-Kutta方法的基本思想 则 由微分中值定理有 y’(x)=f(x,y(x)) 设 称K*为区间[xn,xn+1]上的平均斜率 注:只要对K*提供一种算法，就相应的导出一种计算格式 Euler格式 可改写成 改进的Euler格式 可改写成 这个处理过程提示我们，如果设法在(xn,xn+1)内多预测几个点的斜率值，然后将它们加权平均作为平均斜率K*，则有可能构造出具有更高精度的计算格式 Runge-Kutta方法的基本思想 5.3.3 二阶Runge-Kutta方法 在 上取两点xn和 以这两点处的斜率值K1和K2的加权平均(或称为线性组合)来求取平均斜率K*的近似值K，即 式中:K1为xn点处的切线斜率值， K2为 点处的切线斜率值,比照改进的Euler法,将 视为 ，即可得 这样设计出的计算格式具有如下形式： 上式含有三个待定系数 下面通过适当选取这些系数的值，使得格式具有二阶精度. ( 5.10 ) 将以上结果代入,得 设 则 其中 （在(xn ,yn)Taylor展开） K2 K1 比较系数后可知,只要 成立,格式(5.10)的局部截断误差就等于 二阶Runge-Kutta格式 ( 5.11 ) (1)若取 则 这是无穷多解中的一个解，将以上所解的值代入式(5.10)并改写可得 (2)若取 则 此时二阶Runge-Kutta法的计算公式为 此计算公式称为变形的Euler格式. 5.3.4 三阶Runge-Kutta方法 为了进一步提高精度，设除 外再增加一点 并用三个点 的斜率K1,K2,K3加权平均 得出平均斜率K*的近似值，这时计算格式具有形式: 为了预报点 的斜率值K3,在区间 内有两 个斜率值K1和K2可以用,可将K1,K2加权平均得出 上的平均斜率,从而得到 的预报值 这样设计出的计算格式具有如下形式： 可以通过适当选取 这些系数的值，使得格式具有三阶精度. ( 5.12 ) 从而 和前面的推导类似，只要 成立,格式(5.12)的局部截断误差就等于 三阶Runge-Kutta格式 ( 5.13 ) 下列是其中的一种，称为Kutta公式. 5.3.5 四阶Runge-Kutta方法 如果需要再提高精度，用类似上述的处理方法，只需在区间 上用四个点处的斜率加权平均作为平均斜率K*的近似值，构成一系列四阶Runge-Kutta塔公式. 具有四阶精度，即局部截断误差是 . ( 5.14 ) 由于推导复杂，这里从略，只介绍最常用的一种四阶经典Runge-Kutta公式. 龙格—库塔方法的推导基于Taylor展开方法，因而它要求所求的解具有较