Exponential Lawson integration for nearly Hamiltonian systems arising in optimal control

Abstract: We are concerned with the discretization of optimal control problems when a Runge-Kutta scheme is selected for the related Hamiltonian system. It is known that Lagrangian's first order conditions on the discrete model, require a symplectic partitioned Runge-Kutta scheme for state-costate equations. In the present paper this result is extended to growth models, widely used in Economics studies, where the system is described by a current Hamiltonian. We prove that a correct numerical treatment of the state-current costate system needs Lawson exponential schemes for the costate approximation. In the numerical tests a shooting strategy is employed in order to verify the accuracy, up to the fourth order, of the innovative procedure we propose.

[1]  Fasma Diele,et al.  Steady-state invariance in high-order Runge-Kutta discretization of optimal growth models , 2010 .

[2]  W. Hager.  Rates of Convergence for Discrete Approximations to Unconstrained Control Problems , 1976 .

[3]  K. Judd.  Numerical methods in economics , 1998 .

[4]  William W. Hager,et al.  Runge-Kutta methods in optimal control and the transformed adjoint system , 2000, Numerische Mathematik.

[5]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[6]  Fasma Diele,et al.  Splitting and composition methods for explicit time dependence in separable dynamical systems , 2010, J. Comput. Appl. Math..

[7]  M. Sofroniou,et al.  Sympletic Runge--Kutta Shemes I: Order Conditions , 1997 .

[8]  Philippe Michel,et al.  Discrete-time finite horizon approximation of infinite horizon optimization problems with steady-state invariance , 1994 .

[9]  Monique Chyba,et al.  The role of symplectic integrators in optimal control , 2009 .

[10]  Ngo Van Long,et al.  Optimal control theory and static optimization in economics: Bibliography , 1992 .

[11]  Ander Murua,et al.  On Order Conditions for Partitioned Symplectic Methods , 1997 .

[12]  S. Blanes,et al.  Splitting methods for non-autonomous separable dynamical systems , 2006 .

[13]  J. Frédéric Bonnans,et al.  Computation of order conditions for symplectic partitioned Runge-Kutta schemes with application to optimal control , 2006, Numerische Mathematik.

[14]  Knut Sydsæter,et al.  Optimal control theory with economic applications , 1987 .

[15]  J. D. Lawson.  Generalized Runge-Kutta Processes for Stable Systems with Large Lipschitz Constants , 1967 .

[16]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[17]  Fasma Diele,et al.  SB3A splitting for approximation of invariants in time-dependent Hamiltonian systems , 2010, Appl. Math. Comput..