areafind.asv

来自「可以用黄金分割和斐波那契算法实现优化」· ASV 代码 · 共 31 行

function [a,b]=areafind(t0,h,c)
%利用进退法求出，调用函数为g,t0为初始值，h为步长,c为步长的系数
a=t0;
b=t0+h;

while 1
    
    if g(a)<g(b)
        h=h*c;
        a=a-h;
        b=a+h;
        
        if g(b)<g(b+h/c) && g(b)<g(a) 
            b=b+h/c;
            break;
        end
        
    else
       
        a=b;
        h=h*c;
        b=a+h;
        
        if g(a)<g(b) && g(a)<g(a-h/c) 
           a=a-h/c;
           break;
        end

    end
end