This article discusses a numerical gradient-based method for solving a generalized matrix equation. The iterative method mentioned includes two positive parameters, for which a range is determined to ensure the convergence of the introduced method. It has been shown that the optimal parameters of this method satisfy a constrained optimization problem. Then, specific solutions of this equation, such as the reflexive solution, are examined. Furthermore, to increase the rate of convergence of the proposed method, the idea of momentum is utilized, and a range for the momentum parameter is obtained to ensure the convergence. Finally, the efficiency of this method is investigated through numerical simulations. Additionally, in numerical results section, applications of mentioned matrix equation in image restoration and anti-linear systems are addressed.