Adding Plotting to bisection method

Question

gracias claude 2021-2-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/744297-adding-plotting-to-bisection-method

评论： Alan Stevens 2021-2-13

Not sure what the c is in this bisection method. Also I would like to add plotting of the intervals

function [x,e] = MyBisectFunc(f, a1,b1, number)

format long

c1 = f(a1); d1 = f(b1);

if c1*d1 > 0.0

error ('Func has same sign at both endpoints')

end

for i = 1:number

x = (a1 + b1) /2;

y = f(x);

disp([x,y])

if y == 0.0

e = 0

break

end

if c*y < 0

b = x

else

a = x

end

x = (a1 + b1)/2;

e = (b1 - a1) /2;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Alan Stevens 2021-2-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/744297-adding-plotting-to-bisection-method#answer_622527

在 MATLAB Online 中打开

Like so:

f = @(x) x^2 - 2;  % arbitrary function
a1 = 0; b1 = 1.8;  % initial end-points
c1 = f(a1); d1 = f(b1);
if c1*d1 > 0.0
error ('Func has same sign at both endpoints')
end
number = 20;
for i = 1:number
    
    c1 = f(a1); d1 = f(b1);
    if c1*d1 > 0.0
    error ('Func has same sign at both endpoints')
    end
    x = (a1 + b1) /2;
    y = f(x);
    disp([x,y])
    if abs(y) < 10^-8
       break
    end
    if c1*y < 0    %  c1 not c
        b1 = x;
    else
        a1 = x;
    end
    
end

6 个评论
显示 4更早的评论隐藏 4更早的评论

gracias claude 2021-2-13

not very familar with the syntax

Alan Stevens 2021-2-13

在 MATLAB Online 中打开

This section

if abs(y) < 10^-8
       break
    end

breaks out of the loop if the condition is met. Probably best to remove it if you want to plot the intervals. You'll need to do something like:

f = @(x) x^2 - 2;  % arbitrary function
number = 20;  % number of iterations
lo = zeros(1,number);   % initialise vector
hi = zeros(1,number);
lo(1) = 0;
hi(1) = 1.8;
for i = 1:number-1
    
    flo = f(lo(i)); fhi = f(hi(i));
    if flo*fhi > 0.0
    error ('Func has same sign at both endpoints')
    end
    mid = (lo(i) + hi(i)) /2;
    fmid = f(mid);
   
    if flo*fmid < 0    %  c1 not c
        hi(i+1) = mid;
        lo(i+1) = lo(i);
    else
        lo(i+1) = mid;
        hi(i+1) = hi(i);
    end
    
end
disp([mid fmid])
plot(1:number,abs(hi-lo),'o'),grid
xlabel('iteration number'), ylabel('interval')