f618.m

来自「黄金分割法的matlab程序」· M 代码 · 共 47 行

%F618.m文件

function [t,min]=F618(precision,t0,h)
%初始区间为a,b,相对精度precision 即为原区间的多少分之,t0为初始值，h为步长
%用进退法得到高低高区间，在用0.618法求最值
%precision=0.001;
[a,b]=Fpush(t0,h);
An=a;
Bn=b;
n=1;
while 0.618^(n-1)>=precision
    A=0.618.*(An-Bn)+Bn;
    B=0.618.*(Bn-An)+An;
    if A<B
         if F618Ret(A)<F618Ret(B)
             Bn=B;
         else
             An=A;
         end
     else
         if F618Ret(A)<F618Ret(B)
             An=A;
         else
             Bn=B;
         end
       
    end
    n=n+1;
end
t=An;
min=F618Ret(An);
%=======================================
%测试过程
%   >> [t,min]=F618(0.00001,0,0.8);
%   >> t
%
%   t =
%
%       4.9999
%
%   >> min
%
%   min =
%
%      11.0000
%========================================