ode15s - accept solution
显示 更早的评论
Hello,
I solve an ODE system with the ode15s solver.
The system has very big positiv and negative gradients. If a gradient is negative and the timestep is too big, my solution gets negative.
However, this result is not physicaly correct and my next timestep breaks down because I like to calculate with this negative temperature other values. Normaly the solver "sees" that the timestep was to big and adjusts it self. But in this case (the gradient was to negative) a other function returns an error because the negative temperature and the hole calculatin stops.
How can I modify my code that I don't come to this case? I know that my model is OK, just the automatic time step from the ODE solver is to big.
Thanks for help.
采纳的回答
Torsten
2015-11-30
Use something like
T=y(1);
if T < 0
T = 1e-3;
end
in the routine where you supply your gradients.
Best wishes
Torsten.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Translated by 
