一维线性搜索进退法+黄金分割法

在最优化理论中经常用到一位搜索算法，进退法用来确定搜索区间，黄金分割法用来进行精确搜索

function t = Gold_single（ xPkeps）
[ab]=minJT（xPk）;
d=single（b-a）;
c=single（0.618）;
ta=single（a+c*（b-a））;
tb=single（a+（1-c）*（b-a））;
m=0;
% fx=matlabFunction（f）;
while d>eps 
xa=x+ta*Pk;
xb=x+tb*Pk;
fa=single3_8_1（xa）;
fb=single3_8_1（xb）; 
 disp（‘fa=‘）;disp（fa）;
 disp（‘fb=‘）;disp（fb）;
    if fa>fb
        b=single（ta）;
        ta=single（tb）;
        tb=single（a+（1-c）*（b-a））;
    elseif fa==fb
        a=single（tb）;
        b=single（ta）;
        ta=single（a+c*（b-a））;
        tb=single（a+（1-c）*（b-a））;        
    else
        a=single（tb）;
        tb=single（ta）;
        ta=single（a+c*（b-a））;
    end
    d=abs（b-a）;m=m+1;
%disp（‘d=‘）;disp（d）;
%  disp（‘a=‘）;disp（a）;
%  disp（‘b=‘）;disp（b）;
% disp（‘ta=‘）;disp（ta）;
% disp（‘tb=‘）;disp（tb）;
% if m==100
%     disp（‘???????С???‘）;
%     t = NaN;
% return;
% end
end
    t=single（（a+b）/2）;
%    disp（‘m = ‘）;  disp（m）;
    disp（‘t = ‘）;  disp（t）;
%     disp（‘a = ‘）;  disp（a）;
%     disp（‘b = ‘）;  disp（b）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件          80  2017-07-06 21:48  y.txt
     文件        1002  2017-09-24 20:34  Gold_single.m
     文件         788  2017-09-20 20:15  minJT.m
     文件        1396  2017-11-15 21:40  single3_8_1.m
     文件         249  2017-09-29 10:46  Untitled3.m
     文件          73  2017-07-06 21:48  x1.txt
     文件          81  2017-07-06 21:48  x2.txt
     文件          84  2017-07-06 21:48  x3.txt