running MATLAB ode45 back-to-back
显示 更早的评论
Hi All,
I am curious what happens with the first solution point when we run ode45 to integrate a right hand side function. The first state solution that is returned, is it the same as the initial condition? I would assume so because the associated time point for the first solution always seems to be a zero.
Next, suppose I am to run ode45 back-to-back because a physical phenomena abruptly changes. e.g. closing a valve in gas filled tank. In this case, a final solution for the previous ode integration becomes the initial condition for the next ode integration, i.e. successive iteration. Now, for the sake of continuity in time, I want to combine the solutions from two different ode integrations. What is the best way to do this? Should we remove the initial solution from the second integration and combine with the first solution? Also, for the corresponding time, should we consider the first solution time point of the second iteration to be exactly the same as the ending time point of the first integration?
Thank you so much!
-Taehun
采纳的回答
Star Strider
2021-1-18
1 个投票
‘The first state solution that is returned, is it the same as the initial condition?’
In my experience, yes. Sometimes it is the only solution if there are problems with the code and the other values are all not finite (NaN or ±Inf).
‘Now, for the sake of continuity in time, I want to combine the solutions from two different ode integrations. What is the best way to do this?’
I vertically concatenate the time vectors, and vertiocally concatenate the integrated solution matrices.
‘Should we remove the initial solution from the second integration and combine with the first solution?’
I generally do not, since the last row of the preceding integration likely overplots the first row of the next integration. It probably does not make any difference.
‘Also, for the corresponding time, should we consider the first solution time point of the second iteration to be exactly the same as the ending time point of the first integration?’
It most likely will be. You can always check to be sure.
.
更多回答（1 个）
Bjorn Gustavsson
2021-1-18
Well if you have a problem like this (improvising for simplicity)
ode1 = @(t,y,alpha) alpha*y; % exponential growth
t_span1 = [0,3]
t_span2 = [3,12];
y0 = 1;
[t1,y1] = ode45(@(t,y) ode1(t,y,1),t_span1,y0);
[t2,y2] = ode45(@(t,y) ode1(t,y,-1/2),t_span2,y1(end));
t_both = [t1;t2];
y_both = [y1;y2];
You should be fine, if you want you can use unique to remove identical rows:
t_y = unique([t_both,y_both],'rows');
to get rid of the duplicate points. Or you can just remove either the last of the first or the first of the second outputs when concatenating.
HTH
类别
在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
