微分方程组数值求解四阶显式龙格库塔法自编程实现-龙格-库塔法求解常微分方程-Matalab自编程序实现.pdf 常微分方程（组）的数值解法有多种，其中Runge-Kutta方法是最为流行算法之一， 几乎可以解决所有常微分方程组的数值解问题。 论坛中不少朋友执着于利用MATLAB自置函数solve求解微分方程。其实solve函数存在局限性 （比如系统输入是时变的，solve函数就不易处理），也不利于大家掌握算法。 本人在此提供两个微分方程组的数值算例及程序代码，供大家参考（附件）。算法采用显式4阶Runge-Kutta算法。 关于Runge-Kutta算法，大家可以参考一般 数值方法 的教材即可。

自己用matlab语言编写的用于求解微分方程的.m文件

计算速率方程龙格库塔法C语言程序

数值解与理论解对比可知，四阶龙格-库塔法的精度已经很高，用它来解一般常微分方程已经足够了。 有程序运行说明 数值解与理论解对比可知，四阶龙格-库塔法的精度已经很高，用它来解一般常微分方程已经足够了。 有程序运行说明

输入、输出说明 fcn为f(x,y),a,b分别为x的上下限，stepsize为步长,y0为初值。输出为x,及对应的y。

龙格库塔法matlab（2阶与3阶）算例代码与text文件代码一样

基于龙格—库塔法（R-K）Matlab代码实现

通过本次实验，熟悉求解常微分方程初值问题的有关方法和理论，主要是欧拉法、改进欧拉法、四阶龙格库塔法，学会编制这两种方法的计算程序。体会这两种解法的功能、优缺点及适用场合。 解初值问题，并在屏幕上按适当的比例和位置画出坐标轴及解的函数曲线。 matlab7.0以上版本,图形界面,带实验报告

龙格库塔法求二阶常微分方程，改几下也就可以求解n阶常微分方程

数值解与理论解对比可知，四阶龙格-库塔法的精度已经很高，用它来解一般常微分方程已经足够了。 有程序运行说明