[代码资源] 怎么进行一位精确搜索

[复制链接]

字体大小: 正常放大

silence123

7 主题	8 听众	317 积分

升级 5.67%

TA的每日心情

	慵懒 2016-10-11 09:39

签到天数: 173 天

[LV.7]常住居民III

自我介绍: 想学习数学建模

群组: 数学建模培训课堂2

群组: 2011年第一期数学建模

群组: 数学建模培训课堂1

电梯直达

1^#

发表于 2013-4-6 09:23 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

fminbnd 这里不会写了
%用最速下降法求解问题;
eps=0.1; %给出误差线;
p=[];
syms y z lamd
f=exp(y^2+z^2-y*z)*cos(y+z);
grd=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=');
  disp(x);
  else
for k=1:20;
   p=-X./XX;
   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);
   XX=sqrt(X(1)^2+X(2)^2); %x的模;
   if XX<=eps;
      disp('所求点为x=');
   disp(x);
   end
end
end

zan

搜索