解非线性规划的同伦JFNK方法

来源 :大连理工大学 | 被引量 : 1次 | 上传用户:haolong12345
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
同伦连续方法现已发展成为广泛应用于非线性规划中的大范围收敛算法,通过构造适当的同伦,并经过追踪其数值路径可以很好地获得原始非线性规划问题的KKT点.然而随着问题规模的增加,其数值算法效率亟待提高. JFNK方法已发展成为专门求解大型稀疏非线性方程组的一类高效算法,其显著优点为:一是运用Krylov子空间迭代法非精确求解Newton方程,减少计算量;二是Krylov子空间迭代法只需Jacobian矩阵与向量的乘积,而此乘积可通过向量函数的差分来近似代替,因而无需形成和存储Jacobian矩阵,可大大节省存储容量. 本文分别选取组合同伦、凝聚同伦和平整化凝聚同伦方法为例,将Jacobian-Free Newton-Krylov方法应用于占据数值同伦连续方法主要计算量的校正步,给出了基于JFNK方法的数值同伦方法,以此改进传统Newton-PLU方法在校正阶段所显现的“过解”和计算效率低下的弊端.其中我们分别选取Krylov子空间方法中较有代表性的3类算法GMRES、BICGSTAB和TFQMR;并与传统Newton-PLU方法在Linux平台下的CUTEr-Matlab测试环境中进行了对比实验,数值结果表明改进后算法效率得到明显提升.在Krylov子空间迭代法中,相比BICGSTAB和TFQMR而言,内嵌GMRES迭代法的JFNK算法执行效果最佳,这种嵌入JFNG方法的同伦方法很值得推广至工程应用. 全文篇章结构安排如下:第一章概括介绍了本文研究的背景、选题意义及其同伦方法发展史与当前国内外研究现状;第二章介绍了三种同伦方法,为后续研究提供理论铺垫和算法支撑;第三章,结合不精确牛顿法和Krylov子空间方法给出了嵌入JFNK方法的同伦方法;第四章是全文的数值实验部分,首先详述了CUTEr测试环境的配置过程,并在此环境下进行了同伦方法改进前后的对比实验;结论部分简短地介绍本文所取得的主要成果.
其他文献
区间值模糊集是对Zadeh模糊集的一种形式上的扩充和推广,近来,人们对于区间值模糊集的研究逐渐增多,主要是因为在实际的应用过程中,数据处理的结果用区间值模糊集表示更能反映其模糊性质和不确定性,同时,基于区间值模糊集的处理同样也可以非常有效地减少模糊信息的丢失.在实际的应用方面,该模糊集已应用于决策分析,人工智能等相关领域. 本文给出了区间值模糊集的区间值水平截集的概念,建立了区间值模糊集的
学位
集值映射广义向量变分不等值(GVVTI)是变分不等式的重要推广形式,是研究多目标规划、均衡问题以及其他数学和工程领域中等问题的重要理论基础和工具,对这一问题的研究涉及到集值分析、凸分析、线性与非线性分析、非光滑分析、泛函分析等,有重要的学术价值。本文主要从理论上研究Hausdorff拓扑向量空间上一类集值映射广义向量变分不等式解的存在性问题,通过添加下半连续、C-单调等条件,运用经典的KKM-Fa
学位
本文主要研究了三维欧式空间中曲线所对应的五维Lorentz空间中的共形曲线。n维Mobius群同构于n+2维保定向的Lorentz变换群,三维空间中的球可以对应到五维Lorentz空间的点,利用三维欧氏空间中曲线的密切球串,可以把三维空间中的曲线对应到五维Lorentz空间中的曲线,并由此可以建立一个五维自然标架{T1,T2,T3,T4,T5} 本文利用这个五维自然标架场,得到共形曲线的一
学位
vu-分解算法是2000年Lemarechal, Mifflin, Sagastizabal和Oustry对具有光滑子结构的非光滑优化函数提出的.现在我们将这种方法应用到半无限问题上来,将其不光滑特征集中于v空间中,借助一个中间函数,u-拉格朗日函数,得到切于u-空间的某个光滑轨道上的二阶展式,进而设计算法并证明其收敛性.这种算法是超线性收敛的,因为在vu-分解算法中,它包括v-步和u-步,而u-
学位
随着模糊优化问题在日常生活和现代科技中的普遍应用,国内外许多学者对带有t-norm算子的模糊关系方程以及不等式约束的优化问题已经进行了很多研究,而与t-norm同等重要的s-norm算子的研究相对较少,由于s-norm与t-norm的性质不同,因此约束优化问题的可行域也会有差别。本文主要研究了带有min-Archimedean s-norm算子的方程的求解,以及带有Archimedean s-no
学位
对复发风险感知的起源与发展、评估工具、脑卒中病人复发风险感知的影响因素以及干预措施进行综述,旨在提高脑卒中病人以及临床工作者对脑卒中复发风险感知的认识,为促进脑卒中病人健康行为方式、改善其健康结局提供参考。
期刊
本文主要应用滑动区域方法证明有界区域上的退化p-laplace方程的正解在某一方向上的的单调性及该解的唯一性。我们首先可以通过所研究的方程的弱C1(Ω)解的正则性结果得到一个弱比较原理。然后应用该弱比较原理及滑动区域方法证明了此方程的正解的单调性及唯一性。
学位
混凝土是一种多相复合材料,因其方便取材、适用性强、耐久性好等优点,被广泛应用于实际工程中。由于混凝土材料的复杂性,在制作、施工以及使用等一系列过程中难免会出现一些微小的裂缝或者其他缺陷,这些裂缝或者缺陷会在某种应力状态下逐渐扩展,并最终造成构件的断裂破坏。为更好地掌握混凝土断裂破坏的全过程,国内外学者进行了大量的试验和研究,揭示出很多有关混凝土材料裂缝扩展及能量消耗等方面的规律。但是任何单一的试验
学位
在量子力学的许多实际问题中,由于体系的哈密顿算符的复杂性,大多数的Schrodinger方程和主方程是很难精确求解,甚至是不可解,因此,需要引入近似方法,微扰法在描述量子系统实际问题时是一个很重要求近似解的工具,微扰理论在线性系统中是成立的,而对于非线性系统还是一个未知的问题。 本文首先介绍了非线性二能级系统模型和一种简单的非线性开放二能级系统模型,引入了两种最常用的微扰方法一定态Schr
学位
2006年,于波、商玉凤提出了求解非凸规划问题的动边界组合同伦方法,在较弱条件下证明了同伦路径的存在性和收敛性.与已有的拟法锥条件、伪锥条件下的修正组合同伦方法相比,所需的条件更弱,同伦构造更容易,并且不要求初始点是可行集的内点,因此它更便于应用. 在求解仅带有不等式约束的动边界组合同伦方法中,为了保证全局收敛性,需要构造一类具有下述两个性质的动约束函数:积极约束满足正独立性及初始约束集合
学位