An Iterative Algorithm for the Least Squares Generalized Reflexive Solutions of the Matrix Equations 𝐴𝑋𝐵=𝐸,𝐶𝑋𝐷=𝐹

Abstract

The generalized coupled Sylvester systems play a fundamental role in wide applications in several areas, such as stability theory, control theory, perturbation analysis, and some other fields of pure and applied mathematics. The iterative method is an important way to solve the generalized coupled Sylvester systems. In this paper, an iterative algorithm is constructed to solve the minimum Frobenius norm residual problem: ‖𝐴𝑋𝐵𝐶𝑋𝐷−𝐸𝐹‖= min over generalized reflexive matrix 𝑋. For any initial generalized reflexive matrix 𝑋1, by the iterative algorithm, the generalized reflexive solution 𝑋∗ can be obtained within finite iterative steps in the absence of round-off errors, and the unique least-norm generalized reflexive solution 𝑋∗ can also be derived when an appropriate initial iterative matrix is chosen. Furthermore, the unique optimal approximate solution 𝑋 to a given matrix 𝑋0 in Frobenius norm can be derived by finding the least-norm generalized reflexive solution 𝑋∗ of a new corresponding minimum Frobenius norm residual problem: min‖𝐴𝑋𝐵𝐶𝑋𝐷−𝐸𝐹‖ with 𝐸=𝐸−𝐴𝑋0𝐵, 𝐹=𝐹−𝐶𝑋0𝐷. Finally, several numerical examples are given to illustrate that our iterative algorithm is effective.