Order reduction methods for solving large-scale differential matrix Riccati equations

Abstract

We consider the numerical solution of large-scale symmetric differential matrix Riccati equations. Under certain hypotheses on the data, reduced order methods have recently arisen as a promising class of solution strategies, by forming low-rank approximations to the sought after solution at selected timesteps. We show that great computational and memory savings are obtained by a reduction process onto rational Krylov subspaces, as opposed to current approaches. By specifically addressing the solution of the reduced differential equation and reliable stopping criteria, we are able to obtain accurate final approximations at low memory and computational requirements. This is obtained by employing a two-phase strategy that separately enhances the accuracy of the algebraic approximation and the time integration. The new method allows us to numerically solve much larger problems than in the current literature. Numerical experiments on benchmark problems illustrate the effectiveness of the procedure with respect to existing solvers.

Tridiagonalization of systems of coupled linear differential equations with variable coefficients by a Lanczos-like method

“…Systems of coupled linear ODEs: background. Systems of coupled linear differential equations with variable coefficients naturally arise in a variety of contexts in mathematics [22,15,6,4,3] and beyond, from engineering to quantum physics [14,1,11,13,2,25,16,26]. Yet, determining the solutions of such systems both formally and numerically remains surprisingly difficult, their widespread applicability making these difficulties only more acute.…”

Lanczos-like algorithm for the time-ordered exponential: The $\ast$-inverse problem

“…Further applications are found via differential Lyapunov and Riccati matrix equations, which frequently appear in control theory, filter design, and model reduction problems [41,30,9,5,3]. Indeed, the solutions of such differential equations involve time-ordered exponentials [29,1,23,28].…”