四阶龙格—库塔算法及其C-语言程序设计论文.docVIP

包头师范学院 2012 -2013学年第一学期期末考试课试卷 考试科目 计算物理 成绩 院系 物理学院 专业 10 级 本一 姓名 张艳有 任课教师签名： 院系负责人签名： 四阶龙格—库塔算法及C 语言程序设计 物理科学与技术学院 10级物理系一班 1009320036 张艳有 摘　要　本文叙述了四阶龙格—库塔算法的C 语言实现过程、数据存储及其结果并以具体实例说明了这一过程。 关键词　龙格—库塔算法　数据存储 例题 程序设计 计算结果 引言 在科学技术中常常需要求解常微分方程的定解问题, 这就需要一种合适的数值解法求出常微分方程的解。在诸多数值算法中龙格—库塔算法具有较高的精确度, 是一种优先选取的算法。 具体内容 一、四阶龙格—库塔算法简述 从常微分方程数值解法的几何意义看，欧拉方法取一点处的斜率作为平均斜率，因此欧拉方法近似公式为 向后的欧拉方法则采用点处的斜率作为平均斜率，即 所以这两种方法也称作矩形法。改进的欧拉方法则取点处和点处斜率 和的平均值作为平均斜率，即 因此改进的欧拉方法又称为梯形方法。可以预见，若去多点处斜率的加权平均值作为平均斜率，误差会更小，这就是龙格-库塔方法。 最常用的是四阶龙格库卡近似计算公式，即： 式中 其中: 　 表示计算过程中选取的步长; K1 表示点处的斜率; 　　 K2 表示利用K1 求得的 点处的斜率; 　　 K3 表示利用K2 求得的点处的斜率; 　　 K4 表示利用K3 求得的 点处的斜率; 二、使用标准四阶龙格——库塔法求解初值问题实例 1、用标准四阶龙格——库塔法求解 　???????? 　解：取步长， 2、计算所用程序如下所示 #include stdio.h #include conio.h float func(float x,float y) { return(2*x*y); } float runge_kutta(float x0,float xn,float y0,int n) { float x,y,y1,y2,h,xh; float d1,d2,d3,d4; int i; x=x0; y=y0; h=(xn-x0)/n; for(i=1;i=n;i++) { xh=x+h/2; d1=func(x,y); d2=func(xh,y+h*d1/2.0); d3=func(xh,y+h*d2/2.0); d4=func(xh,y+h*d3); y=y+h*(d1+2*d2+2*d3+d4)/6.0; x=x0+i*h; } return(y); } void main() { float x0,xn,y0,e; int n; printf(\ninput n:\n); scanf(%d,n); printf(input x0,xn:\n); scanf(%f%f,x0,xn); printf(input y0:\n); scanf(%f,y0); e=runge_kutta(x0,xn,y0,n); printf(y(%f)=%6.6f,y0,e); } 3、部分计算结果如图所示 三、 结束语 利用C 语言可方便地实现四阶龙格—库塔算法, 这对于科学技术中所遇到的常微分方程的求解具有重要意义。 参考文献： 物理学与计算物理，物理学与计算机翻译组译，北京：科学出版社，1986 计算机数值计算方法，计九三译，北京清华大学出版社，1987