Biconjugate residual algorithm for solving general Sylvester-transpose matrix equations

Abstract

The present paper is concerned with the solution of the coupled generalized Sylvester-transpose matrix equations {A1XB1 + C1XD1 + E1XTF1 = M1, A2XB2 + C2XD2 + E2XTF2 = M2, including the well-known Lyapunov and Sylvester matrix equations. Based on a variant of biconjugate residual (BCR) algorithm, we construct and analyze an efficient algorithm to find the (least Frobenius norm) solution of the general Sylvester-transpose matrix equations within a finite number of iterations in the absence of round-off errors. Two numerical examples are given to examine the performance of the constructed algorithm.