雅可比在枫树中的应用

Question

如何将下面编写的Jacobian用作(x，y)中的函数？

g := (x, y) -> x - y
u := (x, y) -> x^2 + y^2

J := jacobian([g(x, y), u(x, y)], [x, y]);

我的想法是做这样的趣事

Jf := (u, v) -> subs(x = u, y = v, J(x, y))

但它会返回丑陋的矩阵，里面有虚张声势。

我用枫树17

Answer 1

不建议使用linalg包(它导出jacobian命令)和小写matrix。使用LinearAlgebra和矩阵代替，并使用VectorCalculus:-Jacobian。

另外，注意unapply的使用。

restart:

g := (x, y) -> x - y:

u := (x, y) -> x^2 + y^2:

J:=VectorCalculus:-Jacobian([g(x,y),u(x,y)],[x,y]);

                                [ 1     -1 ]
                           J := [          ]
                                [2 x    2 y]

Jf:=unapply(J,[x,y]):

Jf(1,1);

                               [1    -1]
                               [       ]
                               [2     2]

Jf(s,t);
                             [ 1     -1 ]
                             [          ]
                             [2 s    2 t]