0.618法(黄金分割法)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%本函数用于最优化方法中一维搜索的0.618法（黄金分割法）
function xmin=goldensection(a1,b1,e)
%xmin=a1;%此处xmin为所求极小点;%这里k是迭代次数
%global rackx racky h z1 m m11 m12 KK
a(1)=a1;b(1)=b1;
k=1;
while b(k)-a(k)>e
x=[a(k),b(k)];
y(k)=a(k)+0.382*(b(k)-a(k));
u(k)=a(k)+0.618*(b(k)-a(k));
if F(y(k))>F(u(k))
a(k+1)=y(k);
b(k+1)=b(k);
y(k+1)=u(k);
u(k+1)=a(k+1)+0.618*(b(k+1)-a(k+1));
else a(k+1)=a(k);
b(k+1)=u(k);
u(k+1)=y(k);
y(k+1)=a(k+1)+0.382*(b(k+1)-a(k+1));
end
k=k+1;
end
xmin=(a(k)+b(k))./2;
k=k-1 ; %k-1为迭代次数
function f=F(x)
% global rackx racky h z1 m m11 m12 KK
% for j=1:h
% for i=1:KK
% slope1(j,i)=tan((rackx(j+1,(i-1)*4+2)-rackx(j+1,(i-1)*4+1))/(racky(j+1,(i-1)*4+2)-racky(j+1,(i-1)*4+1)));
% slope2(j,i)=tan((rackx(j+1,(i-1)*4+4)-rackx(j+1,(i-1)*4+3))/(racky(j+1,(i-1)*4+4)-racky(j+1,(i-1)*4+3)));
% end
% end
f=x^2-4*x+4;