怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
%用最速下降法求解问题;
eps=0.1; %给出误差线;
# A; ~, P+ B0 y0 S- ~. Op=[];
syms y z lamd
8 l. Q+ x1 G6 M1 Mf=exp(y^2+z^2-y*z)*cos(y+z);
3 |  E# V0 v  L2 K6 Rgrd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
# @1 @# r* @4 i& g, w6 b3 M2 C& SX=subs(grd,[y,z],x); %计算x点出的梯度;
! q, r* Q# O, q# w  j: m# PXX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
6 |9 r: y. N* X' Z3 K  q* R      disp('所求点为x=');
  disp(x);
0 ^8 S+ u; @5 y  else
: C1 x$ @0 H! [4 i6 U   for k=1:20;
      p=-X./XX;
      H=subs(f,[y,z],x+lamd*p);
) q; u+ ~* I, W+ s- d" U     %求解最优步长lamd是关键问题;
     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
$ f) E& i: ^- s* e       x=x+lamd*p
      X=subs(grd,[y,z],x);
+ s$ W" w$ `: u: z& b      XX=sqrt(X(1)^2+X(2)^2); %x的模;
% y' e5 v/ z3 S! ~      if XX<=eps;
          disp('所求点为x=');
3 I9 Y; s/ d2 [! a6 x      disp(x);
      end
   end
end