找出谁是初始传播者python

Python016

找出谁是初始传播者python,第1张

SIR模型是传染病模型中最经典的一个。SIR是三个单词首字母的缩写,其中S是Susceptible的缩写,表示易感者;I是Infective的缩写,表示感染者;R是Removal的缩写,表示移除者。这个模型本身是在研究这三者的关系。在病毒最开始的时候,所有人都是易感者,也就是所有人都有可能中病毒;当一部分人在接触到病毒以后中病毒了,变成了感染者;感染者会接受各种治疗,最后变成了移除者。

琢磨了下SIR模型,核心是一个微分方程组,怎么解老早忘了

发现各路专家都是用Matlab,Mathematica,python等工具求解

但这种事情怎么少得了Excel大法,请看附件 SIR

绿格子是可调的,其他基本都是公式

基本参数及假设:

12月初发现感染者,8号发病,到月底人均接触数(密切接触者人数)6人,一代病毒感染率4%,康复率2%

1月初开始采取隔离措施,人均接触数快速下降5人,4人,3人,2人,一月底开始二代病毒感染率降到3%,治疗手段也上来了康复率到5%

3月份开始复工率提高,尽管大家隔离意思强,但人均接触数提高到3人,三代病毒感染率2%,加上随着气温升高,感染率4月降到1%

照此推算,峰值出现在2月末,6月中下旬感染人数降到百人

SIR模型是基础模型,提供个框架思路和大致的预测,各类假设加上去就没边了

用excel的好处在于对接触数感染率康复(死亡)率可以作每天调整,以拟合报道数据

还有,门槛低,不需要解微分方程组,中学生就会

再细化一下,可以用规划求解或者水晶球做蒙特卡洛,提高历史数据拟合度,找到更加适合的参数。