The RADI algorithm for solving large-scale algebraic Riccati equations

In recent years, several new approaches for solving the large-scale continuous-time algebraic Riccati equation have appeared in the literature. Amodei and Buchot suggest computing a low-dimensional invariant subspace of the associated Hamiltonian matrix. Simoncini and Lin also target the Hamiltonian matrix, but in a different way: they iterate on the Cayley-transformed matrix with various shifts. Wong and Balakrishnan directly generalize the Lyapunov ADI-method to the Riccati equation. In this paper we introduce another method, inspired by the Cholesky-factored variant of the Lyapunov ADI-method. The advantage of the new algorithm is in its immediate and efficient low-rank formulation, and a simpler implementation compared to the three algorithms mentioned above. We discuss the theoretical properties of the new method, as well as various shift selection strategies. Finally, we show that all of the seemingly different methods listed above in fact produce exactly the same iterates when used with the same parameters: they are algorithmically different descriptions of the same approximation sequence to the Riccati solution.

[1]  Valeria Simoncini,et al.  A New Iterative Method for Solving Large-Scale Lyapunov Matrix Equations , 2007, SIAM J. Sci. Comput..

[2]  Peter Benner,et al.  An improved numerical method for balanced truncation for symmetric second-order systems , 2013 .

[3]  Michael L. Overton,et al.  Large-Scale Optimization of Eigenvalues , 1990, SIAM J. Optim..

[4]  L. Reichel,et al.  A generalized ADI iterative method , 1993 .

[5]  Peter Benner,et al.  On the Solution of Large-Scale Algebraic Riccati Equations by Using Low-Dimensional Invariant Subspaces , 2016 .

[6]  Hung-Yuan Fan,et al.  A structure-preserving doubling algorithm for continuous-time algebraic Riccati equations , 2005 .

[7]  V. Mehrmann,et al.  Defect correction method for the solution of algebraic Riccati equations , 1988 .

[8]  John Sabino,et al.  Solution of Large-Scale Lyapunov Equations via the Block Modified Smith Methods , 2006 .

[9]  Ngai Wong,et al.  Fast Positive-Real Balanced Truncation Via Quadratic Alternating Direction Implicit Iteration , 2007, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[10]  Ngai Wong,et al.  Quadratic alternating direction implicit iteration for the fast solution of algebraic Riccati equations , 2005, 2005 International Symposium on Intelligent Signal Processing and Communication Systems.

[11]  Peter Benner,et al.  Numerical solution of large‐scale Lyapunov equations, Riccati equations, and linear‐quadratic optimal control problems , 2008, Numer. Linear Algebra Appl..

[12]  Gene H. Golub,et al.  Matrix computations , 1983 .

[13]  Peter Benner,et al.  On the ADI method for Sylvester equations , 2009, J. Comput. Appl. Math..

[14]  Peter Benner,et al.  Efficient handling of complex shift parameters in the low-rank Cholesky factor ADI method , 2012, Numerical Algorithms.

[15]  M. Heyouni,et al.  AN EXTENDED BLOCK ARNOLDI ALGORITHM FOR LARGE-SCALE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATI ON ∗ , 2008 .

[16]  Valeria Simoncini,et al.  Computational Methods for Linear Matrix Equations , 2016, SIAM Rev..

[17]  Peter Benner,et al.  A Reformulated Low‐Rank ADI Iteration with Explicit Residual Factors , 2013 .

[18]  Bruno Iannazzo,et al.  Numerical Solution of Algebraic Riccati Equations , 2012, Fundamentals of algorithms.

[19]  E. Wachspress.  Iterative solution of the Lyapunov matrix equation , 1988 .

[20]  Eugene L. Wachspress,et al.  The ADI Model Problem , 2013 .

[21]  Jean-Marie Buchot,et al.  An invariant subspace method for large-scale algebraic Riccati equation , 2010 .

[22]  Mark R. Opmeer,et al.  Analysis of an Iteration Method for the Algebraic Riccati Equation , 2016, SIAM J. Matrix Anal. Appl..

[23]  Jens Saak,et al.  Efficient Numerical Solution of Large Scale Algebraic Matrix Equations in PDE Control and Model Order Reduction , 2009 .

[24]  Jan G. Korvink,et al.  Oberwolfach Benchmark Collection , 2005 .

[25]  Axel Ruhe.  The rational Krylov algorithm for nonsymmetric eigenvalue problems. III: Complex shifts for real matrices , 1994 .

[26]  A. Laub.  A schur method for solving algebraic Riccati equations , 1978, 1978 IEEE Conference on Decision and Control including the 17th Symposium on Adaptive Processes.

[27]  Thomas Wolf,et al.  The ADI iteration for Lyapunov equations implicitly performs pseudo-optimal model order reduction , 2013, Int. J. Control.

[28]  Valeria Simoncini,et al.  On two numerical methods for the solution of large-scale algebraic Riccati equations , 2014 .

[29]  Peter Benner,et al.  Self-Generating and Efficient Shift Parameters in ADI Methods for Large Lyapunov and Sylvester Equations , 2014 .

[30]  Peter Benner,et al.  An inexact low-rank Newton-ADI method for large-scale algebraic Riccati equations , 2016 .

[31]  Peter Benner,et al.  Low-rank Newton-ADI methods for large nonsymmetric algebraic Riccati equations , 2016, J. Frankl. Inst..

[32]  P. Benner,et al.  Numerical solution of large and sparse continuous time algebraic matrix Riccati and Lyapunov equations: a state of the art survey , 2013 .

[33]  Jacob K. White,et al.  Low Rank Solution of Lyapunov Equations , 2002, SIAM J. Matrix Anal. Appl..

[34]  Valeria Simoncini,et al.  Adaptive rational Krylov subspaces for large-scale dynamical systems , 2011, Syst. Control. Lett..

[35]  Peter Benner,et al.  Theory and Numerical Solution of Differential and Algebraic Riccati Equations , 2015 .

[36]  P. Benner,et al.  Efficient solution of large scale Lyapunov and Riccati equations arising in model order reduction problems , 2008 .

[37]  Leiba Rodman,et al.  Algebraic Riccati equations , 1995 .

[38]  Valeria Simoncini,et al.  A new subspace iteration method for the algebraic Riccati equation , 2013, Numer. Linear Algebra Appl..