新建 文本文档 (2).txt

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

function[x,fx,m]=huangjing(fun,a,b,ep,m_max)
    
l=a+0.382*(b-a);
u=a+0.618*(b-a);
k=1;
while k<=m_max
    if feval(fun,l)>feval(fun,u)
       if b-l<=ep
           x=u;
           break;
       else
           a=l;
           l=u;
           u=a+0.618*(b-a);
           k=k+1;
       end
   else
       if u-a<=ep
           x=l;
           break;
       else
           b=u;
           u=l;
           l=a+0.382*(b-a);
           k=k+1;
       end
   end
end
fx=feval(fun,x);
m=k;