怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
8 g1 C! ?4 W. ~6 K  v3 K%用最速下降法求解问题;
" e8 K& t0 w* e7 }& z6 j4 ueps=0.1; %给出误差线;
& }: w" T0 b3 ]! Ap=[];
3 v: P5 e9 T1 E+ k8 f& csyms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
( @' Q8 l2 E& L. h, q) |grd=jacobian(f,[y,z]);
& {+ l5 }1 N& x$ Z. L) J9 ax=[1 1];%给出起始点;
5 ^. u7 P* I& ]2 U2 wX=subs(grd,[y,z],x); %计算x点出的梯度;
XX=sqrt(X(1)^2+X(2)^2); %x的模;
( R6 c& y) m' W7 Z  if XX<=eps;
      disp('所求点为x=');
' ?# s3 O/ X, V4 X. u3 ~$ H  disp(x);
2 Y' T9 A3 j+ e1 Q/ u6 z  else
   for k=1:20;
3 ]5 B; j/ S0 {( ^( B% \5 i; d      p=-X./XX;
      H=subs(f,[y,z],x+lamd*p);
  K/ D3 r8 `, z4 a- F6 ^     %求解最优步长lamd是关键问题;
' m- Q' D- ~' d" T3 h! {9 Z% s     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
0 ~( Q3 s. C) K$ w- \0 ?$ E      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
  W" W: d/ f0 {      if XX<=eps;
% X& ?4 ^/ O( Y/ p7 ]* a  c4 x          disp('所求点为x=');
      disp(x);
4 }; j$ C" Y5 J% H4 d  ~      end
   end
end
: M# E' d- B- u/ o/ t# {9 M) ]