This study presents a computational method for the solution of the fractional optimal control problems subject to fractional systems with equality and inequality constraints. The proposed procedure is based upon Fibonacci wavelets. The fractional derivative is described in the Caputo sense. The Riemann–Liouville operational matrix for Fibonacci wavelets is obtained. Then, we use this operational matrix and the Galerkin method to reduce the given problem into a system of algebraic equations. We discuss the convergence of the algorithm. Several numerical examples are included to observe the validity, effectiveness, and accuracy of the suggested scheme. Moreover, fractional optimal control problems are studied through a bibliometric viewpoint.