How to find Jacobian of a vector field given as a dataset?

18 次查看（过去 30 天）

Mirlan Karimov 2020-6-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/545999-how-to-find-jacobian-of-a-vector-field-given-as-a-dataset

编辑： J. Alex Lee 2020-6-12

I have a dataset given in 3x 1 x m x n x l. 3x1 vector at each point of m x n x l 3D tensor holds the velocity vector [u,v,w]. How do I estimate the Jacobian matrix at each point?

7 个评论
显示 5更早的评论隐藏 5更早的评论

Mirlan Karimov 2020-6-11

在 MATLAB Online 中打开

Hi John,

Thank you for your reply. The velocity field data may be a result of a simulation or experiment. I do not understand what you mean when you say "talking about the Jacobian of a point set is meaningless". Of course, numerical Jacobian is defined at each point to understand how the velocity field maps from point to point.

What I have done initially is to calculate the Jacobian as follows:

%Calculate Gradients
[gradxVx,gradyVx,gradzVx] = gradient(vx,stepx,stepy,stepz);
[gradxVy,gradyVy,gradzVy] = gradient(vy,stepx,stepy,stepz);
[gradxVz,gradyVz,gradzVz] = gradient(vz,stepx,stepy,stepz);
gradV = zeros(3,3,length(yy),length(xx),length(zz));
gradVT = zeros(3,3,length(yy),length(xx),length(zz));
for i=1:length(yy)
    for j=1:length(xx)
        for k=1:length(zz)
        gradV(1,1,i,j,k) = gradxVx(i,j,k);
        gradV(2,1,i,j,k) = gradxVy(i,j,k);
        gradV(3,1,i,j,k) = gradxVz(i,j,k);
        
        gradV(1,2,i,j,k) = gradyVx(i,j,k);
        gradV(2,2,i,j,k) = gradyVy(i,j,k);
        gradV(3,2,i,j,k) = gradyVz(i,j,k);
        
        gradV(1,3,i,j,k) = gradzVx(i,j,k);
        gradV(2,3,i,j,k) = gradzVy(i,j,k);
        gradV(3,3,i,j,k) = gradzVz(i,j,k);
        end
    end
end

Mirlan Karimov 2020-6-12

Alex, the domain is fixed.

Yes, I found that velocity gradient as a term is quite confusing. As you noted, the gradient well suits the situations where there is deformation and we measure at each point how the point moves.

Anyway, my method works correctly as I benchmarked it with a data set generated from an analytical velocity field. I just wanted to know if there is a better way of doing this.

J. Alex Lee 2020-6-12

编辑：J. Alex Lee 2020-6-12

Great, glad that you have a working method!

Just to note that I don't think "velocity gradient" is confusing at all...I think "deformation gradient" could potentially be confusing.

Back to the actual question at hand:

I think your method is as good as can be unless you want to try to compute with higher accuracy finite differences (TMW's gradient() function uses central differences; docs also state that edge gradients are approximated by one-sided differences, which are only 1st order accurate). Added effort for higher "accuracy" would probably go to waste unless your velocities were computed using similarly accurate FDs (if you method is based on FD).

If you're extracting the nodal solutions from an FEM-based code (I've seen it before on Answers), you might be better off computing the gradients within the FEM framework and exporting it, if you have control over that.

请先登录，再进行评论。

请先登录，再回答此问题。