This paper introduces a new numerical approach to solving a system of fractional differential equations (FDEs) using the Legendre wavelet operational matrix method (LWOMM). We first formulated the operational matrix of fractional derivatives in some special conditions using some notable characteristics of Legendre wavelets and shifted Legendre polynomials. Then, the system of fractional differential equations was transformed into a system of algebraic equations by using these operational matrices. At the end of this paper, several examples are presented to illustrate the effectivity and correctness of the proposed approach. Comparing the methodology with several recognized methods demonstrates that the advantages of the Legendre wavelet operational matrix method are its accuracy and the understandability of the calculations.