matlab二分法代码

二分法是一种常用的数值计算方法，可以用来求解函数的根。在MATLAB中，实现二分法的代码如下：

function [x,iter] = bisection(f,a,b,tol,maxiter)
% f: 待求解的函数句柄
% a, b: 求解区间
% tol: 迭代精度控制
% maxiter: 最大迭代次数控制

fa = f(a);
fb = f(b);

if fa*fb > 0
    error('The function has the same sign at both ends of the interval.');
end

for iter = 1:maxiter
    x = (a+b)/2;
    fx = f(x);
    if abs(fx) < tol
        break;
    end
    if fx*fa < 0
        b = x;
        fb = fx;
    else
        a = x;
        fa = fx;
    end
end

if iter == maxiter
    warning('The maximum number of iterations has been reached.');
end
end

其中，f为待求解的函数句柄，a、b为求解区间，tol为迭代精度控制，maxiter为最大迭代次数控制。函数返回值x为求解得到的根，iter为迭代次数。