【源码】二分法及MATLAB实现

二分法是一种求解给定函数根的数值方法。二分法的实质是通过不断地把函数f(x)零点所在的区间一分为二,使区间的两个端点逐步逼近零点,进而得到零点近似值。
在这里插入图片描述

在这里插入图片描述

将m = (a+b)/2作为新的b端点,继续重复以上步骤,直到逼近误差足够小为止。

二分法的求解步骤:
在这里插入图片描述

举例:求解f(x) = x2– 3在区间[0, 4]上的根。
在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

完整MATLAB求解函数:

function p = bisection(f,a,b)

% f为待求解函数,a与b分别为待求解区间的两个端点值

% provide the equation you want to solvewith R.H.S = 0 form.

% Write the L.H.S by using inline function

% Give initial guesses.

% Solves it by method of bisection.

% A very simple code. But may come handy

if f(a)*f(b)>0

disp(‘Wrong choice bro’)

else

p= (a + b)/2;

err = abs(f§);

while err > 1e-7

if f(a)*f§<0

  b = p;

else

  a = p;          

end

p= (a + b)/2;

err = abs(f§);

end

end

英文课件下载地址:

http://page2.dfpan.com/fs/7l6cdjf262211269163/

更多精彩文章请关注微信号:在这里插入图片描述

展开阅读全文

没有更多推荐了,返回首页