擅长:python、mysql、java
<p>要计算向量函数的旋度,还可以使用<a href="https://pypi.python.org/pypi/Numdifftools" rel="nofollow noreferrer">numdifftools</a>进行自动数值微分,而无需绕道符号微分。Numdifftools不提供<code>curl()</code>函数,但它计算一个或多个变量的向量值函数的雅可比矩阵,这提供了向量场的所有分量相对于所有变量的导数;这是计算旋度所必需的。</p>
<pre><code>import import scipy as sp
import numdifftools as nd
def h(x):
return sp.array([3*x[0]**2,4*x[1]*x[2]**3, 2*x[0]])
def curl(f,x):
jac = nd.Jacobian(f)(x)
return sp.array([jac[2,1]-jac[1,2],jac[0,2]-jac[2,0],jac[1,0]-jac[0,1]])
x = sp.array([1,2,3)]
curl(h,x)
</code></pre>
<p>返回旋度在<code>x</code>:<code>array([-216., -2., 0.])</code>的值
绘制如上文所述。</p>