1.官方回答
(1).[t,y] = ode45(odefun,tspan,y0)
(其中 tspan = [t0 tf]
)求微分方程组 y′=f(t,y) 从 t0
到 tf
的积分,初始条件为 y0
。解数组 y
中的每一行都与列向量 t
中返回的值相对应。
所有 MATLAB® ODE 求解器都可以解算 y′=f(t,y) 形式的方程组,或涉及质量矩阵 M(t,y)y′=f(t,y) 的问题。求解器都使用类似的语法。ode23s
求解器只能解算质量矩阵为常量的问题。ode15s
和 ode23t
可以解算具有奇异质量矩阵的问题,称为微分代数方程 (DAE)。使用 odeset 的 Mass
选项指定质量矩阵。
ode45
是一个通用型 ODE 求解器,是您解算大多数问题时的首选。但是,对于刚性问题或需要较高准确性的问题,其他 ODE 求解器可能更适合。