DFP算法及Matlab程序.doc_第1页
DFP算法及Matlab程序.doc_第2页
DFP算法及Matlab程序.doc_第3页
DFP算法及Matlab程序.doc_第4页
DFP算法及Matlab程序.doc_第5页
免费预览已结束，剩余1页可下载查看

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

作业二 用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的梯度函数%初始点的梯度和函数值g0=subs(fi,x1 x2,x0);f0=subs(f,x1 x2,x0);H0=eye(2);%输出x0,f0,g0x0f0g0xk=x0;fk=f0;gk=g0;Hk=H0;k=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次寻优xk = 2.0000 0.5000fk = -5.5000gk = -1 -2Hk = 0.8400 0.3800 0.3800 0.4100*第2次寻优xk = 4 2fk = -8gk = 0 0Hk = 1.0000 0.5000 0.5000 0.5000结果如下：best_x = 4 2best
人人文库> 全部分类> 应用文书 > 事务文书

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论