DFP法BFGS共轭梯度法.docVIP

1 DFP法 给定控制误差ε. Step1,给定初始点x0,初始矩阵H0 通常取单位矩阵 ,计算g0,令k 0. Step2,令pk -Hkgk. Step3,由精确一维搜索确定步长ak. f xk+akpk minf xk+apk . a 0 ; Step4 令xk+1 xk+akpk. Step5 若||gk+1|| ε,则x* xk+1停;否则令 sk xk+1-xk, yk gk+1-gk. Step6,由DFP修正公式 3.39 得Hk+1.令k k+1,转Step2. 返回值函数 double min double X[] return X[0]*X[0]+X[1]*X[1]-X[0]*X[1]-10*X[0]-4*X[1]+60; 定义结构体 typedef struct ARRAY int Row,Col; double *Addr; Array; 初始化结构体变量 int Init_Array Array*pArray,int Row,int Col pArray- Row Row; pArray- Col Col; if ! pArray- Addr double * malloc Row*Col*sizeof double printf Faile in Initial Array! ; exit 1 ; return 1; 2 B