Abstract:In this paper, a structure-preserving algorithm is developed for the computation of a semi-stabilizing solution of a Generalized Algebraic Riccati Equation (GARE). The semi-stabilizing solution of GAREs has been used to characterize the solvability of the (J, J ')-spectral factorization problem in control theory for general rational matrices which may have poles and zeros on the extended imaginary axis. The main difficulty in solving such a GARE lies in the fact that its associated Hamiltonian/skew-Hamiltonian pencil has eigenvalues on the extended imaginary axis. Consequently, it is not clear which eigenspace of the associated Hamiltonian/skew-Hamiltonian pencil can characterize the desired semi-stabilizing solution. That is, it is not clear which eigenvectors and principal vectors corresponding to the eigenvalues on the extended imaginary axis should be contained in the eigenspace that we wish to compute. Hence, the well-known generalized eigenspace approach for the classical algebraic Riccati equations cannot be employed directly. The proposed algorithm consists of a structure-preserving doubling algorithm (SDA) and a postprocessing procedure to determine the desired eigenvectors and principal vectors corresponding to the purely imaginary and infinite eigenvalues. Under mild assumptions, linear convergence of rate 1/2 for the SDA is proved. Numerical experiments illustrate that the proposed algorithm performs efficiently and reliably.

Large-Scale Discrete-Time Algebraic Riccati Equations — Numerically Low-Ranked Solution , SDA and Error Analysis ∗

Solving Large-Scale Continuous-Time Algebraic Riccati Equations by Doubling

Solving Large-scale Discrete-time Algebraic Riccati Equations by Doubling

Solving Large-Scale Nonsymmetric Algebraic Riccati Equations By Doubling

Decoupled Structure-Preserving Doubling Algorithm with Truncation for Large-Scale Algebraic Riccati Equations

A structure-preserving doubling algorithm for nonsymmetric algebraic Riccati equation

Low-rank alternating direction doubling algorithm for solving large-scale continuous time algebraic Riccati equations

Structure-Preserving Doubling Algorithms That Avoid Breakdowns for Algebraic Riccati-Type Matrix Equations

The Structure-preserving Doubling Algorithm for the Discrete Coupled Riccati Matrix Equations

A decoupled form of the structure-preserving doubling algorithm with low-rank structures

On the generalized structure-preserving doubling algorithm for Generalized Discrete-time Algebraic Riccati equations

A STRUCTURE-PRESERVING ALGORITHM FOR SEMI-STABILIZING SOLUTIONS OF GENERALIZED ALGEBRAIC RICCATI EQUATIONS

On a family of low-rank algorithms for large-scale algebraic Riccati equations

The Intrinsic Toeplitz Structure and Its Applications in Algebraic Riccati Equations

An RADI-type method for stochastic continuous-time algebraic Riccati equations

ISC: an RADI-type Method for Stochastic Continuous-Time Algebraic Riccati Equations.

Using $LDL^{T}$ factorizations in Newton's method for solving general large-scale algebraic Riccati equations

An iterative method for sloving stochastic riccati algebralc equations

Systolic algorithms and architectures for algebraic riccati equation

Comments on "upper Matrix Bounds for the Discrete Algebraic Riccati Matrix Equation" [with Author's Reply]