A POD projection method for large-scale algebraic Riccati equations

The solution of large-scale matrix algebraic Riccati equations is important for instance in control design and model reduction and remains an active area of research. We consider a class of matrix algebraic Riccati equations (AREs) arising from a linear system along with a weighted inner product. T...

Full description

Bibliographic Details
Main Authors: Singler, John, Kramer, Boris
Other Authors: Massachusetts Institute of Technology. Department of Aeronautics and Astronautics
Format: Article
Language:en_US
Published: American Institute of Mathematical Sciences 2017
Online Access:http://hdl.handle.net/1721.1/108516
Description
Summary:The solution of large-scale matrix algebraic Riccati equations is important for instance in control design and model reduction and remains an active area of research. We consider a class of matrix algebraic Riccati equations (AREs) arising from a linear system along with a weighted inner product. This problem class often arises from a spatial discretization of a partial differential equation system. We propose a projection method to obtain low rank solutions of AREs based on simulations of linear systems coupled with proper orthogonal decomposition. The method can take advantage of existing (black box) simulation code to generate the projection matrices. We also develop new weighted norm residual computations and error bounds. We present numerical results demonstrating that the proposed approach can produce highly accurate approximate solutions. We also briefly discuss making the proposed approach completely data-based so that one can use existing simulation codes without accessing system matrices