cut.m

来自「单变量寻优基本算法的原码」· M 代码 · 共 25 行

x0=1
h=1
f0=df(x0);
x1=x0;
x2=x0;
while(df(x1)*df(x2)>0)   %寻找解区间
   if(f0<0)
      x2=x1+h;
   else
      x2=x1-h;
   end
end

for i=1:10
x3=x2-(x2-x1)*df(x2)/(df(x2)-df(x1));
fz=df(x3)
   if(fz>0)
       x2=x3;
   else
       x1=x3;
   end
end