怎么用matlab求m文件中的函数的梯度

问题描述：

怎么用matlab求m文件中的函数的梯度
function sh=fun1(x,m)
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
end
这是fun1的函数,m已知（已赋值）,x是长度为24的未知行矩阵,怎么求函数的梯度,即对x中所有未知量的偏导数,结果应该是公式组成的行矩阵而不是数
matlab应该是可以求的吧，但是这个涉及到矩阵，用diff或者jacobian怎么写
写出一个diff，拼凑也行，但是我就是一个的偏导数也算不出来，不知道哪里写错了

1个回答分类：综合 2014-10-03

问题解答：

我来补答

试试,下面的代码,如果xi是实数将eval(['syms ',sx,' ',sm])改为eval(['syms ',sx,' ',sm,' real']),计算结果可能更简洁一些,这里由于不知道你的m值,所以也将其参数化了,你可以将其换为你需要的值.若有问题,可继续追问
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'

展开全文阅读

上一页：解不等式

下一页：我的第一次文学尝试。