function?best_x=golddiv（xd）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%黃金分割法求單峰下降區間的極小點%%%%%%%%%%%%
%輸入：
%[jintui_ajintui_b]=jintuifa（xd）
%function?best_x=golddiv（xd）
%?目標函數：f=mrf_obj（xi）
%先由進退法搜索出單峰下降區間，再在單峰下降區間里面利用黃金分割方法...
%搜索出迭代步長，得到的步長為精確步長
%黃金分割法
%best_x為最優的x值，x為初始向量，d為搜索方向
%精確度：epsilon
%當保留的區間長度|b-a|<=epsilon時停止迭代
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
global?i
[jintui_ajintui_b]=jintuifa（xd）;
epsilon=1e-8;%設置默認的epsilon
z2=jintui_a+0.618*（jintui_b-jintui_a）;%h1左點
f2=mrf_obj（x+z2*di）;
z1=jintui_a+0.382*（jintui_b-jintui_a）;%h2，右點
f1=mrf_obj（x+z1*di）;
while（abs（jintui_b-jintui_a）>epsilon）?%abs?a=0b=1?golden?linesearch如果區間不夠小
????if?f1????????jintui_b=z2;%右點往左縮進
????????z2=z1;
????????f2=f1;
????????z1=jintui_a+0.382*（jintui_b-jintui_a）;%轉h2
????????f1=mrf_obj（x+z1*di）;
????elseif?f1==f2%如果左右相等
????????jintui_a=z1;%一起往內縮進
????????jintui_b=z2;
????????z2=jintui_a+0.618*（jintui_b-jintui_a）;
????????f2=mrf_obj（x+z2*di）;%轉h1
????????z1=jintui_a+0.382*（jintui_b-jintui_a）;
????????f1=mrf_obj（x+z1*di）;
????else
????????jintui_a=z1;
????????z1=z2;
????????f1=f2;
????????z2=jintui_a+0.618*（jintui_b-jintui_a）;%h5
????????f2=mrf_obj（x+z2*di）;
????end
end%while
best_x=（jintui_a+jintui_b）/2;%最優的x值
%best_fx=subs（yxbest_x）;%最優的函數值

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1593??2014-09-26?13:33??無約束0dB\golddiv.m
?????文件????????1439??2014-09-26?13:33??無約束0dB\jintuifa.m
?????文件?????3320304??2014-09-26?13:33??無約束0dB\m.mat
?????文件?????????113??2014-09-26?13:33??無約束0dB\mrf_obj.m
?????文件?????3245338??2014-09-26?13:33??無約束0dB\smooth_lnm.mat
?????文件?????3598403??2014-09-26?13:33??無約束0dB\syn_d.mat
?????文件????????3192??2014-09-26?13:33??無約束0dB\Untitled.m