优化黄金分割程序matlab

我可以为您提供一个基于黄金分割法的最优化算法的matlab代码示例：

function [xmin,fmin] = golden_section_search(f,xL,xU,tol)
phi = (1+sqrt(5))/2; %黄金分割比例
x1 = xU - (xU - xL)/phi; %计算内部点
x2 = xL + (xU - xL)/phi;
f1 = f(x1); %计算函数值
f2 = f(x2);
while abs(xU - xL) > tol %迭代直到满足精度要求
if f1 < f2 %缩小右侧区间
xU = x2;
x2 = x1;
f2 = f1;
x1 = xU - (xU - xL)/phi;
f1 = f(x1);
else %缩小左边区间
xL = x1;
x1 = x2;
f1 = f2;
x2 = xL + (xU - xL)/phi;
f2 = f(x2);
end
end
xmin = (xL + xU)/2; %最优解
fmin = f(xmin); %最优解的函数值

end

请注意，这个算法可以用来优化任意一个单变量函数。