资源介绍
x0=[1;1];%初始点
xk=x0;
g0=zeros(2,1);
g1=zeros(2,1);
g2=zeros(2,1);
d0=zeros(2,1);
d1=zeros(2,1);
syms x1
syms x2
xk1=xk;
% 计算xk点的梯度及梯度值
fun=fun(x1,x2);
fx1=diff(fun,'x1');
fx2=diff(fun,'x2');
fun=inline(fun);
fx1=inline(fx1);
fx2=inline(fx2);
funval=feval(fun,xk1(1),xk1(2));
gradx1=feval(fx1,xk1(1),xk1(2));
gradx2=feval(fx2,xk1(1),xk1(2));
% 计算搜索方向d0
d0(1)=-gradx1;
d0(2)=-gradx2;
g0(1)=gradx1;
g0(2)=gradx2;
% 沿搜索方向d0进行一维搜索