This paper focuses on a linearized fully discrete projection scheme for time-dependent magnetohydrodynamics equations in three-dimensional bounded domain. It is shown that the proposed projection scheme allows for a discrete energy inequality and is unconditionally stable. In addition, we present a rigorous analysis for the rates of convergence.