怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
! ~4 S# K( A$ K2 x%用最速下降法求解问题;
( b+ }4 m5 k, w( ~: {8 Z: |eps=0.1; %给出误差线;
6 Y3 q5 g( O4 M: m  |8 q7 lp=[];
* ?* X+ ?6 _4 m& ~syms y z lamd
2 H5 C& d% \% V+ S: g: pf=exp(y^2+z^2-y*z)*cos(y+z);
# W' \" m! P' F! cgrd=jacobian(f,[y,z]);
9 e$ Y: E" w% M9 I; b# L0 S+ C+ Hx=[1 1];%给出起始点;
8 f) ~" Q: o* t7 `  C6 j; UX=subs(grd,[y,z],x); %计算x点出的梯度;
XX=sqrt(X(1)^2+X(2)^2); %x的模;
& J6 N/ n. q  Y2 Z  if XX<=eps;
      disp('所求点为x=');
  disp(x);
  else
   for k=1:20;
* A. ^! Z$ ?* t; ?; T0 A      p=-X./XX;
) k5 p6 @+ m; Q! Z  g/ j) l      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
% o8 \; q6 a8 f- L2 A( V' z     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
! f( `/ E9 H/ f* m       x=x+lamd*p
! R5 Q& B7 F& R! u, R      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
) ^2 Y$ }1 `% x3 k" Y+ j      if XX<=eps;
$ ]$ n) K9 x5 d6 f9 U6 Q8 V          disp('所求点为x=');
      disp(x);
      end
   end
end
6 w1 J2 i: U0 j  E  v* c! o# c