Premium
An efficient iterative method for solving the matrix equation AXB + CYD = E
Author(s) -
Peng Zhenyun,
Peng Yaxin
Publication year - 2006
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.470
Subject(s) - mathematics , iterative method , matrix norm , matrix (chemical analysis) , norm (philosophy) , symmetric matrix , matrix difference equation , matrix splitting , state transition matrix , mathematical analysis , algebra over a field , pure mathematics , algorithm , differential equation , eigenvalues and eigenvectors , riccati equation , physics , quantum mechanics , materials science , political science , law , composite material
This paper presents an iterative method for solving the matrix equation AXB + CYD = E with real matrices X and Y . By this iterative method, the solvability of the matrix equation can be determined automatically. And when the matrix equation is consistent, then, for any initial matrix pair [ X 0 , Y 0 ], a solution pair can be obtained within finite iteration steps in the absence of round‐off errors, and the least norm solution pair can be obtained by choosing a special kind of initial matrix pair. Furthermore, the optimal approximation solution pair to a given matrix pair [ X̄ , Ȳ ] in a Frobenius norm can be obtained by finding the least norm solution pair of a new matrix equation AX̃B + CỸD = Ẽ , where Ẽ = E − AX̄B − CȲD . The given numerical examples show that the iterative method is efficient. Copyright © 2005 John Wiley & Sons, Ltd.