新建 文本文档.txt

来自「matlab中GUI界面来实现运筹学中常用的几种最优算法」· 文本 代码 · 共 62 行

function c=bb(fun,a,b,ess)
ya=feval(fun,a);
yb=feval(fun,b);
if yb==0
    c=b;
end
if ya*yb>0
    disp('这不是有根区间');
else c=(a+b)/2;
  yc=feval(f,c);
   if   yc==0
       a=c;
       b=c;
       break
   else if yb*yc>0
           b=c;
           yb=yc;
       else a=c;
           ya=yc;
       end
     if (b-a)<ess
         break;
     end
 end
     c=(a+b)/2;
     err=abs(b-a);
     yc=feval(f,c);
 end












function x=bisect(fname,a,b,e)
feval(fname,a);
feval(fname,b);
if feval(fname,a)*feval(fname,b)>0
disp('函数两端的值应为异号');
end
x=(a+b)/2
while(b-a)>(2*e)
    feval(fname,x)
  if feval(fname,a)*feval(fname,x)<0
      b=x;
      feval(fname,b)=  feval(fname,x)
  else

      a=x;
   feval(fname,a)=feval(fname,x);
  end
  x=(a+b)/2
end