Question

如何使用下面写的雅可比作为（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))

但它会返回带有内置刹车的丑陋矩阵。

P上。 S.我使用Maple 17

Answer 1

不推荐使用linalg包（导出jacobian命令）和小写matrix。请改用LinearAlgebra和Matrix，以及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]

Jacobian在Maple中使用

1 个答案: