Simulation makes no sense

1 次查看（过去 30 天）
Juan Carlos Cao Yao
回答： Juan Carlos Cao Yao 2019-6-28
I have some more troubles with Chance1960_Glycolysis_Respiration model. I have set the glucose's trigger concentration to 0.3. When I run it for 1000 s, I get the step increase of glucose and mildly decreases in time, while oxygen, lactate or pyruvate stay at near 0 levels:
1.PNG
When I run it for 3600 s, I get an error, while the graph remains the same:
Error.PNG
Conversely, when I run it for 1 hour which is the equivalent, results are different.
When I run it for 2 days, the step increase in glucose that was seen in the previously when run at shorter periods of time, is no longer seen, and glucose level is flat from the very beginning at near 0:
Glucose is gone.PNG
I know these profiles are not logical and conflicting, which points to me that the simulation is obviously wrong. How could I fix that? Thanks very much!
Juan Carlos
  1 个评论
Juan Carlos Cao Yao
Juan Carlos Cao Yao 2019-6-28
By the way, my model comes from BioModels database's curated files.

请先登录，再进行评论。

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

回答（3 个）

Jeremy Huard
Jeremy Huard 2019-6-28
Hi Juan Carlos,
as Pax suggested in his answer to your previous question, this model does not contain units, so you can't make use of the unit conversion feature in SimBiology. Which also means that changing the unit of the StopTime won't have any effect.
This StopTime will always be in the unit chosen by the model author to define all quantities, which is here 'second'.
If you want to simulate for 1 hour, you will have to set StopTime to 3600 sec as you previously did.
Now, the error you get suggests that the integration tolerances are not met. You can try to decrease them and see it runs.
You can modify the integration tolerances in the simulation settings:
Capture.PNG
Capture2.PNG
Jeremy Huard
Jeremy Huard 2019-6-28
By the way, I think had changed the glucose injection to 0.3 instead of 0.003 mainly to make easier to spot the event.
But if you revert it to 0.003, the model runs and simulate 3600 sec with the default tolerances without issues.
Juan Carlos Cao Yao
Wow, Jeremy, thank you for this amazing skills showcase. It understand my problem perfectly now and solved it. Integration tolerance is something deep about mathematics, but I get that depending on tolerance, there is a limit to the amount of substance the software can handle, say glucose. Thank you very much.

社区

更多回答在  SimBiology Community

类别

在 Help Center 和 File Exchange 中查找有关 Extend Modeling Environment 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by