黄金分割法代码
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()进行输出