Question

MATLAB编程求梯度值

下面的程序该怎么修改，使得该函数可以接收参数传递，求得在某点的梯度值
function [t]=tidu1(x)
n=length(x)
for i=1:n
s=0;
p=1/2*x(i)^2-log(x(i));
s=s+p;
end
t=jacobian(s,x);
t
end

Answer 1

x=[];
m=[];
for i=1:24
sx=['x',num2str(i)];
sm=['m',num2str(i)];
eval(['syms ',sx,' ',sm])
eval(['x=[x,',sx,'];']);
eval(['m=[m,',sm,'];']);
end
sh=sum(x.^2)/2*1e6;
y=0;
for i=1:length(x)
y=y+sum(x(1:i));
sh=sh+(y-sum(x(1:i))/2)^2*m(i);
end
jac=jacobian(sh,x);
jac'

追问

看不太懂啊，尤其sh的表达式不太懂啊

Answer 2

matlab里有一个强大的求梯度的函数jacobian(f)