数值分析ppt第9章-常微分方程初值问题数值解法.pptVIP

第9章常微分方程初值问题数值解法;9.1引言;虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法.;初值问题的数值解法有个根本特点，他们都采取“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进.描述这类算法，只要给出用信息yn,yn-1,yn-2,?计算yn+1的递推公式.;9.2简单的数值方法与根本概念;x0;一般地，设已做出该折线的顶点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).;例(两种方法的精度比较);;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出发，先以h为步长