SEIR模式代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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')

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：SEI function Legend Figure xlabel

这是matlab的代码吧？想请问博主这个模型需要导入什么要的数据，就网络矩阵 x 吗？如何生成这种节点间的网络矩阵呢？另外 t 是时期数吗？谢谢！