This paper presents a new numerical method for solving fractional differential equations based on Bernoulli polynomials. Bernoulli operation matrix of fractional derivative is deduced, and fractional differential equations are simplified to algebraic equations by combining Tau and collocation methods. The validity and applicability of this method are illustrated by examples.