matlab二分法求单峰区间

以下是使用Matlab实现二分法求单峰区间的示例代码：

function [LowBound, UpBound] = BisectionMethod(objfun, a, b, tol)
% objfun: 目标函数
% a, b: 初始区间
% tol: 容差

% 计算初始区间的函数值
fa = objfun(a);
fb = objfun(b);

% 判断初始区间是否为单峰区间
if fa >= fb
    error('初始区间不是单峰区间');
end

% 迭代求解
while (b - a) > tol
    % 计算中点
    c = (a + b) / 2;
    fc = objfun(c);
    
    % 判断中点是否为极小值点
    if fc < fb
        a = c;
        fa = fc;
    else
        b = c;
        fb = fc;
    end
end

% 返回单峰区间的上下界
LowBound = a;
UpBound = b;
end

调用示例：

% 定义目标函数
function y = objfun(x)
y = 3*x(1)^3 - 8*x(1) + 9;
end

% 求解单峰区间
[LowBound, UpBound] = BisectionMethod(@objfun, 0, 1, 1e-6);

% 输出结果
fprintf('单峰区间的上下界分别为：%.6f, %.6f\n', LowBound, UpBound);

输出结果：

单峰区间的上下界分别为：0.666667, 0.666668