龙格库塔法例题.pdfVIP

四阶龙格一库塔法 通常所说的龙格一库塔法是指四阶而言的．我们可以仿二阶、三阶的情形推导出常用的标准四阶 龙格一 库塔公式 (9.22) 公式 (9.22) 的局部截断误差的阶为 ． 龙格一库塔法具有精度高，收敛，稳定 (在一定的条件下 )，计算过程中可以改变步长，不需要计算高阶 导数值等优点．但仍需计算 在一些点上的值，如四阶龙格－库塔法每计算一步需要算四次 的值，这给实际计算带来一定的复杂性，因此，多用来计算 “表头 ”．(即开始几点的近似值 ) ． 例 3.用标准龙格一库塔法求初值问题 在 处的解． 解 因 与 ．若应用标准龙格一库塔方法公式 (9.22)计算， 对于 n=0 时，则有 于是得 这个值与准确解 在 处的值 已十分接近． 再对 n=1 ，2 ，3，4 应用式 (9.22)计算，具体计算结果如表 3 所示： 例 3 写出用四阶龙格――库塔法求解初值问题 的计算公式，取步长 h＝0.2 计算 y (0.4) 的近似值。至少保留四位小数。 f x y y 解 此处 ( ， ) ＝8－3 ，四阶龙格――库塔法公式为 其中 1 ＝f ( xk ，yk ) ； 2 ＝f ( xk +0.5 h，yk+0.5 h 1 ) ； f xk h yk h f xk h yk h 3 ＝ ( +0.5 ， +0.5 2) ； 4 ＝ ( + ， + 3) 本例计算公式为： 1 k 2 k 其中 ＝8－3y ； ＝5.6 －2.1 y ； 3 ＝6.32 －2.37 yk ； 4 ＝4.208 －1.578 yk ＝1.2016 ＋0.5494 yk ( k ＝0，1，2，…) 当 x0 ＝0，y0 ＝2， y (0.2) ≈y 1＝1.2016 ＋0.5494 y0 ＝1.2016 ＋0.5494 ×2 ＝2.3004 y (0.4) ≈y2 ＝1.2016 ＋0.5494 y 1 ＝1.2016 ＋0.5494 ×2.3004 ＝2.4654