怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
; _6 U  Z0 G7 ^%用最速下降法求解问题;
. ]$ d5 P* @% d) C2 neps=0.1; %给出误差线;
' q) b* G, n2 d$ }2 op=[];
( ]# K& j7 p2 t- o8 d# Y: `$ Msyms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
) L) N) U9 S7 U$ Q3 Ogrd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
  s5 v7 ^& m: zX=subs(grd,[y,z],x); %计算x点出的梯度;
% n# B+ P: l3 f  `9 Q/ |XX=sqrt(X(1)^2+X(2)^2); %x的模;
, A0 d# K7 y9 [6 Z  if XX<=eps;
      disp('所求点为x=');
" O2 s5 L! o1 a" w# z  disp(x);
  else
   for k=1:20;
! j5 C1 o1 r+ D# V) A      p=-X./XX;
      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
3 u- `3 l$ u1 o5 K- x. H, D' C      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
: A- l2 A( [/ K3 b      if XX<=eps;
          disp('所求点为x=');
      disp(x);
      end
   end
end