资源描述:
《用共轭梯度法求函数最优值的代码.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、%%用共轭梯度法求函数f=1.5*x1^2+0.5*x2^2-x1*x2-2*x1的极小点symsx1x2x01x02adisp('请输入初始点值x0,它是一个二维向量');x0=input('x0>>');disp('请输入收敛精度');e=input('e>>');ee=1;k=0;f=1.5*x1^2+0.5*x2^2-x1*x2-2*x1;fx1=diff(f,x1);fx2=diff(f,x2);f2x1=diff(diff(f,x1),x1);f2x2=diff(diff(f,x2),x2);fx1x2=diff(diff(f,
2、x1),x2);fx2x1=diff(diff(f,x2),x1);g=[fx1;fx2];H=[f2x1fx1x2;fx2x1f2x2];g0=subs(g,{x1,x2},{x0(1),x0(2)});d0=-g0;sign=0;while(double(ee)>e)xx=x0+a*d0;k=k+1y=subs(f,{x1,x2},{xx(1),xx(2)});aa=solve(diff(y,a),'a');xx=subs(xx,a,aa);g1=subs(g,{x1,x2},{xx(1),xx(2)});b=(norm(double(
3、g1),2))^2/(norm(double(g0),2))^2d1=-g1+b*d0;ee=norm(double(g1),2);x0=xx;d0=d1;G=subs(H,{x1,x2},{x0(1),x0(2)});ifdouble(H(1))>0&double(det(H))>0sign=1;endendifsign==1disp('极值点如下:');disp(x0(1));disp(x0(2));disp('极值如下:');disp(double(subs(f,{x1,x2},{x0(1),x0(2)})));elsedisp('所
4、得海赛矩阵不是正定的,请重新计算');end