Python欧拉法模拟SIR类模型时遇TypeError错误求助

首先咱们先拆解这个错误：TypeError: unsupported operand type(s) for -: 本质就是你做减法运算时，两边变量类型不兼容——比如一边是列表、一边是数值，或者某个变量压根没被正确初始化为数值类型。我结合SIR模型欧拉法的实现逻辑，给你梳理下排查和解决思路：

1. 先确认核心微分方程的函数实现 #

SIR模型的核心微分方程是：

• dS/dt = -βSI/N
• dI/dt = βSI/N - γ*I
• dR/dt = γ*I

这里要注意：所有参数（β、γ、N）都得是int/float数值类型，不能是列表或字符串。给你一个正确的方程函数示例：

def sir_model(S, I, R, beta, gamma, N):
    dS_dt = -beta * S * I / N
    dI_dt = beta * S * I / N - gamma * I
    dR_dt = gamma * I
    return dS_dt, dI_dt, dR_dt

2. 初始化存储列表的正确姿势 #

你提到初始化了空列表存各时间步的值，这里要注意：初始的S0、I0、R0得是数值，先把初始值加入列表，再启动循环。比如：

# 初始化模型参数
total_pop = 1000  # 总人口数
initial_S = 990   # 初始易感者数量
initial_I = 10    # 初始感染者数量
initial_R = 0     # 初始康复者数量
beta = 0.3        # 感染率
gamma = 0.1       # 康复率
max_time = 100    # 模拟总时长
time_step = 1     # 时间步长

# 初始化存储各时间步结果的列表
S_list = [initial_S]
I_list = [initial_I]
R_list = [initial_R]
time_list = [0]

3. 循环计算的正确写法 #

循环里必须取列表的最后一个元素来计算下一个值——这是最容易踩的坑！要是直接用列表本身去做运算（比如S_list = S_list - ...），肯定会触发类型错误。正确写法如下：

for _ in range(max_time):
    # 获取当前时间步的数值（取列表最后一个元素）
    current_S = S_list[-1]
    current_I = I_list[-1]
    current_R = R_list[-1]
    
    # 计算当前时间步的微分值
    dS, dI, dR = sir_model(current_S, current_I, current_R, beta, gamma, total_pop)
    
    # 欧拉法更新下一个时间步的数值
    next_S = current_S + dS * time_step
    next_I = current_I + dI * time_step
    next_R = current_R + dR * time_step
    
    # 将新值追加到对应列表中
    S_list.append(next_S)
    I_list.append(next_I)
    R_list.append(next_R)
    time_list.append(time_list[-1] + time_step)

4. 常见错误场景排查 #

• 直接用列表变量（比如S_list）和数值做运算，比如S_list = S_list - beta*S_list*I_list/total_pop，此时S_list是列表，和数值做减法必然报错；
• 初始化时把初始值设成了列表（比如initial_S = [990]），后续所有基于它的运算都会是列表与数值的混合操作，触发类型错误；
• 参数beta或gamma不小心被赋值成了字符串（比如beta = "0.3"），乘法运算时也会触发类似的类型错误。

你可以对照自己的代码检查这些点，把列表操作和数值运算分开，确保每一步参与运算的都是数值类型，应该就能解决这个TypeError了。

内容的提问来源于stack exchange，提问作者Howsikan