function?[mHdaXg1fv]?=?DFP（GbcXe）
%X：為初始點(diǎn)，列向量表示
%G：對(duì)稱正定矩陣；2*2
%b和c：函數(shù)參數(shù)，列向量表示
%e：精度
%m：最后導(dǎo)數(shù)值（精度）；H：參數(shù)H用于求方向；g1：導(dǎo)數(shù)值數(shù)組
%a：各階段步長；d：迭代方向；X：各階段迭代點(diǎn)；fv：函數(shù)極小值
if?nargin<5
????error（‘輸入?yún)?shù)必須為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;%第一個(gè)搜索方向
%a=-（d‘*g1）/（d‘*G*d）;%第一步步長

%----------中點(diǎn)法求零點(diǎn)----------------------------------
????df=[g1（:1）+lema*G*d（:1）]‘*d（:1）%輸出代入具體函數(shù)值的關(guān)于步長的一階導(dǎo)數(shù)，方便輸入?yún)^(qū)間
????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;
??????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????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??變尺度法劉敏\變尺度法課程設(shè)計(jì).doc

????I.A....?????77824??2013-11-17?23:40??變尺度法劉敏\最優(yōu)化課程設(shè)計(jì)要求.doc

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

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

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