经典Runge-Kutta方法和显式四阶Adams方法 .pdfVIP

微分方程数值解实验报告

姓名丁建伟学号200708020211日期2010．10．23

实验项目经典Runge-Kutta方法和显式四阶Adams指导教师徐强

方法

一、上机实验的问题和要求（需求分析）：

λx

考虑一阶常微分方程初值问题dy/dx=λy,0=x=1,y(0)=1，其精确解为y=e，求初

值问题数值解。

要求：（1）λ=-2时，分别以h=1/16,1/32,1/64,1/128为步长时，使用两种方法求解初值问

题，给出右端点的误差，分析误差阶。



（2）λ=-100时，取h=1/16，用经典Runge-Kutta法求解初值问题，并验证h(2.785,0)

时数值不稳定性。

二、程序设计的基本思想，原理和算法描述：

算法的基本思想：

若已知dy/dx=f（x，y），0=x=1;y(0)=1;

经典Runge-Kutta法公式：

h





yyk2k2kk

k1k61234



kf(x,y)

1kk





1h

kf(xh,yk)

2k2k21

1h

kf(xh,yk)

3kk2

22

kf(xh,yhk)

4kk3



显性四阶Adams公式：

h



yy55f59f37f9f

k1k24kk1k2k3

算法描述：先设置两个数组y[]和Y[]，用来存储结果。

经典Runge-Kutta算法：

（1）置h=1/n；x=0，y=1；f=my

（2）设置一个Runge-Kutta函数，每一步求出对应的x和y值，共循环进行n次，并

每次将结果存储到一个数组里y[]。

（3）将结果输出，并求在右端点处截断误差。

显性四阶Adams：

（1）从y0到y3