This paper focuses on the least squares symmetric problem of matrix equation AXB + CXD = E. The explicit expressions for least squares symmetric solution with the least norm of matrix equation AXB + CXD = E are derived. Numerical algorithms and numerical examples show the feasibility of our methods