0.618法matlab

0.618法
建立golds.m文件
function [s,phis,k,G,E]=golds(phi,a,b,delta,epsilon)
%输入：phi是目标函数，a,b是搜索区间的两个端点
% delta,epsilon分别是自变量和函数值的容许误差
%输出：s,phis分别是近似极小点和极小值，G是nx4矩阵
% 其第k行分贝时a,p,q,b的第k次迭代值[ak,pk,qk,bk]
% E=[ds,dphi],分别是s和phis的误差限
t=(sqrt(5)-1)/2;h=b-a;
phia=feval(phi,a);phib=feval(phi,b);
p=a+(1-t)h;q=a+th;
phip=feval(phi,p);phiq=feval(phi,q);
k=1;G(k,:)=[a,p,q,b];
while(abs(phib-phia)>epsilon)|(h>delta)
if(phip<phiq)
b=q;phib=phiq;q=p;phiq=phip;
h=b-a;p=a+(1-t)h;phip=feval(phi,p);
else
a=p;phia=phip;p=q;phip=phiq;
h=b-a;q=a+th;phiq=feval(phi,q);
end
k=k+1;G(k,:)=[a,p,q,b];
end
ds=abs(b-a);dphi=abs(phib-phia);
if(phip<=phiq)
s=p;phis=phip;
else
s=q;phis=phiq;
end
E=[ds,dphi];