常微分方的差分方法.pptVIP

数值分析简明教程 亚当姆斯格式3 数值分析简明教程 隐式亚当姆斯格式1 同样，我们也可导出如下隐式的二阶、三阶和四阶亚当姆斯格 式： 数值分析简明教程 隐式亚当姆斯格式2 数值分析简明教程 作业 数值分析简明教程 证明3阶Adams格式 数值分析简明教程 亚当姆斯预报-校正系统 仿照改进的欧拉格式的构造方法，将显式和隐式两种亚当姆斯格 式相匹配，可构成下列亚当姆斯预报－校正系统： 预报 校正 数值分析简明教程 改进的亚当姆斯预报-校正系统 我们可以方便地估计出亚当姆斯预报校正系统的截断误差，从而 依据这种估计将该系统 就可改进为如下精度更高的计算方案： 预报 改进 校正 改进 数值分析简明教程 记住的公式 数值分析简明教程 收敛性问题1 在用差分格式求解微分方程时我们要考虑差分格式的收敛性。我 们称差分格式是收敛的，如果对任意固定的 ，数值解 当 （同时 ）时趋于准确解 。 以下我们研究欧拉方法的收敛性。我们记 ，记 为 关于 的李普希兹常数，经反复递推，可得 其中 为常数。 若初始 是准确的，即 ，则当 时，有 。 这说明欧拉方法是收敛的。 数值分析简明教程 收敛性问题2 数值分析简明教程 收敛性问题3 数值分析简明教程 收敛性问题4 数值分析简明教程 稳定性问题1 关于收敛性的讨论有个前提，即必须假定差分方法的每一步计 算都是准确的。然而实际计算中往往由于有舍入误差等原因而产生 扰动，而这些扰动有可能 “淹没” 真解，所以我们还要考虑稳定性问 题。收敛性主要考虑算法的内在误差,而稳定性考虑外在误差。 我们称差分方法是稳定的，如果在节点值 上大小为 的扰动 于以后各节点值 上产生的偏差值均不超过 。 稳定性比较复杂。为简化讨论，我们仅考察下列模型方程 可以验证 ，对于该模型方程，欧拉格式是条件稳定的，而隐式欧拉 格式是恒稳定的。 数值分析简明教程 稳定性问题2 数值分析简明教程 数值分析简明教程 数值分析简明教程 数值分析简明教程 一阶方程组 前面我们研究了单个方程 的差分方法，只要把 和 理解为向量，所提供的各种算法即可推广应用到一阶方程组的情形。 譬如，对于方程组 令 ，以 表示节点 上的近似解， 则其欧拉格式具有形式： 数值分析简明教程 一阶方程组 则其改进的欧拉格式具有形式： 预报 校正 数值分析简明教程 数值分析简明教程 数值分析简明教程 化高阶方程为一阶方程组 关于高阶微分方程（或方程组）的初值问题，原则上总可以归结 为一阶方程组来求解。譬如，对于下列二阶方程的初值问题 若引进新的变量 ，即可化为一阶方程组的初值问题 从而可以运用前面介绍的各种算法求解。 数值分析简明教程 * 数值分析简明教程 安工大2004 4.* 数值分析简明教程 第三章 常微分方程的差分方法 1 欧拉方法 2 改进的欧拉方法 3 龙格-库塔方法 4 亚当姆斯方法 5 收敛性与稳定性 6 方程组与高阶方程的情形 7 边值问题 数值分析简明教程 引言 科学技术当中常常需要求解常微分方程的定解问题。这类问题的 最简单的形式，是本章着重要考察的一阶方程的初值问题： 本章中我们假定右函数适当光滑以保证初值问题解的存在唯一。 虽然求解常微分方程有各种各样的解析方法，但求解从实际问题中归 结出来的微分方程要靠数值解法。 差分法是一类重要的数值方法，这类方法是要寻求 离散节点 上的近似解 ，相邻节点间距 称为步长。 初值问题的各种差分方法都采用“步进式”，即求解过程顺着节点排列 的次序一步一步地向前推进。描述这类算法，只要给