Authors’ abstract: By extending the idea of the conjugate gradient (CG) method, we construct an iterative method to solve the general coupled matrix equations
(including the generalized (coupled) Lyapunov and Sylvester matrix equations as special cases) over generalized bisymmetric matrix group . By using the iterative method, the solvability of the general coupled matrix equations over the generalized bisymmetric matrix group can be determined in the absence of roundoff errors. When the general coupled matrix equations are consistent over generalized bisymmetric matrices, a generalized bisymmetric solution group can be obtained within finite iteration steps in the absence of roundoff errors.
The least Frobenius norm generalized bisymmetric solution group of the general coupled matrix equations can be derived when an appropriate initial iterative matrix group is chosen. In addition, the optimal approximation generalized bisymmetric solution group to a given matrix group in Frobenius norm can be obtained by finding the least Frobenius norm generalized bisymmetric solution group of new general coupled matrix equations. The numerical results indicate that the iterative method works quite well in practice.