局部截断误差.PPTVIP

§3 Convergency and Stability §3 Convergency and Stability §3 Convergency and Stability 例：隐式龙格-库塔法 而显式 1~ 4 阶方法的绝对稳定区域为 其中2阶方法 的绝对稳定区域为 0 Re Img k = 1 k = 2 k = 3 k = 4 - 1 - 2 - 3 - - - 1 2 3 Re Img 无条件稳定 注：一般来说，隐式欧拉法的绝对稳定性比同阶的显式法的好。 小结 尤拉两步公式与尤拉单步公式相比，使用两个节点上的已知信息将精度提高一阶。可以设想，计算y(xn+1)时，充分利用前面已经求出的节点上的 y 及 y’ 值的线性组合来近似y(xn+1)，精度会大大提高。----线性多步法的基本思想。 构造线性多步法有多种途径，这里介绍两种： 基于数值积分的构造方法； 基于Taylor展开的构造方法。 §4 线性多步法 /* Multistep Method */ ) ... ( ... 1 1 0 1 1 1 1 0 1 r n r n n n r n r n n n f f f f h y y y y - - + - - - + + + + + + + + + = b b b b a a a 线性多步法的通式可写为： 当 ??1?0 时，为隐式公式; ??1=0 则为显式公式。 ? 基于数值积分的构造法 将 在 上积分，得到 只要近似地算出右边的积分 ，则可通过 近似y(xn+1) 。而选用不同近似式 I，可得到不同的计算公式。例如利用左矩形积分公式得到尤拉公式；梯形积分公式得到梯形公式。 一般地，利用插值原理所建立的一系列数值积分方法也可以导出解微分方程的一系列计算公式。运用插值方法的关键在于选取合适的插值节点。假设已构造出f(x,y(x))的插值多项式Pr(x),则 §4 Multistep Method ? 亚当姆斯显式公式 /* Adams explicit formulae */ 利用r+1 个节点上的被积函数值 构造 r 阶牛顿后插多项式, 有 Newton 插值余项 /*Adams 显式公式 */ 外推技术 /* extrapolation */ table 5-6, p.181 j ?j 0 1 1 1/2 2 5/12 3 3/8 实际计算时，将差分展开 i j r j=i 局部截断误差为： 注：Br 与yn+1 计算公式中 fn , …, fn?k 各项的系数 均可查表得到 。 见下表。 1 0 1 2 3 r fn fn?1 fn?2 fn?3 … Br … … … … … … … 例： 常用的是 r = 3 的4阶亚当姆斯显式公式 Adams显式公式用 作为插值节点，在求积区间[xn, xn+1]上用插值函数Nr(x)近似f(x,y(x))，而xn+1不在插值节点内，因此是一个外推的过程。虽然公式是显式的，便于计算，但效果并不理想，比如稳定性较差等。 因此改用通过r+1个节点 的插值多项式Nr(x)近似f(x,y(x))，由于xn+1是其中一个插值点，因此是内插多项式，但导出的公式是隐式的。 §4 Multistep Method ? 亚当姆斯隐式公式 /* Adams implicit formulae */ 利用r+1 个节点上的被积函数值 fn+1 , fn , …, fn?r+1 构造r 阶牛顿前插多项式。与显式多项式完全类似地可得到一系列隐式公式。 j 0 1 1 -1/2 2 -1/12 3 -1/24 局部截断误差为： 1 0 1 2 3 k fi+1 fi fi?1 fi?2 … Br … … … … … … … ~ 小于Br 注： 与yn+1 计算公式中 fn+1 , …, fn+1?k 各项的系数 均可查表得到 。 见下表。 常用的是 k = 3 的4阶亚当姆斯隐式公式 较同阶显式稳定 例： §4 Multistep Method ? 亚当姆斯预测-校正系统 /* Adams predictor-corrector system */ Step 1: 用Runge-Kutta 法计算前 k 个初值； Step 2: 用A