How can I put condition in ode solver
显示 更早的评论
I'm trying to put if condition in ode solver.
but it seems that I can't control variables in ode function. (can control just derivative of variables)
for example
function dcdt = exo(t,c)
dcdt = zeros(3,1)
dcdt(1) = -c(1)*exp(-1/c(3));
if c(1) < 0
c(1) = 0; ----->>> this doesn't work!!!
end
dcdt(2) = -c(2)*exp(-1/c(3));
dcdt(3) = - dcdt(1) + dcdt(2);
end
How can I change this function to make c(1) = 0 if c(1) < 0
采纳的回答
Walter Roberson
2020-12-31
0 个投票
The above solution is really only appropriate if the reason that c(1) becomes negative is due to round-off or overshooting (that is, a smaller step size would not have gone negative.)
If you are trying to do the equivalent of bounding off a floor, then it is not appropriate to use this solution, and instead you should use an ode event function to terminate the integration, then change the boundary condition, and then restart the integration.
更多回答（1 个）
Steven Lord
2020-12-31
0 个投票
Use odeset to create an options structure with the NonNegative option set to the elements of the solution that need to remain non-negative and pass that options structure into the ODE solver.
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
