第2章常微分方程初值问题的解法.pptVIP

31 § 6.3 线性多步法 思想 : 用若干节点处的 y 及 y 值的线性组合来近似 y ( x n +1 ) 。 其通式可写为： 1 0 1 1 1 1 0 1 1 ( ) n n n k n k n n n k n k y y y y h f f f f + - - - + - - = a + a + + a + b + b + b + + b L L ) , ( n n n y x f f ? 当 ? ? 1 ? 0 时，为 隐式公式 ; ? ? 1 =0 则为显式公式。 Adams 方法是线性多步法的一个代表，它是利用插值多项式进 行积分得出来的，这样构造线性多步法的方法称为数值求积法， 它是构造线性多步法的一种途径，另外还有 Taylor 法 。 32 ? 构造线性多步法的途径： 1 数值积分法 显式 Adams 方法 隐式 Adams 方法 其它方法 ? Simpson 方法（二步法） ? ? ? ? ? 1 )) ( , ( ) ( ) ( 1 k k x x k k dx x y x f x y x y 方程 (1.1) ) , ( y x f dx dy ? ) , ( 1 ? k k x x 在 上积分，即 由 Simpson 求积公式，得 ) ( 90 1 ))] ( , ( )) ( , ( 4 )) ( , ( [ 3 ) ( ) ( ) 5 ( 5 1 1 1 1 1 k k k k k k k k k y h x y x f x y x f x y x f h x y x y ? ? ? ? ? ? ? ? ? ? ? 用 近似 ，得 Simpson 方法 ) ( ), ( 1 k k x y x y ? k k y y , 1 ? ) 4 ( 3 1 1 1 ? ? ? ? ? ? ? k k k k k f f f h y y 局部截断误差为 ) ( ) ( 90 6 ) 5 ( 5 1 h O x y h k k ? ? ? ? ? 结论： Simpson 方法是四阶隐式方法。 例 ] ) 2 ( 4 ) ( ) ( [ 6 ) ( 1 1 1 ? ? ? ? ? ? ? ? i i x x i i i i x x f x f x f h dx x f 33 ? Milne 法 ? ? ? ? ? ? ? 1 3 )) ( , ( ) ( ) ( 3 1 k k x x k k dx x y x f x y x y 方程 ) , ( 1 3 ? ? k k x x 在 ) , ( y x f dx dy ? 上积分，即 用过点 的插值多项式近似 ) , ( ), , ( ), , ( 2 2 1 1 ? ? ? ? k k k k k k f x f x f x )), ( , ( x y x f 得 Milne 方法 ) 2 2 ( 3 4 2 1 3 1 ? ? ? ? ? ? ? ? k k k k k f f f h y y 局部截断误差为 ) ( ) ( 45 14 6 ) 5 ( 5 1 h O x y h k k ? ? ? ? 结论： Milne 方法是四阶显式线性多步法。 1 计 算 方 法 华中科技大学数学与统计学院 第六章 常微分方程初值 问题的数值解法 计算方法课程组 2 § 6.1 基本离散方法 § 6.2 Runge-Kutta 方法 § 6.3 线性多步法 § 6.4 收敛性与稳定性 § 6 常微分方程数值解法 3 考虑一阶常微分方程的初值问题 : 例如： , [ , ] ( ) ? ? ? ? ? ? 0 1 0 1 y = x + y x y 其解析解为： ? ? ? ? 1 2 x y x e § 6.1 基本离散方法 ? ? ? ? ? ? ? ? 0 ) ( ] , [ ) , ( y a y b a x y x f dx dy 4 [ , ] x dt x ? ? ? ? 0 1 0 1 2 -t y e 但是 , 只有一些特殊类型的微分方程问题能够得到用解 析表达式表示的函数解，而大量的