PDE和RPDE最優(yōu)控制問題的交替方向乘子法
[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.
【學(xué)位授予單位】:吉林大學(xué)
【學(xué)位級(jí)別】:博士
【學(xué)位授予年份】:2017
【分類號(hào)】:O241.82;O232
【相似文獻(xiàn)】
相關(guān)期刊論文 前10條
1 羅金火,潘立平;終端受限的線性-非二次最優(yōu)控制問題[J];復(fù)旦學(xué)報(bào)(自然科學(xué)版);2003年02期
2 邢進(jìn)生,劉人境,李晉玲;一個(gè)有效兩階段最優(yōu)控制問題的算法[J];北京電子科技學(xué)院學(xué)報(bào);2004年04期
3 佟欣;張洪光;;一類生態(tài)系統(tǒng)的最優(yōu)控制問題[J];生物數(shù)學(xué)學(xué)報(bào);2013年03期
4 俞玉森;評(píng)《最優(yōu)控制問題的計(jì)算方法》[J];數(shù)學(xué)研究與評(píng)論;1981年S1期
5 吳鐵軍,呂勇哉;一種求解帶約束最優(yōu)控制問題的算法[J];控制理論與應(yīng)用;1986年04期
6 卪亮壯;醫(yī)學(xué)中的一個(gè)最優(yōu)控制問題[J];北京航空學(xué)院學(xué)報(bào);1988年03期
7 趙寶元;氣-固反應(yīng)中的一個(gè)最優(yōu)控制問題[J];高校應(yīng)用數(shù)學(xué)學(xué)報(bào)A輯(中文版);1990年02期
8 王玲,李建國(guó),斯洛齊克;解決最優(yōu)控制問題的準(zhǔn)梯度方法(英文)[J];控制理論與應(yīng)用;1999年03期
9 楊然,周鋼,許曉鳴;求解最優(yōu)控制問題的改進(jìn)辛幾何算法[J];上海交通大學(xué)學(xué)報(bào);2000年04期
10 楊然,周鋼,許曉鳴;求解最優(yōu)控制問題的改進(jìn)辛幾何算法[J];上海交通大學(xué)學(xué)報(bào);2000年05期
相關(guān)會(huì)議論文 前10條
1 潘立平;周淵;;線性非二次最優(yōu)控制問題的一種解法[A];第二十七屆中國(guó)控制會(huì)議論文集[C];2008年
2 張寶琳;樊銘渠;;一類奇異時(shí)滯系統(tǒng)奇異二次指標(biāo)最優(yōu)控制問題的近似方法[A];第二十七屆中國(guó)控制會(huì)議論文集[C];2008年
3 李春發(fā);陳華;;古地溫度場(chǎng)系統(tǒng)的參數(shù)識(shí)別及最優(yōu)控制問題[A];中國(guó)運(yùn)籌學(xué)會(huì)第六屆學(xué)術(shù)交流會(huì)論文集(上卷)[C];2000年
4 高彩霞;馮恩民;;一類以脈沖系統(tǒng)為約束最優(yōu)控制問題的優(yōu)化算法[A];中國(guó)運(yùn)籌學(xué)會(huì)第八屆學(xué)術(shù)交流會(huì)論文集[C];2006年
5 唐萬生;李光泉;;時(shí)變廣義系統(tǒng)最優(yōu)控制問題[A];全國(guó)青年管理科學(xué)與系統(tǒng)科學(xué)論文集(第1卷)[C];1991年
6 雍炯敏;;具有狀態(tài)約束的二階半線性橢圓型方程的最優(yōu)控制問題[A];1991年控制理論及其應(yīng)用年會(huì)論文集(下)[C];1991年
7 肖華;吳臻;;一類線性二次正倒向隨機(jī)控制系統(tǒng)的最優(yōu)控制問題[A];第二十三屆中國(guó)控制會(huì)議論文集(上冊(cè))[C];2004年
8 陶世明;朱經(jīng)浩;;Canonical對(duì)偶方法與一類最優(yōu)控制問題[A];中國(guó)運(yùn)籌學(xué)會(huì)第九屆學(xué)術(shù)交流會(huì)論文集[C];2008年
9 楊富文;;求一類H~∞最優(yōu)控制問題的非迭代算法[A];1992年中國(guó)控制與決策學(xué)術(shù)年會(huì)論文集[C];1992年
10 王水;朱經(jīng)浩;;線性規(guī)劃在半定二次最優(yōu)控制問題中的應(yīng)用[A];中國(guó)運(yùn)籌學(xué)會(huì)第八屆學(xué)術(shù)交流會(huì)論文集[C];2006年
相關(guān)博士學(xué)位論文 前10條
1 李景詩(shī);PDE和RPDE最優(yōu)控制問題的交替方向乘子法[D];吉林大學(xué);2017年
2 邵殿國(guó);若干正倒向隨機(jī)比例系統(tǒng)的最優(yōu)控制問題[D];吉林大學(xué);2015年
3 鞏本學(xué);具有隨機(jī)場(chǎng)系數(shù)偏微分方程的最優(yōu)控制問題數(shù)值方法[D];山東大學(xué);2016年
4 王海洋;時(shí)間不相容的隨機(jī)控制問題和弱形式的正倒向隨機(jī)微分方程[D];山東大學(xué);2016年
5 張倩;幾類PDE約束最優(yōu)控制問題的數(shù)值方法研究[D];南京師范大學(xué);2016年
6 劉平;控制變量參數(shù)化最優(yōu)控制問題計(jì)算方法研究[D];浙江大學(xué);2017年
7 張穩(wěn);若干微分方程最優(yōu)控制問題的譜方法[D];上海大學(xué);2009年
8 郭磊;混合動(dòng)態(tài)系統(tǒng)建模、穩(wěn)定性及最優(yōu)控制問題研究[D];山東大學(xué);2006年
9 李彬;含狀態(tài)和控制約束的最優(yōu)控制問題和應(yīng)用[D];哈爾濱工業(yè)大學(xué);2011年
10 唐躍龍;兩類最優(yōu)控制問題變分離散方法的研究[D];湘潭大學(xué);2012年
相關(guān)碩士學(xué)位論文 前10條
1 張培勇;時(shí)標(biāo)上一類最優(yōu)控制問題研究[D];貴州大學(xué);2009年
2 管文君;發(fā)展方程的能控性和最優(yōu)控制問題[D];東北師范大學(xué);2015年
3 黃啟燦;數(shù)值天氣預(yù)報(bào)模式誤差項(xiàng)的最優(yōu)控制問題研究[D];蘭州大學(xué);2015年
4 方研;帶有終端角度和攻擊時(shí)間約束的協(xié)同制導(dǎo)律設(shè)計(jì)[D];哈爾濱工業(yè)大學(xué);2015年
5 夏云飛;一類滿足Lotka-Volterra互惠關(guān)系的生物種群最優(yōu)控制問題[D];哈爾濱師范大學(xué);2015年
6 邵志政;帶有非線性干擾補(bǔ)償?shù)腁DP控制方法及在風(fēng)機(jī)變槳控制的應(yīng)用[D];東北大學(xué);2014年
7 李越;基于空間分?jǐn)?shù)階擴(kuò)散方程及點(diǎn)態(tài)受限約束的三維最優(yōu)控制問題的快速算法[D];山東大學(xué);2016年
8 孫肖斌;帶擴(kuò)散的對(duì)偶模型的最優(yōu)分紅與注資[D];曲阜師范大學(xué);2016年
9 劉志博;Navier-Stokes方程約束最優(yōu)控制問題的分裂預(yù)處理迭代方法[D];南京師范大學(xué);2016年
10 蔡超;三類發(fā)展型方程的系數(shù)反演問題[D];蘭州交通大學(xué);2016年
,本文編號(hào):2217309
本文鏈接:http://www.sikaile.net/kejilunwen/yysx/2217309.html