domen0.m

一些优化的源代码

function t=domen0(a,b)                   %[a,b]为初始区间  t:极小值点   
e=0.001;
a0=a;
b0=b;
F0=(b0-a0)/e;
F(1)=1;
F(2)=1;
for i=3:inf
    F(i)=F(i-1)+F(i-2);
    if(F(i)>=F0)
        n=i;
        break;
    end
end

rafa1=a0+(1-F(n-1)/F(n))*(b0-a0);
rafa2=a0+F(n-1)/F(n)*(b0-a0);
y1=fai(rafa1);
y2=fai(rafa2);
for i=1:n
if(y1<=y2)
    if((rafa2-a0)<=e)
        t=rafa1;
        break;
    else
        a0=a0;
        b0=rafa2;
        rafa2=rafa1;
        y2=y1;
        rafa1=a0+(1-F(n-i-1)/F(n-i))*(b0-a0);
        y1=fai(rafa1);
    end
else
    if((b0-rafa1)<=e)
        t=rafa2;
        break;
    else
        a0=rafa1;
        b0=b0;
        rafa1=rafa2;
        y1=y2;
        rafa2=a0+F(n-i-1)/F(n-i)*(b0-a0);
        y2=fai(rafa2);
    end
end
end
y=fai(t)