我想找到这个功能的衍生物:
function [f] = f4(x)
f=0.5*(x(1)^2+10*x(2)^2);
end
我正在使用符号来找到jacobian:
x = sym('x', [2 1]);
f=0.5*(x(1)^2+10*x(2)^2);
grad=jacobian(f,x)
返回
grad = [ x1, 10*x2]
然后我手动更改它看起来像这样我可以使用feval
:
grad = [ x(1,1), 10*x(2,1)];
我发现使用feval
比使用subs
更快。我想找到我可以使用feval
评估的任何函数的派生,并避免手动执行任何操作。