半无限规划的有效数值算法研究

来源 :上海大学 | 被引量 : 5次 | 上传用户:qcolin
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
半无限规戈(Semi-Infinite Programming,简写为SIP)不仅在工程设计、最优控制、信息技术、经济均衡等领域有着广泛而直接的应用,而且对Chebyshev逼近理论、数学物理、模糊集、鲁棒优化等学术方面起着重要作用.因此,研究半无限规划的有效数值算法具有很强的应用价值,在国际上已引起学者们极大的关注和研究.许多学者利用光滑非线性规划的技术提出了求解半无限规划的各种算法,其中很多是基于非线性规划中的罚函数方法进行设计.罚函数方法具有初始点任意、结构简单、在适当的条件下可得到算法的收敛性等优点.然而,其罚参数不易调整,迭代点及最终最优解可能不满足可行性要求.为此,简金宝教授在两阶段方法[Polak,1978]、罚函数方法及可行方向法思想的基础上,提出并研究了一类强次可行方向法.该方法的初始点任意选取,能保证迭代点的可行性单调增加,且一旦迭代点落入可行域即可变成可行方向法.近年来,强次可行方向法和序列二次规划(Sequential Quadratic Programming,简写为SQP)、序列二次约束二次规划(Sequential Quadratically Constrained Quadratic Programming,简写为SQCQP)、模松弛等技术充分结合,得出了一系列成果,如强次可行模松弛SQP、强次可行模松SQCQP等方法.最近也基于罚函数对强次可行方向法中的线搜索进行了改进,提出了新型强次可行SQP算法.这些算法具有收敛性好、数值稳定性好等优点,其中强次可行方向法和模松弛SQP的部分成果已被应用到斜拉桥索力优化和滤波器组的设计.因此,研究如何将上述强次可行系列算法的优秀思想更好地应用到半无限规划的算法构建,以形成半无限规划一个新的求解体系,具有一定的理论意义和应用价值.本学位论文首先针对半无限规划离散化问题提出了三个新的算法:强次可行模松弛SQP算法、强次可行模松弛SQCQP算法和基于罚函数的新型强次可行SQP算法.然后,基于离散化方法提出了一个求解半无限规划的新型两阶段SQP算法,该算法内迭代本身就是一个求解半无限规划离散化问题的新型两阶段SQP算法.最后,分别基于离散化方法和局部约化方法给出求解半无限规划的两个算法结构,并对第一个基于离散化方法求解半无限规划的算法结构进行收敛性分析.第一章给出了本文的研究背景,简要阐述了半无限规划研究的历史与现状,从算法的主要思想、步骤、相关理论等方面介绍了半无限规划的基本算法,如离散化方法、局部约化方法和交换集方法,并对其他算法进行概述.最后,分析了本文算法设计基础,给出本文的主要研究内容和结构安排.第二章结合模松弛可行SQP算法和强次可行方向法的思想,提出了一个求解半无限规划离散化问题的强次可行模松弛SQP算法.算法的初始点任意选取,每次迭代只需求解一个搜索方向子问题即可获得主搜索方向.通过修正其约束指标集来减少约束个数.基于强次可行方向法的性质,改进了现有SIP可行SQP算法中的约束指标集修正技术,不仅能够保证算法的全局收敛性,而且还大大降低了搜索方向子问题的计算成本,避免了半无限规划离散化问题因约束个数过多造成的数值困难.初步的数值试验说明算法是有效的.此外,通过精心构造高阶校正方向对算法进行改进,使其达到超线性收敛.最后,基于约束指标集的修正技术,将算法推广到离散半无限极大极小问题.第三章提出了一个求解半无限规划离散化问题的强次可行模松弛SQCQP算法.该算法适用于求解半无限规划约束函数非线性程度较高的问题.通过对二次约束二次规划(Quadratically Constrained Quadratic Programming,简写为QCQP)子问题中约束指标集的修正及约束函数二阶Hessian阵校正技术的精心设计,算法不仅大大减少了QCQP子问题的约束个数,极大降低了其计算成本,还保持了非线性规划中强次可行模松弛SQCQP算法的优点:从任意初始点开始,能保证有限步之后落入可行域;不需要任何校正方向,在适当的条件下,就能证明算法的全局和超线性收敛性.数值试验表明算法是稳定有效的.第四章提出了一个求解半无限规划离散化问题的基于罚函数的新型强次可行SQP算法.该算法的特点为:初始点任意选取;根据迭代点偏离可行域的程度,构造了两个搜索方向子问题,并将其形式统一.当迭代点远离可行域时,基于罚函数构建线搜索;反之,基于强次可行方向法的思想构建线搜索.每次迭代算法只需求解一个二次规划(QuadraticProgramming,简写为QP)子问题即可获得搜索方向,通过选择适当的约束指标集可降低计算成本.所构建的新型线搜索比强次可行SQP方法的更具有弹性.在较温和的条件下,可以证明该算法具有全局收敛性.数值试验表明算法是有效的.第五章基于离散化方法提出了一个求解半无限规划的新型两阶段SQP算法.首先给出了算法的内迭代算法,该算法本身就是一个求解半无限规划离散化问题的新型两阶段SQP算法,初始点可任意选取,所构建的QP子问题结构简单,通过选择合适的约束可大大降低计算成本,线搜索基于两阶段方法构建.在适当的条件下,可证明算法具有全局收敛性.数值试验表明算法是有效的.最后,基于离散化方法,将上述算法作为内迭代,提出了一个求解半无限规划的新型两阶段SQP算法,并给出收敛性结果.第六章基于本文前面几章提出的求解半无限规划离散化问题的算法构建求解半无限规划算法.首先基于离散化方法提出一个求解半无限规划的算法结构,其内迭代可选取本文前面几章提出的求解半无限规划离散化问题的算法.在较温和的条件下,证明了该算法具有全局收敛性.最后,基于局部约化提出一个求解半无限规划的两阶段算法结构.算法的第一阶段采用前面几章提出的求解半无限规划离散化问题的算法求得SIP的近似解,第二阶段以此近似解作为初始点求解SIP的局部既约问题,从而获得原SIP问题的最优解.第七章对本文的工作进行总结,并提出可进一步研究的问题.最后,需要指出的是,本文提出的求解半无限规划离散化问题的四个算法可获得原半无限规划问题的近似解,在精度要求不高时可直接应用,其理论依据见§6.1.2节.在精度要求高时,可将其获得的近似解作为另外一种方法的初始点进一步构建SIP算法,如§6.2节提出的基于局部约化的求解SIP的两阶段算法.
其他文献
本文首先研究基于Kac-Moody代数s1(2,C[λ-1,λ)与so(3,R)所得新的谱问题.从so(3,R)出发,借助半单Lie代数的理论,从扩展零曲率方程得到相关矩阵谱问题的Dirac孤子方程族.在一类非半单圈代数中,结合变分恒等式构造一族可积耦合的Dirac孤子方程族,给出相应的Hamilton结构及Liouville可积性.新的谱问题可以视为AKNS谱问题的推广,由此引出耦合Burger
线性锥规划是目前优化领域中最热门的研究课题之一. Nemirovskii在2006年国际数学家大会一小时报告[72]中指出,锥规划是近20年凸优化研究工作中的一项重大突破性进展,主要包括线性规划,二阶锥规划和半正定锥规划这三类常见的对称锥规划.这三类规划可以转化成统一的线性锥规划模型.线性锥规划不仅具有良好的结构特征和对偶理论,而且具备很强的数学表达能力和实用性.许多应用领域中的实际问题都可以描述
本论文的研究内容属于凸几何分析,内容涉及仿射等周不等式与OrliczBrunn-Minkowski理论Brunn-Minkowski理论是凸几何分析的核心内容,它的发展起源于对两个凸体与它们的Minkowski和之间的体积关系的研究,由此产生了著名的Brunn-Minkowski不等式,整个Brunn-Minkowski理论的奠基石.仿射等周不等式主要研究凸体的一些仿射量的极值问题,是凸几何分析理
本文主要研究离散可积方程的Lax对在研究可积特征方面的应用.建立利用方程的Lax对获得方程的无穷多守恒律和对称的方法.这些离散的可积方程包括:·Adler-Bobenko-Suris方程,Nijhoff-Quispel-Capel方程;·离散Boussinesq方程,离散非线性Schrodinger方程等多分量方程.论文首先对Adler-Bobenko-Suris链方程、Nijhoff-Quisp
无网格方法是继有限元法之后发展起来的一种新型数值方法,因其形函数构造仅需要问题所在区域或边界的节点信息,不需要形成区域或边界网格,因而具有前处理简单、计算精度高、可消除体积闭锁现象等特点,是目前科学和工程计算的重要数值方法之一.移动最小二乘法是无网格方法中构造形函数的最重要方法之一,很多重要的无网格方法都是基于移动最小二乘法建立的.但是移动最小二乘法的形函数不满足Kronecker函数的性质,使得
本文使用非线性流固耦合偏微分方程组:Cauchy守恒方程和P-T/T方程,描述复杂接触表面的形变和应力分布变化,其中方程的初边值条件根据实验数据分析结果而设置。该方程组采用有限元和有限差分相结合的方法求解,即,使用有限元方法对问题的空间变量进行离散,得到常微分方程组的初值问题后,进一步使用有限差分方法对问题的时间变量进行离散,从而得到原问题的全离散格式。空间变量的离散格式取决于所使用的有限元形函数
期刊
在自然科学和工程技术的实际计算中,最基本的方法是用函数的Taylor展开的部分和作为该函数的近似,而Pade逼近则是一种特殊的有理逼近,它是Taylor多项式逼近的延伸.Pade逼近的研究和发展与数学中的解析函数论,逼近论,矩量理论,连分式以及差分方程等分支有着紧密的联系,并且它在数值分析,量子场论,临界现象和控制论等自然科学领域中已有若干功效卓著的应用.本文主要给出了基于广义逆多元Newton-
声波照射到水中目标会产生散射,不同方向的散射波声压分布与入射波性质以及目标的大小、形状、声学参数密切相关。已知入射声波、目标几何形状和周围介质的声学参数,计算散射声场分布特性是声学正问题。求解散射问题可采用解析法或数值法,但仅有少数形状规则的物体可用解析法获得场的精确解。不规则形状以及复杂材料目标的散射声场可用时域有限差分法(finite difference timedomain, FDTD)、
离散Boussinesq型系统以三元的形式可以视为定义在四方格上满足多维相容性的离散可积系统.本文主要利用双线性化方法,研究离散Boussinesq型系统的双线性结构,并给出其Casorati行列式解.本文的结构安排如下:首先我们回顾所有已知的离散的Boussinesq型方程.这些都可视为由Hi-etarinta所分类的离散Boussinesq型方程[A-2],[B-2]和[c-3]方程或其特殊情