How to use if statement with ode23
5 次查看(过去 30 天)
显示 更早的评论
I have a formula in the form of:
dy/dt = f - a*y
f(t) = x-k*t
I can use the ode23 function to solve this for both:
x = 2;
a = 3;
x1 = 4;
a1 = 5;
Tb = 0;
Te = 0.6;
Ti = 0.15;
F = @(t,y) (f*t)-(a*y);
[t,y] = ode23(F,[0 0.6],0);
F1 = @(t,y) (f*t2)-(k*y2);
[t2,y2] = ode23(F1,[0 0.6],0);
When f(t) is between min en max use formula 1, if f(t) is between max and min use formula 2.
Hope someone can help me out here.
回答(1 个)
Jan
2017-9-15
ode23 can handle smooth functions only. See Answers: Discontinuities in ODE45 . The only reliable solution is an event function, which stops the integration and restarts it with the current values and the new function. See https://www.mathworks.com/help/matlab/math/ode-event-location.html
F1 = @(t,y) (f*t)-(a*y);
F2 = @(t,y) (f*t2)-(k*y2);
opt = odeset('Events', @(t,y) DetectEvent(t, y, minValue, maxValue));
t = Tb;
y0 = 0;
while t < Te
if "y is inside the interval" % (This is not clear in your description)
F = F1;
else
F = F2;
end
[tt, yy] = ode23(F, [t, Tb], y0, opt);
t = tt(end);
y0 = yy(end);
... some code to collect the yy values
end
6 个评论
m4 Chrennikov
2017-9-15
Hmmm, my fault. The only reason I want to use sol struct is that I use deval function to process ode solution. If you inspect its source code you can see that it utilizes not only data mentioned above but some 3d matrix sol.idata.f3d. My idea for next day was to try to concat this this matrix if there's no out-of-box solution. Hope I was clear enough)
m4 Chrennikov
2017-9-20
works for me:
function sol1 = concatSol( sol1, sol2 )
%CONCATSOL Summary of this function goes here
% Detailed explanation goes here
if ~strcmp(sol1.solver,'ode45')
error('ode45 solver required')
end
if ~strcmp(sol1.solver,sol2.solver)
error('sol1 & sol2 sovers are not the same')
end
sol1.x = [sol1.x, sol2.x(2:end)];
sol1.y = [sol1.y, sol2.y(:,2:end)];
sol1.xe = [sol1.xe, sol2.xe];
sol1.ye = [sol1.ye, sol2.ye];
sol1.ie = [sol1.ie, sol2.ie];
sol1.stats.nsteps = sol1.stats.nsteps + sol2.stats.nsteps;
sol1.stats.nfailed = sol1.stats.nfailed + sol2.stats.nfailed;
sol1.stats.nfevals = sol1.stats.nfevals + sol2.stats.nfevals;
sol1.idata.f3d = cat(3,sol1.idata.f3d,sol2.idata.f3d(:,:,2:end));
end
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!