An iterative method for the bisymmetric solutions of the consistent matrix equations
In this paper, an iterative method is presented for finding the bisymmetric solutions of a pair of consistent matrix equations A1XB1=C1, A2XB2=C2, by which a bisymmetric solution can be obtained in finite iteration steps in the absence of round-off errors. Moreover, the solution with least Frobenius norm can be obtained by choosing a special kind of initial matrix. In the solution set of the matrix equations, the optimal approximation bisymmetric solution to a given matrix can also be derived by this iterative method. The efficiency of the proposed algorithm is shown by some numerical examples.