species go to negative values
2 次查看(过去 30 天)
显示 更早的评论
Hi
I developed a model in Simbiology simulating the effect of a drug on a tissue mass. Sometime, I received a message that the integration toerance does not met and other times I had a mass species with negative values.Why does the software allow the value of a species representing the mass of this tissue go to negative values? what is the cause of this issue and how can I solve it? I want the content of this species to go as small as it can go without jumping to negative values.
Thank you
2 个评论
emjey
2024-5-14
You may find this discussion useful https://uk.mathworks.com/matlabcentral/answers/2002822-modeling-dissociation-using-hill-function-provides-misleading-results?s_tid=srchtitle
In general it is not impossible to get negative values as solutions. Imagine for example zero-order degradation in your system, it will lead to negative solutions eventually.
To help you with your specific issue it would be good to have a concrete example.
回答(2 个)
John D'Errico
2024-5-14
编辑:John D'Errico
2024-5-14
Remember that the software does not truly understand your model, and what you are attempting to simulate in that model. It only understands numbers. And numbers have no non-negativity constraint, so in certain circumstances, you could have a problem where negative numbers made sense.
I would suggest you work in terms of the log of the mass, I'll call it m. Wherever you need the mass itself (M) all you need do is use the simple transformation
M = exp(m)
Note that the log of the mass now has no constraints, and the mass will always be a positive value. This simple transformation makes a lot of sense when you have a problem where a value can span orders of magnitude, as otherwise, there will often be numerical problems. This sort of transformation is often appropriate when you have a variable which really only works in terms of positive numbers.
Jeremy Huard
2024-5-15
I would also recommend to have a look at the Model Debugger for SimBiology. This tool will allow you to define breakpoint conditions (e.g. negative species values) to pause a model simulation and investigate which term is problematic.
社区
更多回答在 SimBiology Community
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Scan Parameter Ranges 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!