欧拉方法培训课件.ppt

第一部分：常微分方程的数值解法 研究的对象： 1.问题(*)解的存在唯一性定理 2.求解问题(*)数值方法的基本思想 3. 数值方法的具体操作过程 第一节 欧拉方法 同样地，有 总的， 欧拉折线法名称的由来： 例1. 原方程的精确解为 二、梯形方法 一个比较直接自然的想法， 这样，我们有以下三种数值方法 例2. 作业 * 其中, f(x,y)为已知函数，y0 为给定的初始值， It is the exception rather than a rule when a differential equation of the above general form can be solved exactly and explicitly by elementary methods. 在区间 I 上的解，此处不妨设 (*) 求一阶常微分方程初值问题 设上述问题的精确解为 y(x). 在区域 上连续， 内满足李普希兹(Lipschitz)条件，即 ，使得对于 内任意两点 与 ,恒有 则初值问题(*)的 解存在并且唯一。 证：此定理的证明可在普通的《常微分方程》教材 上找到，故略。 定理 设函数 且在区域 D 存在正数 L 上任取 个节点 在这些节点上采用离散化方法（通常用数值积分、微 分、泰勒公式等）将上述初值问题化成关于离散变量 的相应问题，并且将所得离散问题的解 作为真解 的近似值，这样求得的 在节点 上的数值解。 求解问题(*)数值方法的基本思想是在求解区域 I 一般说来，不同的离散过程将导致不同的方法， 从而得到不同的解。 就是上述初值问题 第一步，网格剖分（在解区域上作网格剖分）： 称为 的步长， 这些 可以不相等，但有时为了计算方便，特别是 编写程序方便，一般取成相等，此时就有等步长 这样得到的剖分网格称为一致网格， 称为网格节点。 在 I 上任取 个节点 第三步，在节点上采用离散化方法（忽略高阶项）, 对连续方程作节点离散，使连续方程仅在 离散节点上成立，从而得到节点离散方程。 并用 作为真解 的近似值， 可得含有 差分格式或差分方程。 第五步，从理论上研究数值格式的局部截断误差（即 相容性）、稳定性以及收敛性与整体误差。 第六步，分析数值结果与理论分析是否一致，考察 有无局限性及可改进的方法。 第四步，将差分方程改写为线性方程组的形式并求 解，解出 第二步， 第一步，对问题(*)的求解区域作一致网格剖分， 得到一致网格节点 第二步，对连续方程作节点离散，得到节点离散方程。 第三步，由泰勒公式 改写上式为 一、欧拉方法(欧拉折线法) 得 由离散节点方程 忽略高阶项, 并用 作为真解 的近似值, 得差分方程 即 以上方法称为欧拉方法或欧拉折线法。 这就是著名的欧拉(Euler)公式。 利用欧拉方法求解常微分方程初值问题 ， 取步长为0.1。 解： 因步长为0.1, 故对[0,1]等距剖分10份， 得节点坐标 利用欧拉方法得差分方程为 数值结果见下表。 （原方程为伯努利方程，常微教材中有精确求解方法） 0111.0 0110.9 0110.8 0110.7 0110.6 0110.5 0110.4 0110.3 0110.2 0110.1 0110.0 误差 精确解 欧拉方法 节点 我们可以直接对(*)中的方程两边从 若利用左矩形公式计算右端积分，即 则得 从而得欧拉方法 欧拉法的精度之所以很低是由于采用低精度的矩形 公式计算右端积分的结果。　 若利用右矩形公式计算右端积分，即 则得 从而得后退（隐式）欧拉方法 如果改用梯形公式计算 右端的积分，可期望得到较高的精度。 即， 从而可以得到差分格式： 这个方