9常微分程数值解法.pptVIP

第九章 常微分方程初值问题数值解法;本章内容;9.1 引 言; 虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法.; 初值问题的数值解法有个基本特点，他们都采取“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进. 描述这类算法，只要给出用已知信息yn,yn-1,yn-2,?计算yn+1的递推公式.;9.2 简单的数值方法与基本概念; 一般地，设已做出该折线的顶点Pn,过Pn(xn, yn)依方向场的方向再推进到Pn+1(xn+1, yn+1)，显然两个顶点Pn,Pn+1的坐标有关系; 例1 用欧拉公式求解初值问题;依次计算下去，部分计算结果见下表. ; 欧拉公式具有明显的几何意义, 就是用折线近似代替方程的解曲线，因而常称公式(2.1)为欧拉折线法.; 为了分析计算公式的精度，通常可用泰勒展开将y(xn+1)在xn处展开，则有; 如果对方程(1.1)从xn到xn+1积分，得; ; 设用欧拉公式;由于f(x, y)对y满足Lipschitz条件(1.3). 由(2.6)减(2.5)得;9.2.2 梯形方法;　为了分析迭代过程的收敛性, 将(2.7)与(2.8)相减, 得;9.2.3 单步法的局部截断误差与阶; 定义1 设y(x)是初值问题(1.1),(1.2)的准确解, 称;即为(2.3)的结果. 这里 称为局部截断误差主项. 显然Tn+1=O(h2). 一般情形的定义如下;　　若将(2.10)展开式写成;　　同样对梯形法(2.7)有;9.2.4 改进的欧拉公式;这样建立的预测—校正系统通常称为改进的欧拉公式:; 例2 用改进的欧拉法解例1中的初值问题(2.2).; 例 (两种方法的精度比较); xn;9.3 龙格—库塔方法;9.3.1 显式龙格—库塔法的一般形式;它比欧拉法的?(xn, yn, h)=f(xn, yn), 增加了计算一个右函数f 的值，可望 p=2.若要使得到的公式阶数p更大, ? 就必须包含更多的f 值. 实际上从方程(1.1)等价的积分形式(2.4) ，即;一般说来，点数r 越多，精度越高，上式右端相当于增量函数?(xn, yn, h)，为得到便于计算的显式方法，可类似于改进欧拉法(3.1),(3.2)，将公式表示为; 当r=1, ?(xn, yn, h)=f(xn, yn)时，就是欧拉法，此时方法的阶为p=1. 当r=2时，改进欧拉法(3.1),(3.1)是其中一种，下面将证明其阶p=2. 要使公式(3.4),(3.5)具有更高的阶p，就要增加点数r. 下面我们只就r=2推导R-K方法. 并给出 r=3,4 时的常用公式，其推导方法与r=2时类似，只是计算较复杂.;9.3.2 二阶显式R-K方法;其中;将以上结果代入(3.7)，则有;要使公式(3.6)具有p=2阶，必须使; 则由此可以看出在改进的欧拉公式中相当于取(xn,yn), (xn+1,yn+1)两点处斜率的平均值，近似代替平均斜率，其精度比欧拉公式提高了.;称为中点公式(变形的欧拉公式)，相当于数值积分的中矩形公式.也可以表示为; 对r=2的R-K公式(3.6)能否使局部误差提高到O(h4)? 为此 需把K2多展开一项，从(3.8)的 看到展开式中的项 是不能通过选择参数削掉的，实际上要使 h3 的项为零，需增加3个方程，要确定4个参数c1, c2, λ2, μ21，这是不可能的. 故r＝2的显式R-K方法的阶只能是p=2，而不能得到三阶公式.;9.3.3 三阶与四阶显式R-K方法;只要K1, K2将按二元泰勒展开，使Tn+1＝O(h4)，可得待定参数满足方程;这是8个未知数6个方程的方程组，解不是唯一的. 可以得到很多公式. 满足条件(3.12)的公式(3.11)统称为三阶R-K公式. 下面只给出其中一个常见的公式. ; 继续上述过程，经过较复杂的数学演算，可以导出各种四阶R-K公式，下列经典公式是其中常用的一个：; 然而值得指出的是，龙格-库塔方法的推导基于泰勒展开方法，因而它要求所求的解具有较好的光滑性质. 反之，如果解的光滑性差，那么，使用龙格-库塔方法求得的数值解，其精度可能反而不如改进的欧拉方法. 实际计算时，我们应当针对问题的具体特点选择合适的算法.;9.3.4 变步长的龙格-库塔方法;　　我们考察四阶R-K公式(3.13) ，从节点xn出发，先