How to solve a system of three ODE

Question

James Beard 2012-6-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41195-how-to-solve-a-system-of-three-ode

I was looking into epidemiology, and I came to a system of three ODEs to solve the spread. For the life of me I can't remember how to do this in Matlab. Here is my system, with the initial conditions.

S(0) = 19,465,197
E(0) = 0
I(0) = 0
dS/dt=-147199.4991*S+10*S*I
dE/dt=10*S*I-147203.5121*E
dI/dt=4*E-147199.5197*I

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

请先登录，再进行评论。

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

Answer 1

Benjamin Schwabe 2012-6-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41195-how-to-solve-a-system-of-three-ode#answer_50793

Hi James,

start writing it in a vectorform with x=(S,E,I); Then you have x(0)=x0 as defined and you can write dx/dt=f(x); create a function doin this. the you might use ode45 oder another ode solver for you problem. The help is quite detailed here.

Hope that helps, Benjamin

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

请先登录，再进行评论。

Answer 2

Walter Roberson 2012-6-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/41195-how-to-solve-a-system-of-three-ode#answer_50794

在 MATLAB Online 中打开

WIth the symbolic toolbox, it would look somewhat like

syms S(t) E(t) I(t)
dsolve(S(0) == 19465197, E(0) == 0, I(0) == 0, diff(S(t), t) == -1.471994991*10^5*S(t)+10*S(t)*I(t), diff(E(t), t) == 10*S(t)*I(t)-1.472035121*10^5*E(t), diff(I(t), t) = 4*E(t)-1.471995197*10^5*I(t))

with solution

E(t) = 0, I(t) = 0, S(t) = 19465197*exp(-147199.4991*t)