Efficient and reliable algorithms for condition estimation of Lyapunov and Riccati equations

Efficient and reliable condition estimators for continuoustime or discrete-time matrix Lyapunov and Riccati equations are discussed. The basic calculation for obtaining these estimates is the evaluation of the 1-norm of the inverse of a Lyapunov operator. This is performed by evaluating this operator and its dual on suitably selected values, or, equivalently, by solving related Lyapunov equations. Various strategies for exploiting the symmetric structure of the problem are investigated and compared. The results show that the proposed strategies are highly accurate and very fast.