怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
4 _" N3 c% `( A7 f; A- C( k%用最速下降法求解问题;
3 z/ A% ]  T+ H% c) `' Eeps=0.1; %给出误差线;
p=[];
  X2 n% O- Z  h; i( N  e5 `6 Jsyms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
# I& V. a" }' ogrd=jacobian(f,[y,z]);
x=[1 1];%给出起始点;
! }* f2 `/ Q( E% B. a' \$ jX=subs(grd,[y,z],x); %计算x点出的梯度;
( d7 ]' J' R2 x, I7 Z- h" u# wXX=sqrt(X(1)^2+X(2)^2); %x的模;
: N# U( Y( J" v9 j$ R9 t2 C  if XX<=eps;
      disp('所求点为x=');
! e) {1 ]& B8 m( x8 s3 t3 F  T  disp(x);
, [0 w% k: C, C7 ~6 W8 {  else
   for k=1:20;
      p=-X./XX;
      H=subs(f,[y,z],x+lamd*p);
  ^: g$ \" t/ S6 S7 I$ t9 _     %求解最优步长lamd是关键问题;
  ?7 O. x/ z3 h( \3 Z0 z/ A     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
+ t) n1 b( Y9 M4 n+ T  a- m+ a       x=x+lamd*p
      X=subs(grd,[y,z],x);
8 ~5 _' b3 M* H2 N      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
0 W5 m  H! e+ m) m/ u7 n          disp('所求点为x=');
4 I& u0 B: D- y3 o0 X5 O" U8 H      disp(x);
+ w2 A' g$ ~" A9 d; s& `( U" h      end
   end
; {( \4 C: T- F2 I2 M7 J end