怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
+ X$ Z4 K/ ~! Q, g%用最速下降法求解问题;
eps=0.1; %给出误差线;
p=[];
syms y z lamd
# e' x; n) _. E2 p" Xf=exp(y^2+z^2-y*z)*cos(y+z);
/ B/ i7 H7 n. g9 n) egrd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
XX=sqrt(X(1)^2+X(2)^2); %x的模;
  if XX<=eps;
      disp('所求点为x=');
0 K  D7 J  O! l7 M7 k: s/ K, \; M- k  disp(x);
  else
   for k=1:20;
      p=-X./XX;
, N0 j7 X) e7 `      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
0 I+ A8 a# b2 W- s" f. D     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
" _0 F9 U& W- i- P* S# A      X=subs(grd,[y,z],x);
* O3 q( L* \' T8 p      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
8 _% V! k, S$ ?8 C2 B! x          disp('所求点为x=');
) k4 \3 X5 W* \4 e8 {8 {      disp(x);
5 w7 e0 H( b5 n- s6 X4 O  `      end
   end
end