常微分方程数值解法第一节欧拉法

第五章常微分方程数值解/*NumericalMethodsforOrdinaryDifferentialEquations*/,待求解的问题：一阶常微分方程的初值问题/*Initial-ValueProblem*/:,解的存在唯一性（“常微分方程”理论）：只要f(x,y)在a,bR1上连续，且关于y满足Lipschitz条件，即存在与x,y无关的常数L使对任意定义在a,b上的y1(x)和y2(x)都成立，则上述IVP存在唯一解。,解析解法：（常微分方程理论）只能求解极少一类常微分方程；实际中给定的问题不一定是解析表达式，而是函数表，无法用解析解法。,如何求解,步进式：根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。,-EulersMethod,1欧拉方法/*EulersMethod*/,1EulersMethod,Taylor展开法,几何意义,亦称为欧拉折线法/*Eulerspolygonalarcmethod*/,几何直观是帮助我们寻找解决一个问题的思路的好办法哦,说明,1EulersMethod,截断误差:实际上，y(xn)yn，yn也有误差，它对yn+1的误差也有影响，见下图。但这里不考虑此误差的影响，仅考虑方法或公式本身带来的误差，因此称为方法误差或截断误差。局部截断误差的分析：由于假设yn=y(xn)，即yn准确，因此分析局部截断误差时将y(xn+1)和yn+1都用点xn上的信息来表示，工具：Taylor展开。,欧拉法的局部截断误差：,Rn+1的主项/*leadingterm*/,1EulersMethod,欧拉法具有1阶精度。,在xn点用一阶向前差商近似一阶导数,在第2章讨论牛顿插值公式时介绍了差商的概念和性质，各阶差商可以近似各阶导数，具有不同的精度，且可以用函数值来表示。上一章中数值微分的方法之一就是用差商近似导数,Eulersmethod,1EulersMethod,1EulersMethod,欧拉公式的改进：,隐式欧拉法或后退欧拉法/*implicitEulermethodorbackwardEulermethod*/,隐式或后退欧拉公式,由于未知数yn+1同时出现在等式的两边，故称为隐式/*implicit*/欧拉公式，而前者称为显式/*explicit*/欧拉公式。隐式公式不能直接求解，一般需要用Euler显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好,收敛性,1EulersMethod,见上图，显然，这种近似也有一定误差，如何估计这种误差y(xn+1)yn+1？方法同上，基于Taylor展开估计局部截断误差。但是注意，隐式公式中右边含有f(xn+1,yn+1),由于yn+1不准确，所以不能直接用y(xn+1)代替f(xn+1,yn+1),设已知曲线上一点Pn(xn,yn),过该点作弦线，斜率为(xn+1,yn+1)点的方向场f(x,y)方向,若步长h充分小，可用弦线和垂线x=xn+1的交点近似曲线与垂线的交点。,几何意义,1EulersMethod,隐式欧拉法的局部截断误差：,1EulersMethod,1EulersMethod,隐式欧拉法的局部截断误差：,即隐式欧拉公式具有1阶精度。,1EulersMethod,比较尤拉显式公式和隐式公式及其局部截断误差,显式公式,隐式公式,1EulersMethod,若将这两种方法进行算术平均，即可消除误差的主要部分/*leadingterm*/而获得更高的精度,称为梯形法,1EulersMethod,梯形公式/*trapezoidformula*/,显、隐式两种算法的平均,注：的确有局部截断误差，即梯形公式具有2阶精度，比欧拉方法有了进步。但注意到该公式是隐式公式，计算时不得不用到迭代法，其迭代收敛性与欧拉公式相似。,梯形法的迭代计算和收敛性,收敛性,1EulersMethod,梯形法的简化计算迭代计算量大，且难以预测迭代次数。为了控制计算量，通常只迭代一次就转入下一点的计算。用显式公式作预测，梯形公式作校正，得到如下预测校正系统，也称为改进尤拉法：,改进欧拉法/*modifiedEulersmethod*/,1EulersMethod,注：此法亦称为预测-校正法/*predictor-correctormethod*/。可以证明该算法具有2阶精度，同时可以看到它是个单步递推格式，比隐式公式的迭代求解过程简单。后面将看到，它的稳定性高于显式欧拉法。,1EulersMethod,其它形式,1EulersMethod,令x=x1,得,Anotherpointofview,对右端积分采用左矩形、右矩形、梯形积分公式，即可得尤拉显式、隐式、梯形公式,1EulersMethod,中点欧拉公式/*midpointformula*/,假设，则可以导出即中点公式也具有2阶精度，且是显式的。,需要2个初值y0和y1来启动递推过程，这样的算法称为双步法/*double-stepmethod*/，而前面的三种算法都是单步法/*single-stepmethod*/。,1EulersMethod,令x=x2,得,Anotherpointofview,对右端积分采用中矩形公式即得中点公式,1EulersMethod,summary,算例分别用显式Euler方法，梯形方法和预估校正Euler方法解初值问题,解：,取h=0.1，,梯形方法为：,续,算例分别用显式Euler方法，梯形方法和预估校正Euler方法解初值问题,解：,取h=0.1，,梯形方法为：,预估校正Euler方法：,续,数值例子表明