SEIR模型是一种常用的传染病模型,用于模拟传染病在人群之间的传播。SEIR模型包含四个状态:
- S(易感):健康,未感染传染病的人。
- E(潜伏):已感染传染病,但还未发病的人。
- I(感染):已感染传染病,并且已经发病的人。
- R(恢复):已恢复的人。
下面是一个使用Python实现SEIR模型的完整代码示例:
```python import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt
设置参数
beta = 0.3 # 传播系数 gamma = 0.1 # 恢复率 sigma = 1.0 / 3.0 # 潜伏期持续时间 N = 1000 # 总人数
定义SEIR模型
def seir_model(y, t, beta, gamma, sigma): S, E, I, R = y dSdt = -beta * S * I / N dEdt = beta * S * I / N - sigma * E dIdt = sigma * E - gamma * I dRdt = gamma * I return dSdt, dEdt, dIdt, dRdt
初始化参数
S0, E0, I0, R0 = 999, 1, 0, 0 # 假设有1个人已感染,其余999个人健康 y0 &#