Question

我对odeint有疑问。我必须解决一阶差分系统，然后是二阶系统，但我对第一阶系统有点困惑。你能解释我标记为错误的吗？谢谢:)）

import scipy.integrate as integrate
import numpy as np
def fun(t,y):
  ys = np.array([y[1], (1-y[0]**2)*y[1]-y[0]])
  return(ys)
N = 3
x0 = np.array([2.00861986087484313650940188,0])
t0tf = [0, 17.0652165601579625588917206249]
T=([0 for i in range (N+1)])
T[0]= t0tf[0]
Pas = (t0tf[1]-t0tf[0])/N
for i in range (1,N+1):
       T[i]= t0tf[0] + i*Pas
X = integrate.odeint(fun, x0,T,Dfun=None, col_deriv=0,full_output=True)
T = np.array(T)
T = T.reshape(N+1,1)
S = np.append(X,T,axis=1)
print(S)

返回的错误是：

ys = np.array（[y [1]，（1-y [0] ** 2）* y [1] -y [0]]）

TypeError：＆＃39; float＆＃39;对象不可订阅

Answer 1

您需要颠倒衍生函数的参数顺序 - 它应该是f(y, t)，而不是f(t, y)。这与scipy.integrate.ode类使用的顺序相反。

此外，串联S = np.append(X,T,axis=1)将失败，因为X是包含您的积分和字典的元组。请改用S = np.append(X[0],T,axis=1)。

用于差分系统的odeint

1 个答案: