问题描述
我的ode求解器出现了一些问题，我试图解决SEIR问题，尽管我的代码基于非常相似的代码，但我仍然遇到相同的错误.我的代码是:
I'm having some issues with my ode solver, I am trying to solve an SEIR problem and I keep getting the same errors dispite the code that i have based my code on being very similar. My code is:
import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt
# Total population, N.
N1 = 55600
# Initial number of infected and recovered individuals, I0 and R0.
I10, R10, E10 = 1, 0, 0
# Everyone else, S0, is susceptible to infection initially.
S10 = N1 - I10 - R10 - E10
# parameters
B = 0.05
a = 0.000001
d = 0.0167
g = 0.0167
z = 0.0167
M = 100000
# A grid of time points (in months)
t = np.linspace(0, 160, 160)
# The SIR model differential equations.
def deriv(y, t, N1, B, a, d, g, z, M):
S1, E1, I1, R1 = y
dS1dt = B*N1 + d*(R1) - S1/N1 * (M*a(I1))
dE1dt = S1/N1 * M*a(I1) - g * E1
dI1dt = g * E1 - z * I1
dR1dt = z * I1 - d * R1
return dS1dt, dE1dt, dI1dt, dR1dt
# Initial conditions vector
y0 = S10, E10, I10, R10
# Integrate the SIR equations over the time grid, t.
ret = odeint(deriv, y0, t, args=[N1, B, a, d, g, z, M])
S1, E1, I1, R1 = ret.T
我不断收到错误消息:
文件"C:/Users/Angus/PycharmProjects/firstAttempt/bugfinder.py"，第44行 在
File "C:/Users/Angus/PycharmProjects/firstAttempt/bugfinder.py", line 44, in
ret = odeint(deriv, y0, t, args=[N1, B, a, d, g, z, M],)
在odeint中，文件"C:\ Python36 \ lib \ site-packages \ scipy \ integrate \ odepack.py"，第215行ixpr，mxstep，mxhnil，mxordn，mxords)odepack.error:额外的参数必须在元组中
File "C:\Python36\lib\site-packages\scipy\integrate\odepack.py", line 215, in odeintixpr, mxstep, mxhnil, mxordn, mxords)odepack.error: Extra arguments must be in a tuple
任何帮助将不胜感激！
推荐答案
尝试替换:
ret = odeint(deriv, y0, t, args=[N1, B, a, d, g, z, M],)
与此:
ret = odeint(deriv, y0, t, args=(N1, B, a, d, g, z, M))
来自scipy 文档:
额外的参数传递给函数.
Extra arguments to pass to function.
此外，Google 差异黑白列表和元组.
Also, google differences b/w list and tuple.
这篇关于odepack.error:额外的参数必须在元组中的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！