如何在Mathematica中计算雅可比行列式

时间:2014-02-20 16:42:10

标签: wolfram-mathematica

说我想以这种方式计算矢量函数的雅可比行列式(在Mathematica 9中):

Clear[Z, Z1, Z2, Z3, p];
phi1[Z2_, Z3_, p_] := p (1 - (1 - Z2)^2) (1 - (1 - Z3)^2);
phi2[Z1_, Z3_, p_] := p (1 - (1 - Z1)^2) (1 - (1 - Z3)^2);
phi3[Z1_, Z2_, p_] := p (1 - (1 - Z1)^2) (1 - (1 - Z2)^2);
JacobianMatrix[{phi1[Z2, Z3, p], phi2[Z1, Z3, p], phi3[Z2, Z3, p]}, {Z1, Z2, Z3}]

为什么不起作用? 我也试过D [{phi1 [Z2,Z3,p],phi2 [Z1,Z3,p],phi3 [Z2,Z3,p]},{Z1,Z2,Z3}],但它没有救命。从这里 https://mathematica.stackexchange.com/questions/5790/how-to-make-jacobian-automatically-in-mathematica 我认为它可能是函数或参数p?

的问题

更新 从其中一条评论中我了解到雅可比矩阵已经过时了,但更重要的是,这个问题似乎是通过在这样的自变量集周围增加两个大括号来解决的:

D[{phi1[Z2, Z3, p], phi2[Z1, Z3, p], phi3[Z2, Z3, p]}, {{Z1, Z2, Z3}}]

结果:

{{0, 2 p (1 - Z2) (1 - (1 - Z3)^2), 
2 p (1 - (1 - Z2)^2) (1 - Z3)}, {2 p (1 - Z1) (1 - (1 - Z3)^2), 0, 
2 p (1 - (1 - Z1)^2) (1 - Z3)}, {0, 2 p (1 - Z2) (1 - (1 - Z3)^2), 
2 p (1 - (1 - Z2)^2) (1 - Z3)}}

0 个答案:

没有答案