%SEIR模型修正
%假设1月15日开始出现第一例确诊,;1月23号武汉政府开始封城,此时其他省市也相应作出响应隔离措施,大约距离15号后的11天政府管控发挥明显作用
% N=13.95*100000000;%人口数
N=1395380000;
load quanguo.mat
%第一列为累计确诊人数,第二列为累计死亡人数,第三列为累计治愈人数
E=0;%潜伏者
D=0;%死亡患者人数
I=1;%感染人数
S=N-I;%易感人数
R=0;%康复者人数
r=17;%感染者接触数量
% r=19;
B=0.602;%传染概率
% a=0.17;%潜伏者转化为感染者的概率
% a=0.175;
a=0.198;%潜伏者转化为感染者的概率
% r2=8;%潜伏者接触人数
r2=15;%潜伏者接触人数
% B2=0.03;%潜伏者传染正常人的概率
B2=0.05;
y=0.05;%康复概率
k=0.0001;%日致死率
B3=0.1;%转阴率
% T=1:200;
T=1:180;
for idx=1:length(T)-1
%若以1月18日为疫情起点,在政府发出管控号召时间以及各个地方响应延迟时间,此处采用11天后为临界点,
%相当于11天后,感染者与潜伏者流动性和医疗配置发生明显变化,具体为接触人数
if idx>=11
r=0.20;%感染者接触数量
r2=1.8;%感染者接触人数量
y=0.15;%康复率上升为0.16
a=0.12;%潜伏者转化为感染者的概率
% a=0.10;
k=0.0001;%日致死率暂无改变
end
if idx<11
S(idx+1)=S(idx)-r*B*S(idx)*I(idx)/N-r2*B2*S(idx)*E(idx)/N;%易感人群迭代
E(idx+1)=E(idx)+r*B*S(idx)*I(idx)/N-a*E(idx)+r2*B2*S(idx)*E(idx)/N;%潜伏者迭代
I(idx+1)=I(idx)+a*E(idx)-(k+y)*I(idx);%感染人数迭代
R(idx+1)=R(idx)+0.05*I(idx);%康复人数迭代
D(idx+1)=R(idx)+k*I(idx);%死亡患者人数迭代
else
S(idx+1)=S(idx)-r*B*S(idx)*I(idx)/N-r2*B2*S(idx)*E(idx)/N+B3*E(idx-10);%易感人群迭代
E(idx+1)=E(idx)+r*B*S(idx)*I(idx)/N-a*E(idx)+r2*B2*S(idx)*E(idx)/N-B3*E(idx-10);%潜伏者迭代
I(idx+1)=I(idx)+a*E(idx)-(k+y)*I(idx);%感染人数迭代
% Y参数有问题
R(idx+1)=R(idx)+0.045*I(idx-9);%康复人数迭代
D(idx+1)=R(idx)+k*I(idx);%死亡患者人数迭代
end
end
B={'01-19','02-08','02-28','03-19','04-08','04-28','05-18','06-07','06-27','07-17','08-06'};
% plot(1:1:102,huibei(:,1)-huibei(:,2)-huibei(:,3),'r*');hold on
plot(1:1:102,quanguo(:,1)-quanguo(:,2)-quanguo(:,3),'g-');hold on
plot(1:1:102,quanguo(:,3),'k-');hold on
% legend('实际患病','实际康复')
% xlabel('天数');
% ylabel('人数');
% legend('实际患病')
plot(T,R,'b',T,I,'r');
grid on;
hold on;
plot([7 7],[0 1000]);
set(gca,'XTickLabel',B)
xlabel('日期');
ylabel('人数');
legend('实际患病','实际康复','预测康复者','预测患者');
title('采取隔离措施的SEIR模型');
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
改进SEIR模型的matlab代码.zip (106个子文件)
SEIR3.m 2KB
quanguo.mat 868B
3.10全国.pdf 577KB
1.30全国.pdf 394KB
1.31全国.pdf 393KB
1.22全国.pdf 392KB
2.1全国.pdf 386KB
1.29全国.pdf 382KB
1.28全国.pdf 377KB
2.12全国.pdf 320KB
2.19全国.pdf 309KB
3.28全国.pdf 307KB
3.25全国.pdf 306KB
3.26全国.pdf 306KB
3.24全国.pdf 306KB
4.16全国.pdf 306KB
4.17全国.pdf 306KB
4.11全国.pdf 306KB
3.31全国.pdf 306KB
4.30全国.pdf 306KB
4.14全国.pdf 306KB
4.10全国.pdf 306KB
4.12全国.pdf 306KB
3.27全国.pdf 306KB
4.15全国.pdf 306KB
3.19全国.pdf 306KB
4.18全国.pdf 305KB
4.19全国.pdf 305KB
4.28全国.pdf 305KB
4.9全国.pdf 305KB
4.13全国.pdf 305KB
3.23全国.pdf 305KB
4.23全国.pdf 305KB
4.2全国.pdf 305KB
4.7全国.pdf 305KB
4.8全国.pdf 305KB
4.27全国.pdf 305KB
4.22全国.pdf 305KB
4.4全国.pdf 305KB
4.20全国.pdf 305KB
4.24全国.pdf 305KB
2.21全国.pdf 305KB
4.5全国.pdf 305KB
3.22全国.pdf 305KB
3.21全国.pdf 305KB
4.21全国.pdf 305KB
4.25全国.pdf 305KB
3.20全国.pdf 305KB
4.3全国.pdf 304KB
3.16全国.pdf 304KB
4.29全国.pdf 304KB
4.26全国.pdf 304KB
4.1全国.pdf 304KB
4.6全国.pdf 304KB
2.13全国.pdf 304KB
3.18全国.pdf 304KB
3.15全国.pdf 304KB
3.14全国.pdf 304KB
3.5全国.pdf 304KB
3.11全国.pdf 304KB
3.12全国.pdf 304KB
3.6全国.pdf 304KB
3.17全国.pdf 303KB
3.29全国.pdf 303KB
3.13全国.pdf 303KB
3.8全国.pdf 303KB
3.7全国.pdf 303KB
3.9全国.pdf 303KB
3.4全国.pdf 303KB
3.30全国.pdf 303KB
2.9全国.pdf 303KB
2.11全国.pdf 303KB
2.26全国.pdf 302KB
2.10全国.pdf 302KB
2.22全国.pdf 302KB
2.23全国.pdf 302KB
2.20全国.pdf 302KB
2.17全国.pdf 302KB
2.18全国.pdf 302KB
3.3全国.pdf 302KB
2.16全国.pdf 302KB
2.28全国.pdf 302KB
2.29全国.pdf 302KB
2.27全国.pdf 302KB
3.2全国.pdf 301KB
3.1全国.pdf 301KB
2.14全国.pdf 301KB
2.25全国.pdf 301KB
2.24全国.pdf 301KB
2.8全国.pdf 301KB
2.15全国.pdf 301KB
2.5全国.pdf 301KB
2.7全国.pdf 300KB
2.4全国.pdf 300KB
2.3全国.pdf 300KB
2.6全国.pdf 299KB
1.24全国.pdf 298KB
1.23全国.pdf 298KB
1.26全国.pdf 297KB
1.21全国.pdf 297KB
共 106 条
- 1
- 2
资源评论
- 夏日のblog2020-11-14代码很好,所以博客在哪,博主博客里没找到
- 妙蛙种子1592022-04-29为什么打不开代码
- 慎独8062021-04-01打不开呀,请问怎么打开,谢谢
- qgguuvio2020-07-03为什么打不开呀
- miao18050562020-06-03可以运行效果很好 请问一下代码中的传染率 接触率 是如何计算的
我是苏~格~拉
- 粉丝: 188
- 资源: 4
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功