以下代码的目标是将功能“ ddW”插入odeint中,以找到给定X的W(以后)。我正在使用打印功能来确保功能运行)。
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,因此应该合适。
假设您对第32行的定义(定义d [0])将替换w
为W
您获得的大小w
,则您需要对要放入d [1]的数量执行相同的操作。
第=
32行的右侧是无量纲的Quantity
对象,而第=
33行的右侧是其Quantity
尺寸standard_atmosphere / kilogram
不能按d
原样放入数组的对象,因为它仅需为大小,因此可以转换为float类型。
如果要保留数组d
类型,请尝试以下操作float
:
val = -a/2 *P0**2 * (1-X)/P #dP/dW
d[1] = val.magnitude
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句