凸约束的非线性方程系统的仿射内点信赖域法

来源 :上海师范大学 | 被引量 : 1次 | 上传用户:lovele
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
最优化问题是在有限种或无限种可行方案中挑选出最优的方案。它在工农业、国防、交通、金融、通信等领域都有着广泛的应用。伴随着计算机的高度发展和软件的完善,解决最优化问题在生产和生活的各个方面变得越来越重要,也成为现实。 求解最优化问题的方法有很多种,其中线搜索技术与信赖域策略是保证算法的全局收敛性的两个重要手段。本文主要针对凸约束,尤其是有界变量约束的非线性方程系统,将其转换成最优化问题,提出了各类结合非单调线搜索技术的内点信赖域法。 现代信赖域方法的基本思想是在当前迭代点的某个邻域(称为信赖域)内极小化目标函数的一个合适的二次模型,并反复校正信赖域半径,得到可接受方向步。对于仅带有线性不等式的约束优化问题,Coleman和Li在[6]提出了“双信赖域方法”,巧妙地构造了一个仿射变换矩阵,以及合理的近似二次函数和信赖域子问题,克服了不等式约束带来的困难。Heinkenschloss et al.在文献[20]中提出了一类仿射矩阵,得到了与此有关的一阶必要条件。而文献[13]中给出的仿射矩阵则帮助我们在求解信赖域子问题时不涉及有界约束。 本文给出了求解有界约束的非线性方程组的非单调信赖域内点算法。在求解子问题的过程中,通过引入一个仿射矩阵巧妙的去掉了有界约束,从而将原子问题转换成一个只具有椭球约束的信赖域子问题。解此子问题即得严格位于可行域内部的点,通过非单调线搜索获得下一迭代点并保证其有足够下降量。在合理的假设条件下,所给出的这类算法具有全局收敛性和超线性收敛速率。数值结果表明了算法的有效性。 共轭梯度法是解优化问题时的一类常用方法,由于具有算法简便,只需一阶信息,易于编程以及存储空间小等优点,共轭梯度法已经成为求解大规模问题的一种主要方法。Bulteau与Vial在[3]中构造了无约束最优化问题的共轭梯度路径,其基本思想是将标准共轭方向法应用于无约束优化目标函数的局部二次近似函数,得到一组共轭方向序列。共轭梯度路径即为该共轭方向序列的线性组合。在本文中,我们引入另一类仿射变化矩阵,构造了共轭梯度路径。考虑将信赖域子问题中的信赖域约束去掉,沿着这条路径搜索得到迭代方向。当该迭代方向步不严格可行时。利用非单调回代线搜索技术得到可接受的步长因子,并且此步长因子保证了新的迭代点有足够的下降量并且位于可行域的内部。证明了当共轭梯度路径中的参数趋于无穷时,产生的搜索方向即为牛顿步或拟牛顿步从而具有超线性收敛速率。数值测试表明算法的可行性与有效性。 Lanczos方法和共轭梯度路径法的思想启迪我们,对优化问题的近似二次模型应用Lanczos方法过程中同时应用共轭梯度法,即对问题三对角化的同时也计算出了共轭方向序列,这样我们可以得到Lanczos方向序列和共轭方向序列,由此生成一条新的路径,这里命名为Lanczos路径。这条路径有类似于共轭梯度路径的一些重要性质,这对于算法的整体收敛性和超线性收敛性的分析很重要。利用仿射Lanczos路径法求解有界变量约束非线性方程组能大大的减少计算量,在合理的假设下,本文也证明了这类算法具有整体收敛性和局部超线性收敛速率。数值测试结果表明了算法的有效性和稳定性。投影梯度法是解决最优化问题的又一类方法[4,5],算法虽然具有较快的收敛速度,但每次迭代可能要计算几次投影,这样大大增加了工作量。文献[9]与[36]给出了每次只需计算一次投影的算法。在求解问题的最优化方法中,在最优点x*处的非奇异假设是一个常用的条件([36]),文献[24,40]中,用Levenberg-Marquardt法求解凸约束非线性方程组时把非奇异假设用一个较弱的条件——局部误差界条件来代替。受此启发,改变[36]中所提供的投影牛顿法中的投影区域类似的得到投影牛顿步并结合局部误差界这一条件,本文给出了局部误差界的有界约束非线性方程组的非单调信赖域算法。证明了算法全局收敛性并且证明了在局部误差界这一较弱条件下算法具有超线性收敛速率。数值实验表明了以上所给算法的可行性和有效性。 对于具有一般凸约束的非线性方程组。[9]给出了求一般凸约束的投影梯度算法,参考此文得到搜索方向,并通过求解一个简单的信赖域子问题得到搜索步长,得到求解凸约束非线性方程组的非单调信赖域算法,证明了算法在合理假设下全局收敛且在满足增长条件(局部误差界的特殊形式)下算法超线性收敛。本文给出的凸约束非线性方程组的投影信赖域算法主要参考了文献[4]与[42]中的算法,避免直接求解信赖域子问题,通过近似求解满足两个与信赖域相关条件的步进而得到搜索方向,再通过一维搜索得到下一迭代点。所提供的投影信赖域算法全局收敛且在局部误差界条件下算法具有1.5阶的收敛速率。算法的数值测试将作为进一步研究工作。 本文最后对所做工作进行总结并提出了进一步研究方向。
其他文献
随着Intemet网络规模的迅速增长、网络技术的不断成熟完善,传统的“尽力而为”网络机制并不能满足QOS(Quality of Service,服务质量)通信的要求。所以,IP网络的OoS路由问题已经
设S是图G的一个顶点子集,若S中任意两个顶点G中均不相邻,则称S是G的一个独立集;空集φ是G的一个独立集.用Sk表示图G中基数为k的独立集数目,则σ(G)=n∑k=0Sk,其中n为G的顶点数.σ(G)
1925年夏天,周恩来与邓颖超这两位志同道合、心心相印的革命战友,经组织同意结为夫妻。可这时周恩来在广州,担任中共广东区委会军事部长兼黄埔军校政治部主任,工作忙得不可开
金融是一种价值流动、资金交易的经济活动,金融在一定程度上促进了经济的进一步发展,因此,全世界都将其置于战略位置,并且其地位不可动摇.作为每个国家的核心竞争力,金融的地
中图分类号:E232.5文献标识码:A 文章编号:  近几年来,全国各地根据本地实际情况积极摸索住房公积金会计核算模式,并取得阶段性成果,本文就公积金模式的现状及规范发展与大家共同探讨,以供商榷。   一、住房公积金会计核算管理模式的现状  现阶段,住房公积金会计核算模式大致分为两类:一类是直接归集模式;另一类是委托归集模式,两种核算模式各有利弊。  (一)直接归集模式  住房公积金管理中心(以下
期刊
本文考察了向海底无限延伸的滨海含水层系统中由海潮及其载荷效率引起的地下水头波动。该系统由承压层和它下面的弱透水层所组成,并且向海底无限延伸。本文综合考虑了弱透水
科学和工程中的许多问题可归结为外部问题,例如:流体力学中大量存在的障碍问题等。求解此类问题的最简单的方法是设定一个人工边界,加上人工边界条件,然后在有限子区域中用通常的
众所周知,积分方程解的性态是积分方程理论中一个重要而又基本的问题,其中关于方程的周期解、概周期型解的存在性问题更是具有重要的理论意义和应用价值。K.Cook和J.Kaplan于19
稳定性理论是微分方程,时滞微分方程理论研究中的一个基本而又重要的研究课题。本文主要研究了几类不等式和微分方程,时滞微分方程稳定性。首先推广了几类积分不等式,并利用积分
对于一类含单参数Hamilton平面向量场的n次多项式扰动系统(x)=-2y+3λy2+∈P(x,y),(y)=-x-x2+∈Q(x,y),其中λ为小参数,∈为扰动参数,0