In this paper we present a new numerical method for solving fractional differential equations (FDEs) based on Genocchi polynomials operational matrix through collocation method. The operational matrix of fractional integration in Riemann-Liouville sense is derived. The upper bound for the error of the operational matrix of fractional integration is also shown. The properties of Genocchi polynomials are utilized to reduce the given problems to a system of algebraic equations. Illustrative examples are finally given to show the simplicity, accuracy and applicability of the method. [ABSTRACT FROM AUTHOR]