function?[xvalk]=bfgs（fungfunx0varargin）
%功能:?用BFGS算法求解無約束問題:??min?f（x）
%輸入:?x0是初始點?fun?gfun分別是目標函數及其梯度;
%?varargin是輸入的可變參數變量?簡單調用bfgs時可以忽略它
%?但若其它程序循環調用該程序時將發揮重要的作用
%輸出:??x?val分別是近似最優點和最優值??k是迭代次數.
maxk=500;???%給出最大迭代次數
rho=0.55;?sigma1=0.4;?epsilon1=1e-5;?
k=0;???n=length（x0）;?
Bk=eye（n）;???%Bk=feval（‘Hess‘x0）;?
while（k????gk=feval（gfunx0varargin{:}）;?%計算梯度
????if（norm（gk）????dk=-Bk\gk;??%解方程組?計算搜索方向
????m=0;?mk=0;
????while（m<20）???%?用Armijo搜索求步長?
????????newf=feval（funx0+rho^m*dkvarargin{:}）;
????????oldf=feval（funx0varargin{:}）;
????????if（newf????????????mk=m;?break;
????????end
????????m=m+1;
????end
????%BFGS校正
????x=x0+rho^mk*dk;??
????sk=x-x0;??yk=feval（gfunxvarargin{:}）-gk;
????if（yk‘*sk>0）
????????Bk=Bk-（Bk*sk*sk‘*Bk）/（sk‘*Bk*sk）+（yk*yk‘）/（yk‘*sk）;
????end
????k=k+1;?????x0=x;
end
val=feval（funx0varargin{:}）;?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-10-29?14:18??擬牛頓法\
?????文件????????1083??2009-10-15?20:56??擬牛頓法\bfgs.m
?????文件????????1036??2009-09-07?22:15??擬牛頓法\broyden.m
?????文件?????????862??2009-09-07?10:39??擬牛頓法\dfp.m
?????文件??????????52??2009-08-29?20:40??擬牛頓法\fun.m
?????文件??????????81??2009-08-29?20:54??擬牛頓法\gfun.m
?????文件?????????785??2009-09-06?20:20??擬牛頓法\sr1.m