数值分析Euler方法.ppt

《数值分析Euler方法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第二节Euler方法5.2.1.Euler方法设节点为xk=x0+kh(h=(b-a)/nk=0,1,…n)方法一泰勒展开法（将y(xk+1)在xk泰勒展开得）则可得：方法二数值微分法（用向前差商近似导数）方法三数值积分法依上述公式逐次计算可得：也称Euler为单步法，又称为显格式的单步法。2欧拉法的几何意义：也称欧拉折线法.从上述几何意义上得知，由Euler法所得的折线明显偏离了积分曲线，可见此方法非常粗糙。3.欧拉法的局部截断误差：定义在假设yi=y(xi)，即第i步计算是精确的前提下，考虑的截断误差Ri=y(xi+1)yi+1称为局部截断误差定义若某算法的局部截断误差为

2、O(hp+1)，则称该算法有p阶精度。欧拉法的局部截断误差：欧拉法具有1阶精度。5.2.2后退的欧拉公式（隐式欧拉公式）向后差商近似导数由于未知数yn+1同时出现在等式的两边，故称为隐式欧拉公式，而前者称为显式欧拉公式。隐式公式不能直接求解，一般需要用Euler显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好。几何意义：向后差商近似导数x0x1))(,()(1101xyxfhyxy+见上图，显然，这种近似也有一定误差，如何估计这种误差y(xn+1)yn+1？方法同上，基于Taylor展开估计局部截断误差。但是注意，隐式公式中右边含

3、有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)隐式欧拉法的局部截断误差：§1Euler’sMethod隐式欧拉法的局部截断误差：即隐式欧拉公式具有1阶精度。§1Euler’sMethod比较欧拉显式公式和隐式公式及其局部截断误差显式公式隐式公式若将这两种方法进行算术平均，即可消除误差的主要部分而获得更高的

4、精度,称为梯形法5.2.3梯形公式在用数值积分的方法推导欧拉公式时，右端的积分用梯形积分公式可得：梯形法的迭代计算和收敛性注：的确有局部截断误差，即梯形公式具有2阶精度，比欧拉方法有了进步。但注意到该公式是隐式公式，计算时不得不用到迭代法，其迭代收敛性与欧拉公式相似。5.2.4改进的欧拉格式欧拉方法容易计算，但精度较低；梯形公式精度高，但是隐式形式，不易求解；若将二者结合，可得到改进的欧拉格式。上述方法也可以表示为下述两种形式：5.2.5欧拉两步公式中心差商近似导数x0x2x1假设，则可以导出即中点公式也具有2阶精度，且是显式的。需要2个初值y0和y1来启动递推过程，这样的算法称

5、为双步法预测-校正系统中点法具有二阶精度，且是显式的，与梯形公式精度相匹配，用中点公式作预测，梯形公式作校正，得到如下预测校正系统：校正误差约为预测误差的1/4预测误差和校正误差的事后误差估计式利用上两式可以估计预测值和校正值与准确值的误差，可以期望，利用这两个误差分别作预测值和校正值的补偿，有可能提高精度。设pn,cn分别为第n步的预测值和校正值，即此时cn+1未知，故用pn-cn代替预测-校正-改进公式注：利用该算法计算yn+1时，需要例：试分别用欧拉格式和改进的欧拉格式求解下列初值问题：解：欧拉格式的具体算式为：改进的欧拉格式为：改进的欧拉格式的具体算式为：结果列表如下：x

6、k欧拉法yk欧拉预校yk准确解y(xk)0.11.0954451150.21.1832159570.31.2649110640.41.3416407870.51.4142135620.61.4832396970.71.5491933390.81.6124515500.91.6733200531.01.7320508081.11.0959090911.1918181821.1840965691.2774378341.2662013611.3582126001.3433601511.4351329191.4164019291.5089662541.4859556021.5803382

7、381.5525140911.6497834311.6164747831.7177793481.6781663641.7847708321.737867401