常微分方程组的数值解.ppt

预报-校正格式 Adams预测-校正公式 修正： 校正： 修正： §3 单步法的收敛性与稳定性 /* Convergency and Stability */ 前面介绍了两大类微分方程数值解法：一类是用差商近似导数得到的尤拉系列公式，另一类是基于平均斜率概念的Runge—Kutta公式。基本思想都是通过某种离散化手续，将微分方程转化为差分方程（代数方程）来求解。 Q1. 这种转化是否合理？要看差分问题的解yn当h?0时是否收敛到微分方程的解y(xn),即是否成立 yn ? y(xn)， h?0. -----收敛性问题 Q2. 实际计算时，由于舍入误差的影响，差分方程的解本身也有误差，这种误差在计算过程中会不会扩大 -----稳定性问题 本例的Euler公式为 由此式递推可得： 定义 ： 若一种数值方法对任意固定的 当 （同时 ） 时 有 ，则称方法是收敛的。 考察Euler方法的收敛性 例如：对初值问题 所谓的单步法(显式)，就是在计算 时只用到它前一步的信息 。 Taylor级数法，Runger-Kutta方法等都是单步法的例子。显然 单步法的共同特征是，他们都是将 加上某种形式的增量得 出 ，其计算公式形如： 式中的 称作增量函数。例如： 对Euler公式，有 对改进的Euler公式有： 关于单步法有下述的收敛定理： (1) (2) 定理：设单步法（1）具有p阶的精度 ，且增量函数 关于y的满足Lipschitz条件 又设初值 是准确的 ，即y0=y(x0),则其的整体截断误差为: (3) 而整体截断误差为 证明: 设 表示当yn =y(xn)时， 由公式(1)求得的结果，即 则局部截断误差为 (2) (2)-(1),得 §3 Convergency and Stability 从而 §3 Convergency and Stability §3 Convergency and Stability Euler法的收敛性 ： ?(x,y)=f(x,y),故当f(x,y)满足Lipschitz条件时，尤拉法收敛； §3 Convergency and Stability 因此改进的Euler方法 也是收敛 的。 类似的可以证明其 它的Runger-Kutta方法的收敛性。 §3 Convergency and Stability ? 稳定性 /* Stability */ 例：考察初值问题 在区间[0, 0.5]上的解。 分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。 0.0 0.1 0.2 0.3 0.4 0.5 精确解 改进欧拉法 欧拉隐式 欧拉显式 节点 xi 1.0000 ?2.0000 4.0000 ?8.0000 1.6000?101 ?3.2000?101 1.0000 2.5000?10?1 6.2500?10?2 1.5625?10?2 3.9063?10?3 9.7656?10?4 1.0000 2.5000 6.2500 1.5626?101 3.9063?101 9.7656?101 1.0000 4.9787?10?2 2.4788?10?3 1.2341?10?4 6.1442?10?6 3.0590?10?7 定义 若某算法在计算过程中任一步产生的误差在以后的计算中都逐步衰减，则称该算法是绝对稳定的。 一般分析某算法的稳定性时，为简单起见，只考虑模型方程或试验方程 当步长取为 h 时，将某算法应用于上式，并假设只在初值产生误差 ，则若此误差以后逐步衰减，就称该算法相对于 绝对稳定， 的全体构成绝对稳定区域。我们称算法A 比算法B 稳定，就是指 A 的绝对稳定区域比 B 的大。 h l h = h §3 Convergency and Stability §3 Convergency and Sta