怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
%用最速下降法求解问题;
eps=0.1; %给出误差线;
  M- t) I1 ]+ @  {8 J. N2 Rp=[];
syms y z lamd
, H: c$ a3 P4 T# Q$ Vf=exp(y^2+z^2-y*z)*cos(y+z);
' K+ ]; a. t$ k' j2 O$ d( [1 ~grd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
$ z, a9 P: D: n; A( _! pX=subs(grd,[y,z],x); %计算x点出的梯度;
! M" l+ w8 I+ b% S# IXX=sqrt(X(1)^2+X(2)^2); %x的模;
% L- ], Q9 L1 Z3 T2 P  if XX<=eps;
      disp('所求点为x=');
. K5 N0 x$ O+ Z  ~+ a7 {; e0 ^  disp(x);
) H2 i/ w& Y- m* u  else
   for k=1:20;
      p=-X./XX;
      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
* q% ^" v. v/ Z1 |. K* N     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
# G. A0 `) O% F# M) {* E8 b# e       x=x+lamd*p
& U8 [' k0 u- \      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
. \& _% J( P3 s% r5 m      if XX<=eps;
# R4 U9 ^3 t/ E2 J          disp('所求点为x=');
      disp(x);
9 x9 Y$ J$ w* V& d: ]4 E. Z7 u5 n      end
   end
end