《DFP算法+wolfe性非线性搜索解决无约束问题的matlab程序》由会员分享,可在线阅读,更多相关《DFP算法+wolfe性非线性搜索解决无约束问题的matlab程序(2页珍藏版)》请在人人文库网上搜索。
1、function x,val,k=dfp(fun,gfun,x0) %功能: 用DFP算法求解无约束问题:min f(x) %输入: x0是初始点,fun, gfun分别是目标函数及其梯度%输出: x, val分别是近似最优点和最优值, k是迭代次数maxk=1e5; rho=0.25;sigma=0.75; epsilon=1e-3; alpha=1; a=0; b=Inf; k=0; n=length(x0); Hk=eye(n); while(k=sigma*gk*dk)a=alpha;alpha=min(2*alpha, (b+alpha)/2);continue;endbreak;e。
2、ndx=x0+alpha*dk;sk=x-x0; yk=feval(gfun,x)-gk; if(sk*yk0) Hk=Hk-(Hk*yk*yk*Hk)/(yk*Hk*yk)+(sk*sk)/(sk*yk); end k=k+1; x0=x; end val=feval(fun,x0); function f=fun(x) f=100*(x(1)2-x(2)2+(x(1)-1)2; function gf=gfun(x) gf=400*x(1)*(x(1)2-x(2)+2*(x(1)-1), -200*(x(1)2-x(2); clearx0=0,0;x,val,k=dfp(fun,gfun,x0。