怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
* ~' I5 C. t( [9 v9 C%用最速下降法求解问题;
0 X) D) `& H! Q2 M3 Weps=0.1; %给出误差线;
6 W7 o) H) m) ?. Fp=[];
/ t- \3 B7 |1 [' X5 bsyms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
grd=jacobian(f,[y,z]);
+ S$ }) M$ U+ R! ]x=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
XX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
/ p+ q0 C; M# m8 ~* w' ~      disp('所求点为x=');
* C. S% A$ p& Y" l  disp(x);
  else
. l$ P+ C. d7 x0 Q   for k=1:20;
      p=-X./XX;
8 x/ o8 c+ g2 P( i6 q; v      H=subs(f,[y,z],x+lamd*p);
1 F% ~4 C. R! s0 M     %求解最优步长lamd是关键问题;
+ ?( [) \. }0 {- w7 w5 s+ s     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
      X=subs(grd,[y,z],x);
' K0 d. t/ E$ ~4 j8 c) E0 E      XX=sqrt(X(1)^2+X(2)^2); %x的模;
- l  Y2 U1 ]- E! E3 X% M      if XX<=eps;
; m& J% ]2 ^# ?; a' H3 i          disp('所求点为x=');
0 Y( A+ u$ I6 X# i* @; a      disp(x);
& M* t2 g* w' L0 V  p      end
   end
4 }& w; p+ x# R; k7 ^" V end