gold.m

来自「黄金分割法 可以解决 单峰函数 的极值问题 比较的准确和快速」· M 代码 · 共 16 行

%《运筹学》清华大学第三版 149页 例5
%单谷函数f=x^2-x+2;
%ab为区间  e0为精度
%徐射雕
a=-1;b=3;t1=0.382*(b-a)+a;t2=0.618*(b-a)+a;e0=1e-9;
f1=t1^2-t1+2;
f2=t2^2-t2+2;
while(b-a>e0)
    if f1<f2
        b=t2;t2=t1;f2=f1;t1=0.382*(b-a)+a;f1=t1^2-t1+2;
    else
        a=t1;t1=t2;f1=f2;t2=0.618*(b-a)+a;f2=t2^2-t2+2;
    end
end
x=(a+b)/2%最后x的取值为他们中点
fval=min(f1,f2)%极小值