PDE和RPDE最优控制问题的交替方向乘子法
[Abstract]:In real life, most physical, medical and financial problems can be described by partial differential equations (PDEs) or stochastic partial differential equations (RPDEs). In many cases, people are not only concerned about the properties of PDEs or RPDEs solutions, but also about whether some variables in the governing equations can make other variables reach the desired state, while ensuring that This is a typical PDE or RPDE optimal control problem with minimal cost. Driven by actual demand, the PDE optimal control problem with deterministic coefficients and stochastic coefficients has attracted extensive research and attention ([44,56,58,131,144]). It is difficult to solve directly and a large number of algebraic equations need to be solved in the process of optimization. For RPDE optimal control problem, stochastic space is introduced on the basis of PDE optimal control problem, so the problem is more complicated. On the one hand, the stochastic space discretization method suitable for the structure of the problem should be considered, on the other hand, the original definite problem should be solved. The difficulty of solving the problem is also increased by the introduction of the random space. The existing algorithms, after the discretization of the random space, usually use iterative method to solve the problem. In each iteration, they need to solve a large number of PDEs, and the amount of memory occupied is very large. In this paper, two core numerical methods are introduced. Alternating Direction Method of Multipliers (ADMM) is an efficient numerical method for solving structural convex optimization problems. Based on the idea of splitting, large-scale optimization problems are decomposed into discrete variables. Several small-scale sub-problems. The method has global convergence and the worst-case convergence rate is O(1/k), where k is the iterative step. Multimode Expansion (MME) is an effective algorithm for solving RPDEs with random coefficients. The expansion coefficients satisfy a series of iteration equations, and the coefficients in the iteration equations are deterministic and the same function. Only the right end of the iteration equation contains random terms, which reduces the computational pressure caused by the inverse of stiffness matrix. In the optimal control problem, several kinds of optimal control problems with elliptic PDE constraints are solved and the complete convergence analysis is given. Secondly, combining MME method with Monte Carlo method, we extend ADMM to RPDE optimal control problem. The convergence of the algorithm is proved theoretically and the effectiveness of the algorithm is verified numerically. The first part is the second chapter of this paper. Taking Poisson equation as an example, the optimal control problem with elliptic PDE constraints is studied. The control variables are distributed control, Dirichlet boundary control and Robin boundary control respectively under unconstrained and box constraints. The problem is described as distributed control problem and boundary control problem. For control problems, the integral interval B is taken as D or (?) D respectively, and DZ is taken as DX or ds. Here, the control constraints are taken as two different cases: Uad = U (unconstrained) and Uad = {u (a) - {u | u a U < u (x) < {UB box constraints. The state equation E (y, u) = 0 is taken as the corresponding variational form of the following three types of control problems: Firstly, the above-mentioned question is considered. The problem is discretized by finite element method, and the discrete schemes of three kinds of state equations are written into a unified linear system of equations. Thus, the model problem is transformed into a finite-dimensional optimization problem in the following forms: Then, for uncontrolled constraints and box control constraints, we use ADMM and two-block ADMM to solve the above discrete problems respectively. Compared with other algorithms, the advantage of the proposed algorithm is that it avoids solving large-scale discrete systems and does not need to solve equations at each iteration. It only needs to find the inverse of the coefficient matrix twice outside the iteration cycle. K is the iteration step of ADMM. Let u* and u H K represent the solution of the original problem and the iterative solution of ADMM, respectively. Then the error estimate _*-Ru_h~k_ (L~2(D))= O(h p) + O (k~(-1/2)) holds, where p is a constant (for different models, see Chapter 2 of this paper). Finally, numerical experiments show that the proposed algorithm can deal with PDE most effectively. The second part is composed of three and four chapters of this paper. We generalize ADMM to the following RPDE optimal control problems. In the third chapter, the state equation s (zeta, y (x, zeta), u (x))= 0 in the above model is taken as the following stochastic Poisson equation. Based on Monte Carlo method and finite element discretization, the discrete optimization problem of the above model can be easily obtained. The Title is: Traditional methods (such as Monte Carlo method to discretize random space and Newton iterative method or SQP algorithm to solve discrete optimization problems) will be coupled to form a very large-scale system. In this chapter, we first use Monte Carlo method and finite element method to solve the discrete optimization problems. Then, according to the global uniformity of the discrete problem, the ADMM splitting method is used to solve the problem. Each sample corresponds to a sub-problem of the same size as the fixed PDE optimal control problem in the iterative process, and based on different samples, we can adopt parallel computing. The advantages of the proposed algorithm in this chapter are still in comparison. In ADMM, the solution of large-scale discrete systems is avoided. Only the corresponding low-dimensional sub-problems are solved for each sample in the iterative process. Particularly, the whole process does not involve solving PDEs, only the inverse of the coefficient matrix is obtained for each sample point outside the iterative cycle of ADMM, which is saved and used directly in the iterative process. The results are as follows: Theorem M is the number of Monte Carlo samples used, h is the mesh step of the finite element method, and K is the ADMM iterative step. If U * and U < are the solution of the original problem and the ADDM iterative solution respectively, then the following discrete error estimates -Ruhk L2 (D) O (M1/2) + O (h2) + O (k-1/2) are obtained. In Chapter 4, we continue to study RPDE optimal control. This paper mainly considers the stochastic Helmholtz equation in which the state equation s (zeta, y (x, zeta).U (x) = 0 is taken as the following. Three algorithms are used to solve this problem. The first algorithm is to extend the method in Chapter 3 directly to the optimal control problem with stochastic Helmholtz equation constraints. Because the model problem in this chapter needs to be solved in complex space. It needs to solve the inverse of M+1 coefficient matrix and M is the number of samples; (2) The memory cost of the algorithm is O (MN~2), and N is the degree of freedom of physical space. When M and N are very large, the amount of computation and memory required by the algorithm is unbearable. The general methods of DE optimal control problems, such as Stochastic Collocation combined with CG method, Monte Carlo method combined with SQP method, also have the above two problems. To solve the above two difficulties, MME method is used to preprocess the optimal control problem. In algorithm 2, based on MME, finite element discretization and Monte Carlo method, the original problem is solved. This method only needs to solve the inverse of the coefficient matrix twice in the whole process. It has a remarkable advantage over the general method in computation, but the memory cost of this algorithm is still O (MN~2). Then, according to some iterative characteristics of MME, we propose a new method. Using this equation, an algorithm is designed on the basis of algorithm 2. This method can transfer all the random fields in the RPDE optimal control problem to the target functional with expected coefficients, making the stochastic process easier. After simple calculation, the original stochastic problem can be changed into a PDE optimal control problem. The algorithm only needs the inverse of the coefficient matrix twice in the whole process, and does not need to solve the equations. The memory cost is only O (N~2), and it is independent of M. It essentially solves the above two difficulties and is the best of the three algorithms. In order to use the sample number of Monte Carlo method, e is the order of perturbation in random refractive index, Q is the expansion term of MME method, h is the mesh step of finite element method and K is the iteration step, then the error estimation between the optimal solution u* and the numerical solution u h Q, K obtained by the proposed algorithm is finally verified by numerical simulation. The efficiency of the algorithm.
【学位授予单位】:吉林大学
【学位级别】:博士
【学位授予年份】:2017
【分类号】:O241.82;O232
【相似文献】
相关期刊论文 前10条
1 罗金火,潘立平;终端受限的线性-非二次最优控制问题[J];复旦学报(自然科学版);2003年02期
2 邢进生,刘人境,李晋玲;一个有效两阶段最优控制问题的算法[J];北京电子科技学院学报;2004年04期
3 佟欣;张洪光;;一类生态系统的最优控制问题[J];生物数学学报;2013年03期
4 俞玉森;评《最优控制问题的计算方法》[J];数学研究与评论;1981年S1期
5 吴铁军,吕勇哉;一种求解带约束最优控制问题的算法[J];控制理论与应用;1986年04期
6 卪亮壮;医学中的一个最优控制问题[J];北京航空学院学报;1988年03期
7 赵宝元;气-固反应中的一个最优控制问题[J];高校应用数学学报A辑(中文版);1990年02期
8 王玲,李建国,斯洛齐克;解决最优控制问题的准梯度方法(英文)[J];控制理论与应用;1999年03期
9 杨然,周钢,许晓鸣;求解最优控制问题的改进辛几何算法[J];上海交通大学学报;2000年04期
10 杨然,周钢,许晓鸣;求解最优控制问题的改进辛几何算法[J];上海交通大学学报;2000年05期
相关会议论文 前10条
1 潘立平;周渊;;线性非二次最优控制问题的一种解法[A];第二十七届中国控制会议论文集[C];2008年
2 张宝琳;樊铭渠;;一类奇异时滞系统奇异二次指标最优控制问题的近似方法[A];第二十七届中国控制会议论文集[C];2008年
3 李春发;陈华;;古地温度场系统的参数识别及最优控制问题[A];中国运筹学会第六届学术交流会论文集(上卷)[C];2000年
4 高彩霞;冯恩民;;一类以脉冲系统为约束最优控制问题的优化算法[A];中国运筹学会第八届学术交流会论文集[C];2006年
5 唐万生;李光泉;;时变广义系统最优控制问题[A];全国青年管理科学与系统科学论文集(第1卷)[C];1991年
6 雍炯敏;;具有状态约束的二阶半线性椭圆型方程的最优控制问题[A];1991年控制理论及其应用年会论文集(下)[C];1991年
7 肖华;吴臻;;一类线性二次正倒向随机控制系统的最优控制问题[A];第二十三届中国控制会议论文集(上册)[C];2004年
8 陶世明;朱经浩;;Canonical对偶方法与一类最优控制问题[A];中国运筹学会第九届学术交流会论文集[C];2008年
9 杨富文;;求一类H~∞最优控制问题的非迭代算法[A];1992年中国控制与决策学术年会论文集[C];1992年
10 王水;朱经浩;;线性规划在半定二次最优控制问题中的应用[A];中国运筹学会第八届学术交流会论文集[C];2006年
相关博士学位论文 前10条
1 李景诗;PDE和RPDE最优控制问题的交替方向乘子法[D];吉林大学;2017年
2 邵殿国;若干正倒向随机比例系统的最优控制问题[D];吉林大学;2015年
3 巩本学;具有随机场系数偏微分方程的最优控制问题数值方法[D];山东大学;2016年
4 王海洋;时间不相容的随机控制问题和弱形式的正倒向随机微分方程[D];山东大学;2016年
5 张倩;几类PDE约束最优控制问题的数值方法研究[D];南京师范大学;2016年
6 刘平;控制变量参数化最优控制问题计算方法研究[D];浙江大学;2017年
7 张稳;若干微分方程最优控制问题的谱方法[D];上海大学;2009年
8 郭磊;混合动态系统建模、稳定性及最优控制问题研究[D];山东大学;2006年
9 李彬;含状态和控制约束的最优控制问题和应用[D];哈尔滨工业大学;2011年
10 唐跃龙;两类最优控制问题变分离散方法的研究[D];湘潭大学;2012年
相关硕士学位论文 前10条
1 张培勇;时标上一类最优控制问题研究[D];贵州大学;2009年
2 管文君;发展方程的能控性和最优控制问题[D];东北师范大学;2015年
3 黄启灿;数值天气预报模式误差项的最优控制问题研究[D];兰州大学;2015年
4 方研;带有终端角度和攻击时间约束的协同制导律设计[D];哈尔滨工业大学;2015年
5 夏云飞;一类满足Lotka-Volterra互惠关系的生物种群最优控制问题[D];哈尔滨师范大学;2015年
6 邵志政;带有非线性干扰补偿的ADP控制方法及在风机变桨控制的应用[D];东北大学;2014年
7 李越;基于空间分数阶扩散方程及点态受限约束的三维最优控制问题的快速算法[D];山东大学;2016年
8 孙肖斌;带扩散的对偶模型的最优分红与注资[D];曲阜师范大学;2016年
9 刘志博;Navier-Stokes方程约束最优控制问题的分裂预处理迭代方法[D];南京师范大学;2016年
10 蔡超;三类发展型方程的系数反演问题[D];兰州交通大学;2016年
,本文编号:2217309
本文链接:https://www.wllwen.com/kejilunwen/yysx/2217309.html