The numerical solution of a fractional optimal control problem having a quadratic performance index is proposed and analyzed. The performance index of the fractional optimal control problem is considered as a function of both the state and the control variables. The dynamic constraint is expressed as a fractional differential equation that includes an integer derivative in addition to the fractional derivative. The order of the fractional derivative is taken as less than one and described in the Caputo sense. Based on the shifted Legendre orthonormal polynomials, we employ the operational matrix of fractional derivatives, the Legendre-Gauss quadrature formula and the Lagrange multiplier method for reducing such a problem into a problem consisting of solving a system of algebraic equations. The convergence of the proposed method is analyzed. For confirming the validity and accuracy of the proposed numerical method, a numerical example is presented along with a comparison between our numerical results and those obtained using the Legendre spectral-collocation method.
