解微分方程
2001年3月焦作大学学报
JOURNALOFJIAOZUOUNIVERSITY№ 1Mar.2001第1期
四阶龙格一库塔算法的C语言实现
毋玉芝
(焦作财会学校)
摘要本文叙述了四阶龙格一库塔算法的C语言实现过程、数据存储及其结果的曲线显
示,并以具体实例说明了这一过程。
关键词龙格库塔算法数据存储曲线显示
在科学技术中常常需要求解常微分方程的定解问题,这就需要一种合适的数值解法求出常微分方程的解。在诸多数值算法中龙格一库塔算法具有较高的精确度,是一种优先选取的算法。
1.四阶龙格一库塔算法简述
龙格一库塔方法实际上是间接地使用台劳级数法的一种技术。
龙格一库塔算法的数学描述如下:
Y。+l=y。+h (Kl+2K2+2K3+K4)/6;
Kl=厂(z。,Y。);
K2=f(z。+hi2,Y。+K1 h/2);
K3=f(z。+h12,Y。+K2 h/2);
K4=f(X。+h,y。+K3 h);
其中:h表示计算过程中选取的步长;K。表示X。点处的斜率;
K:表示利用K。求得的(z。+h/2)点处的斜率;
K,表示利用K:求得的(.717。+hi2)点处的斜率;
K。表示利用K,求得的(z。+h)点处的斜率;
2.C语言的实现过程
2.1算法实现
龙格一库塔算法关键是选择步长h,必须根据题目的要求选出合适的步长,这对龙格一库塔算法结果的精确度及其平滑性尤为重要。在选择了恰当的步长后,利用上述迭代表达式,并根据题目要求的迭代次数,或求解的精度,利用C语言加以实现。
2.2数据存储
由于此计算结果数据庞大,程序运行后数据不可能一屏显示,鉴于此首先利用fopen()函数创建并打开一文本文件,利用fprintf()函数将数据存储到数据文件,可将此文件打印输出,以检验结果的正确性。实现过程如下:
if(fp==NULL)
{
printf(”Can’t
exit(O);
}
for(i=0;i<=J;i++)
{t=ts+i*h:openthisfile\n,,);