English

An iterative method to solve Lyapunov equations

Numerical Analysis 2023-09-26 v2 Numerical Analysis

Abstract

We present here a new splitting method to solve Lyapunov equations of the type AP+PAT=BBTAP + PA^T=-BB^T in a Kronecker product form. Although that resulting matrix is of order n2n^2, each iteration of the method demands only two operations with the matrix AA: a multiplication of the form (AσI)B^(A-\sigma I) \hat{B} and a inversion of the form (AσI)1B^(A-\sigma I)^{-1}\hat{B}. We see that for some choice of a parameter the iteration matrix is such that all their eigenvalues are in absolute value less than 1, which means that it should converge without depending on the starting vector. Nevertheless we present a theorem that enables us how to get a good starting vector for the method.

Cite

@article{arxiv.2309.12143,
  title  = {An iterative method to solve Lyapunov equations},
  author = {Licio Hernanes Bezerra and Felipe Wisniewski},
  journal= {arXiv preprint arXiv:2309.12143},
  year   = {2023}
}

Comments

7 pages

R2 v1 2026-06-28T12:28:26.274Z