怎么用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);endend这是fun1的函数,m已知（已赋值）,x是长度为24的未知行矩阵,怎么求函
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/27 02:30:26
怎么用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);endend这是fun1的函数,m已知（已赋值）,x是长度为24的未知行矩阵,怎么求函
怎么用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，拼凑也行，但是我就是一个的偏导数也算不出来，不知道哪里写错了
怎么用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);endend这是fun1的函数,m已知（已赋值）,x是长度为24的未知行矩阵,怎么求函
试试,下面的代码,如果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'
举个例子
syms x y z
f=x^2+x*y+z;
gradient=jacobian(f,[x,y,z])%求梯度
%gradient =
%[ 2*x + y, x, 1]
x=-1;y=2;z=3;
tiduzhi=eval(gradient) %求在（-1，2，3）的梯度值
%
%tiduzhi =
%<...
全部展开
举个例子
syms x y z
f=x^2+x*y+z;
gradient=jacobian(f,[x,y,z])%求梯度
%gradient =
%[ 2*x + y, x, 1]
x=-1;y=2;z=3;
tiduzhi=eval(gradient) %求在（-1，2，3）的梯度值
%
%tiduzhi =
%
% 0 -1 1
对于补充的问题，那就没什么函数，你直接用diff求微分算了
gradient=[diff(f,x),diff(f,y),diff(f,z)]
收起