species go to negative values

Question

Day 2024-5-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2118776-species-go-to-negative-values

评论： 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 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.

Day 2024-5-30

Thank you

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

John D'Errico 2024-5-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2118776-species-go-to-negative-values#answer_1457541

编辑：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.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Day 2024-5-30

Thank you

请先登录，再进行评论。

Answer 2

Jeremy Huard 2024-5-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2118776-species-go-to-negative-values#answer_1457906

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.