fpush.m

黄金分割法的matlab程序

%Fpush.m文件
function [x1,x2]=Fpush(t0,h)
%进退法，调用函数为,F618Ret,t0为初始值，h为步长
a=t0;
b=t0+h;
while 1
    if F618Ret(a)<F618Ret(b)
        h=0-h;
        a=a+h;
        b=a-h;
        if F618Ret(b)<F618Ret(b-h) &&F618(b)<F618Ret(a) 
            break;
        end
    else
        a=b;
        b=a+h;
        if F618Ret(a)<F618Ret(b) && F618Ret(a)<F618Ret(a-h) 
            break;
        end
    end
end
x1=a;
x2=b;
%============================================
%测试数据
%   >> [x1,x2]=Fpush(0,0.8);
%   >> x1
%
%   x1 =
%
%       4.8000
%
%   >> x2
%
%   x2 =
%
%       5.6000
%
%   >>