matlab实现二分法

二分法是一种简单有效的数值计算方法，可以用于求解函数的零点问题。在 MATLAB 中实现二分法，可以按照以下步骤进行：

1. 定义待求解的函数，例如：

function y = myfun(x)
    y = x^3 - 2*x - 5;
end

2. 定义求解区间的下界和上界，例如：

a = 1;
b = 2;

3. 定义迭代终止的条件，例如：

maxiter = 100; % 设置最大迭代次数
tol = 1e-6; % 设置收敛精度

4. 开始进行迭代，直到满足终止条件。具体实现代码如下：

iter = 0;
while (iter < maxiter) && (abs(b-a) > tol)
    c = (a+b)/2; % 计算区间的中点
    if myfun(c) == 0 % 如果中点就是零点，则直接返回
        break;
    elseif myfun(a)*myfun(c) < 0 % 如果中点和下界函数值的符号不同，则将中点作为新的上界
        b = c;
    else % 如果中点和上界函数值的符号不同，则将中点作为新的下界
        a = c;
    end
    iter = iter + 1;
end

5. 最后输出结果，例如：

fprintf('迭代次数：%d，近似解：%f\n', iter, (a+b)/2);

完整代码如下：

function y = myfun(x)
    y = x^3 - 2*x - 5;
end

a = 1;
b = 2;
maxiter = 100;
tol = 1e-6;

iter = 0;
while (iter < maxiter) && (abs(b-a) > tol)
    c = (a+b)/2;
    if myfun(c) == 0
        break;
    elseif myfun(a)*myfun(c) < 0
        b = c;
    else
        a = c;
    end
    iter = iter + 1;
end

fprintf('迭代次数：%d，近似解：%f\n', iter, (a+b)/2);