Python“ValueError:使用序列设置数组元素”。我看不出哪个部分不起作用

Question

下面代码的目标是将函数“ddW”插入到odeint中，以查找给定X的W（稍后）。我正在使用print函数来确保函数运行）

import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import odeint
import pint
u = pint.UnitRegistry()

R = 8.31446261815324  * u.J/(u.mol*u.K)
Rgas = R.to            (u.atm*u.L/(u.mol*u.K))

P0  = 10    * u.atm
T   = 400   * u.K
v0  = 2     * u.L/u.min
α   = 0.02  * u.kg**-1

k   = 1.4   * u.L**2/(u.mol*u.kg*u.min)

FA0 = 0.5*P0*v0/(Rgas*T) #Assuming I.G.

def ddW(param,w):                               #param: [X,P]
    X       = param[0]
    P       = u.Quantity(param[1]).magnitude    * u.atm
    W       = u.Quantity(w).magnitude           * u.kg

    d       = np.zeros(2)
    d[0]    = k*FA0/(v0**2*P0)  * P*W*(1-X)     #dX/dW
    d[1]    = -α/2 *P0**2       * (1-X)/P       #dP/dW
    return d

param0 = [0,P0]
Wrange = np.linspace(0,100) *u.kg

#PBR     = odeint(ddW,param0,Wrange)


#plt.plot(Wrange,PBR)

print(ddW(param0,0))

Traceback (most recent call last):

  File "<ipython-input-1-83022bc3b5da>", line 1, in <module>
    runfile(REDACTED, wdir=REDACTED)

  File "C:\Users\Spencer\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 705, in runfile
    execfile(filename, namespace)

  File "C:\Users\Spencer\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 102, in execfile
    exec(compile(f.read(), filename, 'exec'), namespace)

  File "REDACTED", line 43, in <module>
    print(ddW(param0,0))

  File "REDACTED", line 33, in ddW
    d[1]    = -α/2 *P0**2       * (1-X)/P       #dP/dW

ValueError: setting an array element with a sequence.

在我在第29行定义了W之前，我之前对现在的第32行（定义d[0]）有相同的错误消息。困扰我的是，d[1]（据我所知）的输入中没有一个具有大小>；1，所以它应该合适

Answer 1

使用以下命令创建数组。您正在创建的数组是float类型，并且您在d[0]和d[1]中获得的值不是float

d = np.array(np.zeros(2), dtype=np.object)

Answer 2

假设您对第32行（定义d[0]）的修复是将w替换为W，其中取了w的大小，那么您需要对尝试放入d[1]的数量执行相同的操作

第32行的=右侧是一个无量纲的Quantity对象，而第33行的=右侧是一个具有standard_atmosphere / kilogram维度的Quantity对象，它不能像现在一样放入数组d，因为它只需要是数量级，就可以转换为浮点类型

如果要将数组d保持为类型float，请尝试以下操作：

val     = -a/2 *P0**2       * (1-X)/P       #dP/dW
d[1]    = val.magnitude