seir模型启发式算法_SEIR模式代码

function[y]=seir(t,x)

a=1/14;b=0.03;c=1/14;r=10;

DS=-r*b*x(1)*x(3);

DE=r*b*x(1)*x(3)-a*x(2);

DI=a*x(2)-c*x(3);

DR=c*x(3);

y=[DS;DE;DI;DR];

end

[t,x]=ode45('seir',[0:2:150],[0.95 0.01 0.03 0.03]);

figure(1)

plot(t,x(:,1),'-.r',t,x(:,2),'-.b',t,x(:,3),'-.k',t,x(:,4),'-.g')

title('SEIR模型')

xlabel('时间')

ylabel('人群占比');

legend('易感者','潜伏期','病人','康复者','4')

2.改变参数传染者的传染概率

，将

由0.03降低到0.02时，绘制的SEIR模型的程序代码：

function[y]=seir(t,x)

a=1/14;b=0.02;c=1/14;r=10;

DS=-r*b*x(1)*x(3);

DE=r*b*x(1)*x(3)-a*x(2);

DI=a*x(2)-c*x(3);

DR=c*x(3);

y=[DS;DE;DI;DR];

end

[t,x]=ode45('seir',[0:2:150],[0.95 0.01 0.03 0.03]);

figure(1)

plot(t,x(:,1),'-.r',t,x(:,2),'-.b',t,x(:,3),'-.k',t,x(:,4),'-.g')

title('SEIR传染率减小模型')

xlabel('时间')

ylabel('人群占比');

legend('易感者','潜伏期','病人','康复者','4')