基于外梯度的变分不等式求解算法

来源 :南京大学 | 被引量 : 0次 | 上传用户:kingxing
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
变分不等式问题及其衍生出来的逆变分不等式问题在数学规划、交通控制以及经济平衡等诸多领域有着广泛的应用。在过去的几十年里,对于单调变分不等式,学者们给出了很多切实有效的数值算法。现实生活中的很多问题并没有显式的函数表达式,只能根据给定的自变量,观测到相应的函数值。只用函数值的方法(例如外梯度方法)在单调变分不等式的求解中仍然起着某种不可替代的重要作用。由于函数值的观测往往是代价不菲的,在应用只用函数值的算法的同时,考虑尽量少用函数值就特别有意义。   一般的外梯度算法都包含两步:一步预测和一步校正,其中预测步给出下降的方向,校正步中确定新的迭代点。投影收缩算法的贡献则相当于在外梯度算法校正步中运用了最优步长而提高收敛速度,最终达到求解问题少用函数值的目的。本文针对结构型的变分不等式,将约束优化理论与方法的研究中已经得到公认的一些思想与技术,用到距离函数下降方向的构造中来,达到提高算法收敛速度的效果。   论文分为三个部分,第一章的内容是变分不等式和逆变分不等式的一些介绍,包括基本定义、性质以及一些成熟的算法;第二章和第三章着重讨论变分不等式的一些具体的算法和应用;第四章和第五章主要研究逆变分不等式的一些性质、算法及其应用。   第二章的主要工作是结合了增广拉格朗日乘子法(ALM)和交替方向法(ADM)来求解结构型变分不等式。在预测步中,我们首先利用ALM的思想,将约束条件考虑在内,得到了Lagrange乘子的一个初步预测点。在此基础上,根据ADM,得到一个自变量和Lagrange乘子的共同预测点,这个过程实际上是采用了一个Gauss-Seidel迭代,这样可以更加及时的利用已有的迭代信息,有效的减少调用函数值的次数,而计算量却没有增加。另外我们也注意到[131]中提出了一种新的步长选取准则,在本章中,我们根据这种准则对校正步的最优步长做了进一步的研究。注意到分析要考虑所有的情形,不可避免地会用到一些不等式,据此得到的最优步长是保守的。我们对此做了更精细的分析,以确定新的步长,虽然这需要一些额外的计算,但是在每步迭代中并不增加函数值的调用次数,而在总体上达到减少计算函数值次数的目的。数值试验表明,基于这两点改进,我们的算法有较好的表现,比原先的方法减少约10%的调用函数值的次数。   第三章中,对单调变分不等式,经典的邻近点算法一般采用形如H(μ-μκ)的邻近项进行求解,这里的H(μ-μκ)是某个距离函数的导数,因此日必须是对称正定的。在本章中,我们对每个子问题的预测步采用了一个线性邻近项M(μ-μκ),这里的M只需正定而无需对称性。因此,M(μ-μκ)也不再是一个函数的导数。这种邻近项的修改使得算法具有更广泛的适用性,可以帮助我们解决多种类型的结构型变分不等式。另外在校正步中,我们提出了两种不同的下降方向和相应的最优步长选取准则。第一种下降方向是一个线性函数,因此迭代公式里只用到简单计算而无需投影。数值试验说明该方法不仅可行,而且有效的降低了计算量。第二种下降方向在采用相同的最优步长时,理论与实践上有更好的收敛性。我们同时证明了采用这两种下降方向的方法可以归结到一个统一框架里[82],不仅在表达上更加简洁直观,而且证明了第二种下降方向的优越性。在这个统一框架下,校正步采用第二种下降方向,并且利用第一种下降方向来计算最优步长。数值试验表明,算法的校正步经过这样的修改后表现更好。   在第四章中,我们首先从一些实际例子中引入逆变分不等式的概念,然后通过分析得出,它实际上可以扩展为高维的变分不等式。在此基础上,我们讨论了逆变分不等式的解的一些性质,并且对一类线性逆变分不等式提出一种简单的投影算法,我们给出了下降方向和步长选取准则,同时证明了算法的收敛性。数值试验也表明,该方法是切实可行的。   在第五章中,我们将逆变分不等式推广到一个连续化模型,对其采用逆变分不等式的算法来求解,在相当宽松的条件下,这个模型是全局收敛且渐近稳定的。
其他文献
风险测度是一个广义的概念,它可以用不同的测量方式表达。20世纪90年代发展起来的风险价值(VaR)方法就是一种衡量风险的新测量方法。它是金融史上风险可量化的一个重要标志,它
近年来,随着国民经济的飞速发展,一维下料问题在建筑、电力、水利等领域获得了越来越广泛的应用。寻找一种最优的下料方案,不仅可以节省原材料,降低生产成本,而且能够为企业带来直
学位
矩阵广义逆和算子广义逆在理论和应用方面都有十分重要的地位,因此也就得到了许多学者的关注。现在由于非交换微分几何的发展,人们需要对Fredholm模的陈指标进行深入的研究,比如
众所周知,扩散和时滞现象在事物的演化过程中往往是不可避免的,因此时滞反应扩散方程引起了众多学者的关注,其中最关注的就是行波解的存在性问题.然而,自上个世纪九十年代以来,人
冠心病是一类由遗传与环境因素相互作用引起的复杂疾病,是世界范围内死亡和致残的一个重要原因。对冠心病的全基因组关联研究是近年来的研究热点。  以往的冠心病全基因组关
本文主要研究了半导体Doping Profile的反演。我们对瞬时Drift-diffusion系统建立了全局Carleman估计,并且用此估计得到了可反演性的结果。我们得到的结果是Doping Profile是L
Rees矩阵半群是一种重要的半群.许多人研究过Rees矩阵半群的结构和性质,得到了许多很好的结论.其中最著名的结论是群上的在Rees矩阵半群和完全单半群是等价的.在Rees矩阵半群中,幺
近年来用遗传算法解决神经网络的优化设计问题受到广泛重视,尤其在预测领域的应用。本文将就这一课题进行进一步地研究和探索,充分利用遗传算法的全局搜索特性,达到对前馈神
本论文研究了几类在时间尺度上的具有反周期边值条件的脉冲微分方程的最值解的存在性及其相关问题,并得到了一系列新的结果。 本论文的结构如下. 第一章,应用单调迭代技巧
分红问题一直是保险公司研究的主要问题。从最初De Finitti开始提出分红策略,到Gerber首次研究了经典风险模型中的最优分红问题。分红问题一直不停的在进行深入研宄。随着时间