怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
) v* o6 I+ U- Z8 o- }7 l%用最速下降法求解问题;
eps=0.1; %给出误差线;
p=[];
- T2 b* ?3 o2 j4 e+ [1 `: ^syms y z lamd
7 R" c2 c# y0 k7 M% |f=exp(y^2+z^2-y*z)*cos(y+z);
) B+ W% n( i  p9 P7 ]grd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
8 f( I% p! c: BX=subs(grd,[y,z],x); %计算x点出的梯度;
9 [5 J9 |. Y7 B* Q- BXX=sqrt(X(1)^2+X(2)^2); %x的模;
) T* P) o- g7 ?7 B' t, s  if XX<=eps;
      disp('所求点为x=');
# \, M$ h& `, \$ n$ H7 u  disp(x);
  else
, m3 Z% x0 I0 d   for k=1:20;
      p=-X./XX;
) Q- C% y+ Z5 O2 A      H=subs(f,[y,z],x+lamd*p);
     %求解最优步长lamd是关键问题;
     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
      X=subs(grd,[y,z],x);
% g5 l- _  V* _( b# U      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
          disp('所求点为x=');
' {4 N$ B8 p8 d" M1 z' o% c' T      disp(x);
1 o! d/ C2 Z/ Z0 l      end
   end
3 w8 b$ P" g) ?8 i' ^0 k9 ^ end