新冠肺炎SEIR模型(修正版)

初始模型 seir模型

clear;clc;
N=1400000000;
I=1;
R=0;
D=0; 
E=0;
S=N-I;
r=5;
r11=0;
a=0.11;
B=0.05;
B1=0.04;
count=0;
y=0.02;
r2=15;
k=0.025373;
T=20:1000;
for idx =1:length(T)-1
     if(idx>=25)
         r=1;
         r2=1;
         B=0.01;
         y=0.6;
         k=0.03;
         B1=0.01;
     end
     if(idx>=100)
         r=0;
         r2=0;
        y=0.9;
       k=0.05;
     end
   S(idx+1)=S(idx)-r*B*I(idx)*S(idx)/N(idx)-r2*B*E(idx)*S(idx)/N(idx);   
   E(idx+1)=E(idx)-a*E(idx)+r*B*I(idx)*S(idx)/N(idx)+r2*B*E(idx)*S(idx)/N(idx); 
   I(idx+1)=I(idx)+a*E(idx)-(y+k)*I(idx); 
   R(idx+1)=R(idx)+y*I(idx);  
   D(idx+1)=D(idx)+k*I(idx);
   N(idx+1)=N(idx)-D(idx);
end
plot(T,S,T,E,T,I,T,R,T,D);
hold on
grid on;
xlabel('日期 ');
ylabel(' 人数');
legend('易感者','潜伏者','传染者','康复者','死亡者');
title('模型')

相比与之前的模型,这个模型考虑到了随着疫情的发展,各项参数都会发生变化.现在只能说是再次对整个疫情的一次数学模拟,并没有考虑地区的自然死亡率,以及出生率.
最终结果是模拟康复人数为8.4万人.
在这里插入图片描述参考文档1

参考文档2
(致敬所有的抗疫英雄)

  • 10
    点赞
  • 99
    收藏
    觉得还不错? 一键收藏
  • 9
    评论
评论 9
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值