0.618法求解凸函数极小值(内附matlab代码)

0.618法求解凸函数极小值(内附matlab代码)

%%%%      0.618法      %%%%%
%%%%%    Liu Deping    %%%%%
%%%%%    2020.09.21    %%%%%
clc;
clear all
format short
s=input('请输入函数表达式：f = ','s');
f=inline(s);
a = input('请输入区间左端点a的值:');
b = input('请输入区间左端点b的值:');
eps= input('请输入运算精度:');    %%“|b-a|<=eps”
k=0;
x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
fprintf(' k          [a,b]          x1        x2        f(x1)       f(x2)\n ');
fprintf('%d      [%.3f,%.3f]    %.3f     %.3f     %.3f      %.3f\n', k,a,b,x1,x2,f(x1),f(x2));
while abs(a-b)>eps
    k=k+1;
    x1=a+0.382*(b-a);
    x2=a+0.618*(b-a);
    if  f(x1)>f(x2)
        a=x1;
        b=b;
    elseif  f(x1)<=f(x2)
        a=a;
        b=x2;
    end
   fprintf(' %d      [%.3f,%.3f]    %.3f     %.3f     %.3f      %.3f\n', k,a,b,x1,x2,f(x1),f(x2));
end
fprintf('经过%d次迭代，可取%.3f作为近似最优解',k,(a+b)/2);