2012数值分析 -lec24--线性多步法.pptVIP

多步法的绝对稳定性和绝对稳定域 将线性多步法的公式应用到试验方程 进行考察。这里假定Re 0,即试验方程本身是稳定的。 记 得 是常系数差分方程，其特征方程为 记它的 k个特征根为 并设它们是互异的。显然根与 有关，不妨记为 注意到当 时 这时由特征方程得 由线性多步法的相容性条件得 是一个根。不妨设, 差分方程的解为 其中系数由线性多步法的出发值确定。 另一方面，y(0)=1的试验方程的精确解为 ， 设多步法截断误差为 ，由此可得 ， 我们称 为主根，其它根都为增根。 定义：线性多步法的绝对稳定区域 对给定的 ，如果特征方程的特征根 皆按模小于1，则线性多步法关于u是绝对稳定的。使得 成立的 构成绝对稳定区域。 注：从误差角度来看绝对稳定区域的方法是一个理想的方法。这样，绝对稳定区域越大，方法适用性越广，因而越优越。 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期三下午2:00－5:00 答疑地点：主216 * 第二十四讲线性多步法 第七章 常微分方程初值问题的数值解法 * 线性多步法 数值格式的构造 误差估计（收敛性分析） 稳定性分析 * 插值法构造数值格式 * * Adams 外插法 (k=2) 3阶3步 显式 表 4 外插系数bki值 图 3阶3步外插法 * * Adams 内插法 (k=2) 4阶3步 图 7 4步3阶Adams内插公式 * * * 图 8 一般化插值形式 * * * * 定义 　　　若某算法在计算过程中任一步产生的误差在以后的计算中都逐步衰减，则称该算法是绝对稳定的 /*absolutely stable */。 一般分析时为简单起见，只考虑试验方程 /* test equation */ 常数，可以是复数 当步长取为 h 时，将某算法应用于上式，并假设只在初值产生误差 ，则若此误差以后逐步衰减，就称该算法相对于 绝对稳定， 的全体构成绝对稳定区域。我们称算法A 比算法B 稳定，就是指 A 的绝对稳定区域比 B 的大。 h l h = h 例：考察显式欧拉法 由此可见，要保证初始误差?0 以后逐步衰减， 必须满足： 0 - 1 - 2 Re Img 例：考察隐式欧拉法 可见绝对稳定区域为： 2 1 0 Re Img 注：一般来说，隐式欧拉法的绝对稳定性比同阶的显式法的好。 * 线性多步法的绝对稳定性： * 定义： 绝对稳定。 绝对稳定区域。 * Milne * 表 6 计算结果 * 预估--校正方法 （Predictor--Corrector Method） * 预估--校正方法 （Predictor--Corrector Method） * 注意：一步校正的计算量 ? 预估计算量。 所以要适当选取h才能发挥PC的优点。 设绝对稳定区域： 达到精度的校正次数为N，则h的选取， 应满足： 否则可用N步显式算法稳定达到目的。 h 作业 课后习题13、15、17、18、20、21、22 * * Thanks for your attention! Any questions?! i 0 1 2 3 4 5 b0i 2 b1i 12 b2i 24 b3i 720 b4i 1440 b5i 1 3 23 55 1901 4277 -1 -16 -59 -2774 -7923 5 37 2616 9982 -