四阶龙格库塔法c语言程序,四阶龙格_库塔算法的C语言实现_毋玉芝

2001年3月焦作大学学报 1 第1期JOURNALOFJIAOZUOUNIVERSITYMar.2001

四阶龙格 库塔算法的C语言实现

毋玉芝

(焦作财会学校)

摘 要 本文叙述了四阶龙格 库塔算法的C语言实现过程、数据存储及其结果的曲线显示,并以具体实例说明了这一过程。

关键词 龙格 库塔算法 数据存储 曲线显示

在科学技术中常常需要求解常微分方程的定解问题,这就需要一种合适的数值解法求出常微分方程的解。在诸多数值算法中龙格 库塔算法具有较高的精确度,是一种优先选取的算法。

1.四阶龙格 库塔算法简述

龙格 库塔方法实际上是间接地使用台劳级数法的一种技术。

龙格 库塔算法的数学描述如下:

yn+1=yn+h (K1+2K2+2K3+K4)/6;

K1=f(xn,yn);

K2=f(xn+h/2,yn+K1 h/2);

K3=f(xn+h/2,yn+K2 h/2);

K4=f(xn+h,yn+K3 h);

其中: h表示计算过程中选取的步长;K1表示xn点处的斜率;

K2表示利用K1求得的(xn+h/2)点处的斜率;

K3表示利用K2求得的(xn+h/2)点处的斜率;

K4表示利用K3求得的(xn+h)点处的斜率;

2.C语言的实现过程

2 1 算法实现

龙格 库塔算法关键是选择步长h,必须根据题目的要求选出合适的步长,这对龙格 库塔算法结果的精确度及其平滑性尤为重要。在选择了恰当的步长后,利用上述迭代表达式,并根据题目要求的迭代次数,或求解的精度,利用C语言加以实现。

2 2 数据存储

由于此计算结果数据庞大,程序运行后数据不可能一屏显示,鉴于此首先利用fopen()函数创建并打开一文本文件,利用fprintf()函数将数据存储到数据文件,可将此文件打印输出,以检验结果的正确性。实现过程如下:

if(fp==NULL)

{

printf( Can topenthisfile\n );

exit(0);

}

for(i=0;i<=j;i++)

{