修正DFP和Broyden族校正公式及其正定遗传性分析_第1页
修正DFP和Broyden族校正公式及其正定遗传性分析_第2页
修正DFP和Broyden族校正公式及其正定遗传性分析_第3页
修正DFP和Broyden族校正公式及其正定遗传性分析_第4页
全文预览已结束

下载本文档

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

文档简介

本文格式为Word版,下载可任意编辑——修正DFP和Broyden族校正公式及其正定遗传性分析分析拟牛顿算法正定遗传性问题,对DFP和Broyden族校正公式修正,得出相应的正定性一般性结论和正定遗传性的等价条件,并证明白等价条件的结论.

校正公式;修正;拟牛顿算法;正定遗传性

1.引言

我们知道,对于DFP校正公式,由Hk的正定性要求对称矩阵Hk+1正定的等价条件是sTkyk>0.

这里记yk=gk+1-gk,sk=xk+1-xk,gk=f(xk),Hk+1yk=sk.

在Broyden族校正公式中:

设Hk+1=Hk+asksTk+b(HkyksTk+skyTkHk)+cHkykyTkHk.

由拟牛顿条件Hk+1yk=sk,假定Hkyk,sk线性无关,引入一个参数,则得到关于的校正公式:

Hk+1=Hk+sksTksTkyk-HkykyTkHkyTkHkyk+vkvTk

=HDFPk+1+vkvTk

=HBFGSk+1+(-1)vkvTk.

在文献[1]中定理5.5.2.

定理设Hk正定,对Broyden族校正公式,Hk+1正定的充分必要条件是sTkyk>0且>(sTkyk)2(sTkyk)2-yTkHkyksTkBksk.

探讨了关于Hk+1正定遗传性的等价问题.事实上,假定Hk正定,关于Broyden族校正公式的Hk+1正定性问题有更一般性结论.

2.问题分析与几个结论

条件修正DFP校正公式:

Hk+1=Hk+sgn(sk,yk)sksTksTkyk-HkykyTkHkyTkHkyk.

(1)

结论1若=0,则对于公式(1)Hk+1正定的充分必要条件是sTkyk≠0.

证明由于Hk正定,故存在可逆下三角矩阵Lk∈Rn×Rn,满足Hk=LkLTk,对任意非零向量z∈Rn,由公式(1),

则zTHk+1z=zTLkLTkz+sgn(sk,yk)zTsksTkzsTkyk-

zTLkLTkykyTkLkLTkzyTkHkyk.

记akLTkz,bkLTkyk,有

zTHk+1z=‖ak‖2-〈ak,bk〉2‖bk‖2+sgn(sk,yk)‖zTsk‖sTkyk.

利用CauchySchwarz不等式以及z的任意性,简单证得‖zTksk‖≠0.

所以,Hk+1正定zTHk+1z正定

‖ak‖2-〈ak,bk〉2‖bk‖2+

sgn(sk,yk)‖zTsk‖sTkyk>0

sTkyk≠0.

由(1)式中引入参数=-bsTkyk,

则a=1sTkyk+yTkHkyk(sTkyk)2,c=-1yTkHkyk.

由Hk+1正定性条件,可推得关于a,b,c的取值有如下结果:

(1)当=0时,有Hk+1=HDFPk+1,这时Hk+1正定等价于sTkyk>0,

则a=1sTkyk>0,b=0,c=-1yTkHkyk0,b=-1sTkyk1时,易知Hk+1正定.

综合上述(1)~(4),我们可以得到如下结论:

结论2若Hk正定,则Hk+1正定sTkyk>0且≥0.

(5)当0且

>(sTkyk)2(sTkyk)2-yTkHkyksTkH-1ksk.

证明可参见文献[1],这个结论指出在sTkyk>0时,Hk+1保持正定性参数取得的最小值.

3.修正Broyden校正公式

条件对Broyden校正公式进行修正:

Hk+1=Hk+sgn(sk,yk)sksTksTkyk-HkykyTkHkyTkHkyk+vkvTk.

(2)

结论3对于公式

温馨提示

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

评论

0/150

提交评论