第7章-2收敛性稳定性R-K方法资料.ppt

（三）m=3 比较（7-16）和（7-17），令 f,h,h2 的系数 相等，并注意 的任意性，得 四个方程不能完全确定六个系数，因此这是含两个参数的三级 三阶方法类。 常见方案有： Heun三阶方法。 此时取 * * * * 常微分方程的数值解法的 收敛性、稳定性 第7章--2 以上我们讨论了求解问题（7-1）,（7-2）的单步法 和多步法。 应关注三个问题： 、数值方法的局部截断误差和阶 二、在离散点tn处的数值解un是否收敛到精确解u(tn) 三、数值方法的稳定性 具体说， 对于上述两类方法求近似解（数值解）还 误差估计、收敛性和稳定性。 对于第一个问题前面我们已经讨论过，而关于数值 方法收敛性问题我们在这里不详细讨论，只给出一些基 本结论性的结果，即： 对单步法，当方法的阶p≥1时，有整体误差 故有 ，因此方法是收敛的。 对于多步法，若方法是k 步p 阶法，那么（7-24）是 一个k阶差分方程，引入多步法（7-24）的第一特征多项 式和第二特征多项式： 定义7.1 若（7-24）的第一特征多项式ρ(λ)的所有 根在单位圆内或圆上（︱λ︱≤1），且位于单位圆周上 的根都是单根，称多步法（7-24）满足根条件。 第二特征多项式 第一特征多项式 定理7.2 若线性多步法(7-24)的阶p≥1，且满足 根条件，则方法是收敛的。 对于常用的数值方法都是满足收敛性条件的。 下面我们着重讨论第三个问题，即数值方法的稳 是有误差的，且这些误差将在计算中传递下去。 定性问题。 误差积累无限增长，则会歪曲真解，这样的算法是不 如果 能用的。 用多步法计算时，各种因素如初值 精确解为 考虑二步三阶显式法： 例如 初值问题 取步长h=0.1，初值u0=1，附加值： -6.96×108 25.0000000 2.0 … … … +367.26392 4.8841000 1.0 -68.639804 4.0000000 1.0 … … … 1.5929935 1.5625000 0.5 1.3388660 1.3456000 0.4 1.1892385 1.1881000 0.3 1.0812000 1.0816000 0.2 1.0201000 1.0201000 0.1 1.0000000 1.0000000 0 数值解 精确解 数值结果表 在开始几步数值解与精确解符合，但在再往后算，数值 解的误差则急剧增长，完全歪曲了真解. 通常人们都是通过模型方程来讨论方法的数值稳定性。 (7-32) 而一般形式的一阶微分方程总能化成（7-32）的形式。 。因为实际计算时，h是固定的。 当某一步un有舍入误差时， 若以后的计算中不会逐步扩大，称这种稳定性为绝对稳定性。 此后，若不做特殊说明，都是指绝对稳定性 。 模型方程为： 本书中数值方法的稳定性也是如此。前提是求解好条件问题， 其中Re(μ)＜0。另外，我们也不考虑h→0时方法的渐近稳定性 例如，对最简单的Euler法 (7-33) 用其求解模型方程（7-32）得到 当un有舍入误差时，其近似解为 ，从而有 取 ，得到误差传播方程 记 ，只要 都不会恶性发展，此时方法绝对稳定。 ，则显式Euler方法的解和误差 从 可得 即 时， (-1,0)为圆心，1为半径的单位圆。 又由于实数μ＜0， （7-33）绝对稳定， 若μ为复数，在 的复平面上，则 表示为以 绝对稳定 区域 绝对稳定区间 定义7.2 一个数值方法用于求解模型问题（7-32），若在 平面中的某一区域D中方法都是绝对稳定的，而在区域D外，方法 是不稳定的，则称D是方法的 绝对稳定区域； 绝对稳定区间。 它与实轴的交称为 例如，显式Euler方法的 绝对稳定区域、区间。如图 现在考察多步法(7-24)，将它用于解模型方程(7-32) 得到k阶线性差分方程 (7-34) 若取 ，则记（7-34）的特征方程为 (7-35) 其中 由k阶线性差分方程的性质我们可以得到如下结论， 区域： 例如，对于k=1时，考虑隐式方法中最简单的后退Euler法 方程（7-35）的根都在单位圆内（︱λ︱＜1) ，则线性多步法 （7-4）关于 绝对稳定,其绝对稳定域是复平面 上的 其特征方程为： 若特征 得 当 时， 故 就是 隐式Euler法的绝对稳定区域。 当μ＜0为实数时，绝对稳定区间为 (-∞,0)。 平面上以（1.0）为圆心的单位圆外区域。 它是 当Re μ＜0时，它位于