In the numerical solution of the algebraic Riccati equation A∗X + XA - XBB∗X + C∗C = 0, where A is large, sparse, and stable, and B, C have low rank, projection methods have recently emerged as a possible alternative to the more established Newton-Kleinman iteration. In spite of convincing numerical experiments, a systematic matrix analysis of this class of methods is still lacking. We derive new relations for the approximate solution, the residual, and the error matrices, giving new insights into the role of the matrix A - BB∗X and of its approximations in the numerical procedure. In the context of linear-quadratic regulator problems, we show that the Riccati approximate solution is related to the optimal value of the reduced cost functional, thus completely justifying the projection method from a model order reduction point of view. Finally, the new results provide theoretical ground for recently proposed modifications of projection methods onto rational Krylov subspaces.

Simoncini, V. (2016). Analysis of the rational Krylov subspace projection method for large-scale algebraic riccati equations. SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 37(4), 1655-1674 [10.1137/16M1059382].

Analysis of the rational Krylov subspace projection method for large-scale algebraic riccati equations

SIMONCINI, VALERIA
2016

Abstract

In the numerical solution of the algebraic Riccati equation A∗X + XA - XBB∗X + C∗C = 0, where A is large, sparse, and stable, and B, C have low rank, projection methods have recently emerged as a possible alternative to the more established Newton-Kleinman iteration. In spite of convincing numerical experiments, a systematic matrix analysis of this class of methods is still lacking. We derive new relations for the approximate solution, the residual, and the error matrices, giving new insights into the role of the matrix A - BB∗X and of its approximations in the numerical procedure. In the context of linear-quadratic regulator problems, we show that the Riccati approximate solution is related to the optimal value of the reduced cost functional, thus completely justifying the projection method from a model order reduction point of view. Finally, the new results provide theoretical ground for recently proposed modifications of projection methods onto rational Krylov subspaces.
2016
Simoncini, V. (2016). Analysis of the rational Krylov subspace projection method for large-scale algebraic riccati equations. SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 37(4), 1655-1674 [10.1137/16M1059382].
Simoncini, V.
File in questo prodotto:
File Dimensione Formato  
16m1059382.pdf

accesso aperto

Tipo: Versione (PDF) editoriale
Licenza: Licenza per accesso libero gratuito
Dimensione 731.07 kB
Formato Adobe PDF
731.07 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/586007
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 38
  • ???jsp.display-item.citation.isi??? 33
social impact