DFP算法及Matlab程序(共6页).docx

精选优质文档-倾情为你奉上作业二 用DFP算法求解，取，。一、求解：（1） 求迭代点x1令，得的极小值点，所以得：于是，由DFP修正公式有下一个搜索方向为（2） 求迭代点x2令，得的极小值点于是得：，所以：，因Hesse阵为正定阵，为严格凸函数，所以为整体 极小点。二、DFP算法迭代步骤如下：（1）给定初始点，初始矩阵（通常取单位阵），计算，令k=0，给定控制误差。（2）令。（3）由精确一维搜索确定步长，（4）令。（5）若，则停； 否则令 ， 。（6）由DFP修正公式得。令k=k+1，转步骤（2）3、 DFP算法matlab程序实现function best_x,best_fx,count=DFP(x0,ess)syms x1 x2 t;f=x1*x1+2*x2*x2-2*x1*x2-4*x1;fx=diff(f,x1);%求表达式f对x1的一阶求导 fy=diff(f,x2);%求表达式f对x2的一阶求导fi=fx fy;%构造函数f