golden.m

线性规划 包括单纯形法

%-------------用黄金分割法求最佳lamuda--------------%
function [f_min,xinla] = Golden(region,x,s)

    sigma = 1/10^10;
    GoldRate = (sqrt(5)-1)/2;
    rate = GoldRate*( region(2)-region(1));
    d1 = region(2)-rate;
    f1 = funla(d1,x,s);
    d2 = region(1)+rate;
    f2 = funla(d2,x,s);
    
    while((region(2)-region(1))> sigma )
    if (f1 <= f2)
        region(2) = d2;
        d2 = d1;
        f2 = f1;
        d1 = region(1)+region(2)-d2;
        f1 = funla(d1,x,s);
    else
        region(1) = d1;
        d1 = d2;
        f1 = f2;
        d2 = region(1)+region(2)-d1;
        f2 = funla(d2,x,s);
    end
    end
    xinla = (region(1)+region(2))/2;
    f_min = funla(xinla,x,s);