因此,我试图解决一个由三个ode组成的系统,并使用ODEint开发了以下代码来解决它们。但是当我运行时,ODEint在调用我的函数来处理方程组时遇到了一个问题。在
from scipy.integrate import odeint
#initial parameters
X0 = 75.
M0 = 150.
G0 = 105.
N0 = 80.
T0 = 7.
u10 = 0.0301231859
u20 = 0.0078947020
u30 = 0.0010708464
Rxg = 1.92
Rxm = 3.84
Rxma = 5.76
#define system of equations
def sugar(t,y,u10,u20,u30,Rxg,Rxm,Rxma):
ddt = [-u10*(X0+(Rxg*(G0-y[0]))+(Rxm*(M0-y[1]))+(Rxma*(N0-y[2]))),
-u20*(X0+(Rxg*(G0-y[0]))+(Rxm*(M0-y[1]))+(Rxma*(N0-y[2]))),
-u30*(X0+(Rxg*(G0-y[0]))+(Rxm*(M0-y[1]))+(Rxma*(N0-y[2])))]
return(ddt)
#time parameter
tinit = 0.
tend = 10.
h = 2.
t = arange(tinit,tend+h,h)
#initial guess
y0 = [G0,M0,N0]
#run odeint
y = odeint(sugar,y0,t,args = (u10,u20,u30,Rxg,Rxm,Rxma))
print 't =',t
print 'y =',y
在我的输出中我得到:
^{2}$改变y0和t的顺序仍然会导致相同的错误。我不理解这个错误,也不能指导我解决这个问题。任何帮助都将不胜感激。在
根据the docs,
sugar
的签名应为而不是
^{pr2}$收益率
相关问题 更多 >
编程相关推荐