matlab_常微分方程数值解法【爆款】.pptVIP

欧拉方法（前、后）和改进的欧拉方法优点是算法简单，但计算精度较低，不能满足实际问题求解对精度的要求，所以使用较少。 龙格-库塔(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,y]=ode45(′ fun′, [x0,xn], y0,option] 说明： 适用于求解一阶常微分方程组 fun定义微分方程组的函数文件名 [x0,xn]求解区域 y0初始条件向量 option可选参数，由ODESET函数设置，比较复杂 x输出自变量向量，y输出[y, y ′, y ″,..] .精品课件. * 没有一种算法可以有效地解决所有的 ODE 问题，因此MATLAB 提供了多种ODE函数。