syms x1 x2;f=4*(x1)2-4*(x1)*(x2)+3*(x2)2+x1;%函数表达式%f=8*(x1)2+4*(x1)*(x2)+5*(x2)2;X=-10:0.1:10; Y=-10:0.1:10;X,Y=meshgrid(X,Y);Z=4.*X.2-4.*X.*Y+3.*Y.2+X;mesh(X,Y,Z)contour(X,Y,Z)x0=10 10;%初始点v=x1,x2; mu=0.01;%最小误差gradf=gradient(f);%函数的梯度H=jacobian(gradf,v);h0=subs(H,x1;x2,x0);%在点x0处的Hessiang0=subs(gradf,x1;x2,x0);%在点x0处的梯度值if g0*g0mu%检验是否满足精度条件 minf=subs(f,x1;x2,x0);%函数的最小值 return;endd0=-g0;%搜索方向alpha=-(g0*d0)/(d0*h0*d0);%步长xk
