常微分方程的差分方法()培训教程文件.pptVIP

第三章 常微分方程的差分方法 * 重庆大学数理学院 数 值 分 析 第八讲 主讲教师： 谭 宏 1．教学内容： Euler方法：Euler公式，单步显式公式极其局部截断误差；后退Euler公式，单步隐式公式极其局部截断误差；梯形公式，预测校正公式与改进Euler公式。 2．重点难点： Euler公式，预测校正公式与改进Euler公式 3．教学目标： 了解欧拉方法的几何意义、对给出的初值问题，能利用Euler公式，改进Euler公式进行数值求解 科学技术当中常常需要求解常微分方程的定解问题。这类问题的最简单的形式，是本章着重要考察的一阶方程的初值问题： （1） （2） 本章中我们假定右函数适当光滑以保证初值问题解的存在唯一。虽然求解常微分方程有各种各样的解析方法，但求解从实际问题中归结出来的微分方程要靠数值解法。 初值问题（1）、（2）局部解的唯一存在条件：若 连 续且满足Lipschitz条件，即存在常数L，对一切 有 差分法是一类重要的数值方法，这类方法是要寻求 离散节点 上的近似解 ，相邻节点间距 称为步长。 初值问题的各种差分方法都采用“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进。描述这类算法，只要给出从已知信息 计算 的递推公式，这类计算格式统称为差分格 式。 微分方程初值问题（1）、（2）的数值解法，就是求它 的解y(x) 在一系列节点 上的近似值 ，用 。 称为步长，一般总取h为常数。 3、1 欧拉方法 1、欧拉格式 微分方程的本质特征是方程中含有导数项，这也是它难于求解的症结所在。数值解法的第一步就是设法消除其导数项，这项手续称为离散化。实现离散化的基本途径就用差商代替导数。 譬如，若在点 处列出方程 并用差商 代替 ，结果有 设用 的近似值 代入上式右端，记所求结果为 ，这样导出的计算公式 （3） 已先期算出 已知节点 步长 这就是众所周知的欧拉（Euler）格式，若初值 是已知的，则依据上式即可逐步算出数值解 y x 0 仿此不断地作下去…… 欧 拉 方 法 的 几 何 解 释 Y1 例1 求解初值问题 （其解析解为） 解： 设步长 h=0.1，由欧拉公式（3）有： 所以， …… 计算结果表 1.0 0.9 0.8 0.7 0.6 xn 1.7321 1.7848 1.4142 1.4351 0.5 1.6733 1.7178 1.3416 1.3582 0.4 1.6125 1.6498 1.2649 1.2774 0.3 1.5492 1.5803 1.1832 1.1918 0.2 1.4832 1.5090 1.0945 1.1000 0.1 y(xn ) yn y(xn ) yn xn 解析解 数值解 为简化分析，人们常假设在第n步求得的 为准确即 的前提下估计误差 这种误差称为局部截断误差。 误差估计为： y(xn+1)-[y(xn)+hf（xn, y(xn)）] 如果不作这一假定，累积了n 步的误差，称为 整体截断误差。其表达式为 y(xn+1)-yn+1=y(xn+1)-［yn+hf(xn, yn)］ 如果一种数值方法的局部截断误差为 则称它的的精度是 p 阶的，或称之为 p 阶方法。 对于欧拉格式（3），假定 则有： 由此我们可知欧拉格式仅为一阶方法。 将 在 点泰勒展开： 因此有： 虽然欧拉公式（3）的精确度很差，但却体现了 数值方法的基本思想。 2、隐式欧拉格式 设改用向后差商 替代方程 中的导数项 再离散化，即可导出下列格式 （5） 该格式右端含有未知的 它实际上是个关于 的函数方程。故称该格式为隐式欧拉格式。 由于向前差商和向后差商具有同等精度，故隐式欧拉 格式也是一阶方法，精度与欧拉格式相当。