怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
%用最速下降法求解问题;
6 h$ E. Z2 v/ m5 b+ y% `- R* Deps=0.1; %给出误差线;
p=[];
syms y z lamd
5 ?2 S/ L% k6 F% f- {. mf=exp(y^2+z^2-y*z)*cos(y+z);
grd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
; R  w" p; F% h3 T1 p2 a3 MX=subs(grd,[y,z],x); %计算x点出的梯度;
' m* |$ `" @- H, fXX=sqrt(X(1)^2+X(2)^2); %x的模;
% R& t& ?7 Z( c$ n. p# h* U  f  if XX<=eps;
2 |5 s# [5 z; W/ T1 J, s  x      disp('所求点为x=');
  disp(x);
) T; ]+ I" I0 J' Q  else
   for k=1:20;
      p=-X./XX;
1 U  m1 e$ w& P      H=subs(f,[y,z],x+lamd*p);
6 M2 L4 x3 b9 ~; v' C     %求解最优步长lamd是关键问题;
& d' L' V0 G0 w. B2 @     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
. B, t! p; j6 K- `5 y1 X       x=x+lamd*p
      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
& G( j0 g6 ]/ ?( B0 e4 L          disp('所求点为x=');
      disp(x);
      end
   end
end