怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
  ?6 T" p! L3 z& H) a%用最速下降法求解问题;
3 {6 z! b) R7 P/ z. H! Xeps=0.1; %给出误差线;
! {4 a! b" y0 e$ a; D/ Up=[];
. {" Q' R# s. ?) p$ t7 Zsyms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
3 p2 ~, t' g+ _grd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
) M! S) v& p7 m3 oXX=sqrt(X(1)^2+X(2)^2); %x的模;
7 H" r3 J- p( U  if XX<=eps;
5 D6 u. B$ j7 t      disp('所求点为x=');
  disp(x);
  else
  _$ Q8 U9 }, {: U, k   for k=1:20;
5 C6 l' ]) Z2 v7 Z0 L/ {# H      p=-X./XX;
4 t/ ~2 o$ `! H/ v+ w' ]      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
% b  X! D% |7 R! A2 v      X=subs(grd,[y,z],x);
3 M$ x/ L( i" j      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
! o$ O& N  W( e- n6 L  v/ z          disp('所求点为x=');
      disp(x);
1 P& L+ K* M! l  d- O5 W) s+ J6 m      end
   end
3 u8 r! W3 c. V  i2 d8 \ end