冠状病毒SEIR模型仿真

通过已获得武汉疫情的感染者、死亡人数等数据，按照病毒传播的SEIR数学模型，使用python进行仿真，病毒传染、人员流动、医疗救治等外在因素太多，仿真结果与实际疫情肯定有较大出入，仿真结果仅供参考学习。

# -*- 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]）    #易感染者变成潜伏者系数
lamda[0]=1
gamma = 0.08 #患者治愈的几率
sigma = 0.25  #潜伏者变为感染者的几率
#delta = 0.1 #治愈人员重新变为以感染者的几率


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] =