非线性共轭梯度法_共轭梯度法_使用非线性共轭梯度法求解优化问题_共轭梯度

function s=golds(xk,dk) %phi为目标函数 % 进退法搜索单峰区间 [a,b] alpha0=0; %初始点 h0=0.1; %h为搜索单峰的区间步长 k=1; %k为迭代次数 ph0=fun(xk+alpha0*dk);%初始值 while 1 alpha1=alpha0+h0; ph1=fun(xk+alpha1*dk); if ph1<ph0 %若ph1<ph0,则扩大步长继续搜索 h1=2*h0; alpha=alpha0; alpha0=alpha1; ph0=ph1; k=k+1; else if k==1 h1=h0; %否则反向搜索 alpha=alpha1; alpha1=alpha0; ph1=ph0; k=2; else break; end end end a=min(alpha,alpha1); b=max(alpha,alpha1); %黄金分割法求最优步长 s epsilon=1e-2; %epsilon为自变量的误差 t=(sqrt(5)-1)/2; %黄金分割法 t=0.618 p=a+(1-t)*(b-a); %p,q为左右试探点 q=a+t*(b-a); i=1; phia(i)=fun(xk+p*dk); phib(i)=fun(xk+q*dk); while 1 if phia<=phib %计算左试探点 if abs(a-q)<epsilon %停算条件，左右端点误差小于epsilon break else b=q; phib=phia; q=p; p=a+(1-t)*(b-a); phia=fun(xk+p*dk); end else %计算右试探点 if abs(b-p)<epsilon %停算条件，左右端点误差小于epsilon break else a=p; phia=phib; p=q; q=a+t*(b-a); phib=fun(xk+q*dk); end end i=i+1; end if (phia<=phib) s = p; else s = q; end