割线法解方程matlab代码

割线法是求解非线性方程的一种迭代方法。下面是使用Matlab实现割线法的代码。

function x = secant_method(f, x0, x1, tol, max_iter)
% 使用割线法求解方程 f(x)=0 的根

    % 初始化参数
    iter = 0;  % 迭代次数
    x = x0;  % x 的初始值
    delta = tol + 1;  % 误差初始化为tol+1

    % 迭代求解
    while delta > tol && iter < max_iter
        % 计算割线的斜率
        slope = (f(x1) - f(x0)) / (x1 - x0);
        
        % 计算割线与x轴的交点
        x = x1 - f(x1) / slope;

        % 更新x0和x1的值
        x0 = x1;
        x1 = x;

        % 计算误差
        delta = abs(f(x));

        % 更新迭代次数
        iter = iter + 1;
    end

    % 判断是否找到了根
    if iter == max_iter && delta > tol
        disp('未找到满足要求的根');
    else
        disp(['迭代了 ', num2str(iter), ' 次']);
        disp(['找到的根为 ', num2str(x)]);
    end
end

使用时，首先需要定义一个函数 `f`，这个函数代表了要求解的非线性方程。然后选择一个初始的迭代点 `x0` 和 `x1`，还需要设置一个迭代精度 `tol` 和最大迭代次数 `max_iter`。最后调用 `secant_method` 函数，传入这些参数即可获得方程的近似解 `x`。

例如，我们要求解方程 `x^3 - 2*x - 5 = 0`，可以定义如下的函数：

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

然后使用割线法解方程：

x0 = 1;
x1 = 2;
tol = 1e-6;
max_iter = 100;

x = secant_method(@f, x0, x1, tol, max_iter);

运行上述代码后，将会得到近似解 `x` 的值。