怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
4 l8 A/ C) D- H$ ?: ^1 h1 Y%用最速下降法求解问题;
! \' J# q# ^+ `eps=0.1; %给出误差线;
2 E/ e* Y  O' \7 N* m4 Qp=[];
5 C2 g7 d( ~1 N# K6 g" z! }syms y z lamd
) E: q* \: `& K" i/ K( P9 ~: [f=exp(y^2+z^2-y*z)*cos(y+z);
grd=jacobian(f,[y,z]);
2 A  }& l$ U9 E) D" ^! T" qx=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
% g* v& w" M, n$ a0 h8 DXX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
      disp('所求点为x=');
  G$ i) H/ e- [8 W  disp(x);
2 a. \2 d. `! n1 m; D4 V, y  else
   for k=1:20;
      p=-X./XX;
      H=subs(f,[y,z],x+lamd*p);
+ d2 @% i- r: K3 P. I/ n     %求解最优步长lamd是关键问题;
     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
# m/ [, ~& W7 c  S' Z      X=subs(grd,[y,z],x);
$ B/ j( s$ s) l6 _      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
! w+ {7 F) g$ O9 `0 {          disp('所求点为x=');
5 p  J2 M6 X  d# {, C! z9 D      disp(x);
; @8 i3 k; Q% h' D4 y      end
   end
end