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 另一方面，按自然的方法对步长进行减半，即一旦误差大于6， 就对步长进行减半 ． 第二种氖略： 在每一步计算后， 采用如下通 用公 式确定下一步的步长： h 仁 0.9h信）志 其中p是一对Runge-K ut ta方法中笫一个公 式的阶。这个公式既可能培加步长， 也可徒减小步长 四、实验题目： ． 应用RK F45 或RKF54方法分 设 计实现 自适应方法， 求解如下常微分方程 初值问题： { y ＝ 沪 ＋ cos(y — x) y{l)::::: 3 ． 初值步长 取为h 0.01. 在自适 应方法中步长的 选 取采用笫 二种策咯＆ ． 在郁溢出前 终 止口 ． 程序的 输出： O 解的范围：［1 ？ ］ ＠ 提示捡八一个介千上述范围的值， 应用 茼 单的两点 线性插值计算出对拉的函敖值。 这里，我采用了RKF54 方法，容限§=0.000001，步长为 100。 五、问题的解： 编写程序（程序见后面附录），输出结果如下： 111213123456? 89 11 12 13 1 2 3 4 5 6 ? 8 9 i0 ￡(1.010000000000)=3201 ￡(1.024125?11243)=3.6808409138?1 f (1. 036465 797553) iiii4. 468515751528 13 f(i.045608919668)=10.016828439199 44 f(i.045608919668)=10.016828439199 45 f(i.045608919668)=10.016828439199 16 f (i . 04560 891966 8 ) =1” 一” 1 6828 4391 99 f (i. 04225??66 264 =5. 40:809:246!:45:1 ￡1.0444B9922224)=6.435956013431 4? ￡(1.045608919668)=10.016828439199 48 ￡(1.045608919668)=10.016828439199 f(1.04525f.48?952)=?.491203044544 f (1. 045501181985 =8. 48076562576? 19 f(i.045608919668)=10.016828439199 50 f(1 -04 560 891966 8 ) ＝1 0 _ 016828439199 f (1. 045576241496)iiii9.259112207008 51 ￡(1.045608919668)=10.016828439199 f(i.0455990i2868)=9.7i82 9004246 ￡1.045605916318)=9.916093715807 52 f(i.045608919668)=10.016828439199 53 f(i.045608919668)=10.016828439199 ￡(1.0456080091?1)=9.9951!;3394164 f(i.045608643641)=10.00?114498196 54 ￡(1.045608919668)=10.016828439199 55 ￡( 1 . 04560 891966 8 )1 0