How to get calculated variables from ode function?

Question

Alex 2022-11-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1866713-how-to-get-calculated-variables-from-ode-function

评论： Jan 2022-11-30

Hi to all

I solve ODE for location and velocity of fluid, using ode23s embeded function. The solution is good. My question is how can I extract variables, which are calculated inside the ode23s? For example, how can I exytract fluid acceleration, a?

Thanks in advance.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2022-11-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1866713-how-to-get-calculated-variables-from-ode-function#answer_1115463

编辑：Jan 2022-11-30

在 MATLAB Online 中打开

The function to be integrated computes these values. You have all required inputs as output of ODE23S already. Use them as inputs to your function. Convert the function such, that it accepts matrices as inputs.

tspan = [0, 10];

y0 = [0, 0];

[T, Y] = ode23s(@fcn, tspan, y0);

Z = fcn(T.', Y.').';

plot(T, Y); hold('on');

plot(T, Z(:, 2), 'co');

function dy = fcn(t, y)

dy = [y(2, :); ...

sin(t)]; % [EDITED, T -> t]

end

3 个评论
显示 1更早的评论隐藏 1更早的评论

Bjorn Gustavsson 2022-11-30

Because @Jan gave you an "example differential equation" in fcn, implying (according to my interpretation) that you should modify the system of ODEs you have that you use ode23s to integrate. One point to keep track of here is that you need to take proper care of the orientation of the inputs such that you get a column-array for dy when calling the function at one instance in time with a 1-D array for y, but an [ny x nt] array when you call it with an array of times and a 2-D array for y.

Jan 2022-11-30

@Alex: I've fixed the typo "T" to "t" and my example does run. So if your code (not mine!) produces an error, post this code. Then the readers can suggest a solution.

请先登录，再进行评论。