Variational problems of variable fractional order involving arbitrary kernels

: The aim of this work is to study several problems of the calculus of variations, where the dynamics of the state function is given by a generalized fractional derivative. This derivative combines two well-known concepts: fractional derivative with respect to another function and fractional derivative of variable order. We present the Euler–Lagrange equation, which is a necessary condition that every optimal solution of the problem must satisfy. Other problems are also studied: with integral and holonomic constraints, with higher order derivatives, and the Herglotz variational problem


Introduction
The concept of fractional calculus, or arbitrary order calculus, is an extension of the standard calculus, where derivatives and integrals of non-integer order are used (see e.g., [20]). This theory was originated from a question formulated in an exchange of correspondence between Leibniz and l'Hopital, where the interpretation of a derivative of order 1/2 was questioned. Therefore, we can say that the birth of fractional calculus was simultaneous with ordinary calculus, although the first one only had a great development in the last decades. During this period, many famous mathematicians devoted some time to the study of fractional calculus, such as Lagrange, Laplace, Lacroix, Fourier, Abel, Liouville, Riemann, and Grünwald. In the end of the XX century, it was observed that the use of fractional calculus makes it possible to express natural phenomena more precisely when compared to ordinary calculus and, therefore, it can be useful when applied to real world systems. For example, applications in physics [17,23], chemistry [3,26], engineering [13,14,21], biology [18,32], economics [36], and control theory [24,29,37,38], have been found.
There are several definitions for fractional derivatives, although the most common are the Riemann-Liouville and the Caputo ones. However, due to the high number of different concepts, we find several works studying similar problems. One way to overcome this issue is to consider more general definitions with respect to fractional operators. In this work we intend to combine two types of existing generalizations, the fractional derivative with respect to another function [4,27] and fractional derivatives of variable order [16,30,31].
One of the areas where fractional calculus has been applied is in the calculus of variations. The classic problem of the calculus of variations is to find the minimum or maximum value of functionals, usually in the form F(u) := b a L(t, u(t), u (t)) dt, possibly subject to some boundary conditions u(a) = U a , u(b) = U b for some fixed U a , U b ∈ R. In the fractional calculus of variations, this first order derivative u (t) is replace by some kind of fractional derivative D γ u(t). With Riewe's pioneering work in 1996 [28], where he formulated the problem of calculus of variations and obtained the respective Euler-Lagrange equation, numerous works have emerged in this area since then. For example, in [8], the authors considered the isoperimetric problem dealing with the left and right Riemann-Liouville fractional derivatives. In [12], some variational problems were formulated, with dependence on a term, and taking its limit, we obtain the total derivative at the classical level. In the book [22] and in the paper [25], several fractional calculus of variations problems were studied in a general form, where the kernel of the fractional operators is an arbitrary function, for the Riemann-Liouville and Caputo fractional derivatives. Again, due to the large number of definitions for fractional derivatives, we find numerous works in the area of the fractional calculus of variations for different derivatives, but studying similar problems. The aim of this work is to unify some previous works, when considering this new generalized fractional derivative. With this, we generalize some previous works on fractional calculus of variations. In fact, if g(t) = t, then we obtain the usual variable-order fractional operators and such variational problems have been studied extensively e.g., [33,35]. If we fix the order, that is, γ n (·, ·) = γ ∈ R + , then the problem was considered in [5]. In addition, if g(t) = t, then this situation was studied in [1,10] for the Riemann-Liouville fractional derivative and in [2,9,11] for the Caputo fractional derivative. If g(t) = ln t or g(t) = t σ (σ > 0), then the respective variational problems were considered in [6,7,15,19]. Thus, with this paper, we intend to generalize these previous works, and for other choices of the fractional order γ n (·, ·) or the kernel g(·), new results can be obtained.
We start by fixing some notation. For what follows, n is a positive integer, γ n : [a, b] 2 → (n − 1, n) is a function, and u, g : [a, b] → R are two functions with g ∈ C n [a, b] and g (t) > 0, for all t ∈ [a, b]. Definition 1. The generalized variable-order left and right Riemann-Liouville fractional integrals of u, with respect to g and with order γ n , are defined as respectively.
For what concerns the derivatives, two different types are presented.
Definition 2. The generalized variable-order left and right Riemann-Liouville fractional derivatives of u, with respect to g and with order γ n , are defined as respectively.
Definition 3. The generalized variable-order left and right Caputo fractional derivatives of u, with respect to g and with order γ n , are defined as respectively.
We remark that, when g(t) = t, the previous definitions reduce to the classical variable-order fractional operators.
Proof. First observe that Thus, With the change of variable τ = g(s)−g(a) g(t)−g(a) and recalling the definition of the Beta function B(·, ·), we get proving the desired formula.
In an analogous way, we have the following: Lemma 2. Suppose that the fractional order γ n is of form γ n (t, s) = γ n (s), where γ n : [a, b] → (n − 1, n) is a function. Then, for the function u(t) = (g(b) − g(t)) β , with β > n − 1, The paper is structured as follows: in Section 2 we present the integration by parts formulae, dealing with the previous presented fractional derivatives. These formulas will be crucial for the rest of the paper. The main result is given in Section 3, where we prove the fractional Euler-Lagrange equation, which is an important formula to determine if a given curve is a minimizer or a maximizer of a functional. Then, we extended this result by considering additional constraints in the formulation of the problem (Section 4) or in presence of higher order fractional derivatives (Section 5). The Herglotz problem will be considered in Section 6. We end with a conclusion section.

