第9章节常微分方程数值解法.pptVIP

第9章 常微分方程数值解 /* Numerical Methods for Ordinary Differential Equations */ ? 待求解的问题：一阶常微分方程的初值问题 /* Initial-Value Problem */: Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 解析解法：（常微分方程理论） 只能求解极少一类常微分方程；实际中给定的问题不一定是解析表达式，而是函数表，无法用解析解法。 步进式：根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. --------Euler’s Method §1 欧拉方法 /* Euler’s Method */ §1 Euler’s Method Taylor展开法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 几何意义 亦称为欧拉折线法 /* Euler’s polygonal arc method*/ 几何直观是帮助我们寻找解决一个问题的思路的好办法哦 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 说明 §1 Euler’s Method Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 截断误差: 实际上，y(xn) ? yn， yn 也有误差，它对yn+1的误差也有影响，见下图。但这里不考虑此误差的影响，仅考虑方法或公式本身带来的误差，因此称为方法误差或截断误差。 局部截断误差的分析：由于假设yn = y(xn) ，即yn准确，因此分析局部截断误差时将y(xn+1) 和 yn+1都用点xn上的信息来表示，工具：Taylor展开。 ? 欧拉法的局部截断误差： Rn+1 的主项 /* leading term */ §1 Euler’s Method Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. ? 欧拉法具有 1 阶精度。 在xn点用一阶向前差商近似一阶导数 在第6章讨论牛顿插值公式时 介绍了差商的概念和性质， 各阶差商可以近似各阶导数，具有不同的精度， 且可以用函数值来表示。 上一章中数值微分的方法之一 就是用差商近似导数 Euler’s method §1 Euler’s Method Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. §1 Euler’s Method ? 欧拉公式的改进： 隐式欧拉法或后退欧拉法 /* implicit Euler method or backward Euler method*/