常微分方程初值问题数值解法 一、问题提出 科学计算中经常遇到微分方程(组)初值问题需要利用Euler法,改进Euler法Rung-Kutta方法求其数值解,诸如以下问题:
常微分方程初值问题數值解法
科学计算中经常遇到微分方程(组)初值问题需要利用Euler法,改进Euler法Rung-Kutta方法求其数值解,诸如以下问题:
(4)利用四阶标准R- K方法求二阶方程初值问题的数值解
使用Euler法求解运算程序简单,但是计算结果准确度不高使用改进的Euler法求解过程相对复杂,但是准确度会更高准确度最高的是四阶龙格库塔法,求解步骤也是最复杂的问题(1)使用Euler求解,并与准确解对比问题(3)使用改进的Euler法求解。问题(4)(I)(IV)使用四届标准龙格库塔法求解
6.解问题(4)(I):
7.解问题(4)(IV):
四、实验数据结果及分析
对比看,结果相差不大参考结果應该是使用四阶龙格库塔方法求解的。
3.解问题(4)(I):
3.解问题(4)(IV):
使用c语言写代码虽然比Matlab的代码长但是还是c语言适合用来做练習,因为用c写更有结构Matlab隐藏了许多细节部分,不适合对方法的深入理解但是写代码速度快,代码更简练
版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。