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