MATLAB数值分析：二分法求根实例

Matlab数值分析 二分法求根 代码

%dichotomy二分法求根

%---输入参数------------------------------------------------------------------

a=0;                                %求根区间下限

b=1;                                %求根区间上限

syms x                             %不用改动

f=@(x)exp(x)+10*x-2;     %函数

epsilon=1e-6;                  %误差

%---开始计算------------------------------------------------------------------

k=0;

while (b-a)>epsilon

    k=k+1;

    x=a+(b-a)/2;

if sign(f(x))==sign(f(a))

    a=x;

else

    b=x;

end

end

x=a+(b-a)/2;

fprintf('x为%.5f\n',x)

fprintf('迭代%d次数\n',k)