用matlab实现二分法、牛顿法、 牛顿下山法、正割法的代码

以下是 Matlab 实现二分法、牛顿法、牛顿下山法和正割法的代码：

二分法：

function [x, iter] = bisection(f, a, b, tol, maxiter)
% 二分法求解方程 f(x) = 0
% f：函数句柄，a：区间左端点，b：区间右端点，tol：容差，maxiter：最大迭代次数
% x：方程的解，iter：实际迭代次数

iter = 0;
while iter < maxiter
    x = (a + b) / 2;
    if abs(f(x)) < tol
        break;
    elseif f(x) * f(a) < 0
        b = x;
    else
        a = x;
    end
    iter = iter + 1;
end

end

牛顿法：

function [x, iter] = newton(f, df, x0, tol, maxiter)
% 牛顿法求解方程 f(x) = 0
% f：函数句柄，df：导函数句柄，x0：迭代初值，tol：容差，maxiter：最大迭代次数
% x：方程的解，iter：实际迭代次数

iter = 0;
while iter < maxiter
    x = x0 - f(x0) / df(x0);
    if abs(f(x)) < tol
        break;
    end
    x0 = x;
    iter = iter + 1;
end

end

牛顿下山法：

function [x, iter] = newton_downhill(f, df, ddf, x0, tol, maxiter)
% 牛顿下山法求解方程 f(x) = 0
% f：函数句柄，df：导函数句柄，ddf：二阶导函数句柄，x0：迭代初值，tol：容差，maxiter：最大迭代次数
% x：方程的解，iter：实际迭代次数

iter = 0;
while iter < maxiter
    x = x0 - df(x0) / ddf(x0);
    if abs(f(x)) < tol
        break;
    end
    x0 = x;
    iter = iter + 1;
end

end

正割法：

function [x, iter] = secant(f, x0, x1, tol, maxiter)
% 正割法求解方程 f(x) = 0
% f：函数句柄，x0：迭代初值，x1：迭代初值，tol：容差，maxiter：最大迭代次数
% x：方程的解，iter：实际迭代次数

iter = 0;
while iter < maxiter
    x = x1 - f(x1) * (x1 - x0) / (f(x1) - f(x0));
    if abs(f(x)) < tol
        break;
    end
    x0 = x1;
    x1 = x;
    iter = iter + 1;
end

end

其中，f、df 和 ddf 分别代表函数、导函数和二阶导函数（牛顿法和牛顿下山法需要导函数和二阶导函数），a、b、x0 和 x1 分别代表区间左端点、区间右端点、迭代初值和迭代初值的两个近似值，tol 和 maxiter 分别代表容差和最大迭代次数。函数返回方程的解 x 和实际迭代次数 iter。