A fractional integration by parts formula
As a first result, we present two integration by parts formulae for the two Caputo fractional derivatives (left and right). These formulae are important in the follow-up of the work, and will be used in the proofs of the results to be presented.
, then the following fractional integration by parts formulae hold: Proof. Changing the order of integration, we obtain the following double integral: (2.1) If we integrate by parts, (2.1) becomes Integrating again by parts, and performing one last time integration by parts, we get proving the first formula. The second one is obtained using similar techniques.

The fractional Euler-Lagrange equation
The purpose of this section is to present the basic problem of the fractional calculus of variations, involving the fractional derivatives presented in Definition 3. To find the candidates for minimizing or maximizing a given functional, we will have to solve a fractional differential equation, known as the Euler-Lagrange equation (see Eq (3.3)).
We will consider the following fractional calculus of variation problem: minimize or maximize the functional The boundary conditions may be imposed on the problem and, for abbreviation, we introduce the operator [·] defined by functional (3.1) reduces to the one studied in [5]. If g(t) = t, that is, we are in presence of the usual variable order fractional operators, then the variational problem was already considered in [33][34][35].
Remark 3. We say that u ∈ Ω is a local minimizer of F is there exists > 0 such that, whenever u ∈ Ω with u − u < , then F(u ) ≤ F(u). If F(u ) ≥ F(u), then we say that u is a local maximizer of F. In such cases, we say that u is a local extremizer of F.
Theorem 2. Let u ∈ Ω be a local extremizer of F as in (3.1). If the maps , then the following fractional Euler-Lagrange equation is satisfied: If u(a) my take any value, then the following fractional transversality condition Proof. Defining function f ( ) := F(u (t) + δ(t)) in a neighbourhood of zero, then f (0) = 0, where δ ∈ Ω is a perturbing curve. If the boundary conditions (3.2) are imposed on the problem, then δ(a) and δ(b) must be both zero so that the curve u (t) + δ(t) is an admissible variation for the problem.
Integrating by parts (Theorem 1), we prove that If, in the set of admissible functions, the boundary conditions (3.2) and since δ may take any value in (a, b), we conclude that for all t ∈ [a, b], proving (3.3). Otherwise, δ is also arbitrary at t = a and t = b. Replacing (3.3) into (3.5), we have and depending if u(a) or u(b) is arbitrary, we deduce the two transversality conditions (3.4).
Observe that, although the functional only depends on the Caputo fractional derivative, the Euler-Lagrange equation (3.3) also involves the Riemann-Liouville fractional derivative. So, this equation deals with four types of fractional derivatives: the left and right Caputo fractional derivatives, and the left and right Riemann-Liouville fractional derivatives. Therefore, in many situations, it is not possible to determine the exact solution of this equation and numerical methods are usually used to determine an approximation of the solution. Such fractional differential equations are useful to check if a given function may or not be a solution of the variational problem. In some particular situations, using some properties of the fractional operators, we may solve the Euler-Lagrange equation and thus produce the optimal solutions. When such a situation is not possible, then using appropriate numerical methods (for example, discretize the equation and then solve a finite dimensional system), an approximation of the solution is obtained. Then, using some sufficient conditions of optimality (e.g., convexity assumptions) we can prove that the obtained solution is indeed a minimizer or maximizer of the functional.

