Based on chebyshev approximation, an integer-order spectral operator differential matrix is derive, and a recurrence scheme for computing fractional spectral differential operator matrix is also obtained by using the three terms recurrence of chebyshev polynomials and the derivatives of chebyshev polynomials. Numerical examples are given to demonstrate the accuracy and effectiveness of the scheme.