黄金分割法的MATLAB实现

黄金分割法的MATLAB实现

%黄金分割法的matlab实现
function [xmin]=goldSearch(f,a,b,eps)
f=input('请输入所求的函数');
a=input('请输入区间的下界');
b=input('请输入区间的上界');
eps=input('请输入区间的最小阀值长度');
t1=b-0.618*(b-a);%插入点的值
t2=a+0.618*(b-a);
f1=f(t1);
f2=f(t2);
while f1<f2%比较插入点的函数值的大小 
    if abs(t2-a)>eps%循环条件
        b=t2;%进行换名
        t2=t1;
        t1=a+0.382*(b-a);
        f2=f1;
        f1=f(t1);
    else abs(t2-a)<=eps;
        xmin=t1;
        return
    end
        while f1>f2
            if abs(b-t1)>eps
        a=t1;
        t1=t2;
        t2=a+0.618*(b-a);
        f1=f2;
        f2=f(t2);
            else abs(b-t1)<=eps;
            xmin=t2;
            return
            end
        end
end %迭代到满足条件为止就停止迭代