In the present paper, a direct numerical technique for solving fractional optimal control problems based on an orthonormal wavelet, is introduced. First we approximate the involved functions by Sine-Cosine wavelet basis; then, an operational matrix is used to transform the given problem into a linear system of algebraic equations, which is easier. In fact operational matrix of Riemann-Liouville fractional integration and derivative of Sine-Cosine wavelet are employed to achieve a linear algebraic equation. The mentioned matrices are derived via hat functions. The solution of transformed system, gives us the solution of original problem. Two numerical examples are also given. Finally, the paper is ended with conclusion.