An efficient real representation method for least squares problem of the quaternion constrained matrix equation AXB + CY D = E