%********************************************************************
%?王俊波?2017-6-23?于?Canterbury@Kent
%********************************************************************
function?[vcm]=backtracking（fungfunxdkt）
%功能：Backtracking?Line?Search（回溯直線搜索）
%?fun?????????基于障礙法的等價優化目標函數名
%?gfun????????等價優化目標的梯度函數名
%?x???????????搜索的起點
%?dk??????????x處的搜索方向，由Newton法給出
%?t???????????與內點法有關的參數
%輸出：?
%?v???????????步長
%?cm??????????循環次數???????
%********************************************************************

alpha=0.2;%取值范圍：0-0.5
beta=0.5;?%取值范圍：0-1
v=1;?%步長初始值
cmmax=10000;%最大迭代次數

cm=0;
%非常重要，由于在搜索過程中，有可能初始步長過大引起問題的約束條件不滿足，
%此時，障礙函數的對數取值非法，因此必須有此預處理，確保在可行區間內搜索
while（~isreal（feval（funx+v*dkt）））%若步長過大引起問題非法，縮小步長
??v=beta*v;
end;

while?（feval（funx+v*dkt）>feval（funxt）+alpha*v*feval（gfunxt）‘*dk）
????v=beta*v;
????cm=cm+1;
????if?cm>cmmax?
????????break;
????end;
end;

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

?????文件???????1100??2017-06-23?15:04??IPM_SampleV1.1\backtracking.m

?????文件????????283??2017-06-23?15:43??IPM_SampleV1.1\fun1.m

?????文件????????311??2017-06-23?15:50??IPM_SampleV1.1\fun2.m

?????文件????????278??2017-06-23?15:43??IPM_SampleV1.1\gfun1.m

?????文件????????316??2017-06-23?16:09??IPM_SampleV1.1\gfun2.m

?????文件????????289??2017-06-23?15:43??IPM_SampleV1.1\Hess1.m

?????文件????????364??2017-06-23?16:10??IPM_SampleV1.1\Hess2.m

?????文件???????1118??2017-06-23?15:04??IPM_SampleV1.1\ipm.m

?????文件???????1261??2017-06-23?15:04??IPM_SampleV1.1\newton.m

?????文件???????1376??2017-06-23?16:33??IPM_SampleV1.1\test1.m

?????文件???????1422??2017-06-23?16:33??IPM_SampleV1.1\test2.m

?????目錄??????????0??2017-06-23?16:52??IPM_SampleV1.1

?????文件??????50153??2017-06-23?16:51??IPM_SampleV1.1\maple.mw

?????文件????????160??2017-06-23?16:54??IPM_SampleV1.1\readme.txt

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

????????????????58431????????????????????14