怎么进行一位精确搜索

silence123

fminbnd 这里不会写了
! A, ?7 @+ w" @3 J* O%用最速下降法求解问题;
4 O* v7 y6 h$ _) F1 d2 J0 `& d% heps=0.1; %给出误差线;
# J: i2 g  u7 v0 {* T/ {  {p=[];
syms y z lamd
+ H: _- |) Q+ V% ~f=exp(y^2+z^2-y*z)*cos(y+z);
+ o+ ~  y4 l- k. T/ e9 E$ Q! Vgrd=jacobian(f,[y,z]);
2 y* |$ O/ \" A' m& D4 r! }0 Kx=[1 1];%给出起始点;
X=subs(grd,[y,z],x); %计算x点出的梯度;
XX=sqrt(X(1)^2+X(2)^2); %x的模;
4 M# G7 v+ K" b) Z  if XX<=eps;
      disp('所求点为x=');
  disp(x);
% r6 {! j) p# i2 q2 t  else
: i1 S) l3 M) M  V2 u) O( c) s   for k=1:20;
4 ?" d$ w0 |0 B3 J  H+ b      p=-X./XX;
! ^! x* Z  B. G; ?9 I      H=subs(f,[y,z],x+lamd*p);
( H% l' D# c0 G& W8 C. D: K     %求解最优步长lamd是关键问题;
3 Y+ A- X+ `- U" y1 i     lamd=fminbnd(@(lamd),myfun(lamd,x),0,1);
       x=x+lamd*p
" K) }; J9 q2 R7 ]2 J% \      X=subs(grd,[y,z],x);
      XX=sqrt(X(1)^2+X(2)^2); %x的模;
      if XX<=eps;
2 d2 ?6 w" Z0 \5 ?# R          disp('所求点为x=');
      disp(x);
  P. H2 ?% ^7 y  E      end
   end
; ?% J& i. g' g# N end