有没有一种方法可以在Matlab中评估Jacobian矩阵？

Question

I have this code:

syms x y z
f =[x^2+y^2+z^2-100 ,x*y*z-1 ,x-y-sin(z)];
v = [x, y, z];
fp = jacobian(f,v)

这将会得到：

fp =
[ 2*x, 2*y,     2*z]
[ y*z, x*z,     x*y]
[   1,  -1, -cos(z)]

现在我想在 x=y=z=1 的情况下评估此函数。

fp([1,1,1])

给了我

[ 2*x, 2*x, 2*x]

有没有一种方法可以做到这一点？

- ario

Answer 1

你可以使用subs函数：

subs(fp, [x y z], [1 1 1])

- dumbmatter