species go to negative values

2 次查看（过去 30 天）
Day
Day 2024-5-14
评论： Day 2024-5-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
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
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
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 中查找有关 Extend Modeling Environment 的更多信息

标签

产品


版本

R2024a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by