怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
& D4 L. W; O9 ^% W1 J7 W8 I, Z& l# o%用最速下降法求解问题;
$ A' B& J% I; s5 Zeps=0.1; %给出误差线;
( n+ ^" m6 c2 r2 K+ \" H$ J& `# ip=[];
syms y z lamd
% d) ?6 C+ d' K% t8 O5 f( ~f=exp(y^2+z^2-y*z)*cos(y+z);
8 s5 z- ]1 h- Ygrd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
4 K: S+ S6 {' A; s4 g6 Z, HXX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
; Z" `; D" P  w) Y& `8 E      disp('所求点为x=');
2 k" _  n# o1 e, [5 M; c  disp(x);
  else
$ u( W0 o7 E+ x; o, A/ m   for k=1:20;
. E4 ]4 Q& X+ D; t$ {      p=-X./XX;
0 @7 B* P9 ~. n7 v) ^      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
          disp('所求点为x=');
, |& @: k. T# N! Q3 n! \( X0 {1 g      disp(x);
      end
& B6 E0 t5 f/ m: `, D   end
( l- F. ]+ ]! u end