We expose a theoretical analysis of a least-squares mixed finite elements method for the linear elasticity problem in two-and three-dimensional domains. The coerciveness of the weak form of the problem is proved. It is shown that the finite element approximation yields a symmetric positive definite linear system with condition number $O(h^{-2})$. The error estimate is obtained.