Variational problems under additional constraints
Suppose now that, in the formulation of the variational problem, an integral constraint is imposed on the set of admissible functions (what is called in the literature as an isoperimetric problem). For simplicity of the computations, we will assume from now on that the boundary conditions (3.2) are imposed when formulating the problem (if not, transversality conditions similar to Eq (3.4) are derived). The fractional isoperimetric problem is formulated in the following way: minimize or maximize functional F (as in (3.1)), subject to the boundary conditions (3.2) and to the integral constraint where M : [a, b] × R 3 → R is a C 1 function and Υ ∈ R a fixed number.
In our next problem we add a holonomic constraint, that is, an equation that involves the spatial coordinates of the system and time as well. It is described in the following way. Let The goal is to minimize or maximize the functional where L H : [a, b] × R 6 → R is a function of class C 1 , subject to the boundary conditions (4.4) and to the holonomic constrain then there exists a continuous function λ : Proof. Condition (4.6) is obviously meet for i = 2, if we define .
The case i = 1 is proven in the following way. The variation curve of u is given by u (t) + δ(t), where δ ∈ Ω H and δ(a) = δ(b) = (0, 0). Since any variation must be admissible for the problem, condition (4.5) must be verified for this curve and so the equation Applying Theorem 1, and since δ( Observing that proving the case i = 1 in Eq (4.6).
The necessary condition that every extremizer of this problem must satisfy is given in the next result.
Theorem 5. If u ∈ Ω n is a local minimizer or maximizer of F n (5.1), subject to the boundary conditions and if, for i = 1, . . . , n, the maps Proof. A variation of the optimal curve will be given by u (t) + δ(t), where δ ∈ Ω n and δ (i) (a) = δ (i) (b) = 0, for each i = 0, . . . , n − 1, so that the variation curve satisfies the boundary conditions. Since its first variation must vanish, we obtain From Eq (5.3), the desired result (5.2) follows.
Remark 5. Observe that, if n = 1, Theorem 5 reduces to Theorem 2. Also, additional constraints like the ones presented in Section 5 could be added and similar results as those ones are derived.

The Herglotz variational problem
The Herglotz variational problem is an extension of the previous problems. Instead of finding the extremals for the functional (3.1), we are interested in finding a pair (u , z ) for which function z(·) attains its maximum or minimum value, where functions u and z are related by the ODE where L z : [a, b] × R 4 → R is a function of class C 1 , u ∈ Ω and z ∈ C 1 [a, b]. This problem formulation is an extension of the one presented in Section 3. In fact, if L z does not depend on z, then integrating both sides of Eq (6.1), we get that

Let
[u, z](t) := (t, u(t), C D γ 1 a+ u(t), C D γ 1 b− u(t), z(t)). Theorem 6. Let (u , z ) ∈ Ω × C 1 [a, b] be a solution of problem (6.1). Define function λ : If the maps Proof. We begin by remarking that function z not only depends on time t, but also on the state function u and so we will write z(t, u) instead of z(t) when we need to emphasize this dependence. A variation of the curve u will be still denoted by u (t) + δ(t) (δ ∈ Ω with δ(a) = δ(b) = 0) and the associate variation curve of z is given by The first derivative of Z is then given by Solving this ODE, we prove that Using the fractional integration by parts formulae, and since Z(a) = 0 (z(a) is fixed) and Z(b) = 0 (z(b) attains its extremum), we get that b a λ(t) By the arbitrariness of function δ, we obtain the desired formula.
The previous theorem can be generalized for functions of several independent variables. We denote them by t ∈ [a, b] (time coordinate) and s = (s 1 , . . . , s n ) ∈ S (spatial coordinates), where . . , C D γ 1 a n + u(t)) and , where C D γ 1 a+ u and C D γ 1 b− u are to be understood as the left and right partial fractional derivatives of u with respect to variable t, respectively, and for i = 1, . . . , n, C D γ 1 a i + u and C D γ 1 b i − u are to be understood as the left and right partial fractional derivatives of u with respect to variable s i , respectively.
The new problem is formulated in the following way: find a pair (u , z ) for which z (b) is maximum or minimum value, where u and z are related by the system .

Conclusions
In this paper we investigated several fundamental problems of the calculus of variations, involving a fractional derivative of variable order, and with the kernel depending on an arbitrary function g. More specifically, the functional to minimize or maximize depends on time, the state function, and the left and right Caputo fractional derivatives. We have considered the fixed and free endpoint problems, as well as with additional constraints. Then the problem was generalized, first by considering fractional derivatives of any order and then the generalized Herglotz problem. Since our fractional derivative depends on an arbitrary kernel g(·) and the fractional order is not constant, we obtain numerous works already known in the fractional calculus of variations as particular cases of ours. Also, new ones can be produced by the arbitrariness of those functions. We believe that this is a path of research to be followed, to avoid the multiplication of works dealing with similar problems.
A question that deserves study is how to solve the fractional differential equations presented in this work. As is recognized, in most cases there is no method for analytically solving these equations and so numerical methods are used to find approximations to the optimal solution. For this type of fractional derivative, there is still no numerical method developed and this topic will be studied in a future work.