How do I declare and iterate a variable in an ODE driver file?

Question

Miraculous Mongoose 2019-9-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479497-how-do-i-declare-and-iterate-a-variable-in-an-ode-driver-file

评论： Miraculous Mongoose 2019-9-9

I have a pair of function and ode driver files. Currently, all values except initial conditions are declared as variables in the function file.

I would like to iterate on one variable, and check if the result is better than the previous value.

For example,

start with Snh_in_1=1,
run the ODE,
print y(end,3).
+0.1 to Ssnh (i.e. Snh_in_1 = 1.1)
run the ODE,
print new y(end,3).
if new y(end,3), > original y(end, 3), repeat from step 1.
stop when new y(end,3) less than or equal to previous y(end,3), and report this value of Snh_in_1.

These are my code so far

Function file: https://pastebin.com/Mp43ea1v

Driver file: https://pastebin.com/ct8iNL8W

I can imagine an if loop that might work:

result2 = 2;
result1 = 1;
Ssnh_in_1 =  1;
if result2 >= result1
    [t,y]=ode15s(@pha_cont,[0,15],y0);
    result1 = y(end,3);
    Ssnh_in_1 = Ssnh_in_1 + 0.1;
    [t,y]=ode15s(@pha_cont,[0,15],y0);
    result2 = y(end,3);
else
    disp(Ssnh_in_1)
end

But I have no idea how to pass this variable from the driver into the function file.

How can I do what I want in Matlab? Thanks in advance!

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

请先登录，再进行评论。

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

Answer 1

Torsten 2019-9-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479497-how-do-i-declare-and-iterate-a-variable-in-an-ode-driver-file#answer_391004

yold = -Inf;
Ssnh_in_1 =  1;
[t,y]=ode15s(@(t,y)pha_cont(t,y,Ssnh_in_1),[0,15],y0);
ynew = y(end,3);
while yold < ynew
 Ssnh_in_1 = Ssnh_in_1 + 0.1;
 [t,y]=ode15s(@(t,y)pha_cont(t,y,Ssnh_in_1),[0,15],y0);
 yold = ynew;
 ynew = y(end,3);
end
disp(Ssnh_in_1) 

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

Miraculous Mongoose 2019-9-9

Thank you! learned 2 new things today :)

请先登录，再进行评论。

Answer 2

Bjorn Gustavsson 2019-9-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479497-how-do-i-declare-and-iterate-a-variable-in-an-ode-driver-file#answer_390969

在 MATLAB Online 中打开

The way I read your ODE-function your Snh_in_1 is the initial value of the second component of y. If so it should be no more difficult than something like this:

Snh_in_1 = 1;
dSnh_in_1 = 0.1;
for i_iter = 1:10
  y0 = [a0,Snh_in_1+(i_iter-1)*dSnh_in_1,c0,d0,e0]; % I lost track of how many components you used
  [t,y] = ode15s(@pha_cont,[0,15],y0);
  results(i_iter,:) = y(end,:);
  Results{i_iter,1} = t; % If you've integrated all your ODE's you might just
  Results{i_iter,2} = y; % as well save the entire solution - if it doesn't become too big...
end
% Then do the rest of the analysis after all the ODE-integrations