二分法求函数方程根的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:
在这里插入图片描述
希望对大家有所帮助,喜欢的朋友点个赞,加个关注,谢谢。

Logo

为开发者提供学习成长、分享交流、生态实践、资源工具等服务,帮助开发者快速成长。

更多推荐