Optimization with equality and inequality constraints using parameter continuation

We generalize the successive continuation paradigm introduced by Kernevez and Doedel [16] for locating locally optimal solutions of constrained optimization problems to the case of simultaneous equality and inequality constraints. The analysis shows that potential optima may be found at the end of a sequence of easily-initialized separate stages of continuation, without the need to seed the first stage of continuation with nonzero values for the corresponding Lagrange multipliers. A key enabler of the proposed generalization is the use of complementarity functions to define relaxed complementary conditions, followed by the use of continuation to arrive at the limit required by the Karush-Kuhn-Tucker theory. As a result, a successful search for optima is found to be possible also from an infeasible initial solution guess. The discussion shows that the proposed paradigm is compatible with the staged construction approach of the COCO software package. This is evidenced by a modified form of the COCO core used to produce the numerical results reported here. These illustrate the efficacy of the continuation approach in locating stationary solutions of an objective function along families of two-point boundary value problems and in optimal control problems.

[1]  H. B. Keller,et al.  NUMERICAL ANALYSIS AND CONTROL OF BIFURCATION PROBLEMS (II): BIFURCATION IN INFINITE DIMENSIONS , 1991 .

[2]  Hari Om Gupta,et al.  Optimal control of nonlinear inverted pendulum dynamical system with disturbance input using PID controller & LQR , 2011, 2011 IEEE International Conference on Control System, Computing and Engineering.

[3]  Harry Dankowicz,et al.  Staged Construction of Adjoints for Constrained Optimization of Integro-Differential Boundary-Value Problems , 2018, SIAM J. Appl. Dyn. Syst..

[4]  Matthias Gerdts,et al.  Global Convergence of a Nonsmooth Newton Method for Control-State Constrained Optimal Control Problems , 2008, SIAM J. Optim..

[5]  Brian C. Fabien,et al.  Indirect Solution of Inequality Constrained and Singular Optimal Control Problems Via a Simple Continuation Method , 2014 .

[6]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[7]  K. G. Murty,et al.  Complementarity problems , 2000 .

[8]  Thomas F. Fairgrieve,et al.  AUTO 2000 : CONTINUATION AND BIFURCATION SOFTWARE FOR ORDINARY DIFFERENTIAL EQUATIONS (with HomCont) , 1997 .

[9]  A. Ben-Tal,et al.  A unified theory of first and second order conditions for extremum problems in topological vector spaces , 1982 .

[10]  Andrew J Szeri,et al.  Optimization of acoustic scattering from dual-frequency driven microbubbles at the difference frequency. , 2003, The Journal of the Acoustical Society of America.

[11]  A. Mayne Parametric Optimization: Singularities, Pathfollowing and Jumps , 1990 .

[12]  Hexi Baoyin,et al.  Practical Techniques for Low-Thrust Trajectory Optimization with Homotopic Approach , 2012 .

[13]  J. Wilkening,et al.  A fully discrete adjoint method for optimization of flow problems on deforming domains with time-periodicity constraints , 2015, 1512.00616.

[14]  Alexandre Goldsztejn,et al.  On continuation methods for non-linear bi-objective optimization: towards a certified interval-based approach , 2014, Journal of Global Optimization.

[15]  Y. Kuznetsov,et al.  New features of the software MatCont for bifurcation analysis of dynamical systems , 2008 .

[16]  A. Szeri,et al.  Optimized translation of microbubbles driven by acoustic fields. , 2007, The Journal of the Acoustical Society of America.

[17]  J. Meditch,et al.  Applied optimal control , 1972, IEEE Transactions on Automatic Control.

[18]  S. Crescitelli,et al.  On the choice of the optimal periodic operation for a continuous fermentation process , 2010, Biotechnology progress.

[19]  Defeng Sun,et al.  On NCP-Functions , 1999, Comput. Optim. Appl..

[20]  R. Epenoy,et al.  New smoothing techniques for solving bang–bang optimal control problems—numerical results and statistical interpretation , 2002 .

[21]  Thomas Lauß,et al.  The Discrete Adjoint Gradient Computation for Optimization Problems in Multibody Dynamics , 2017 .

[22]  Frank Schilder,et al.  Recipes for Continuation , 2013, Computational science and engineering.

[23]  E. Doedel,et al.  Optimization in Bifurcation Problems using a Continuation Method , 1987 .

[24]  Michael C. Ferris,et al.  Engineering and Economic Applications of Complementarity Problems , 1997, SIAM Rev..

[25]  C. Hillermeier Generalized Homotopy Approach to Multiobjective Optimization , 2001 .

[26]  Pier Luca Maffettone,et al.  A critical appraisal of the ? -criterion through continuation/optimization , 2006 .

[27]  David A. Ham,et al.  Automated Derivation of the Adjoint of High-Level Transient Finite Element Programs , 2012, SIAM J. Sci. Comput..