对分法求方程的根MATLAB完整代码

数值分析中，用对分法求方程根的完整MATLAB代码。首先是对分法的完整代码，，文件中举了一个函数的例子函数可以更改。是数值分析课程很好的实验案例

function bisec_g（f_nameacxminxmaxn_points）
%	ac:end points of  initial  interval
%	tolerance	: tolerance
%	it_limit:limit of iteration number
%	Y_a Y_c : y values of the current end points 
clfhold off
clear Y_aclear Y_c
wid_x=xmax-xmin;
dx=（xmax-xmin）/n_points;
xp=xmin:dx:xmax;
yp=feval（f_namexp）;
plot（xpyp‘r‘）;                        % y= f（x）
xlabel（‘x‘）;ylabel（‘f（x）‘）; title（‘Bisection Method‘）hold on
ymin=min（yp）;ymax=max（yp）;wid_y=ymax-ymin;
yp=0.0*xp;
plot（xpyp）                             % x axes 
fprintf（‘Bisection Scheme \n\n‘）;
fprintf（‘It	 a	 b	 c		fa= f（a）‘）;
fprintf（‘	fc= f（c） abs（fc-fa）	abs（c-a）/2\n‘）;
tolerance=0.000001;
it_limit = 30;
it=0;
Y_a=feval（f_namea）;
Y_c=feval（f_namec）;
plot（[aa][Y_a0]‘black‘）;text（a-0.1*wid_y‘x=a‘）   %line
plot（[cc][Y_c0]‘black‘）;text（c-0.3*wid_y‘x=c‘）   %line
if（Y_a*Y_c>0）
    fprintf（‘f（a）f（c）>0\n‘）;
else
    while 1
        it=it+1;
        b=（a+c）/2;
        Y_b=feval（f_nameb）;
        plot（[bb][Y_b0]‘:‘）;   %line
        plot（b0‘o‘）               %o point
 if it<4
     text（bwid_y/20[num2str（it）]）
 end
    fprintf（‘%3.0f %10.6f %10.6‘itab）;
    fprintf（‘%10.6f %10.6f %10.6‘cY_aY_c）;
    fprintf（‘%12.3e %12.3e\n‘abs（Y_c-Y_a）abs（c-a）/2）;
    if（abs（c-a）/2<=tolerance）
        fprintf（‘Tolerance is satisfied.\n‘）;
        break
    end
   if（it>it_limit）
       fprintf（‘Iteration limit exceeded.\n‘）;
       break
   end
   if  （Y_a*Y_b<=0）
       c=b;Y_c=Y_b;
   else
       a=b;Y_a=Y_b;
   end
    end
fprintf（‘Final result;Root=%12.6f\n‘b）
end
x=b;
plot（[xx][0.05*wid_y0.2*wid_y]‘g‘）;text（x0.25*wid_y‘Final solutin‘）
plot（[x（x-wid_x*0.004）][0.05*wid_y0.09*wid_y]‘r‘）    % arrow line
plot（[x（x+wid_x*0.004）][0.05*wid_y0.09*wid_y]‘r‘）   % arrow line
            
        
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1872  2017-12-06 10:36  对分法求根\bisec_g.m

     文件         41  2017-12-06 10:01  对分法求根\ex_5x2.m

     文件         29  2017-12-06 09:59  对分法求根\main.m

     目录          0  2017-12-06 10:37  对分法求根

----------- ---------  ---------- -----  ----

                 1942                    4