How to use ODE 45 to generate a SIR model

Question

Zhukun Wang 2021-2-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/753964-how-to-use-ode-45-to-generate-a-sir-model

评论： darova 2021-2-24

function Math462hw2Q6parte

alpha=1.99;

beta=1;

gamma=1/7;

S=0.99;

I=0.01;

R=0.1;

%tstart=0;

%tend=100;

y=[S,I,R];

tspan=0:1:100;

[t,soln]=ode45(@(t,y)dotfunctions,tspan,y);

%plot function

plot(t,soln);

xlabel('time');

ylabel('value');

title('double pendulum');

function ddt=dotfunctions(~,sir)

S=sir(1);

I=sir(2);

R=sir(3);

Sdot=-beta*S*I+alpha*R;

Idot=beta*S*I-gamma*I;

Rdot=gamma*I-alpha*R;

ddt=[S,I,R,Sdot,Idot,Rdot];

end

The above is my code and MATLAB keeps saying the line with the function ode45 is not right. I do not know where the mistakes are. Could someone help me figure it out? THANKS!

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

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2021-2-23

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/753964-how-to-use-ode-45-to-generate-a-sir-model#answer_631279

在 MATLAB Online 中打开

Like so:

alpha=1.99;
beta=1;
gamma=1/7;
S=0.99;
I=0.01;
R=0.1;
%tstart=0;
%tend=100;
y=[S,I,R];
tspan=0:1:100;
[t,soln]=ode45(@(t,y)dotfunctions(t,y,alpha,beta,gamma),tspan,y);  % Pass constants to function
%plot function
plot(t,soln);
xlabel('time');
ylabel('value');
title('double pendulum');
    function ddt=dotfunctions(~,sir,alpha,beta,gamma)
        S=sir(1);
        I=sir(2);
        R=sir(3);
        
        Sdot=-beta*S*I+alpha*R;
        Idot=beta*S*I-gamma*I;
        Rdot=gamma*I-alpha*R;
        ddt=[Sdot;Idot;Rdot];  % Just the gradients.  Must be column vector
    end