Is this a right way of computing Jacobian matrix of a velocity field?

Hi there,

I am computing the Jacobian of velocity field V = [vx; vy; vz].

I take a numerical gradient of a vector field

[gradxVx,gradyVx,gradzVx] = gradient(vx,stepx,stepy,stepz);
[gradxVy,gradyVy,gradzVy] = gradient(vy,stepx,stepy,stepz);
[gradxVz,gradyVz,gradzVz] = gradient(vz,stepx,stepy,stepz);

Step sizes are spacings in corresponding directions because: without declaring the spacing, The spacing between points in each direction is assumed to be 1.

In the following step, I just assemble the partial derivatives in each point to get the Jacobian at that point.

Is this is a right way of computing Jacobian matrix?

[Licheng Feng]

Hi, I have the same question, how to compute the jacobian of vector field, given a series discrete points. Do you solve this problem?