我做了一个求解器,可以在
我遇到的问题如下。如果我决定使用
如果我这样做,那么
编辑:
总结问题:
如果将函数f定义为
编辑:
Scipy--版本0.16.0
scipy.integrate.ode
和scipy.integrate.odeint
之间进行转换。下面是代码。def f(y,s,C,u,v):
y0 = y[0] # u
y1 = y[1] # u'
y2 = y[2] # v
y3 = y[3] # v'
dy = np.zeros_like(y)
dy[0] = y1
dy[2] = y3
C = C.subs({u:y0,v:y2})
dy[1] = -C[0,0][0]*dy[0]**2\
-2*C[0,0][1]*dy[0]*dy[2]\
-C[0,1][1]*dy[2]**2
dy[3] = -C[1,0][0]*dy[0]**2\
-2*C[1,0][1]*dy[0]*dy[2]\
-C[1,1][1]*dy[2]**2
return dy
def solve(C,u0,s0,s1,ds,solver=None):
from sympy.abc import u,v
if solver == None: # use lsoda from scipy.integrate.odeint
s = np.arange(s0,s1+ds,ds)
print 'Running solver ...'
return sc.odeint(f,u0,s,args=(C,u,v))
else: # use any other solver from scipy.integrate.ode
r = sc.ode(f).set_integrator(solver) # vode,zvode,lsoda,dopri5,dop853
r.set_f_params(C,u,v)
r.set_initial_value(u0)
#t = []
y = []
print 'Running solver ...'
while r.successful() and r.t <= s1:
r.integrate(r.t + ds)
y.append(r.y)#; t.append(r.t)
return np.array(y)
我遇到的问题如下。如果我决定使用
scipy.integrate.odeint
中的求解器,则必须按照代码中的顺序指定f
的参数。但是,如果我决定使用scipy.integrate.ode
中的求解器,则必须将函数f(y,s,C,u,v)
的参数顺序更改为f(s,y,C,u,v)
,否则会出现错误。TypeError: 'float' object has no attribute '__getitem__'
如果我这样做,那么
scipy.integrate.odeint
针对定义为f(s,y,C,u,v)
的f
会生成相同的错误。无论参数的顺序如何,我该如何操作以使用统一的f
?编辑:
总结问题:
如果将函数f定义为
f(s,y,C,u,v)
,则scipy.integrate.ode
求解器可以工作,如果将函数f定义为f(y,s,C,u,v)
,则scipy.integrate.odeint
求解器可以工作。为什么会发生这种情况?我该怎么修复?编辑:
Scipy--版本0.16.0
odeint
呢?odeint
使用的求解器lsoda
在ode
中也是可用的。 - Igor Markelov