三阶微分方程的数值解法

matlab的旋律

ode45是通常情况下最好的单步解法，次解法基于Dormand-Prince 4-5阶的龙格-库塔公式，是simulink中默认的解法器，通常也是第一选择，是解决是解决连续柔性问题的最理想解法器。 ode23是基于Bogacki-shampine 2-3阶龙格-库塔公式，有时对轻度的刚性方程，它比ode45更有效，对于相同的精度，它需要比ode45更小的步长。 ode113是变阶次Adams-Bashforth-Moulton解法，此揭发使用前几次节点上的值来计算当前节点处的解，因此在相同的精度下，它可能比ode45或ode23更快一些，但此解法不适合用于不连续的系统。 根据需要可以将代码中的ode113换成其他的函数。