值错误：无法用序列设置数组元素

Question

我的代码基本上都能正常工作，但当我把r1放进下面这个函数的公式里时，就出错了。

def u(Substrate):         
    return((u_max*ys[:,0])/(Ks+ys[:,0]))


biomass = ys[:,1]    
u = u(ys[:,0])

def r1(u,biomass):
    r1 = u*biomass*YieldCO2_1
    return r1

r1 = r1(u,biomass)

def F(y,t):    
    Ptot = 710
    Vgas = 2
    D = 0.00826*(273.15+Temp)    
    Cstar_CO2 = KH_CO2 * y[2]
    Cstar_CH4 = KH_CH4 * y[3]
    TG_CO2 = KLa_CO2*(Cstar_CO2-y[0])
    TG_CH4 = KLa_CH4*(Cstar_CH4-y[1])
    Q_CO2 = -D*V*TG_CO2
    Q_CH4 = -D*V*TG_CH4
    Qgas = (Q_CO2+Q_CH4)+Q
    F=np.zeros(4)
    F[0] = Q/V * (CO2_To-y[0]) + TG_CO2 + r1
    F[1] = Q/V * (CH4_Do-y[1]) + TG_CH4 
    F[2] = -Ptot*D*TG_CO2*(V/Vgas)-y[2]*(Qgas/Vgas)
    F[3] = -Ptot*D*TG_CH4*(V/Vgas)-y[3]*(Qgas/Vgas)
    return F 

yinit = np.array([4,3,250,200])             
ts = np.arange(0,4,0.4)  
y = odeint(F,yinit,ts)

当在公式F[0]中看到r1时，我收到了以下错误信息：

F[0] = Q/V * (CO2_To-y[0]) + TG_CO2 + r1
ValueError: setting an array element with a sequence.
odepack.error: Error occurred while calling the Python function named F

不过，当我在没有r1数组的情况下运行这个函数时，就没有错误。所以我在想，可能是把r1数组放进函数里出了问题。如果有人能帮我看看这个问题，我会很感激。

Answer 1

F[0] = expression 这里期望 expression 是一个数字，而不是一个数组。不过 Q/V * (CO2_To-y[0]) + TG_CO2 + r1 是一个和 r1 一样大小的数组。要理解这一点，可以尝试计算下面这一行：

>>> 1 + numpy.array([1,2])
array([2, 3])

为了消除这个异常，你需要根据你想要实现的目标，把这个表达式转换成一个数字。