部分可分非线性方程组与优化问题的稀疏拟牛顿法

来源 :湖南大学 | 被引量 : 4次 | 上传用户:ethel_baby
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
拟牛顿法是求解非线性方程组和最优化问题颇受欢迎的一类算法.该类算法具有超线性收敛速度,而且,若采用适当线性搜索或信赖域技术,算法可具有全局收敛性.然而,拟牛顿法产生的拟牛顿矩阵往往是稠密的,因此不能直接用于求解大规模问题.来自科学计算和工程等领域的许多实际问题一般都具有一些特殊的稀疏结构.比如微分方程采用有限元法或者有限差分法离散化后得到的非线性方程组,其Jacobian矩阵具有一定的稀疏结构.很多无约束优化问题可写成若干个分量函数的和,其中每个分量函数只跟少数几个变量有关.因此,根据问题的稀疏特征或者目标函数的特殊结构,设计高效的稀疏拟牛顿算法,具有重要的理论意义和实际应用价值,也是优化界关注的一个重要课题,并取得了一系列重要成果.目前已提出了多种稀疏拟牛顿法,这些算法充分利用了问题的稀疏特征,并且保留了传统拟牛顿法的超线性收敛性等重要性质.稀疏拟牛顿法已成为求解大规模非线性方程组和最优化问题的一类重要算法.迄今为止,关于稀疏拟牛顿法的研究主要集中于对算法的局部收敛性质的研究,对于算法的全局收敛性质的研究工作很少.这是由于稀疏拟牛顿法产生的拟牛顿矩阵要保持问题的稀疏结构,使得传统拟牛顿法的某些重要性质如对称正定性、最小变化性等发生了改变,从而大大增加了研究稀疏拟牛顿法的全局收敛性的难度.本文在对稀疏拟牛顿法的性质进行深入分析的基础上,采用线性搜索技术,研究求解非线性方程组和最优化问题的一些重要的稀疏拟牛顿法的全局收敛性.我们首先研究求解大规模非线性方程组的稀疏拟牛顿法.Schubert方法是较早被提出的用于求解方程组的稀疏拟牛顿法,作为Broyden秩一修正公式的稀疏推广,Schubert修正公式能够精确保持方程组的Jacobian矩阵的稀疏结构.但在实际计算中,由于Schubert修正公式过于严格地保持矩阵的稀疏结构,数值表现有时不如Broyden秩一方法好.在本文中我们重点关注部分可分离的非线性方程组,借助于分块BFGS算法的思想,我们提出两种分块拟牛顿算法.当非线性方程组的导数信息不可用时,我们提出一种分块Broyden秩一算法,算法中采用分块Broyden秩一修正.该算法无需计算导数,且能够近似保持Jacobian矩阵的稀疏结构.当非线性方程组的导数信息可通过自动微分间接计算时,我们提出一种分块伴随Broyden算法,算法中采用分块伴随Broyden修正.该修正公式也可近似保持Jacobian矩阵的稀疏结构,并保留了伴随Broyden修正公式的线性不变性.我们采用一种无导数非单调线性搜索技术研究算法的全局收敛性,在适当条件下,建立了上述两种稀疏拟牛顿算法的全局收敛性.我们还证明,经过一定的迭代步后,单位步长可以取到.此时,算法在方程组解的局部还原为单位步长稀疏拟牛顿法,从而具有超线性收敛速度.我们还对算法进行数值检验,并将该两种稀疏拟牛顿算法与Broyden秩一方法和伴随Broyden方法进行了比较.结果表明,稀疏算法在迭代次数,函数值计算次数和计算时间方面均有出色表现.我们也将这两种稀疏算法与Schubert方法进行比较,数值结果进一步验证了这两种稀疏拟牛顿算法的优越性.其次,针对目标函数具有部分可分离结构的无约束优化问题,我们提出一种稀疏的投影分块PSB算法.算法针对目标函数的Hessian矩阵的稀疏结构,采用分块PSB修正.由于简单的分块PSB修正公式不能保正Hessian矩阵的近似矩阵的正定性,从而不能保证拟牛顿方向是目标函数的下降方向.因此我们对拟牛顿方向进行投影,提出一种投影的分块PSB算法,该算法可以保证产生目标函数的一个充分下降方向.在适当条件下,我们证明了采用Armijo或者Wolfe搜索的算法用于求解一致凸函数极小化问题时具有全局收敛性和超线性收敛速度.此外我们还通过数值计算,将分块PSB算法与已有的求解大规模最优化问题的著名的分块BFGS算法以及有限内存BFGS(L-BFGS)算法在30个测试问题上进行数值比较.结果表明,本文提出的分块PSB算法在迭代次数,函数值计算次数,梯度值计算次数和计算时间方面均有较好的表现.最后,我们研究求解对称非线性方程组的拟牛顿算法,基于自动微分技术,我们提出两类拟牛顿算法.首先,类似于Powell对称化技术,我们将伴随Broyden修正公式对称化,提出一种对称伴随Broyden修正.该修正公式保持了原伴随Broyden修正公式的最小变化性质和仿射不变性.此外,我们还提出一类新的伴随秩二拟牛顿修正公式,该修正公式不仅具有和BFGS修正公式类似的正定性和最小变化性质,而且能够精确保证拟牛顿矩阵与方程组的Jacobian矩阵沿拟牛顿方向一致.我们证明,在适当条件下这两种算法均具有全局收敛性和超线性收敛速度.数值结果显示,当用于对称非线性方程组求解时,这两类新算法相比于BFGS方法均具有优势.
其他文献
自从开展农村道路交通安全“两站两员”(乡镇交通安全管理站、乡镇交通安全劝导服务站,乡镇专职交通安全管理员、义务交通协管员)建设以来,交通管理触角延伸到农村地区,农村
报纸
<正>大庆石化公司党委高度重视发挥党内监督巡察作用,结合企业面临的新形势、新任务和新要求,在实践中不断探索巡察工作规律,创新适应公司实际的巡察体制机制、思维观念、方
"孝"是中国传统道德文化中的核心元素。"百善孝为先"是先人们教育后辈齐家治国、为人处世的道德标杆。但由于历史的原因,这种孝道思想被当作封建糟粕曾一度受到严重的批判和
随着大理市社会经济的高速发展,道路交通路网覆盖面的不断增加,尤其是“村村通”工程使大理市农村地区道路交通状况发生了翻天覆地的变化,从而带来农村地区各类车辆、驾驶人数量
报纸
秋胡故事,最初以不同面貌频频出现于史传、变文、绘画、诗歌、小说、戏曲等领域,最终成为中国文学史上的一个文学母题。在流传过程中,不同文体对它有着不同的处理,这种处理不
国家职业教育改革实施方案指出,随着我国进入新的发展阶段,各行各业对技术技能人才的需求越来越紧迫,职业教育重要地位和作用越来越凸显。社会对职业教育的关注度越来越高,对
以三亚蜈支洲岛为例,以马蜂窝旅游网站上相关游记文本为数据来源,应用ROST content mining 6软件对相关游记文本数据进行词频分析从而得出结论,综合来看,旅游者对景区的评价
随着我国经济的快速发展,我国的现代信息技术也得到了巨大的进步。信息技术在我国的广泛应用也给我国的各个行业带来了巨大的挑战,医院也要加强信息化建设,从而有利于医院的