-
大小: 2KB文件類型: .py金幣: 1下載: 0 次發(fā)布日期: 2021-01-11
- 語(yǔ)言: Python
- 標(biāo)簽: python??預(yù)測(cè)仿真??SEIR模型??
資源簡(jiǎn)介
通過(guò)已獲得武漢疫情的感染者、死亡人數(shù)等數(shù)據(jù)，按照病毒傳播的SEIR數(shù)學(xué)模型，使用python進(jìn)行仿真，病毒傳染、人員流動(dòng)、醫(yī)療救治等外在因素太多，仿真結(jié)果與實(shí)際疫情肯定有較大出入，仿真結(jié)果僅供參考學(xué)習(xí)。
代碼片段和文件信息
#?-*-?coding:?utf-8?-*-
“““
Created?on?Sat?Feb?15?17:25:25?2020
@author:?Administrator
“““
#?-*-?coding:?utf-8?-*-
“““
Created?on?Sat?Feb?15?12:06:53?2020
新型冠狀病毒感染SEIR模型
@author:?by?Zhang?Haitao
“““
import?numpy?as?np
import?matplotlib.pyplot?as?plt
N?=?1.5e7
T?=?180
susceptible?=?np.zeros（[T]）
exposed?=?np.zeros（[T]）
infected?=?np.zeros（[T]）
recovered?=?np.zeros（[T]）
lamda?=?np.zeros（[T]）????#易感染者變成潛伏者系數(shù)
lamda[0]=1
gamma?=?0.08?#患者治愈的幾率
sigma?=?0.25??#潛伏者變?yōu)楦腥菊叩膸茁?br/>#delta?=?0.1?#治愈人員重新變?yōu)橐愿腥菊叩膸茁?br/>
infected[0]?=?10.0?/?N
susceptible[0]?=?1.5e7?/?N
exposed[0]?=?500.0?/?N
recovered[0]?=?0
for?t?in?range（T?-?1）:
????lamda[t+1]?=?lamda[t]?*?0.98
????susceptible[t?+?1]?=?評(píng)論
共有 條評(píng)論