声明：该研究是一个2022年的一个比赛题目，这个比赛当时我做过一次。本篇文章主要是优秀论文中的介绍实用的方法。未经允许，不得转载。——川川菜鸟 题目
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
已知条件
在这里插入图片描述
模型假设
在这里插入图片描述
第一问分析
原文：
在这里插入图片描述
文字： 试分析长春市、上海市及北京市疫情传播的基本规律及其关键性因素，并预测出上海疫情和北京疫情预计社会面清零的时间节点？
数据描述（基本规律）
在这里插入图片描述
对数据分析描述（也就是基本规律）：
在这里插入图片描述
长春：
在这里插入图片描述
上海：
在这里插入图片描述
北京：
在这里插入图片描述
也可以这样的图形：
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
对上述的柱形图分析得到基本规律：
在这里插入图片描述
接着处理描述以外，还可以拟合方程：
在这里插入图片描述
在这里插入图片描述
关键性因素
在这里插入图片描述
清零预测
建立模型：
在这里插入图片描述
带入数据预测：
在这里插入图片描述
在这里插入图片描述
主要参考描述
参考：
https://blog.tigerxly.com/1202.html
在这里插入图片描述
由于接触、隔离和治疗，这四类人群是可以相互转化的。为了简化模型，我们以固定的概率来体现转化。
符号说明：
在这里插入图片描述
代码：
% seir_model_simulate.m
function seir_model_simulate()
clear();
N = 10000;
I = 1; %infectious
S = N - I; %susceptible
R = 0; %recovered
E = 0; %exposed
r = 20; %接触数
beta = 0.03; %被感染者传染
beta_1 = 0.02; %被潜伏着传染
alpha = 0.1; %潜伏期10天
gamma = 0.1; %康复概率
T = 1:150;
for i = 1:length(T) - 1
S(i + 1) = S(i) - r * (beta * I(i) + beta_1 * E(i)) * S(i) / N;
E(i + 1) = E(i) + r * (beta * I(i) + beta_1 * E(i)) * S(i) / N - alpha * E(i);
I(i + 1) = I(i) + alpha * E(i) - gamma * I(i);
R(i + 1) = R(i) + gamma * I(i);
end
plot(T,S,T,E,T,I,T,R);
grid on;
xlabel('天'); ylabel('人数')
legend('易感者','潜伏者','传染者','康复者')
end
如下：
在这里插入图片描述
如果人们减少出行，做好防护，数据修改如下:
r = 3; %接触数
beta = 0.02; %传染概率
beta_1 = 0.01; %潜伏者传染概率
总结
在这里插入图片描述
第二问
原文：
在这里插入图片描述
文字： 在疫情爆发期间医疗资源的合理分配与调度是至关重要的，你能否提出 充分考虑不同区域人口、面积与经济发展水平的医疗物资的合理指派方案，同时在此方案中你需要充分考虑因未能得到及时治疗所导致的非确诊病例的伤亡损失
分析
在这里插入图片描述
模型
符号：
在这里插入图片描述
假设：
在这里插入图片描述
模型建立：
在这里插入图片描述
建议文：
https://doc.mbalib.com/view/685f36269bc4120ad97524c6941d1518.html OLR算法求解：
在这里插入图片描述
第三问
原文：
在这里插入图片描述
文字： 在疫情爆发期间合理的核酸检测方案的设计是至关重要的，某大型小区具有 103 栋楼，封闭人数 18326 人。其中最高 6 层、3 个单元、36 户及平均入住率 70%的楼 60 栋，最高 11 层、2 个单元、44 户及平均入住率 75%的楼 20 栋，最高 18 层、3 个单元、126 户及平均入住率 80%的楼 23 栋。你能否结合上述数据提供考虑不同类型楼栋特征的最优核算监测分组方案、时间间隔及其具体的时间节点？
分析
我们将不同层数，不同入住率的住楼划分为三个不同的住楼区域类型 i：其中最高 6 层、3 单元、36 户及平均入住率 70%的 60 栋楼记为区域 1；最高 11 层、2 个单元、44 户及平均入住率 75%的 20 栋楼记为区域 2；最高 18 层、3 个单元、126 户及平均入住率 80%的 23 栋楼记为区域 3，三个区域共有 18326 人.
主要参考
新冠感染池检测选取最优池检测人数参考：
https://blog.csdn.net/weixin_42845306/article/detAIls/110426116
在这里插入图片描述
代码：
%蒙特卡洛仿真新冠感染池检测选取最优池检测人数
%仿真时长约十几秒
N=100000;%仿真总人数
% N=30000;%仿真总人数
T_single=10;%单独检测时间
T_all=40;%群体检测时间
P_infect=0.005;%感染率
% P_infect=1/30;%感染率
N_infect=N*P_infect;%感染人数，需要是个整数
T0=N*T_single;%单独检测总时间
x=sparse(zeros(N,1));%模拟每个人，感染为1，没感染为0，用稀疏矩阵储存
for i=1:N_infect
while(1)%循环，随机让一个人感染
index=randi(N);%产生一个随机数
if x(index)==1%该人已感染，换一个
continue;
end
break;
end
x(index)=1;%感染
end
all_count=1:1000;%一次混合检测样本数量
T=[];
for j=1:size(all_count,2)%对每一个混合检测样本研究一下
i=1;
T(j)=0;%初始化时间
while(1)
iend=i+all_count(j)-1;%一段区间的最后
if iend>N%超出总人数了
iend=N;
end
infectnum=sum(x(i:iend));%池检测一次感染人数
T(j)=T(j)+T_all;%加上一次群体检测时间
if infectnum>0%有一人感染
T(j)=T(j)+T_single*(iend-i+1);%如果有人感染，需要对组内全部做单独检测
end
if iend==N%检测完了，结束
break;
end
i=i+all_count(j);%区间向后
end
end
F_T_theory=@(n) (T_all.*(1-P_infect).^n+(T_single*n+T_all).*(1-(1-P_infect).^n))./n;%理论计算公式
T_theory=N.*F_T_theory(all_count);%计算理论时间
plot(all_count,T,'linewidth',2);
hold on
plot(all_count,ones(size(all_count,2),1).*T0,'--','linewidth',2);
plot(all_count,T_theory,'--','linewidth',2);
xlabel("单次池检测人数/个")
ylabel("检测总时间/分钟")
legend("池检测","单独检测","池检测理论时间")
grid on
mincount=all_count(T==min(T))
mincount_theory=all_count(T_theory==min(T_theory))
title(sprintf("理论最优单次池检测人数%d，仿真人数%d，感染率%.4f",mincount_theory,N,P_infect))
注意：蒙特卡洛算法具有随机性，每次仿真的结果可能不一样。
模型一（池检模型）
如下：
在这里插入图片描述
由上述运行结果可得，若要对 18326 人均进行核酸检测，则理论最优单次池检人数为 30 人。
模型二（规划）
在这里插入图片描述
结果：
在这里插入图片描述
在这里插入图片描述
第四问
在疫情管控期间居民生活物资的订购与发放也同样非常关键，在基本生 活物资供给有限的情况下你能否提供一个更为科学的生活物资订购与指派方案？在模型中请你需要明确考虑因生活物资交互而产生的疫情传播风险。
模型（规划）
在这里插入图片描述
在这里插入图片描述
求解
在这里插入图片描述
在这里插入图片描述
第五问
在疫情得到有效的控制情况下，逐步的复工复产方案的设计也是非常有 必要的。你能否提供考虑不同行业和每日新增确诊人数的复工复产计划？
模型一： 疫情复工风险模型
在这里插入图片描述
参考：
https://baike.baidu.com/item/%E7%9B%B8%E5%AF%B9%E9%A3%8E%E9%99%A9%E7%8E%87/5556710 模型 二（空间自相关检验模型）
在这里插入图片描述
指定计划如下：
在这里插入图片描述
以上PDF已发数模群701163024，自行下载即可。如果不在群里，可以私聊我。
公众号：玩转大数据 知乎关注一下哦！免费分享数模~
上一篇：如何看待核酸检测企业发布半年报，4 家净利润增长 50% 以上 ...
下一篇：全球79国取消所有新冠入境限制 |
窥视卡
雷达卡
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
千斤顶
照妖镜
