function?[mHdaXg1fv]?=?DFP（GbcXe）
%X：為初始點，列向量表示
%G：對稱正定矩陣；2*2
%b和c：函數參數，列向量表示
%e：精度
%m：最后導數值（精度）；H：參數H用于求方向；g1：導數值數組
%a：各階段步長；d：迭代方向；X：各階段迭代點；fv：函數極小值
if?nargin<5
????error（‘輸入參數必須為5‘）;
end
n=length（G）;
if?n==2
format?long
syms?x1?x2?lema
f=1/2*[x1x2]*G*[x1;x2]+b‘*[x1;x2]+c;
f=f（11）;
g=[diff（fx1）;diff（fx2）];

g1=subs（subs（gx1X（11））x2X（21））;%代入x0
d=-g1;%第一個搜索方向
%a=-（d‘*g1）/（d‘*G*d）;%第一步步長

%----------中點法求零點----------------------------------
????df=[g1（:1）+lema*G*d（:1）]‘*d（:1）%輸出代入具體函數值的關于步長的一階導數，方便輸入區間
????df=inline（df）;
????lf=-1;
????rt=1;
????while?abs（rt-lf）>e
????????mid=（lf+rt）/2;
????????if?df（mid）==0
????????????break;
????????else?if?df（mid）*df（rt）<0
????????????????lf=mid;
????????????else
????????????????rt=mid;
??????????

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

?????文件???????2438??2013-11-18?23:47??變尺度法劉敏\DFP.m

?????文件???????1376??2013-11-17?22:19??變尺度法劉敏\DFP2.m

?????文件?????814080??2013-11-14?23:21??變尺度法劉敏\例題是唐老師教材_1025_擬牛頓法（DFP算法）.ppt

?????文件?????309738??2013-11-19?00:41??變尺度法劉敏\變尺度法課程設計.doc

????I.A....?????77824??2013-11-17?23:40??變尺度法劉敏\最優化課程設計要求.doc

?????目錄??????????0??2013-11-19?01:19??變尺度法劉敏

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

??????????????1205456????????????????????6