By applying the canonical correlation decomposition (CCD) of matrix pairs,we obtain a general expression of the least-squares solutions of the matrix equation ATXA = D under the restriction that the solution matrix X Rn×n is bisymmetric, where A Rn×m and D Rm×m are given matrices.