matlab-常微分方程数值解法.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 3 改进的欧拉方法 取两点斜率平均值，即： 整体截断误差： O(h2) 第三十页，编辑于星期一：二十三点 三十分。 欧拉方法（前、后）和改进的欧拉方法优点是算法简单，但计算精度较低，不能满足实际问题求解对精度的要求，所以使用较少。 龙格-库塔(R-K)方法就是一种精度较高的较为实用计算常微分方程的方法。 第三十一页，编辑于星期一：二十三点 三十分。 从以上三种方法的比较，可以推测：若取多点处斜率的加权平均作为平均斜率，精度会更高，这就是龙格—库塔法的基本思想。 第三十二页，编辑于星期一：二十三点 三十分。 四阶龙格—库塔法计算公式为： 整体截断误差： O(h4) 第三十三页，编辑于星期一：二十三点 三十分。 第三十四页，编辑于星期一：二十三点 三十分。 例题： time_Euler = 0.2500 time_EulerPro = 0.5000 time_RK4 = 1.0150 第三十五页，编辑于星期一：二十三点 三十分。 4、 MATLAB的常微分方程函数 ode45 ode23 ode113 ode15s ode23s ode23t ode23tb 格式 [x