170 research outputs found
A Subspace Shift Technique for Nonsymmetric Algebraic Riccati Equations
The worst situation in computing the minimal nonnegative solution of a
nonsymmetric algebraic Riccati equation associated with an M-matrix occurs when
the corresponding linearizing matrix has two very small eigenvalues, one with
positive and one with negative real part. When both these eigenvalues are
exactly zero, the problem is called critical or null recurrent. While in this
case the problem is ill-conditioned and the convergence of the algorithms based
on matrix iterations is slow, there exist some techniques to remove the
singularity and transform the problem to a well-behaved one. Ill-conditioning
and slow convergence appear also in close-to-critical problems, but when none
of the eigenvalues is exactly zero the techniques used for the critical case
cannot be applied.
In this paper, we introduce a new method to accelerate the convergence
properties of the iterations also in close-to-critical cases, by working on the
invariant subspace associated with the problematic eigenvalues as a whole. We
present a theoretical analysis and several numerical experiments which confirm
the efficiency of the new method
A numerical comparison of solvers for large-scale, continuous-time algebraic Riccati equations and LQR problems
In this paper, we discuss numerical methods for solving large-scale
continuous-time algebraic Riccati equations. These methods have been the focus
of intensive research in recent years, and significant progress has been made
in both the theoretical understanding and efficient implementation of various
competing algorithms. There are several goals of this manuscript: first, to
gather in one place an overview of different approaches for solving large-scale
Riccati equations, and to point to the recent advances in each of them. Second,
to analyze and compare the main computational ingredients of these algorithms,
to detect their strong points and their potential bottlenecks. And finally, to
compare the effective implementations of all methods on a set of relevant
benchmark examples, giving an indication of their relative performance
Nonsymmetric algebraic Riccati equations associated with an M-matrix: recent advances and algorithms
We survey on theoretical properties and algorithms concerning the problem of solving a nonsymmetric algebraic Riccati equation, and we report on some known methods and new algorithmic advances. In particular, some results on the number of positive solutions are proved and a careful convergence analysis of Newton\u27s iteration is carried out in the cases of interest where some singularity conditions are encountered. From this analysis we determine initial approximations which still guarantee the quadratic convergence
Iterative and doubling algorithms for Riccati-type matrix equations: a comparative introduction
We review a family of algorithms for Lyapunov- and Riccati-type equations
which are all related to each other by the idea of \emph{doubling}: they
construct the iterate of another naturally-arising fixed-point
iteration via a sort of repeated squaring.
The equations we consider are Stein equations , Lyapunov
equations , discrete-time algebraic Riccati equations
, continuous-time algebraic Riccati equations
, palindromic quadratic matrix equations , and
nonlinear matrix equations . We draw comparisons among these
algorithms, highlight the connections between them and to other algorithms such
as subspace iteration, and discuss open issues in their theory.Comment: Review article for GAMM Mitteilunge
- …