数值分析Euler方法讲解.pptVIP

1 第二节 Euler 方法 2 ? 2009, Henan Polytechnic University 2 § 2 Euler 方法 第五章 常微分方程数值解法 5.2.1. Euler 方法 设节点为 x k = x 0 + kh ( h=(b-a)/n k =0,1,… n ) 方法一 泰勒展开法 （将 y(x k+ 1 ) 在 x k 泰勒展开得 ） 2 1 1 1 ) ( 2 ) ( ) )( ( ) ( ) ( k k k k k k k x x y x x x y x y x y ? ? ? ? ? ? ? ? ? ? ? ? 2 ) ( )) ( , ( ) ( 2 ? y h x y x hf x y k k k ? ? ? ? ? n k x y x hf x y x y k k k k ..... 2 , 1 , 0 )) ( , ( ) ( ) ( 1 ? ? ? ? 则可得： 3 ? 2009, Henan Polytechnic University 3 § 2 Euler 方法 第五章 常微分方程数值解法 方法二 数值微分法（用向前差商近似导数） )) ( , ( ) ( k k k k x y x f x y x ? ? 处有 在点 利用数值微分公式 h x y x y x y k k k ) ( ) ( ) ( 1 ? ? ? ? )) ( , ( ) ( ) ( 1 k k k k x y x hf x y x y ? ? ? 得： 4 ? 2009, Henan Polytechnic University 4 § 2 Euler 方法 第五章 常微分方程数值解法 方法三 数值积分法 dx x y x f x y x y k k x x k k ? ? ? ? ? 1 )) ( , ( ) ( ) ( 1 在积分可得： 两端同时 由 ] , [ )), ( , ( ) ( 1 ? ? ? k k x x x y x f x y )) ( , ( )) ( , ( ) ( 1 k k k k k k x y x hf x y x f x x ? ? ? ? )) ( , ( ) ( ) ( 1 k k k k x y x hf x y x y ? ? ? 即： 5 ? 2009, Henan Polytechnic University 5 § 2 Euler 方法 第五章 常微分方程数值解法 ，则可得： 近似结果为 代入上式右端，记所得 的近似值 用 1 ) ( ? k k k y y x y 1 ..... 2 , 1 , 0 ) , ( 1 ? ? ? ? ? n k y x hf y y k k k k 欧拉格式。 此即为欧拉公式，又称 依上述公式逐次计算可得： n y y y , , , 2 1 ? 也称 Euler 为单步法， 又 称为 显格式的单步法 。 )) ( , ( ) ( ) ( 1 k k k k x y x hf x y x y ? ? ? 6 ? 2009, Henan Polytechnic University 6 § 2 Euler 方法 第五章 常微分方程数值解法 2 欧拉法的几何意义： ) ( ) , ( )) ( , ( 0 0 0 0 0 0 x x k y y y x x y x ? ? ? ? 作切线 过 2 1 0 x x x ) , ( )) ( , ( ) ( 0 0 0 0 0 y x f x y x f x y k ? ? ? ? 斜率 求交点， 与 1 0 0 0 0 ) )( , ( x x x x y x f y y ? ? ? ? ， 纵坐标记为 1 y ) , ( 0 0 0 1 y x hf y y ? ? 则 为斜率作直线， 以 过 ) , ( ) , ( 1 1 1 1 y x f y x ) )( , ( 1 1 1 1 x x y x f y y ? ? ? ， 求交点，纵坐标记为 与 2 2 y x x ? ) , ( 1 1 1 2 y x hf y y ? ? 则 ........ .......... .......... ) , ( 1 1 y x ) , ( 2 2 y x 也称 欧拉折线法 . 从上述几何意义上得知，由 Euler 法所得的折线 明显偏离了积分曲线，可见此方法 非常粗糙。 7 ? 2009, Henan Polytechnic University 7 § 2 Eu