二分法求方程的根(Matlab)

二分法在很多地方应该都会见到,这里是通过迭代逼近的方法求出一个方程的根。

function xc = bisection(f,a,b,tol)
% use the bisection method to find the root of the function
% Page 30,computer problem 7(Bisection method)
% input:
% f:the function that transform from the equation
% a,b:the left and right value of the interval which the root is in
% tol:the accuracy
% output:
% xc:the solution of the equation
if sign(f(a)) * sign(f(b)) >=0
    error('f(a)f(b)<0 not satisfied!')
end
if nargin < 3
    disp('The function should at least include 3 parameters');
end
if nargin == 3
    tol = 10^-6;
end
while (b-a)/2 > tol
    c = (a + b)/2;     
    if f(c) == 0         % when f(c) == 0,c is a root of the function
        break 
    end
    if f(a) * f(c) < 0    % a and c form a new interval
        b = c;
    else                  % c and b form a new interval
        a = c;
    end
end
xc = (a+b)/2;             % the mid_rang is the root that we find
©️2020 CSDN 皮肤主题: 大白 设计师: CSDN官方博客 返回首页
实付0元
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值