matlab二分法程序,二分法及其matlab程序-经典.ppt

二分法及其matlab程序-经典

* 三. 二分法及其MATLAB程序 1. 用二分法求方程 f(x)=0的根 x*的近似值 xk 的步骤 步骤1. 若对于a0, disp('注意：ya*yb>0,请重新调整区间端点a和b.'), return end max1=-1+ceil((log(b-a)- log(abtol))/ log(2)); for k=1: max1+1 a;ya=fun(a); b;yb=fun(b); x=(a+b)/2; yx=fun(x); wuca=abs(b-a)/2; k=k-1; [k,a,b,x,wuca,ya,yb,yx] if yx==0 a=x; b=x; elseif yb*yx>0 b=x;yb=yx; else a=x; ya=yx; end if b-a< abtol , return, end end k=max1; x; wuca; yx=fun(x); % ceil是上取整 例: 确定方程x3-x+4=0的实根的分布情况，并用二分法求在开区间 (-2,-1)内的实根的近似值，要求精度为0.001. 0.003 7 0.012 1 -0.004 8 0.001 0 -1.795 9 -1.794 9 -1.796 9 9 0.012 1 0.029 0 -0.004 8 0.002 0 -1.794 9 -1.793 0 -1.796 9 8 0.029 0 0.062 7 -0.004 8 0.003 9 -1.793 0 -1.789 1 -1.796 9 7 0.062 7 0.129 6 -0.004 8 0.007 8 -1.789 1 -1.781 3 -1.796 9 6 -0.004 8 0.129 6 -0.141 8 0.015 6 -1.796 9 -1.781 3 -1.812 5 5 0.129 6 0.390 6 -0.141 8 0.031 3 -1.781 3 -1.750 0 -1.812 5 4 -0.141 8 0.390 6 -0.716 8 0.062 5 -1.812 5 -1.750 0 -1.875 0 3 -0.716 8 0.390 6 -2.000 0 0.125 0 -1.875 0 -1.750 0 -2.000 0 2 0.390 6 2.125 0 -2.000 0 0.250 0 -1.750 0 -1.500 0 -2.000 0 1 2.125 0 4.000 0 -2.000 0 0.500 0 -1.500 0 -1.000 0 -2.000 0 0 函数值f(xk) 函数值f(bk) 函数值f(ak) 中点xk 右端点bk 左端点ak 次数k 3. 用