二分法求函数方程根的matlab实现(内附例题和代码)

二分法求函数方程根的matlab实现(内附例题和代码)

摘要:本文主要给出二分法求函数方程根的matlab代码，并给出例题对代码进行验证。并且代码具有很好的交互性，按照提示操作即可。

对于二分法求方程的根的思想我们并不陌生，这里便不在赘述(看下图或百度即可):

重点来啦，下面给出matlab代码:

%%%% 二分法求函数的实根%%%%%
%%%%%    Liu Deping    %%%%%
%%%%%  2020.06.09      %%%%%
clc;
clear all
format short
s=input('请输入函数表达式：f = ','s');
f=inline(s);
a = input('请输入区间左端点:a=');
b = input('请输入区间右端点:b=');
eps=input('请输入停止精度要求:eps=');  %%“|b-x|<=eps”
k=1;
x=(a+b)/2;
fprintf('        k        a        f(a)       b        f(b)       xk        f(xk)\n ');
T=[k,a,f(a),b,f(b),x,f(x)];
while abs(T(k,4)-T(k,6))>eps/2
    k=k+1;
    if  f(x)*f(a)==0
        a=a;
        b=x;
        x=(a+b)/2;
        T=[T;k,a,f(a),b,f(b),x,f(x)];
        break
    elseif  f(x)*f(a)>0
        a=x;
        b=b;
        x=(a+b)/2;
        T=[T;k,a,f(a),b,f(b),x,f(x)];
    elseif  f(x)*f(a)<0
        a=a;
        b=x;
        x=(a+b)/2;
        T=[T;k,a,f(a),b,f(b),x,f(x)];
    end
 
end
disp(T);
fprintf('经过%d次迭代，函数方程根的近似解为：x=%.8f\n',k-1,T(k-1,6))

下面给出例题
例题1:

例题2:

希望对大家有所帮助，喜欢的朋友点个赞，加个关注，谢谢。