Matlab 利用二分法求解发f(x)=x-sinx在[-1,sqrt2]上的解

最新推荐文章于 2023-03-12 21:54:24 发布

原创最新推荐文章于 2023-03-12 21:54:24 发布 · 2.5k 阅读

6 ·

CC 4.0 BY-SA版权

Matlab 专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种使用MATLAB或Octave等环境下的二分法求解方程x-sin(x)=0的方法，该方程定义在区间[-1,sqrt(2)]内，并给出了具体的实现代码。

利用二分法求解发f(x)=x-sinx在[-1,sqrt2]上的解

function erfenfa
clc,clear all
a=-1;b=sqrt(2);
epsilon=1e-05;
if fun(a)*fun(b)<=0
    disp('方程的解是: ')
    x0=(a+b)/2;
    while abs(fun(x0))>=epsilon
        c=(a+b)/2;
        if fun(a)*fun(c)<=0
            b=c;
        else
            a=c;
        end
        x0=(a+b)/2;
    end
else
    disp('方程无解')
end
x0
function y=fun(x)
y=x-sin(x)