你可以解决它DAE(微分代数方程),而不是一首颂歌。
M (t, y) y′= f (t, y)
这种形式的例子是一个特例。它可以解决
ode15s
和
ode23t
解决者。使用这些解决,您可以直接指定
米
矩阵不需要转化它。
基本上,你定义一个函数
函数dx = mySys (t, x)
dx = A2 * x;
结束
然后,使用ode解决之前解决它
ode15s
和
ode23t
,指定
选择= odeset (“质量”,A1);
接下来,应用在使用解算器的选项如下
[t、y] = ode15s (@mySys tspan, y0,选项);
