编写正割法计算方程根近似值的Matlab程序

以下是使用正割法计算方程根近似值的Matlab程序：

function [root, iter] = secant(f, x0, x1, tol, maxiter)
% f: 函数句柄，表示要计算根的函数
% x0, x1: 初始点
% tol: 容忍误差
% maxiter: 最大迭代次数

iter = 0;
while abs(f(x1)) > tol && iter < maxiter
    iter = iter + 1;
    % 计算斜率
    slope = (f(x1) - f(x0)) / (x1 - x0);
    % 计算根的近似值
    root = x1 - f(x1) / slope;
    % 更新迭代点
    x0 = x1;
    x1 = root;
end

if abs(f(x1)) <= tol
    disp(['找到根的近似值：' num2str(root)]);
else
    disp('超过最大迭代次数，未能找到根的近似值。');
end

end

你可以将自己的函数作为输入参数`f`传递给该函数，例如：

f = @(x) x^3 - 2*x - 5;
[root, iter] = secant(f, 2, 3, 1e-6, 100);

其中，`f`表示要计算根的函数，`2`和`3`是初始点，`1e-6`是容忍误差，`100`是最大迭代次数。函数的输出是根的近似值和实际迭代次数。

相关问题

用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。