Python 步长(单位:numdifftools-jacobian)
我想得到非线性方程组的数值雅可比矩阵。我在Python3.4上使用模块numdifftools 0.9.16。 我构建了一个函数“fun”,它接受一个数组X并返回一个数组Y。 我通过命令获得雅可比矩阵:Python 步长(单位:numdifftools-jacobian),python,numpy,scipy,derivative,numerical-analysis,Python,Numpy,Scipy,Derivative,Numerical Analysis,我想得到非线性方程组的数值雅可比矩阵。我在Python3.4上使用模块numdifftools 0.9.16。 我构建了一个函数“fun”,它接受一个数组X并返回一个数组Y。 我通过命令获得雅可比矩阵:jac=nd.jacobian(fun)(X),效果很好 我还使用命令控制步长:jac=nd.Jacobian(fun,step=1e-1)(X),这也很好 现在,我希望根据方程有两个不同的步长,例如:jac=nd.Jacobian(fun,step=[1,1,1e-1,1e-1])(X),假设有
jac=nd.jacobian(fun)(X)
,效果很好
我还使用命令控制步长:jac=nd.Jacobian(fun,step=1e-1)(X)
,这也很好
现在,我希望根据方程有两个不同的步长,例如:jac=nd.Jacobian(fun,step=[1,1,1e-1,1e-1])(X)
,假设有4个方程。我想用1表示前两个方程,用1e-1表示其他方程
不幸的是,这个命令不起作用。我找不到任何强制不同步长的示例,但阅读手册这应该是可行的。没有人有这方面的经验?谢谢。我将分别对函数的每个分量应用梯度法(使用自己的大小),然后堆叠结果。我将分别对函数的每个分量应用梯度法(使用自己的大小),然后堆叠结果。