Solving large-scale continuous-time algebraic Riccati equations by doubling

We consider the solution of large-scale algebraic Riccati equations with numerically low-ranked solutions. For the discrete-time case, the structure-preserving doubling algorithm has been adapted, with the iterates for A not explicitly computed but in the recursive form A"k=A"k"-"1^2-D"k^(^1^)S"k^-^1[D"k^(^2^)]^@?, with D"k^(^1^) and D"k^(^2^) being low-ranked and S"k^-^1 being small in dimension. For the continuous-time case, the algebraic Riccati equation will be first treated with the Cayley transform before doubling is applied. With n being the dimension of the algebraic equations, the resulting algorithms are of an efficient O(n) computational complexity per iteration, without the need for any inner iterations, and essentially converge quadratically. Some numerical results will be presented. For instance in Section 5.2, Example 3, of dimension n=20209 with 204 million variables in the solution X, was solved using MATLAB on a MacBook Pro within 45 s to a machine accuracy of O(10^-^1^6).

[1]  Ji-guang Sun,et al.  Condition numbers of algebraic Riccati equations in the Frobenius norm , 2002 .

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

[3]  Peter Benner,et al.  Large-Scale Matrix Equations of Special Type , 2022 .

[4]  Hung-Yuan Fan,et al.  Structure-Preserving Algorithms for Periodic Discrete-Time Algebraic Riccati Equations , 2004 .

[5]  Khalide Jbilou,et al.  ADI preconditioned Krylov methods for large Lyapunov matrix equations , 2010 .

[6]  Peter Benner,et al.  A Semi-Discretized Heat Transfer Model for Optimal Cooling of Steel Profiles , 2005 .

[7]  Alain Bensoussan,et al.  Representation and Control of Infinite Dimensional Systems, 2nd Edition , 2007, Systems and control.

[8]  E. Ostertag.  Linear Matrix Inequalities , 2011 .

[9]  E. Chu,et al.  Solving Large-Scale Discrete-Time Algebraic Riccati Equations by Doubling Tiexiang , 2012 .

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

[11]  Peter Benner,et al.  On the numerical solution of large-scale sparse discrete-time Riccati equations , 2011, Adv. Comput. Math..

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

[13]  K. Jbilou,et al.  Projection methods for large Lyapunov matrix equations , 2006 .

[14]  Khalide Jbilou,et al.  Low rank approximate solutions to large Sylvester matrix equations , 2006, Appl. Math. Comput..

[15]  P. Benner,et al.  Solving large-scale control problems , 2004, IEEE Control Systems.

[16]  E. Chu,et al.  PALINDROMIC EIGENVALUE PROBLEMS: A BRIEF SURVEY , 2010 .

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

[18]  G. Hewer,et al.  The sensitivity of the algebraic and differential riccati equations , 1990 .

[19]  Wen-Wei Lin,et al.  Convergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations , 2006, SIAM J. Matrix Anal. Appl..

[20]  Athanasios C. Antoulas,et al.  Approximation of Large-Scale Dynamical Systems (Advances in Design and Control) (Advances in Design and Control) , 2005 .

[21]  A. Laub,et al.  Benchmarks for the numerical solution of algebraic Riccati equations , 1997 .

[22]  Peter Benner,et al.  Dimension Reduction of Large-Scale Systems , 2005 .

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

[24]  Y. Danieli.  Guide , 2005 .

[25]  Athanasios C. Antoulas,et al.  Approximation of Large-Scale Dynamical Systems , 2005, Advances in Design and Control.

[26]  R. Triggiani,et al.  Control Theory for Partial Differential Equations: Continuous and Approximation Theories , 2000 .

[27]  D. Auslander.  Research needs in the control of mechanical systems , 1987, IEEE Control Systems Magazine.

[28]  V. Mehrmann.  The Autonomous Linear Quadratic Control Problem , 1991 .

[29]  Jacob K. White,et al.  Low-Rank Solution of Lyapunov Equations , 2004, SIAM Rev..

[30]  Xun Yu Zhou,et al.  Linear matrix inequalities, Riccati equations, and indefinite stochastic linear quadratic controls , 2000, IEEE Trans. Autom. Control..

[31]  A. Iserles.  A First Course in the Numerical Analysis of Differential Equations: Stiff equations , 2008 .

[32]  Ji-guang Sun.  Perturbation Theory for Algebraic Riccati Equations , 1998, SIAM J. Matrix Anal. Appl..

[33]  Daniel Boley,et al.  Numerical Methods for Linear Control Systems , 1994 .

[34]  Khalide Jbilou,et al.  Block Krylov Subspace Methods for Large Algebraic Riccati Equations , 2003, Numerical Algorithms.

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

[36]  Alan J. Laub,et al.  Backward error, sensitivity, and refinement of computed solutions of algebraic Riccati equations , 1995, Numer. Linear Algebra Appl..

[37]  Ulrike Baur,et al.  Low rank solution of data‐sparse Sylvester equations , 2008, Numer. Linear Algebra Appl..

[38]  Khalide Jbilou,et al.  An Arnoldi based algorithm for large algebraic Riccati equations , 2006, Appl. Math. Lett..