%精確一維搜索法的子函數0.618法gold.m
%輸入的變量x為初始迭代點是四維的向量d為初始迭代方向是四維的向量
%輸出變量是在[010]區間上使函數取得極小值點的步長因子
function?alfa=gold（xd）
f=‘（（x（1）+alfa*d（1））+10*（x（2）+alfa*d（2）））^2+5*（（x（3）+alfa*d（3））-（x（4）+alfa*d（4）））^2+（（x（2）+alfa*d（2））-2*（x（3）+alfa*d（3）））^4+10*（（x（1）+alfa*d（1））-（x（4）+alfa*d（4）））^4‘;
a=0;
b=10;
tao=0.618;
lanmda=a+（1-tao）*（b-a）;
mu=a+tao*（b-a）;
alfa=lanmda;
m=eval（f）;
alfa=mu;
n=eval（f）;
while?1
????if?m>n
????????if?abs（lanmda-b）<1e-4
????????????alfa=mu;
????????????return
????????else
????????????a=lanmda;
????????????lanmda=mu;
????????????m=n;
????????????mu=a+tao*（b-a）;
????????????alfa=mu;
????????????n=eval（f）;
????????end
????else
????????if?abs（mu-a）<1e-4
????????????alfa=lanmda;
????????????return
????????else
????????????b=mu;
????????????mu=lanmda;
????????????n=m;
????????????lanmda=a+（1-tao）*（b-a）;
????????????alfa=lanmda;
????????????m=eval（f）;
????????end
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????563??2005-12-22?19:12??zuisu.m

?????文件???????1000??2005-12-22?19:01??gold.m

?????文件????????765??2005-12-22?19:20??gonge.m

?????文件????????374??2005-12-22?18:55??hesse.m

?????文件????????597??2005-12-22?19:13??newton.m

?????文件????????789??2005-12-22?20:00??ninewton.m

?????文件????????298??2005-12-22?18:48??tidu.m

-----------?---------??----------?-----??----

?????????????????4604????????????????????8