最优控制问题和随机优化问题的数值方法研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:sisi200713
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文考虑两类优化问题的数值方法:最优控制问题和随机优化问题.最优控制问题在过去几十年内迅速地发展为应用和计算数学领域中的重要分支.最优控制问题是科学研究,工业与工程,金融投资等众多学科或过程中涉及到具体实现的关键问题,例如地球物理,气候科学,材料科学,医学成像,形状设计,机械制造,期权定价等等.层出不穷的实际问题,持续地推动着最优控制问题的理论和方法的发展,不停地激发着研究者们设计和分析求解最优控制问题的数值方法的热情与动力.最优控制问题的目标是:以期望的方式影响或者控制由偏微分方程描述的系统,找到系统内满足要求的最优的“控制”,从而使得系统内的“状态”尽可能地接近或者达到理想的状态.我们主要考虑带Possion方程约束最优控制问题的数值求解方法,其形式可以抽象地表述为以下Banach空间中的线性约束凸优化问题:(?)J(y,u)=Φ(y)+Ψ(u)s.t.Ay+Bu=p,u∈uad.这里函数u表示系统中的“控制”,函数y表示系统中的“状态”,线性约束Ay+Bu=p由不同类型的Possion方程的变分形式描述.我们借助有限元方法对连续形式的最优控制问题进行离散,得到Rn中的线性约束优化问题:(?)Jh(yh,uh)=Φh(yh)+Ψh(uh)s.t.Ayh+Buh=ph,uh∈Uadh,这里矩阵A,B为有限元离散后形成的刚度矩阵,yh,uh,ph为Rn中的向量.可以明显地看出,上述问题具有可分离结构,一个自然的想法是借助交替方向乘子法求解这类问题.交替方向乘子法是一类广为人知的优化算法,能够高效地求解带线性约束的可分离结构型凸优化问题.对于不同控制类型的Possion方程约束最优控制问题的离散形式,本文设计了不同的交替方向乘子法迭代格式进行求解.本文数值方法的误差分析,包括有限元离散的逼近误差和交替方向乘子法的迭代误差,同时以目标泛函的形式|J(y*,u*)-Jh(yhK,uhK)|给出总体的误差估计,其中(y*,u*)表示Possion方程约束最优控制问题的最优解,(yhK,yhK)表示离散问题的迭代解,符号K表示迭代总次数.针对本文设计的数值方法,我们进行不同的数值试验,试验结果表明我们的数值方法能够有效地数值求解该类最优控制问题.随机优化问题的身影在人工智能领域随处可见.统计学习,机器学习和深度学习等技术的实现依赖于随机优化问题的求解,由此产生的智能系统:搜索引擎、推荐平台、语音和图像识别软件等等,已经成为我们日常生活不可或缺的一部分.随机优化方法在不同学习模型中的成功应用和实践,始终推动着研究人员寻求解决不同随机优化问题的新方法.本文提出分块镜像随机梯度方法用以求解随机优化问题(?)f(x)=E[F(x,ξ)],和随机复合优化问题#12与一般的问题不同,本文考虑问题的可行域χ具有可分结构,而且变量x相应地能够分为若干独立的子块.本文设计的分块镜像随机梯度方法吸收了经典的镜像随机逼近方法和分块坐标梯度下降方法的特点,前者能够高效地求解随机优化问题或者大规模优化问题,后者处理变量拥有多个分块结构的优化问题的效果显著.基于随机一阶策略构造的随机次梯度,分块镜像随机梯度方法在每一步迭代中以Gauss-Seidel类的更新方式循环地更新每个子块,并且可以任意的选择各个子块的迭代顺序.我们分析了该算法对于随机优化问题在一般的凸情形和在强凸情形下的收敛性,同时也建立了对于随机复合优化问题在一般凸情形,强凸情形以及非凸情形下的收敛性分析.对于一般的凸情形以及强凸情形下的随机优化问题和随机复合优化问题,分块镜像随机梯度方法的收敛率与经典的随机梯度法完全一致.对于随机复合非凸优化问题,我们借助复合投影梯度的概念分析了算法的收敛性质,对迭代点收敛到稳定点的收敛率进行了估计.不论是凸情形还是非凸情形,本文的收敛性分析都不是平凡的.Gauss-Seidel类型的子块更新方式使得传统的梯度误差无偏假设不再成立,因此我们需要提出更有针对性的假设,同时进行相应的分析.本文的数值试验表明,分块镜像随机梯度方法能够有效地用于求解随机优化问题.与已有的随机优化方法相比较,分块镜像随机梯度方法在不同方面的表现更好.
其他文献
同核双原子分子及离子借助高对称性的阿贝尔点群D∞h下的D2h子群来精确的描述基态和激发态的电子结构。这类分子及离子的电子态波函数具有明确的奇偶性,使得电子态之间的跃迁及自旋-轨道耦合效应(SOC)具有较为严格的选择定则。碳族和氧族的双原子分子及离子的激发态具有较高的电子态密度,密集分布的电子态之间的自旋-轨道耦合作用可以导致一些束缚态的预解离,造成谱带的断裂或者弥散,增加了光谱分析和标识的难度。另
生命系统在应对自然环境变化时,通过借助自身的传感反馈机制和驱动功能来调整其复杂结构以适应新的环境。科学界试图通过参考自然的方式开发新材料和应用。然而,摆在研究者面前的难题是如何模仿自然中的活跃微环境以创造出具有动态和可调节特性的响应材料。刺激响应性聚合物,也称为“智能材料”在科学界引起了极大的兴趣。智能材料在受到外界环境条件的微小变化,如温度、光、压力、电等改变时,材料内部会做出响应的微观结构或构
本文以中亚造山带东段额尔古纳地块和兴安地块为研究区,选取了区内不同时代和不同出露位置的代表性花岗岩作为研究对象,并对其中不同单矿物进行了元素和多重同位素分析测试工作,包括长石原位Sr-Pb同位素、磷灰石原位Nd同位素和锆石原位Hf同位素。依据单矿物Sr-Pb-Nd-Hf同位素组成的时空变异,并结合花岗岩全岩地球化学特征,明确了额尔古纳地块和兴安地块的地壳属性;通过锆石Hf同位素数据,建立了额尔古纳
格子Boltzmann方法自20世纪80年代末诞生以来,迅速成为一种处理计算流体力学、计算传热学乃至偏微分方程数值求解领域的一种新型方法。由于格子Boltzmann方法相较于传统的数值求解方法具有算法并行性好、程序易于实现、边界适应性强等优势。它受到了来自各领域科技工作者的持续关注,并在进入21世纪后得到了长足的发展。本文使用格子Boltzmann方法,在非线性Schr(?)dinger方程的体系
时间序列分析在保险精算、环境污染控制、企业经营管理、流行病学等领域都有着广泛的应用.它所研究的实际数据,往往反映着某一现象的统计指标,因而,时间序列背后是某一现象的变化规律.这是一门关心动态数据的学科,比如可以应用于监控医保卡的违规使用情况、理解金融产品的动态发展机制、或是预测某一地区未来各月的降雨量.这也是一门擅长解释自然与社会现象的学科,比如令地球物理学家困惑的问题,随着陆地和海洋温度的升高,
气候变化、环境演变等重大环境问题与人类活动密切相关,是近几十年来科学家们一直研究的热点。尽管对流层的数据可以作为季节性天气预测的数据来源,但中层大气的气象数据在气候变化预测中也变得愈发重要。从21世纪以来,越来越多的数据模型增加了对平流层甚至是中间层的大气动力学研究。与对流层相比,中层大气演化速度更慢,来自中层大气的向下传输影响可能导致地表的持续且可预测性的变化。所以对中层大气特性进行观测并分析其
天然离子通道在生命体的正常生理活动中扮演着非常重要的角色,参与了许多基本的生理活动。高度专一的选择性是天然离子通道的一个吸引人的特点。天然通道由于存在不易提纯、体外容易变性等问题给直接研究蛋白离子通道的选择性带来了很多困难。然而仿生离子通道则具有稳定、结构相对简单等优点,对仿生离子通道的研究有助于我们理解天然蛋白离子通道。仿生离子通道主要分为仿生阳离子通道和仿生阴离子通道。仿生阴离子通道中大多数为
贝叶斯网络(Bayesian network,BN)是一种基于概率论和图论的数学模型,为不确定性条件下的知识表示和推理提供了一种明确清晰、图形化和可解释的表达形式。然而学习最优的BN是一个NP难问题,因此研究人员从概率推理、信息论、结构学习、参数估计、特征选择以及集成学习等角度结合启发式等学习策略提出了许多方法来进行BN的研究。但就目前已经发表的成果来看,BN还存在以下几个问题需要解决:1)目前大
李超代数,作为李代数的自然推广,是李理论的一类重要研究对象.由于李超代数在物理中的超对称问题方面有着重要的应用,故其研究成为了现代数学中的一个重要课题.在李超代数的研究中,李超代数上同调和形变理论是近年来许多学者关心的重要研究课题.本文旨在研究几类李超代数的上同调与形变理论.给定一个李超代数的模,我们可以定义李超代数的系数取自于这个模的上链,上循环,上边缘和上同调.特别地,上同调空间的维数称为Be
波方程是一类重要的微分方程,用于描述自然界中的各种波动现象,例如声波、光波、电磁波和水波等.本文主要对几类非线性波型方程,包括Camass-Holm方程,Schr(?)dinger方程及相关方程组进行定性分析,研究其行波解的存在性、解的适定性和波裂现象等.首先,考虑Camassa-Holm-Kadomtsev-Petviashvili(Camassa-Holm-KP)方程的孤波解.通过相空间分析方