我有一个微分方程系统,需要计算雅可比矩阵.下面的代码抛出AttributeError: 'list' object has no attribute 'ravel'.我错过了什么?
import numpy as np
import numdifftools as ndt
def rhs(z, t=0):
x,y = z
xdot = (x/5 + y)*(-x**2+1)
ydot = -x*(-y**2+1)
return [xdot, ydot]
Jfun = ndt.Jacobian(rhs)
Jfun([1,1])
Run Code Online (Sandbox Code Playgroud)
做就是了:
return np.array([xdot, ydot])
Run Code Online (Sandbox Code Playgroud)
代替.这应该工作......
| 归档时间: |
|
| 查看次数: |
6322 次 |
| 最近记录: |