SIR传染病预测模型的MATLAB代码

1、一、SIR模型简介

(1)、易染状态S(Susceptible),也称易感者

(2)、感染状态I(Infected),也称感染者

(3)、移除状态R(Removed,Refractory or Recovered)，也称移除者

二、模型中涉及的方程、

dS/dt=-aS(t)I(t)

dI/dt=aS(t)I(t)-bI(t)

dR/dt=bI(t)

三、具体模型matlab实现

建立sir函数

function y=sir(t,x)

a=0.75;b=0.25;

y=[a*x(1)*x(2)-b*x(1),-a*x(1),b*x(1)]';

建立运行main函数

[t,x]=ode45(@sir,[0,40],[0.04,0.95,0.01]);

plot(t,x(:,1),t,x(:,2),t,x(:,3),'lineWidth',2.5)

legend('感染人群 (Infectious)','易感人群 (Susceptible)','康复人群 (Recovered)')

运行主函数main。

四、小结

网上说的不是很清楚，可以看这个我已经加入了legend函数。