A numerical algorithm for singular optimal LQ control systems
It provides a practical numerical method for a class of optimal control problems that are difficult to solve due to singular arcs, but the contribution is incremental as it adapts existing geometric and numerical techniques.
The paper presents a numerical algorithm for solving singular linear-quadratic optimal control problems by extending the presymplectic constraint algorithm with singular value decomposition, demonstrating stable behavior on large-scale examples with high index.
A numerical algorithm to obtain the consistent conditions satisfied by singular arcs for singular linear-quadratic optimal control problems is presented. The algorithm is based on the presymplectic constraint algorithm (PCA) by Gotay-Nester \cite{Go78,Vo99} that allows to solve presymplectic hamiltonian systems and that provides a geometrical framework to the Dirac-Bergmann theory of constraints for singular Lagrangian systems \cite{Di49}. The numerical implementation of the algorithm is based on the singular value decomposition that, on each step allows to construct a semi-explicit system. Several examples and experiments are discussed, among them a family of arbitrary large singular LQ systems with index 3 and a family of examples of arbitrary large index, all of them exhibiting stable behaviour.