Skip to content

Commit 37d9bd0

Browse files
committed
updates
1 parent 8e0c84c commit 37d9bd0

File tree

1 file changed

+6
-10
lines changed

1 file changed

+6
-10
lines changed

lectures/sir_model.md

Lines changed: 6 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -31,9 +31,7 @@ kernelspec:
3131
* [NBER工作论文第26867号](https://www.nber.org/papers/w26867)
3232
* [COVID-19工作论文和代码](https://sites.google.com/site/andyatkeson/home?authuser=0)
3333

34-
他的这些笔记主要是介绍了定量建模
35-
36-
传染病动态研究。
34+
他的这些笔记主要是介绍了定量建模传染病动态研究。
3735

3836
疾病传播使用标准SIR(易感者-感染者-移出者)模型进行建模。
3937

@@ -72,11 +70,11 @@ from scipy.integrate import odeint
7270
- 易感者(S):尚未感染,可能被感染的人群
7371
- 潜伏者(E):已感染但尚未具有传染性的人群
7472
- 感染者(I):已感染且具有传染性的人群
75-
- 移出者(R):已经康复或死亡的人群
73+
- 移出者($R$):已经康复或死亡的人群
7674

7775
需要注意的是:
7876
- 一旦康复,就会获得免疫力,不会再次感染
79-
- 处于移出状态(R)的人包括康复者和死亡者
77+
- 处于移出状态($R$)的人包括康复者和死亡者
8078
- 潜伏期的人虽然已感染,但还不能传染给他人
8179

8280
### 时间路径
@@ -171,11 +169,11 @@ def F(x, t, R0=1.6):
171169
"""
172170
s, e, i = x
173171
174-
# 易感人群的新暴露
172+
# 计算新增感染人数
175173
β = R0(t) * γ if callable(R0) else R0 * γ
176174
ne = β * s * i
177175
178-
# 时间导数
176+
# 导数
179177
ds = - ne
180178
de = ne - σ * e
181179
di = σ * e - γ * i
@@ -205,9 +203,7 @@ x_0 = s_0, e_0, i_0
205203
```{code-cell} ipython3
206204
def solve_path(R0, t_vec, x_init=x_0):
207205
"""
208-
通过数值积分求解i(t)和c(t),
209-
给定R0的时间路径。
210-
206+
给定R0的时间路径,计算感染人数i(t)和累计病例c(t)的演变轨迹。
211207
"""
212208
G = lambda x, t: F(x, t, R0)
213209
s_path, e_path, i_path = odeint(G, x_init, t_vec).transpose()

0 commit comments

Comments
 (0)