0

我的代码工作正常，除非我在下面的函数方程中输入 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 数组放入函数中有问题如果有人可以为我的问题提供输入，我会

4

1 回答 1

1

F[0] = expression期望expression在这里是一个数字，而不是一个数组。然而Q/V * (CO2_To-y[0]) + TG_CO2 + r1是一个r1维度数组。要看到这一点，请尝试评估以下行：

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

要消除异常，您应该根据您要实现的目标以某种方式将此表达式转换为一个数字。

于 2013-03-22T13:39:35.540 回答