Matlab ODE求解器的多个输出

Question

我有以下Matlab ODE代码：

[t,y,~,~,ie] = ode23tb(@(t,y) RHSODE(t,y),[0,t_end], [i0;v0],options);

我希望ODE求解器也可以给我结果z，它是y和dy / dt的函数，因此z = f（y，dy / dt）。

有谁知道如何将这样的z加到求解器的输出中？

Answer 1

有两种方法可以做到这一点。 最常见且通常最快的方法是利用集成功能（在您的情况下为RHSODE ）并在执行集成后评估功能f 。 您的代码中并未提供许多详细信息，但可能看起来像这样：

ydot = RHSODE(t,y);
z = f(y,ydot);

其中t和y是ode23tb的输出。 这要求将RHSODE和f进行向量化（或者可以将以上内容包装在for循环中）。

另一种方法要求您在积分函数RHSODE内创建一个附加方程式（如果z是矢量，则创建方程式）。 通常， ode23tb在此函数中集成了任何内容，因此必须将f乘以t才能取消此功能。 同样，您的代码可能看起来像这样：

function ydot = RHSODE(t,y)
ydot0 = ... % Your original ODE(s)
z = f(y,ydot);
ydot = [ydot0;z*t]; % Make column vector