RKF54方法分析和总结.docxVIP

实验七 一、实验名称 RKF54 方法 二、实验目的与要求： 实验目的: 掌握RKF54 方法和RKF45 方法。 实验要求：1.应用 RKF54 方法或者 RKF45 方法，设计实现自适应方法，求解常微分方程； 2.用 C 语言实现算法，运行环境为 Microsoft Visual C++。 三、算法思路： ． 五阶方法 6沁 ＋ h) = y(x) +L a·; Fi 6 i= l = y(x) 16 _ + 6656 _ —28—S61 _ -—9 _ 2 F6 ＋—13.5F1 注这其中F2项的系数为 零 ． 西阶方 法 12825 F3 + 5643.0F 4 5.0压＋ —55· 6沁 ＋ h) = y(x)，+ 区炳 6 i = 1 = y(x) ＋ 25 E_ 1408 _ 219了 + F3+ F4 1 ----=- － F 5 216 注意其 中巳和F6 项 的 系数 为，·，_, 2565 4104 5 自适应补偿的选择策略： 假设要构造一个自这应算法， 试 囡使局部截斩误差e的位限定在预先给定的容限妒之内 ． 第一种策略： 假设y(x) — y (x) 估计的局部截断误差为 Ch 气当加 倍 步长导致C(2hi)5 5/ 4 时， 对 步长进行加倍是合理的日 因此当局部截断误差小于5/ 128 时 对步长进行加倍g 另一方面，按自