Abstract: Aiming at the least-squares solutions of the matrix equation AX=B under the norm inequality constraint, an iterative method is proposed. This method uses the general Lanczos trust region algorithm as the frame method and remedies its defect on solving the matrix equation. Numerical experiments illustrate that the algorithm is effective.