数值分析方法课件.ppt
数值分析方法数值分析方法第1页，此课件共28页哦5.2.1.5.2.1.Euler方法方法设节点为设节点为xk=x0+kh (h=(b-a)/n k=0,1,n)方法一方法一 泰勒展开法泰勒展开法 （将（将y(xk+1)在在xk泰勒展开得泰勒展开得）则可得：则可得：第2页，此课件共28页哦 方法二方法二 数值微分法（用向前差商近似导数）数值微分法（用向前差商近似导数）第3页，此课件共28页哦方法三方法三 数值积分法数值积分法 第4页，此课件共28页哦依上述公式逐次计算可得：依上述公式逐次计算可得：也称也称Euler为单步法，为单步法，又称为又称为显格式的单步法显格式的单步法。第5页，此课件共28页哦 2 2 欧拉法的几何意义：欧拉法的几何意义：也称也称欧拉折线法欧拉折线法.从上述几何意义上得知，由从上述几何意义上得知，由EulerEuler法所得的折线明显法所得的折线明显偏离了积分曲线，可见此方法偏离了积分曲线，可见此方法非常粗糙。非常粗糙。第6页，此课件共28页哦3.3.欧拉法的局部截断误差：欧拉法的局部截断误差：定义定义在在假假设设 yi=y(xi)，即即第第 i 步步计计算算是是精精确确的的前前提提下下，考考虑的截断误差虑的截断误差 Ri=y(xi+1)yi+1 称为局部截断误差称为局部截断误差定义定义若若某某算算法法的的局局部部截截断断误误差差为为O(hp+1)，则则称称该该算算法法有有p 阶精度。阶精度。第7页，此课件共28页哦 欧拉法的局部截断误差：欧拉法的局部截断误差：欧拉法具有欧拉法具有 1 1 阶精度。阶精度。第8页，此课件共28页哦5.2.2 后退的后退的 欧拉公式欧拉公式（隐式欧拉公式）（隐式欧拉公式）向后差商近似导数向后差商近似导数由于未知数由于未知数 yn+1 同时出现在等式的两边，故称为同时出现在等式的两边，故称为隐式隐式 欧欧拉公式，而前者称为拉公式，而前者称为显式显式 欧拉公式。隐式公式不能直接欧拉公式。隐式公式不能直接求解，一般需要用求解，一般需要用Euler显式公式得到初值，然后用显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好。计算复杂，但稳定性好。第9页，此课件共28页哦第10页，此课件共28页哦几何意义几何意义：向后差商近似导数向后差商近似导数x0 x1)(,()(1101xyxfhyxy+第11页，此课件共28页哦 见上图，见上图，显然，这种近似也有一定误差，显然，这种近似也有一定误差，如何估计这种误差如何估计这种误差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的交点近似曲线与垂线的的交点近似曲线与垂线的交点。交点。几何意义几何意义xnxn+1PnPn+1xyy(x)第12页，此课件共28页哦隐式隐式欧拉法的局部截断误差：欧拉法的局部截断误差：第13页，此课件共28页哦1 Eulers Method第14页，此课件共28页哦隐式隐式欧拉法的局部截断误差：欧拉法的局部截断误差：即隐式欧拉公式具有即隐式欧拉公式具有 1 阶精度。阶精度。1 Eulers Method第15页，此课件共28页哦比较欧拉显式公式和隐式公式及其局部截断误差比较欧拉显式公式和隐式公式及其局部截断误差显式公式隐式公式第16页，此课件共28页哦 若将这两种方法进行算术平均，若将这两种方法进行算术平均，即可消除误差即可消除误差的主要部分而获得更高的精度的主要部分而获得更高的精度,称为梯形法称为梯形法5.2.3 梯形公式梯形公式第17页，此课件共28页哦 在用数值积分的方法推导欧拉公式时，右端的积分在用数值积分的方法推导欧拉公式时，右端的积分用梯形积分公式可得：用梯形积分公式可得：第18页，此课件共28页哦梯形法的迭代计算和收敛性梯形法的迭代计算和收敛性注：注：的确有局部截断误差的确有局部截断误差 ，即梯形公式具有即梯形公式具有2 阶精度，比欧拉方法有了进步。但阶精度，比欧拉方法有了进步。但注意到该公式是注意到该公式是隐式隐式公式，计算时不得不用到迭代法，公式，计算时不得不用到迭代法，其迭代收敛性与欧拉公式相似。其迭代收敛性与欧拉公式相似。第19页，此课件共28页哦第20页，此课件共28页哦5.2.4 改进的欧拉格式改进的欧拉格式 欧拉方法容易计算，但精度较低；梯形公式精度高，但欧拉方法容易计算，但精度较低；梯形公式精度高，但是隐式形式，不易求解；若将二者结合，可得到改进的欧拉是隐式形式，不易求解；若将二者结合，可得到改进的欧拉格式。格式。第21页，此课件共28页哦 上述方法也可以表示为下述两种形式：上述方法也可以表示为下述两种形式：第22页，此课件共28页哦5.2.5 欧拉两步公式欧拉两步公式中心差商近似导数中心差商近似导数x0 x2x1假设假设 ，则可以导出，则可以导出即中点公式也具有即中点公式也具有 2 阶精度，且是显式的。阶精度，且是显式的。需要需要2 2个初值个初值 y0和和 y1来启动递推过程，这样的算法称为来启动递推过程，这样的算法称为双步双步法法第23页，此课件共28页哦预测预测-校正系统校正系统中点法具有二阶精度，且是显式的，与梯形公式精度相匹配，用中点公式作预中点法具有二阶精度，且是显式的，与梯形公式精度相匹配，用中点公式作预测，梯形公式作校正，得到如下预测校正系统：测，梯形公式作校正，得到如下预测校正系统：校正误差约为预测误差校正误差约为预测误差的的1/4第24页，此课件共28页哦预测误差和校正误差的预测误差和校正误差的事后误差估计式事后误差估计式利用上两式可以估计预测值和校正值与准确值的误差，可以期望，利用上两式可以估计预测值和校正值与准确值的误差，可以期望，利用这两个误差分别作预测值和校正值的补偿，有可能提高精度。利用这两个误差分别作预测值和校正值的补偿，有可能提高精度。设设pn,cn分别为第分别为第n步的预测值和校正值，即步的预测值和校正值，即此时此时cn+1未知，故未知，故用用pn-cn代替代替第25页，此课件共28页哦预测预测-校正校正-改进公式改进公式注：利用该算法计算注：利用该算法计算yn+1时，需要时，需要第26页，此课件共28页哦例：试分别用欧拉格式和改进的欧拉格式求解下例：试分别用欧拉格式和改进的欧拉格式求解下列初值问题：列初值问题：解：欧拉格式的具体算式为：解：欧拉格式的具体算式为：第27页，此课件共28页哦 改进的欧拉格式为：改进的欧拉格式为：第28页，此课件共28页哦
