我想通過使用 odeint 之類的東西來求解具有可變系數的二階微分方程。這個的問題是如果初始條件很復雜（現在就是這種情況），它就不起作用。你知道一種用類似于 odeint 的東西來解決上述方程的方法嗎？
2 回答
茅侃侃
TA貢獻1842條經驗 獲得超22個贊
你總是可以使用真正的組件（odeint約定）
def odesys(u,t):
z = u[0]+1j*u[1]
dz = u[2]+1j*u[3]
d2z = f(t,z,dz)
return [ dz.real, dz.imag, d2z.real, d2z.imag ]
其中f代表二階 ODE 的顯式形式。
如果我沒記錯的話，您可以使用的其中一種方法（“vzode”？）scipy.integrate.ode直接適用于復雜的狀態變量。
添加回答
舉報
0/150
提交
取消