A quadratic numerical scheme for fractional optimal control problems

This paper presents a quadratic numerical scheme for a class of fractional optimal control problems (FOCPs). The fractional derivative is described in the Caputo sense. The performance index of a FOCP is considered as a function of both the state and the control variables, and the dynamic constraints are expressed by a set of fractional differential equation. The calculus of variation, the Lagrange multiplier, and the formula for fractional integration by parts are used to obtain Euler-Lagrange equations for the FOCP. The formulation presented and the resulting equations are very similar to those that appear in the classical optimal control theory. Thus, the present formulation essentially extends the classical control theory to fractional dynamic systems. The formulation is used to derive the control equations for a quadratic linear fractional control problem. For a linear system, this method results into a set of linear simultaneous equations, which can be solved using a direct or an iterative scheme. Numerical results for a FOCP are presented to demonstrate the feasibility of the method. It is shown that the solutions converge as the number of grid points increases, and the solutions approach to classical solutions as the order of the fractional derivatives approach to 1. The formulation presented is simple and can be extended to other FOCPs.