数值方法 常微分方程初值问题的数值解法.doc

(Runge-Kutta 第三节 龙格-库塔(Runge-Kutta方法 二、四阶经典的R-K公式 h ? ? y i + 1 = y i + ( K 1 + 2 K 2 + 2 K 3 + K 4 , 6 ? ? K 1 = f ( x i , y i , ? h h ? ? K 2 = f ( x i + , y i + K 1 ， i = 0 ,1 , ? , n ? 1 . 2 2 ? ? h h ? K 3 = f ( x i + 2 , y i + 2 K 2 ， ? ? K 4 = f ( x i + h , y i + hK 3 ， ? 21/24 (Euler 第二节 欧拉(Euler (Euler方法 例3 用经典的四阶R-K公式计算: 2x ? ′= y? , 0 ≤ x ≤ 1, ?y y ? ? y ( 0 = 1 ? 的数值解，取步长 h = 0.2, 并与精确值比较 . 解: 0.2 ? ? yi +1 = yi + 6 ( K 1 + 2 K 2 + 2 K 3 + K 4 , ? ? ? K 1 = yi ? 2 x i , ? yi ? 22/24 (Euler 第二节 欧拉(Euler (Euler方法 ? ?K2 ? ? ?K3 ? ? ?K4 ? = = = x i + 0 .1 y i + 0 .1 K 1 ? 2 , y i + 0 .1 K 1 x i + 0 .1 y i + 0 .1 K 2 ? 2 , i = 0,1, ? ,4. y i + 0 .1 K 2 x i + 0 .2 y i + 0 .2 K 3 ? 2 , y i + 0 .2 K 3 2 x0 2× 0 (1 i = 0 : K 1 = y0 ? = 1? = 1, ? .见P174. y0 1 23/24 第九章 学习要求及作业 1、掌握求解一阶初值问题的欧拉公式, 隐式欧拉公式和欧拉中点公式 ； 2、掌握求解一阶初值问题的梯形公式 及其预估—校正法;　 ３、理解构造龙格—库塔法的思想; 掌握二阶R-K公式的中点公式; 4、掌握四阶经典的R-K公式. 作业：P326：1，3, 4. 24/24