最优化作业||黄金分割法

黄金分割法代码

function [x,f_x]=goldfun(f,a0,b0)
%黄金分割法
a = a0;b = b0; %初始区间
index=0;%迭代数
eps = 0.01; %区间精度
while((b - a) >= eps)
    a;
    b;
    index=index+1;
    x1 = a + 0.382 * (b - a);
    x2 = a + 0.618 * (b - a); %黄金分割法主要步骤
    if f(x1) < f(x2) %两种情形的判断
        b = x2;
    else
        a = x1;
    end
end
x = (a + b) / 2; %得到满足条件的最优解
f_x=f(x);
disp(['黄金分割法最优解： x = ',num2str(x)]);
disp(['此时: f(x) = ',num2str(f(x))]);%使用disp函数和num2str()进行输出