怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
%用最速下降法求解问题;
- \) l  w& ~1 w9 veps=0.1; %给出误差线;
p=[];
& i6 J0 j' j- Y2 g0 n" i( dsyms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
/ D% M3 u: C) _grd=jacobian(f,[y,z]);
% x" T) L  ~* g) ]2 g  A4 ]: f5 ?x=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
5 P/ H) A/ d0 q6 i2 u# aXX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
1 c0 ?1 M$ B8 {9 r5 Y0 M. J: F0 `      disp('所求点为x=');
  disp(x);
  else
   for k=1:20;
      p=-X./XX;
  w3 A/ B' `: c; g# I$ u6 N3 ^      H=subs(f,[y,z],x+lamd*p);
4 d/ |) Y' s5 U, I# T, w1 [     %求解最优步长lamd是关键问题;
2 T  I! m4 A# i% I     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
- z2 b7 E5 z; h$ e* \: c       x=x+lamd*p
      X=subs(grd,[y,z],x);
4 `: ^5 ?/ k1 B4 {, ]% [      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
7 Q& h- ^2 G) m1 }! u          disp('所求点为x=');
- I, v- G+ P) ?1 \5 w+ D* V/ A, H+ _      disp(x);
      end
   end
end