传染病模型python_传染病动力学模型 SI => SIS => SIR => SEIR（python）

Parameter

# parameter

N = 10000 # Headcount

I = 1 # The Infected

S = N - I # Susceptible

E = 0 # Exposed

R = 0 # Recover

r = 10 # Number of people per infected person per day.

r2 = 10

B = 0.01 # Transmission probability

B2 = 0.01

Time = 350 # Day

v = 0.02 # Cure Rate

a = 0.1 # Conversion Rate

SI

# SI Model

# dI / dt = r * B * I * S / N

# dS / dt = - r * B * I * S / N

# I(n) = I(n - 1) + r * B * I(n - 1) * S(n - 1) / N

# S(n) = S(n - 1) - r * B * I(n - 1) * S(n - 1) / N

TimeOfVirus = [1]

s_SI = [S]

i_SI = [I]

for day in range(1, Time):

TimeOfVirus.append(day + 1)

s_SI.append(s_SI[day - 1] - r * B * i_SI[day - 1] * s_SI[day - 1] / N)

i_SI.ap