Using loop for calculating ode with different parameter

2 次查看（过去 30 天）

Martin Pecha 2017-11-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366169-using-loop-for-calculating-ode-with-different-parameter

评论： Steven Lord 2017-11-10

How can I create a cycle, where I would use ode45 and be able to get a solution for 3 different values of a constant k so I don't have to change it manually.

example k=1 --> solution, then k=10 --> solution etc.

At the end, I would need to put all solutions into the same plot to compare.

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

请先登录，再进行评论。

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

回答（1 个）

Eric 2017-11-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/366169-using-loop-for-calculating-ode-with-different-parameter#answer_290289

在 MATLAB Online 中打开

Use a for loop. If k is non-linear, just put your values of interest into a vector.

 loop_i = 0;
 for k = [1 10 ....Rest of your numbers];
     loop_i = loop_i+1;
     [t{loop_i}, y{loop_i}] = ode45(...Whatever involving k);
 end

I have to admit I don't know what size outputs the function ode45() will return, so sticking them in cells was simply a safe bet for me. If you always expect the same size vector to be returned, you could use y(:,loop_i) instead, or y(:,:,loop_i), etc. If y(:,loop_i) works for you, plotting will be as easy as plot(t,y).

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

Martin Pecha 2017-11-10

Thanks, this works perfectly, however new problem occured. I have defined k=1 in function. when I use loop it give 3 same results and not taking in account the k = [1 10 100] in loop. How can I change it?

Steven Lord 2017-11-10

Pass k into your ODE function (the first input to the ODE solver) as an additional parameter. The description of the odefun input argument in the documentation for ode45 contains a link "Parameterizing Functions" that describes various techniques you can use to do that.

请先登录，再进行评论。

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