goldsect.m

来自「优化设计中的一维搜索方法」· M 代码 · 共 25 行

function [fx,fy,iter]=GoldSect(func,limt)
if length(limt)~=2
    error('error input limt!');
end
IniL = limt(1);
IniH = limt(2);
lamda = IniL + 0.382*(IniH-IniL);
mu = IniL + 0.618*(IniH-IniL);
iter = 0;
while IniH-IniL>10^4*eps
    iter = iter+1;
    if subs(func,lamda)>subs(func,mu)
        IniL = lamda;
        IniH = IniH;
        lamda = mu;
        mu = IniL + 0.618*(IniH-IniL);
    else
        IniL = IniL;
        IniH = mu;
        mu = lamda;
        lamda = IniL + 0.382*(IniH-IniL);
    end
end
fx = IniL;
fy = subs(func,IniL);