SIS传染病数学模型代码

SIS模型是一种经典的传染病数学模型，其代码实现如下（使用Python语言）：

import numpy as np
import matplotlib.pyplot as plt

# 定义参数
beta = 0.2  # 传染率
gamma = 0.1  # 恢复率
N = 1000  # 总人口数
I0 = 1  # 初始感染人数
T = 200  # 模拟时间

# 初始化数组
S = np.zeros(T+1)
I = np.zeros(T+1)

# 初始值
S[0] = N - I0
I[0] = I0

# 模拟SIS模型
for t in range(T):
    S[t+1] = S[t] - beta * S[t] * I[t] / N + gamma * I[t]
    I[t+1] = I[t] + beta * S[t] * I[t] / N - gamma * I[t]

# 绘图
plt.plot(range(T+1), S, label='Susceptible')
plt.plot(range(T+1), I, label='Infected')
plt.xlabel('Time')
plt.ylabel('Population')
plt.title('SIS Model')
plt.legend()
plt.show()

运行以上代码，可以得到SIS模型的人口感染情况随时间的变化图。