第3章节常微分方程的差分方法[15—16].pptVIP

第三章 常微分方程的差分方法 * 重庆大学数理学院 数 值 分 析 第八讲 主讲教师： 谭 宏 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 1．教学内容： Euler方法：Euler公式，单步显式公式极其局部截断误差；后退Euler公式，单步隐式公式极其局部截断误差；梯形公式，预测校正公式与改进Euler公式。 2．重点难点： Euler公式，预测校正公式与改进Euler公式 3．教学目标： 了解欧拉方法的几何意义、对给出的初值问题，能利用Euler公式，改进Euler公式进行数值求解 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 科学技术当中常常需要求解常微分方程的定解问题。这类问题的最简单的形式，是本章着重要考察的一阶方程的初值问题： （1） （2） 本章中我们假定右函数适当光滑以保证初值问题解的存在唯一。虽然求解常微分方程有各种各样的解析方法，但求解从实际问题中归结出来的微分方程要靠数值解法。 初值问题（1）、（2）局部解的唯一存在条件：若 连 续且满足Lipschitz条件，即存在常数L，对一切 有 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 差分法是一类重要的数值方法，这类方法是要寻求 离散节点 上的近似解 ，相邻节点间距 称为步长。 初值问题的各种差分方法都采用“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进。描述这类算法，只要给出从已知信息 计算 的递推公式，这类计算格式统称为差分格 式。 微分方程初值问题（1）、（2）的数值解法，就是求它 的解y(x) 在一系列节点 上的近似值 ，用 。 称为步长，一般总取h为常数。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 3、1 欧拉方法 1、欧拉格式 微分方程的本质特征是方程中含有导数项，这也是它难于求解的症结所在。数值解法的第一步就是设法消除其导数项，这项手续称为离散化。实现离散化的基本途径就用差商代替导数。 譬如，若在点 处列出方程 并用差商 代替 ，结果有 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 设用 的近似值 代入上式右端，记所求结果为 ，这样导出的计算公式 （3） 已先期算出 已知节点 步长 这就是众所周知的欧拉（Euler）格式，若初值 是已知的，则依据上式即可逐步算出数值解 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. y x 0 仿此不断地作下去…… 欧 拉 方 法 的 几 何 解 释 Y1 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0