怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
%用最速下降法求解问题;
eps=0.1; %给出误差线;
p=[];
syms y z lamd
" M7 R4 K; ]+ ~& ]/ if=exp(y^2+z^2-y*z)*cos(y+z);
0 u4 ?2 w' Y0 f0 v  zgrd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
  O& o+ I  S- \( _7 [/ x! [XX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
/ i' L' @. Z1 y  e2 m- q/ A      disp('所求点为x=');
  disp(x);
! \/ M2 {. @. N; `* M3 ~. g3 F& ?6 h  else
. l5 D: X4 ~% q5 Z: u# B6 n   for k=1:20;
! I3 [) b& A# U9 C% Y3 r      p=-X./XX;
3 p: {+ O8 a. m: [; \      H=subs(f,[y,z],x+lamd*p);
0 u, c+ a0 w& m# \     %求解最优步长lamd是关键问题;
$ u* e+ y* S5 `& o& F3 F% T3 p     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
9 L* b' ^" C0 j2 z& @6 S: D$ j       x=x+lamd*p
      X=subs(grd,[y,z],x);
" F9 N4 I! d, T# u" q7 }% k      XX=sqrt(X(1)^2+X(2)^2); %x的模;
* q/ S2 D% h. K+ r      if XX<=eps;
. |1 f2 I& t; D- H! n          disp('所求点为x=');
- _$ ^' f7 X4 d: l0 }      disp(x);
  O, v9 n! p; h. n6 [7 [      end
   end
end