matlab dfp法,DFP算法及Matlab程序Word版

《DFP算法及Matlab程序Word版》由会员分享，可在线阅读，更多相关《DFP算法及Matlab程序Word版(6页珍藏版)》请在人人文库网上搜索。

1、传播优秀Word版文档 ，希望对您有帮助，可双击去除！作业二 用DFP算法求解，取，。一、求解：(1) 求迭代点x1令，得的极小值点，所以得：于是，由DFP修正公式有下一个搜索方向为(2) 求迭代点x2令，得的极小值点于是得：，所以：，因Hesse阵为正定阵，为严格凸函数，所以为整体 极小点。二、DFP算法迭代步骤如下：(1)给定初始点，初始矩阵(通常取单位阵)，计算，令k=0，给定控制误差。传播优秀Word版文档 ，希望对您有帮助，可双击去除！(2)令。(3)由精确一维搜索确定步长，(4)令。(5)若，则停；否则令 ， 。(6)由DFP修正公式得。令k=k+1，转步骤(2)3、 DFP算法m。

2、atlab程序实现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的梯度函数%初始点的梯度和函数值g0=subs(fi,x1 x2,x0);f0=subs(f,x1 x2,x0);H0=eye(2);%输出x0,f0,g0x0f0g0xk=x0;传播优秀Word版文档 ，希望对您有帮助，可双击去除！fk=f0;gk=g0;Hk=H0;k。

3、=1;while(norm(gk)ess)%迭代终止条件|gk| x0=1 1;ess=1e-6;best_x,best_fx,count=DFP(x0,ess)程序运行结果：x0 =1 1f0 =-3g0 =-4 2*第1次寻优传播优秀Word版文档 ，希望对您有帮助，可双击去除！xk =2.0000 0.5000fk =-5.5000gk =-1 -2Hk =0.8400 0.38000.3800 0.4100*第2次寻优xk =4 2fk =-8gk =0 0Hk =1.0000 0.50000.5000 0.5000结果如下：best_x =4 2best_fx =-8传播优秀Word版文档 ，希望对您有帮助，可双击去除！count =2可以看到，最优点，迭代次数2次，与前面结果一致。