近似黄金分割法MATLAB代码

以下是近似黄金分割法的MATLAB代码：

function [x_min, f_min] = golden_section(func, a, b, tol)
% func: 目标函数
% a, b: 初始区间
% tol: 精度要求
phi = (1+sqrt(5))/2; % 黄金分割比例
x1 = b - (b - a) / phi;
x2 = a + (b - a) / phi;
f1 = feval(func, x1);
f2 = feval(func, x2);
while abs(b - a) > tol
    if f1 > f2
        a = x1;
        x1 = x2;
        f1 = f2;
        x2 = a + (b - a) / phi;
        f2 = feval(func, x2);
    else
        b = x2;
        x2 = x1;
        f2 = f1;
        x1 = b - (b - a) / phi;
        f1 = feval(func, x1);
    end
end
x_min = (a + b) / 2;
f_min = feval(func, x_min);
end

其中，`func`是目标函数的句柄，`a`和`b`是初始区间，`tol`是精度要求。函数返回最小值点`x_min`和最小值`f_min`。在函数内部，先计算出黄金分割比例`phi`，然后循环直到区间长度小于精度要求。在每次循环中，根据目标函数在两个新的点上的取值，更新区间的左右端点，并计算新的两个点。最后返回区间的中点作为最小值点，并计算最小值。