Enrique Ponsoda

Philipp Bader, Sergio Blanes, Enrique Ponsoda

We present structure preserving integrators for solving linear quadratic optimal control problems. This problem requires the numerical integration of matrix Riccati differential equations whose exact solution is a symmetric positive definite time-dependent matrix which controls the stability of the equation for the state. This property is not preserved, in general, by the numerical methods. We propose second order exponential methods based on the Magnus series expansion which unconditionally preserve positivity for this problem and analyze higher order Magnus integrators. This method can also be used for the integration of nonlinear problems if they are previously linearized. The performance of the algorithms is illustrated with the stabilization of a quadrotor which is an unmanned aerial vehicle.

Philipp Bader, Sergio Blanes, Fernando Casas et al.

Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in significant aspects. The first family is addressed to problems with low to moderate dimension, whereas the second is more appropriate when the dimension is large, in particular when the system corresponds to a linear wave equation previously discretised in space. Several numerical experiments illustrate the main features of the new schemes.

Philipp Bader, Sergio Blanes, Enrique Ponsoda et al.

We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian function and the fundamental matrix solution is a symplectic matrix. This is a very important property to be preserved by the numerical integrators. In this work we present new sixth-and eighth-order symplectic exponential integrators that are tailored to the Hill's equation. The methods are based on an efficient symplectic approximation to the exponential of high dimensional coupled autonomous harmonic oscillators and yield accurate results for oscillatory problems at a low computational cost. Several numerical examples illustrate the performance of the new methods.