Есть ли способ оценить якобиана в Matlab?

Question

У меня есть этот код:

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]

Есть ли способ сделать это?

Answer 1

Вы можете использовать функцию subs:

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