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])

要摆脱异常，你应该根据你想要达到的目标将这个表达式转换为数字。

ValueError：请设置一个带序列的数组元素

1 个答案: