数字分析.ppt

第九章 常微分方程初值问题数值解法 9.1 引言 9.2 简单的数值方法与基本概念 9.3 龙格－库塔方法 9.4 单步法的收敛性与稳定性 9.5 线性多步法 9.6 方程组和高阶方程 9.1 引言 本章讨论一阶常微分方程的初值问题： 只要函数 适当光滑—如满足利普希茨条件： 理论上就能保证初值问题的解 存在并且唯一。 所谓数值解法，就是寻求解 在一系列离散点 上的近似值 ，相邻两个点间的距离 称为步长。一般情况下我们取 为常数，这是节点为： 9.1 引言 初值问题的求解有一个基本特点，它们都是采取 “步进式”求解的，即，一步一步地求函数的值。 求解的主要方法是：先对方程进行离散化，建立 求数值解的递推公式，一类在计算 时只用到前面一 步的 ，称为单步法。另一类在计算 时除了用到 理论上就能保证初值问题的解 存在并且唯一。 所谓数值解法，就是寻求解 在一系列离散点 上的近似值 ，相邻两个点间的距离 称为步长。一般情况下我们取 为常数，这是节点为： 初值问题的求解有一个基本特点，它们都是采取 “步进式”求解的，即，一步一步地求函数的值。 求解的主要方法是：先对方程进行离散化，建立 求数值解的递推公式，一类在计算 时只用到前面一 步的 ，称为单步法。另一类在计算 时除了用到 前利用前面一步的 ，还要用到前面的 等，这种 方称为多步法。 其次，要研究公式的局部截断误差和阶。数值解和 精确解的误差估计和收敛性，还有递推迭代公式的数 值稳定性问题。 显式欧拉法 在 平面上，微分方程初值问题的解 称作方 程的积分曲线。积分曲线上每一点 的斜率等于该 点的函数值 。如果按函数 在 平面上建 立一个方向场，那么，积分曲线上每一点的切线方向 均与方向场在该点的方向一致。 方称为多步法。 其次，要研究公式的局部截断误差和阶。数值解和 精确解的误差估计和收敛性，还有递推迭代公式的数 值稳定性问题。 9.2 简单的数值方法与基本概念 显式欧拉法 在 平面上，微分方程初值问题的解 称作方 程的积分曲线。积分曲线上每一点 的斜率等于该 点的函数值 。如果按函数 在 平面上建 立一个方向场，那么，积分曲线上每一点的切线方向 均与方向场在该点的方向一致。 基于上面的几何解释，我们从初始点 出发， 先依方向场在该点的方向推进到 上一点 ，然后 再从 依方向场的方向推进到 上一点 ，循此前 进做出一条折线 。 9.2 简单的数值方法与基本概念 一般地，设已做出该折线上的顶点 ，过 依方向场的方向再推进到 ，显然两个顶点 的坐标有如下关系 即 基于上面的几何解释，我们从初始点 出发， 先依方向场在该点的方向推进到 上一点 ，然后 再从 依方向场的方向推进到 上一点 ，循此前 进做出一条折线 。 一般地，设已做出该折线上的顶点 ，过 依方向场的方向再推进到 ，显然两个顶点 的坐标有如下关系 即 这就是著名的欧拉公式。若初值 已知，则依该公式 可逐步算出： 例1 求解初值问题：(其解为 ) 解：根据欧拉方法，得到： 这就是著名的欧拉公式。若初值 已知，则依该公式 可逐步算出： 例1 求解初值问题：(其解为 ) 解：根据欧拉方法，得到： 取步长 ，计算得到如下结果： 隐式欧拉法 在前面的讨论中，近似计算公式 可以看成是由 在区间 上积分得到，而 右边的积分是利用左矩形公式 近似，再以 代替 得到，现在右端的积分用右矩形公式，则得 到： 取步长 ，计算得到如