_常微分方程数值解法解决方案.pptVIP

第6章  常微分方程数值解法 引言 在科学技术和工程实际问题中，常常需要求解常微分方程。我们学习过几种类型的常微分方程的解析解求解方法，在更多情况下，无法求得解析解。而在多数工程应用问题中，往往不一定要求解析解，只需知道解在若干点上的函数值，即求数值解。 常微分方程初值问题 求其数值解，就是计算出解函数ｙ(ｘ)在离散点 处的近似值 。 6.1 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 初值问题的Euler方法 6.1.2 误差概述 误差概述 误差概述 误差概述 6.1.3 数值稳定性分析 数值稳定性分析 定义6.1.3 若某数值算法的绝对稳定性区域包含hλ平面上的左半平面Re(hλ)＜0，则称该方法是A稳定的。 隐式Euler法是A稳定的。 6.2 Runge-Kutta方法 Runge-Kutta方法 Runge-Kutta方法 Runge-Kutta方法 6.2.2 四阶Runge-Kutta方法 四阶Runge-Kutta方法 6.2.3 R-K法的稳定性 R-K法的稳定性 R-K法的稳定性 6.2.5 隐式R-K法 隐式R-K法 隐式R-K法 隐式R-K法 隐式R-K法 6.3 线形多步法 单步法主要依据yn的信息去计算yn+1。线性多步法是想依据yn，yn-1…，yn-r(r≥1)的信息去计算yn+1。 考虑到线性组合较为方便，因此，线性多步法一般形式可设为 6.3.1 基于数值积分的方法 基于数值积分的方法 基于数值积分的方法 基于数值积分的方法 基于数值积分的方法 基于数值积分的方法 基于数值积分的方法 Adams预估—校正法 预估 校正 并取 6.3.2 基于Taylar展开式的方法 基于Taylar展开式的方法 基于Taylar展开式的方法 6.4 数值解法得收敛性与稳定性 在学习了初值问题 (6.1) 的各种数值解法以后，为了能够正确地运用这些方法，有必要简单了解一下关于常微分方程初值问题数值解法的收敛性及稳定性。收敛性讨论的是当步长ｈ趋于零时，方法的整体截断误差是否趋于零的问题；稳定性则是讨论计算过程中的扰动 (舍入误差) 对计算结果的影响。 数值解法的收敛性 对于给定的初值问题，如果所采用的数值解法对任一固定的节点 ，当步长 数值解收敛于精确解 ，则称该数值方法是收敛得。 前面几节所介绍的显式单步法可以统一写成 其中 称为增量函数， 的具体形式依赖于方程中的f(x,y)以及离散方式。 关于显式单步法有以下收敛性定理： 定理 设初值问题 (6.1) 的数值解计算公式为 (6.25) ，且满足 （1）局部截断误差 （2）增量函数 关于变量y满足Lipschitz条件，则单步法的整体误差是P阶的，即当 时， 单步法收敛。 数值解法的稳定性 一个收敛的数值解法，截断误差的影响随步长ｈ的减小而减小。但另一方面，舍入误差的影响会随步长ｈ的减小而增大。在使用某种数值方法计算的过程中，如果某步长产生的舍入误差以后不能逐步减弱，累积起来势必给结果造成难以估量的影响，这样的数值方法就不宜采用。 如果某种数值方法，在节点 处的 值有大小为 δ的扰动 (舍入误差) ，而在其后的各节点 值 的扰动都不超过δ，则称该数值方法是 (绝对)稳定的。 各种数值方法的稳定性，依赖于算法过程以及方程的形式，这取决于方程右端的表达式，给讨论数值方法的稳定性带来困难。为简便起见，我们以一个简单的微分方程为例，来说明讨论数值解法稳定性的过程。以下针对微分方程 进行讨论，其中λ为常数 (可以是复常数) 。  Euler法显式格式的稳定区域  当 时，第i步运算产生的扰动，在以后的计算中逐步减弱。从而方程 (6.27) 的Euler法显式格式的稳定区域为 。 在复平面上表示以 (-1，0) 为圆心的单位圆内部。特别当λ为负实数 (λ＜0) 时，稳定域为 ，计算时步长应满足 Euler法隐式格式的稳定区域 设在节点 处值 有扰动 ，由此所