How can I define a time-dependent reaction rate for an SIR model in Simbiology

Question

Nicholas DePietro 2015-2-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180764-how-can-i-define-a-time-dependent-reaction-rate-for-an-sir-model-in-simbiology

评论： Arthur Goldsipe 2016-9-21

I am trying to create an SIR model in Simbiology. I input my SIR as species and then define parameters. But the problem is one of my reaction rates is periodic and depends on time. I keep getting the error "Error evaluating reaction rate (from reaction Infectm) at t=0.0000." I'm also a little unsure of the kinetic laws I should be applying for the reactions. I believe I->R should be mass action, but am unsure of the other ones. If anyone has any insight or suggestions that they think would be helpful I would greatly appreciate it.

Thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Ingrid Tigges 2015-3-2

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180764-how-can-i-define-a-time-dependent-reaction-rate-for-an-sir-model-in-simbiology#answer_169849

ExponentialGrowth_parameter.zip

Hello Nicolas,

You can use rate rules to define parameters that are dependent on time. This link of the documentation explains you how to use rate rules: Create a Rate Rule for a Constant Rate of Change While this documentation page talks about species being rate dependent you can also use the same approach for parameters. It is important that the tick for constant is not set in these cases. Attached you can find a very basic example with a parameter being dependent on time which is based on the above documentation page. Does this help?

What are your problems with understanding the kinetic laws to use?

2 个评论
显示无隐藏无

Alberto Merchante 2016-9-20

Dear Ingrid, I have a similar problem I cannot fix: I want to make a function dependent on the values of another function at previous time, for example f2(t) = f1(t-5)

That is a very simple example which could be solved otherwise, but the idea is, to have as variables at a simulation time t2, the values of other variables at earlier time t1 Thanks so much!

Best wishes, Alberto

Arthur Goldsipe 2016-9-21