Next Article in Journal
Inverse Optimization Method for Safety Resource Allocation and Inferring Cost Coefficient Based on a Benchmark
Next Article in Special Issue
New Stability Results for Abstract Fractional Differential Equations with Delay and Non-Instantaneous Impulses
Previous Article in Journal
Oversampling Application of Identifying 3D Selective Laser Sintering Yield by Hybrid Mathematical Classification Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Euler–Lagrange-Type Equations for Functionals Involving Fractional Operators and Antiderivatives

Center for Research and Development in Mathematics and Applications, Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal
Mathematics 2023, 11(14), 3208; https://doi.org/10.3390/math11143208
Submission received: 29 June 2023 / Revised: 19 July 2023 / Accepted: 19 July 2023 / Published: 21 July 2023
(This article belongs to the Special Issue Recent Research on Fractional Calculus: Theory and Applications)

Abstract

:
The goal of this paper is to present the necessary and sufficient conditions that every extremizer of a given class of functionals, defined on the set C 1 [ a , b ] , must satisfy. The Lagrange function depends on a generalized fractional derivative, on a generalized fractional integral, and on an antiderivative involving the previous fractional operators. We begin by obtaining the fractional Euler–Lagrange equation, which is a necessary condition to optimize a given functional. By imposing convexity conditions over the Lagrange function, we prove that it is also a sufficient condition for optimization. After this, we consider variational problems with additional constraints on the set of admissible functions, such as the isoperimetric and the holonomic problems. We end by considering a generalization of the fundamental problem, where the fractional order is not restricted to real values between 0 and 1, but may take any positive real value. We also present some examples to illustrate our results.

1. Introduction

Solving an optimization problem means, as the name implies, looking for the best result according to some pre-established requests. In mathematics optimization problems, we are interested in maximum and minimum problems such as maximum profit, minimum cost, minimum time, or shortest path. An area of mathematics that is very useful in solving optimization problems is calculus of variations (see, e.g., [1,2,3]), which generalizes the theory of maxima and minima of differential calculus for functions whose domain consists of a set of admissible curves. Although the history of the calculus of variations dates back to Ancient Greece, it was not until the 17th century in Western Europe that substantial progress was made. In 1696, Isaac Newton used variational principles to determine the shape of a body moving in the air, so that resistance is minimal. The brothers Johann and Jacob Bernoulli are often considered the inventors of the calculus of variations for having proposed the brachistochrone problem. Lagrange and Euler developed a new theory using variational techniques to determine a differential equation to solve such optimization problems. Such methods are still used nowadays.
Fractional calculus deals with integrals and derivatives of non-integer order [4,5,6]. Unlike integer calculus (where derivatives and integrals are of order n = 1 (first-order derivative and single integral), n = 2 (second-order derivative and double integral), etc.), fractional calculus (with arbitrary fractional order α R + ) has a specific emergence date. The beginnings of fractional calculus refer to an exchange of correspondence between Leibniz and l’Hôpital in 1695, in which Leibniz questioned the generalization from a usual derivative to a derivative with fractional order. Leibniz replied that the possibility of the concept of such a derivative “will lead to a paradox, from which one day useful consequences will be drawn”. Later, numerous important mathematicians tried to develop these concepts, like Laplace, Lacroix, Fourier, Letnikov, etc. It was in the works of Riemann and Liouville that the important definitions of fractional derivatives and integrals appeared, and are still used nowadays. Starting with Cauchy’s formula to evaluate an n-tuple integral, a generalization of integrals for a real order α is given, and then with this new concept, a definition for fractional derivatives was introduced. Later, several definitions were presented, most of them motivated by the works of Riemann and Liouville. More recently, with the concepts of fractional operators with respect to another function, we can generalize some of those operators into a single one and present results valid for a wide class of fractional operators [5,7,8,9,10]. We remark that these fractional operators depend on an arbitrary function g, and for particular choices of such function, we can recover some well-known fractional operators like the Riemann–Liouvile, the Caputo, the Hadamard, or the Erdelyi–Kober fractional operators. With this new approach, we generalize the results presented in [11,12], where the calculus of variation problems was addressed for (classical) fractional operators. In [13], the authors considered functionals defined on a set of non-differentiable functions, but instead of considering the dynamics modeled by a fractional derivative, a quantum derivative was considered. In [14], functionals involving these generalized operators are considered, but only depending on the fractional derivative. Here, we add a fractional integral and an antiderivative of these operators. Applications have been found in different fields. For example, in [4], a number of works are collected with various applications in physics. To mention a few, applications in biophysics, anomalous diffusion, Markovian chains, polymer dynamics, rheological equations, Hamiltonian chaotic systems, and fractal time series were detailed. In [15], a study of a forced mass-spring-damper in a vertical position was carried out, with the derivation of the fractional Euler–Lagrange equation. For further reading on this topic, we suggest [16,17,18,19].
One area where fractional calculus has been applied with success is in the calculus of variations. Typically, the goal is to minimize/maximize a functional, depending on time, the state function, and its first order derivative [20]. With the pioneering work of Riewe [21,22], where the integer-order derivative was replaced by a derivative of order 1 / 2 , these new models proved to be more realistic in modeling real-world phenomena. Since then, this field has attracted the attention of several researchers and, consequently, a lot of work has been carried out (see, e.g., [15,23,24,25,26,27,28]). One of the topics that can be studied is the case where the Lagrange function depends not only on a state function and its fractional derivative, but also on an antiderivative involving all the variables of the Lagrange function [29,30,31]. In this present paper, we extended the work performed with fractional operators by considering a generalized form of fractional integrals and fractional derivatives.
Our objective is to consider different calculus of variation problems, involving the generalized fractional operators and with the presence of an antiderivative. Applying a variational technique and fractional integration by parts formulae, we are able to deduce the necessary and sufficient optimization conditions for the problems to be studied.
The work is organized in the following way: In Section 2, we present a brief description of fractional calculus, needed in for the following work. In Section 3, we obtain the necessary and sufficient conditions for finding the extremum of a given functional. Instead of being described by an ordinary differential equation, it is given by a fractional differential equation, known as the fractional Euler–Lagrange equation. The sufficiency condition is obtained by assuming some convexity condition over the Lagrange function. The case where the lower bound of the integral is greater than the lower bound of the fractional operators is also considered. In Section 4, we consider variational problems subject to additional constraints. Specifically, isoperimetric and holonomic constraints are imposed on the set of admissible functions. We end with Section 5, where the variational problem is defined for arbitrary fractional orders, in contrast to the previous sections, where the order of the derivative belongs to the interval ( 0 , 1 ) .

2. Preliminary Results

First, we introduce some definitions and preliminary facts of fractional calculus theory. Let α , β R + \ N and n be an integer number with α ( n 1 , n ) . Furthermore, consider two functions x , g : [ a , b ] R , with g differentiable and g ( t ) > 0 , for all t [ a , b ] .
Definition 1
([5,7]). The (left) Riemann–Liouville fractional integral of x, order β, with respect to function g, is defined as
I a + β , g x ( t ) = 1 Γ ( β ) a t g ( τ ) ( g ( t ) g ( τ ) ) β 1 x ( τ ) d τ .
If x , g are of C n , the (left) Caputo fractional derivative of the function x is given by the formula
C D a + α , g x ( t ) = 1 Γ ( n α ) a t g ( τ ) ( g ( t ) g ( τ ) ) n α 1 1 g ( τ ) d d τ n x ( τ ) d τ .
For the following, we also will need the two following notions: the right Riemann–Liouville fractional integral of x:
I b β , g x ( t ) = 1 Γ ( β ) t b g ( τ ) ( g ( τ ) g ( t ) ) β 1 x ( τ ) d τ ,
and the right Riemann–Liouville fractional derivative of x:
D b α , g x ( t ) = 1 Γ ( n α ) 1 g ( t ) d d t n t b g ( τ ) ( g ( τ ) g ( t ) ) n α 1 x ( τ ) d τ .
Fractional integration by parts formulae are an important tool to obtain the variational necessary conditions for optimizing functionals. Here, we recall two of them for fractional integrals and for fractional derivatives.
Theorem 1
([32]). Given x , y C [ a , b ] , we have that, for all β > 0 ,
a b x ( t ) · I a + β , g y ( t ) d t = a b I b β , g x ( t ) g ( t ) · y ( t ) g ( t ) d t .
Theorem 2
([7]). Given x , y C n [ a , b ] , we have that, for all α > 0 ,
a b x ( t ) · C D a + α , g y ( t ) d t = a b D b α , g x ( t ) g ( t ) · y ( t ) g ( t ) d t + k = 0 n 1 1 g ( t ) d d t k I b n α , g x ( t ) g ( t ) · 1 g ( t ) d d t n k 1 y ( t ) t = a t = b .
For the particular case when the fractional order α is between 0 and 1, the integration by parts formula given in Theorem 2 reads as
a b x ( t ) · C D a + α , g y ( t ) d t = a b D b α , g x ( t ) g ( t ) · y ( t ) g ( t ) d t + I b 1 α , g x ( t ) g ( t ) · y ( t ) t = a t = b .

3. Fractional Calculus of Variations

Consider the following functional
F ( x ) = a b L ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) d t , x C 1 [ a , b ] ,
with the antiderivative
z ( t ) = a t Z ( τ , x ( τ ) , C D a + α , g x ( τ ) , I a + β , g x ( τ ) ) d τ ,
where α ( 0 , 1 ) , β > 0 , L : [ a , b ] × R 4 R and Z : [ a , b ] × R 3 R are two continuously differentiable functions. To simplify the notation, by [ x ] L and [ x ] Z we mean
[ x ] L ( t ) = ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) a n d [ x ] Z ( t ) = ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) ) ,
for t [ a , b ] . The goal is to find the curves x ¯ for which functional (1) attains an extremum value, possibly under some boundary conditions. For that purpose, consider an arbitrary curve v C 1 [ a , b ] , ϵ an arbitrary real number, and consider the admissible variation [ a , b ] t x ¯ ( t ) + ϵ v ( t ) . We remark that, in the problem’s formulation, if boundary conditions are imposed at t = a or t = b , then function v must satisfy the conditions v ( a ) = 0 or v ( b ) = 0 , so that the curve x ¯ + ϵ v is an admissible variation of the problem. Since functional F attains an extremum value at x ¯ , its first variation must vanish when evaluated on that curve. Starting with
F ( x ¯ + ϵ v ) = a b L ( t , x ¯ ( t ) + ϵ v ( t ) , C D a + α , g x ¯ ( t ) + ϵ C D a + α , g v ( t ) , I a + β , g x ¯ ( t ) + ϵ I a + β , g v ( t ) , z ( t ) ) d t ,
with
z ( t ) = a t Z ( τ , x ¯ ( τ ) + ϵ v ( τ ) , C D a + α , g x ¯ ( τ ) + ϵ C D a + α , g v ( τ ) , I a + β , g x ¯ ( τ ) + ϵ I a + β , g v ( τ ) ) d τ ,
differentiating with respect to ϵ , we obtain that
a b L x [ x ¯ ] L ( t ) v ( t ) + L C D a + α , g x [ x ¯ ] L ( t ) C D a + α , g v ( t ) + L I a + β , g x [ x ¯ ] L ( t ) I a + β , g v ( t ) + L z [ x ¯ ] L ( t ) × a t Z x [ x ¯ ] Z ( τ ) v ( τ ) + Z C D a + α , g x [ x ¯ ] Z ( τ ) C D a + α , g v ( τ ) + Z I a + β , g x [ x ¯ ] Z ( τ ) I a + β , g v ( τ ) d τ d t = 0 .
Using fractional integration by parts (cf. Theorems 1 and 2), we get
a b L C D a + α , g x [ x ¯ ] L ( t ) C D a + α , g v ( t ) d t = a b D b α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) · v ( t ) g ( t ) d t + I b 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) · v ( t ) t = a t = b
and
a b L I a + β , g x [ x ¯ ] L ( t ) I a + β , g v ( t ) d t = a b I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) · v ( t ) g ( t ) d t .
Now, using standard integration by parts, we obtain
a b L z [ x ¯ ] L ( t ) a t Z x [ x ¯ ] Z ( τ ) v ( τ ) d τ d t = a b d d t t b L z [ x ¯ ] L ( τ ) d τ · a t Z x [ x ¯ ] Z ( τ ) v ( τ ) d τ d t = a b t b L z [ x ¯ ] L ( τ ) d τ · Z x [ x ¯ ] Z ( t ) v ( t ) d t .
Finally, combining standard and fractional integration by parts, we get
a b L z [ x ¯ ] L ( t ) a t Z C D a + α , g x [ x ¯ ] Z ( τ ) C D a + α , g v ( τ ) d τ d t = a b D b α , g t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) · v ( t ) g ( t ) d t + I b 1 α , g t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) · v ( t ) t = a t = b
and
a b L z [ x ¯ ] L ( t ) a t Z I a + β , g x [ x ¯ ] Z ( τ ) I a + β , g v ( τ ) d τ d t = a b I b β , g t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) · v ( t ) g ( t ) d t .
Therefore, we conclude from the previous relations that
a b [ L x [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) ] · v ( t ) d t + I b 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) · v ( t ) t = a t = b = 0 .
From the last relation, and from the arbitrariness of function v, we conclude the following:
Theorem 3.
Assume that x ¯ minimizes or maximizes functional F given by (1). Then, for all t [ a , b ] , the equation
L x [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0
holds. Moreover, if x ( a ) or x ( b ) is arbitrary, then the natural boundary condition
I b 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0
holds at t = a or t = b , respectively.
The fractional differential Equation (2) is usually referred to in the literature as the fractional Euler–Lagrange equation associated with the variational problem. We remark that, when g ( t ) = t , we obtain the result proven in [29].
Remark 1.
If x ¯ = ( x ¯ 1 , , x ¯ n ) is an extremizer for functional
F ( x ) = a b L ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) d t ,
with the antiderivative
z ( t ) = a t Z ( τ , x ( τ ) , C D a + α , g x ( τ ) , I a + β , g x ( τ ) ) d τ ,
where x C 1 [ a , b ] × × C 1 [ a , b ] , C D a + α , g x ( t ) = C D a + α , g x 1 ( t ) , , C D a + α , g x n ( t ) , I a + β , g x ( t ) = I a + β , g x 1 ( t ) , , I a + β , g x n ( t ) , L : [ a , b ] × R 3 n + 1 R , and Z : [ a , b ] × R 3 n R , then the associated fractional Euler–Lagrange equations are obtained as in (2), replacing x with x i , for i { 1 , , n } , obtaining in this way the n fractional differential equations:
L x i [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x i [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x i [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x i [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x i [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x i [ x ¯ ] Z ( t ) 1 g ( t ) = 0 .
If some convexity condition over the Lagrange function is assumed, then we can obtain a sufficient condition for optimization. First, we recall the concepts of convex and concave functions.
Definition 2.
Let k { 1 , , m } and Φ : D R m R be a function such that Φ t i exists and is continuous for all i { k , , m } . We say that Φ is convex in ( t k , , t m ) if
Φ ( t 1 + v 1 , , t m + v m ) Φ ( t 1 , , t m ) Φ t k ( t 1 , , t m ) v k + + Φ t m ( t 1 , , t m ) v m ,
for all ( t 1 + v 1 , , t m + v m ) , ( t 1 , , t m ) D . We say that Φ is concave in ( t k , , t m ) if Equation (4) holds, replacing “≥” with “≤”.
Using the same notation as used before, the result reads as follows:
Theorem 4.
Assume that x ¯ satisfies Equation (2). If L is convex in ( x , C D a + α , g x , I a + β , g x , z ) and one of the two following conditions hold:
1. 
Z is convex in ( x , C D a + α , g x , I a + β , g x ) and L z [ x ¯ ] L ( t ) 0 for all t [ a , b ] ;
2. 
Z is concave in ( x , C D a + α , g x , I a + β , g x ) and L z [ x ¯ ] L ( t ) 0 for all t [ a , b ] ,
then x ¯ is a solution of the following variational problem: minimize F given by (1), subject to the boundary conditions x ( a ) = x a and x ( b ) = x b , with x a , x b R .
Proof. 
Let v C 1 [ a , b ] , with v ( a ) = 0 and v ( b ) = 0 . Then, since
F ( x ¯ + v ) F ( x ¯ ) = a b L ( t , x ¯ ( t ) + v ( t ) , C D a + α , g x ¯ ( t ) + C D a + α , g v ( t ) , I a + β , g x ¯ ( t ) + I a + β , g v ( t ) , a t Z ( τ , x ¯ ( τ ) + v ( τ ) , C D a + α , g x ¯ ( τ ) + C D a + α , g v ( τ ) , I a + β , g x ¯ ( τ ) + I a + β , g v ( τ ) ) d τ ) L t , x ¯ ( t ) , C D a + α , g x ¯ ( t ) , I a + β , g x ¯ ( t ) , a t Z ( τ , x ¯ ( τ ) , C D a + α , g x ¯ ( τ ) , I a + β , g x ¯ ( τ ) ) d τ d t ,
using the assumptions of the theorem, we conclude that
F ( x ¯ + v ) F ( x ¯ ) a b [ L x [ x ¯ ] L ( t ) v ( t ) + L C D a + α , g x [ x ¯ ] L ( t ) C D a + α , g v ( t ) + L I a + β , g x [ x ¯ ] L ( t ) I a + β , g v ( t ) + L z [ x ¯ ] L ( t ) a t Z ( τ , x ¯ ( τ ) + v ( τ ) , C D a + α , g x ¯ ( τ ) + C D a + α , g v ( τ ) , I a + β , g x ¯ ( τ ) + I a + β , g v ( τ ) ) d τ Z ( τ , x ¯ ( τ ) , C D a + α , g x ¯ ( τ ) , I a + β , g x ¯ ( τ ) ) d τ ) ] d τ ] d t a b [ L x [ x ¯ ] L ( t ) v ( t ) + L C D a + α , g x [ x ¯ ] L ( t ) C D a + α , g v ( t ) + L I a + β , g x [ x ¯ ] L ( t ) I a + β , g v ( t ) + L z [ x ¯ ] L ( t ) a t Z x [ x ¯ ] Z ( τ ) v ( τ ) + Z C D a + α , g x [ x ¯ ] Z ( τ ) C D a + α , g v ( τ ) + Z I a + β , g x [ x ¯ ] Z ( τ ) I a + β , g v ( τ ) d τ ] d t .
Applying standard and fractional integration by parts, we get
F ( x ¯ + v ) F ( x ¯ ) a b [ L x [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) ] · v ( t ) d t .
Since x ¯ satisfies (2), we conclude that F ( x ¯ + v ) F ( x ¯ ) 0 , and by the arbitrariness of v, we prove that x ¯ is in fact a solution of the variational problem. □
Remark 2.
If, in Theorem 4, we replace “convex” by “concave”, and vice versa, we prove that x ¯ maximizes F.
For example, let x ¯ ( t ) = ( g ( t ) g ( 0 ) ) 2 , where g C 1 [ 0 , 1 ] with g ( t ) > 0 for all t [ 0 , 1 ] . Then, given α ( 0 , 1 ) and β > 0 , we have the following:
C D 0 + α , g x ¯ ( t ) = 2 Γ ( 3 α ) ( g ( t ) g ( 0 ) ) 2 α a n d I 0 + β , g x ¯ ( t ) = 2 Γ ( 3 + β ) ( g ( t ) g ( 0 ) ) 2 + β .
Consider the following functional:
F ( x ) = 0 1 ( C D 0 + α , g x ( t ) 2 Γ ( 3 α ) ( g ( t ) g ( 0 ) ) 2 α 2 + I 0 + β , g x ( t ) 2 Γ ( 3 + β ) ( g ( t ) g ( 0 ) ) 2 + β 2 + z ( t ) ) d t ,
defined on C 1 [ 0 , 1 ] , where the antiderivative is defined as
z ( t ) = 0 t x ( τ ) ( g ( τ ) g ( 0 ) ) 2 2 d τ .
Furthermore, the boundary conditions
x ( 0 ) = 0 a n d x ( 1 ) = ( g ( 1 ) g ( 0 ) ) 2
are imposed in the formulation of the problem. The fractional Euler–Lagrange Equation (2) reads as
( 1 t ) x ( t ) ( g ( t ) g ( 0 ) ) 2 + g ( t ) D 1 α , g C D 0 + α , g x ( t ) 2 Γ ( 3 α ) ( g ( t ) g ( 0 ) ) 2 α 1 g ( t ) + g ( t ) I 1 β , g I 0 + β , g x ( t ) 2 Γ ( 3 + β ) ( g ( t ) g ( 0 ) ) 2 + β 1 g ( t ) = 0 ,
for t [ 0 , 1 ] . Taking into consideration the formulae presented in (5), we conclude that x ¯ satisfies the Euler–Lagrange equation. Furthermore, since the Lagrange function is convex and condition 1. of Theorem 4 is satisfied, then we conclude that x ¯ is a minimizer of the given functional. Observe that for all x C 1 [ 0 , 1 ] , F ( x ) 0 and F ( x ¯ ) = 0 .
In the previous theorems, the lower bounds of the fractional operators are the same as the lower bound of the integral of the Lagrange function ( t = a ). We can extend this result by considering a different lower bound for the integral. Let A ( a , b ) and consider the functional
F A ( x ) = A b L [ x ¯ ] L ( t ) d t , w h e r e z ( t ) = A t Z [ x ¯ ] Z ( τ ) d τ .
In this case, the necessary conditions are given in the next theorem.
Theorem 5.
Assume that x ¯ minimizes or maximizes functional F A . Then, the following equations hold:
g ( t ) · D b α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) g ( t ) · D A α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) g ( t ) · I A β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0 ,
and
L x [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0 .
for all t [ a , A ] and for all t [ A , b ] , respectively.
Furthermore, we can obtain the following natural boundary conditions:
1. 
If x ( b ) is free, then
I b 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0 ,
at t = b .
2. 
If x ( A ) is free, then
I A 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0 ,
at t = A .
3. 
If x ( a ) is also free, then
I b 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) = I A 1 α , g L C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) ,
at t = a .
Proof. 
The proof is similar to the one presented in Theorem 3, but before applying the fractional integration by parts formulae, the integral defined on the interval [ A , b ] is first split into two in the following way
A b d t = a b d t a A d t ,
and only after applying the fractional integration by parts formulae to each integral. □

4. Optimization Problems Subject to Constraints

In the previous section, we already considered variational problems under some additional constraint (boundary conditions). In this section, we impose two type of constraints on the set of admissible functions to the problem. In the first one, an integral constraint is imposed (knows as an isoperimetric problem) and in the second one, a relation between the position variables and time (known as a holonomic constraint) is established.
For the first case, we assume that we intend to extremize functional (1) subject to the constraint
G ( x ) = a b L G ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) d t = K , K R ,
with the antiderivative
z ( t ) = a t Z G ( τ , x ( τ ) , C D a + α , g x ( τ ) , I a + β , g x ( τ ) ) d τ ,
with α , β , L G and Z G as before. Furthermore, we say that a curve x is an extremal of functional (6) if is satisfies the Euler–Lagrange Equation (2); that is, if we replace L by L G and Z by Z G in (2), then x is a solution of such equation.
Theorem 6.
Let x ¯ be a solution of the variational problem:
extremize F in (1), subject to the constraint G ( x ) = K , ( K R ) and the boundary conditions
x ( a ) = x a , x ( b ) = x b ( x a , x b R ) .
If x ¯ is not an extremal for functional G, then there exists a Lagrange multiplier λ R such that, if we define functions H : [ a , b ] × R 4 R , H ¯ : [ a , b ] × R 4 R 2 , and h ¯ : [ a , b ] × R 3 R 2 by H = L + λ L G , H ¯ = ( L , λ L G ) , and h ¯ = ( Z , Z G ) , then for all t [ a , b ] ,
H x [ x ¯ ] L ( t ) + t b H ¯ z [ x ¯ ] L ( τ ) , h ¯ x [ x ¯ ] Z ( t ) d τ + g ( t ) · D b α , g H C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b H ¯ z [ x ¯ ] L ( τ ) , h ¯ C D a + α , g x [ x ¯ ] Z ( t ) d τ 1 g ( t ) + g ( t ) · I b β , g H I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b H ¯ z [ x ¯ ] L ( τ ) , h ¯ I a + β , g x [ x ¯ ] Z ( t ) d τ 1 g ( t ) = 0 .
where · , · denotes the usual inner product of R 2 .
Proof. 
Consider a variation of x ¯ of the form [ a , b ] t x ¯ ( t ) + ϵ 1 v 1 ( t ) + ϵ 2 v 2 ( t ) , where v 1 , v 2 C 1 [ a , b ] with v i ( a ) = v i ( b ) = 0 , for i = 1 , 2 , and ϵ 1 , ϵ 2 R . Next, define the functions f , g , in a neighborhood of ( 0 , 0 ) , as
f ( ϵ 1 , ϵ 2 ) = F ( x ¯ + ϵ 1 v 1 + ϵ 2 v 2 ) and g ( ϵ 1 , ϵ 2 ) = G ( x ¯ + ϵ 1 v 1 + ϵ 2 v 2 ) K .
Since
g ϵ 2 ( 0 , 0 ) = a b [ L G x [ x ¯ ] L ( t ) + t b L G z [ x ¯ ] L ( τ ) d τ · Z G x [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L G C D a + α , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L G z [ x ¯ ] L ( τ ) d τ · Z G C D a + α , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L G I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L G z [ x ¯ ] L ( τ ) d τ · Z G I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) ] · v 2 ( t ) d t ,
and x ¯ is not an extremal for functional G, then there exists a function v 2 for which g ϵ 2 ( 0 , 0 ) 0 . Furthermore, g ( 0 , 0 ) = 0 and so, by the Implicit Function Theorem, there exists a neighborhood of 0 and a differentiable function ϕ , defined on that neighborhood, such that g ( ϵ 1 , ϕ ( ϵ 1 ) ) = 0 .
Now, we prove the desired necessary condition. For this, observe that the variational problem is equivalent to the following finite dimensional one:
extremize   f ,   subject   to   the   constraint   g = 0 ,
and since ( 0 , 0 ) is a solution of this problem, and we proved that g ( 0 , 0 ) ( 0 , 0 ) , by the Lagrange multiplier rule, there exists a real λ such that ( f + λ g ) ( 0 , 0 ) = ( 0 , 0 ) . Computing ( f + λ g ) ϵ 1 ( 0 , 0 ) = 0 , we prove the result. □
Remark 3.
In Theorem 6, we can drop the assumption that “ x ¯ is not an extremal for functional G”. In such a case, we can prove that there exist Lagrange multipliers λ 0 , λ (both not zero) such that, if we define functions by H = λ 0 L + λ L G and H ¯ = ( λ 0 L , λ L G ) (and h ¯ as before), the result is still valid. In fact, if x ¯ is not an extremal for functional G, then we consider λ 0 = 1 and, in the other case, we take λ 0 = 0 and λ = 1 .
In the next problem that we will study, a holonomic constraint is imposed in the problem’s formulation. Furthermore, we will consider a two-dimensional state function, that is, x = ( x 1 , x 2 ) , where its fractional derivative and integral are given by C D a + α , g x = C D a + α , g x 1 , C D a + α , g x 2 and I a + β , g x = I a + β , g x 1 , I a + β , g x 2 , respectively. The functional to be considered is the following:
F ( x ) = a b L ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) d t ,
with
z ( t ) = a t Z ( τ , x ( τ ) , C D a + α , g x ( τ ) , I a + β , g x ( τ ) ) d τ ,
where L : [ a , b ] × R 7 R and Z : [ a , b ] × R 6 R are two continuously differentiable functions, and x = ( x 1 , x 2 ) C 1 [ a , b ] × C 1 [ a , b ] . Furthermore, we will assume some boundary conditions x ( a ) = X a and x ( b ) = X b for some fixed X a , X b R 2 . Again, we will use the notation
[ x ] L ( t ) = ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) a n d [ x ] Z ( t ) = ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) ) ,
for t [ a , b ] , where x is a two-dimensional vector.
The holonomic constraint that all admissible functions to the problem must satisfy is given by
Λ ( t , x 1 ( t ) , x 2 ( t ) ) = 0 , t [ a , b ] ,
where Λ : [ a , b ] × R 2 R is a C 1 function. The next result provides a necessary condition for optimizing functional (8).
Theorem 7.
Assume that x ¯ = ( x ¯ 1 , x ¯ 2 ) minimizes or maximizes functional (8), and that
Λ x 2 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) 0 , t [ a , b ] .
Then, there exists a continuous function λ : [ a , b ] R such that, for all t [ a , b ] , equations
L x 1 [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x 1 [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x 1 [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x 1 [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x 1 [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x 1 [ x ¯ ] Z ( t ) 1 g ( t ) + λ ( t ) Λ x 1 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) = 0
and
L x 2 [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x 2 [ x ¯ ] Z ( t ) + g ( t ) · D b α , g L C D a + α , g x 2 [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α , g x 2 [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x 2 [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x 2 [ x ¯ ] Z ( t ) 1 g ( t ) + λ ( t ) Λ x 2 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) = 0
hold.
Proof. 
To simplify the notation, by E L E 1 and E L E 2 we mean the Euler–Lagrange Equation (3) with respect to the functions x 1 and x 2 , respectively. If we define function λ as
λ ( t ) = E L E 2 Λ x 2 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) ,
then (11) is proved immediately. To prove (10), a variational technique is used. Consider a variation of the optimal curve, given by [ a , b ] t x ¯ ( t ) + ϵ v ( t ) , where v = ( v 1 , v 2 ) C 1 [ a , b ] × C 1 [ a , b ] and ϵ R . Since we are assuming boundary conditions at t = a and t = b , we need to assume that v ( a ) = ( 0 , 0 ) and v ( b ) = ( 0 , 0 ) . Furthermore, since the variation curve must satisfy the holonomic constraint (9), that is,
Λ ( t , x ¯ 1 ( t ) + ϵ v 1 ( t ) , x ¯ 2 ( t ) + ϵ v 2 ( t ) ) = 0 , t [ a , b ] ,
we conclude that
Λ x 1 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) v 1 ( t ) + Λ x 2 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) v 2 ( t ) = 0 , t [ a , b ] .
On the other hand, the first variation of functional (8) is zero when evaluated at the optimal curve, and so after performing integration by parts, as explained in the previous section, we arrive at
a b [ E L E 1 ] v 1 ( t ) + [ E L E 2 ] v 2 ( t ) d t = 0 .
By (12) and the definition of function λ , we conclude that
[ E L E 2 ] v 2 ( t ) = λ ( t ) Λ x 2 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) v 2 ( t ) = λ ( t ) Λ x 1 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) v 1 ( t ) ,
and so (13) becomes
a b [ E L E 1 ] + λ ( t ) Λ x 1 ( t , x ¯ 1 ( t ) , x ¯ 2 ( t ) ) v 1 ( t ) d t = 0 ,
proving formula (10) from the arbitrariness of function v 1 . □
As an example, again consider the curve x ¯ ( t ) = ( g ( t ) g ( 0 ) ) 2 , where g C 1 [ 0 , 1 ] with g ( t ) > 0 for all t [ 0 , 1 ] . Consider the following isoperimetric problem: extremize the functional
F ( x ) = 0 1 ( C D 0 + α , g x ( t ) 2 + 2 Γ ( 3 α ) ( g ( t ) g ( 0 ) ) 2 α 2 + I 0 + β , g x ( t ) 2 Γ ( 3 + β ) ( g ( t ) g ( 0 ) ) 2 + β 2 + z ( t ) ) d t ,
defined on C 1 [ 0 , 1 ] , where
z ( t ) = 0 t x ( τ ) ( g ( τ ) g ( 0 ) ) 2 2 d τ ,
subject to the integral constraint
0 1 C D 0 + α , g x ( t ) · 2 Γ ( 3 α ) ( g ( t ) g ( 0 ) ) 2 α d t = K
and the boundary conditions
x ( 0 ) = 0 a n d x ( 1 ) = ( g ( 1 ) g ( 0 ) ) 2 .
Here,
K = 0 1 2 Γ ( 3 α ) ( g ( t ) g ( 0 ) ) 2 α 2 d t .
If we take λ = 2 , we can show that x ¯ is a solution of the fractional differential equation given in Theorem 6.

5. Optimization Problems with Arbitrary Fractional Orders

In this section, we allow the fractional order of the derivative to take any positive real value. More precisely, let n N , α k R with α k ( k 1 , k ) , for k { 1 , , n } , and define α = ( α 1 , , α n ) and
C D a + α , g x = C D a + α 1 , g x , , C D a + α n , g x .
Consider the following functional with dependence on higher-order fractional derivatives:
F ( x ) = a b L ( t , x ( t ) , C D a + α , g x ( t ) , I a + β , g x ( t ) , z ( t ) ) d t ,
with the antiderivative
z ( t ) = a t Z ( τ , x ( τ ) , C D a + α , g x ( τ ) , I a + β , g x ( τ ) ) d τ ,
where C D a + α , g x is understood as being the vector defined in (14), and functions L : [ a , b ] × R n + 3 R and Z : [ a , b ] × R n + 2 R are of class C 1 . Furthermore, [ x ] L and [ x ] Z denote the same vectors as before, with C D a + α , g x as defined in (14). We will also assume that x ( k ) ( a ) and x ( k ) ( b ) are fixed for k { 0 , 1 , , n 1 } .
Theorem 8.
If x ¯ is a minimizer of maximizer curve of the functional defined in (15), then
L x [ x ¯ ] L ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z x [ x ¯ ] Z ( t ) + k = 1 n g ( t ) · D b α k , g L C D a + α k , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α k , g x [ x ¯ ] Z ( t ) 1 g ( t ) + g ( t ) · I b β , g L I a + β , g x [ x ¯ ] L ( t ) 1 g ( t ) + t b L z [ x ¯ ] L ( τ ) d τ · Z I a + β , g x [ x ¯ ] Z ( t ) 1 g ( t ) = 0 ,
for all t [ a , b ] .
Proof. 
Let v C n [ a , b ] be a function with v ( k ) ( a ) = 0 and v ( k ) ( b ) = 0 for k { 0 , 1 , , n 1 } , and consider the variation of the optimal curve as being [ a , b ] t x ¯ ( t ) + ϵ v ( t ) . The first variation of the functional, evaluated along this variation, is the following:
a b L x [ x ¯ ] L ( t ) v ( t ) + k = 1 n L C D a + α k , g x [ x ¯ ] L ( t ) C D a + α k , g v ( t ) + L I a + β , g x [ x ¯ ] L ( t ) I a + β , g v ( t ) + L z [ x ¯ ] L ( t ) × a t Z x [ x ¯ ] Z ( τ ) v ( τ ) + k = 1 n Z C D a + α k , g x [ x ¯ ] Z ( τ ) C D a + α k , g v ( τ ) + Z I a + β , g x [ x ¯ ] Z ( τ ) I a + β , g v ( τ ) d τ d t = 0 .
Since v ( k ) ( a ) = 0 and v ( k ) ( b ) = 0 for k { 0 , 1 , , n 1 } , then
1 g ( t ) d d t k v ( t ) = 0 , k { 0 , 1 , , n 1 } ,
when we evaluate it at t = a or at t = b , and so the fractional integration by parts reads as
a b L C D a + α k , g x [ x ¯ ] L ( t ) C D a + α k , g v ( t ) d t = a b D b α k , g L C D a + α k , g x [ x ¯ ] L ( t ) 1 g ( t ) · v ( t ) g ( t ) d t
and
a b L z [ x ¯ ] L ( t ) a t Z C D a + α k , g x [ x ¯ ] Z ( τ ) C D a + α k , g v ( τ ) d τ d t = a b D b α k , g t b L z [ x ¯ ] L ( τ ) d τ · Z C D a + α k , g x [ x ¯ ] Z ( t ) 1 g ( t ) · v ( t ) g ( t ) d t ,
for all k { 0 , 1 , , n 1 } . The rest of the proof is similar to the one given in Theorem 3. □
For example, let n N , α ( n 1 , n ) and β > 0 . Consider the functional
F ( x ) = 0 1 C D 0 + α , g x ( t ) Γ ( α + 2 ) ( g ( t ) g ( 0 ) ) 2 + z ( t ) d t ,
defined on C n [ 0 , 1 ] , subject to the boundary conditions x ( 0 ) = 0 and x ( 1 ) = ( g ( 1 ) g ( 0 ) ) α + 1 , where
z ( t ) = 0 t x ( τ ) ( g ( τ ) g ( 0 ) ) α + 1 2 d τ .
If we define function x ¯ C n [ 0 , 1 ] as x ¯ ( t ) = ( g ( t ) g ( 0 ) ) α + 1 , and since
C D 0 + α , g x ¯ ( t ) = Γ ( α + 2 ) ( g ( t ) g ( 0 ) ) ,
we conclude that x ¯ is a solution for this arbitrary fractional order problem.

6. Conclusions and Future Work

In this paper, we studied some variational problems for the case where the Lagrange functions depend on some fractional operators and on an antiderivative. To make the work as general as possible, we included the fractional integral and fractional derivative of the state function with dependence on an arbitrary kernel. We proved the fractional Euler–Lagrange equation, and under convexity conditions, we showed that in fact it is also a sufficient condition. Then, we imposed some additional constraints on the set of state functions (isoperimetric and holonomic constraints) and deduced the respective necessary conditions for optimization. We ended by studying the case where fractional orders may take any positive real value. With the help of some examples, we show how the variational problem can be solved by computing the solutions of a certain fractional differential equation. For future work, it will be interesting to consider an optimal control approach to these variational problems, where the Lagrange function as well as the state equation depend on the same fractional operators. The goal will be to deduce the fractional Pontryagin’s minimum principle for such a situation. Another path for research would be to develop numerical methods to determine an approximation of the solution to the variational problem. Although there are already numerous works dealing with different fractional operators, we have not yet found many studies dealing with these fractional integrals and derivatives with respect to another function, as presented in the present paper. Usually, such methods consist of discretizing the fractional operators and the cost integral, and thus converting the fractional problem into a (finite) system of difference equations.

Funding

This work was supported by Portuguese funds through the CIDMA—Center for Research and Development in Mathematics and Applications, and the Portuguese Foundation for Science and Technology (FCT-Fundação para a Ciência e a Tecnologia), reference UIDB/04106/2020.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the four anonymous reviewers for their valuable comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Gelfand, I.M.; Fomin, S.V. Calculus of Variations (Revised English Edition. Transl. from Russian); Silverman, R.A., Ed.; Prentice-Hall: Hoboken, NJ, USA, 1963. [Google Scholar]
  2. Ioffe, A.D.; Tihomirov, V.M. Theory of Extremal Problems (Transl. from Russian); Elsevier: Amsterdam, The Netherlands, 1979. [Google Scholar]
  3. Mesterton-Gibbons, M. A Primer on the Calculus of Variations and Optimal Control Theory; American Mathematical Society: Providence, RI, USA, 2009; Volume 50. [Google Scholar]
  4. Hilfer, R. (Ed.) Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  5. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies, 204; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives (Translated from the 1987 Russian Original); Gordon and Breach: New York, NY, USA, 1993. [Google Scholar]
  7. Almeida, R. A Caputo fractional derivative of a function with respect to another function. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 460–481. [Google Scholar] [CrossRef] [Green Version]
  8. Osler, T.J. Leibniz Rule for Fractional Derivatives and an Application to Infinite Series. SIAM J. Appl. Math. 1970, 18, 658–674. [Google Scholar] [CrossRef]
  9. Yang, Y.; Ji, D. Properties of positive solutions for a fractional boundary value problem involving fractional derivative with respect to another function. AIMS Math. 2020, 5, 7359–7371. [Google Scholar] [CrossRef]
  10. Seemab, A.; Rehman, M.; Alzabut, J.; Adjabi, Y.; Abdo, M.S. Langevin equation with nonlocal boundary conditions involving a ψ-Caputo fractional operators of different orders. AIMS Math. 2021, 6, 6749–6780. [Google Scholar] [CrossRef]
  11. Agrawal, O.P. Fractional variational calculus and the transversality conditions. J. Phys. A 2006, 39, 10375. [Google Scholar] [CrossRef]
  12. Almeida, R.; Torres, D.F.M. Calculus of variations with fractional derivatives and fractional integrals. Appl. Math. Lett. 2009, 22, 1816–1820. [Google Scholar] [CrossRef] [Green Version]
  13. Almeida, R.; Torres, D.F.M. Fractional variational calculus for nondifferentiable functions. Comput. Math. Appl. 2011, 61, 3097–3104. [Google Scholar] [CrossRef] [Green Version]
  14. Almeida, R. Optimality conditions for fractional variational problems with free terminal time. Discret. Contin. Dyn. Syst. 2018, 11, 1–19. [Google Scholar] [CrossRef]
  15. Baleanu, D.; Ullah, M.Z.; Mallawi, F.; Saleh Alshomrani, A. A new generalization of the fractional Euler–Lagrange equation for a vertical mass-spring-damper. J. Vib. Control 2021, 27, 2513–2522. [Google Scholar] [CrossRef]
  16. Butkovskii, A.G.; Postnov, S.S.; Postnova, E.A. Fractional integro-differential calculus and its control-theoretical applications. I. Mathematical fundamentals and the problem of interpretation. Autom. Remote Control 2013, 74, 543–574. [Google Scholar] [CrossRef]
  17. Butkovskii, A.G.; Postnov, S.S.; Postnova, E.A. Fractional integro-differential calculus and its control-theoretical applications. II. Fractional dynamic systems: Modeling and hardware implementation. Autom. Remote Control 2013, 74, 725–749. [Google Scholar] [CrossRef]
  18. Baleanu, D.; Tenreiro Machado, J.A.; Luo, A.C.J. Fractional Dynamics and Control; Springer: New York, NY, USA, 2012. [Google Scholar]
  19. Zhou, Y. Fractional Evolution Equations and Inclusions: Analysis and Control; Academic Press: Amsterdam, The Netherlands, 2016. [Google Scholar]
  20. van Brunt, B. The Calculus of Variations; Universitext Springer: New York, NY, USA, 2004. [Google Scholar]
  21. Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E 1997, 55, 3581–3592. [Google Scholar] [CrossRef]
  22. Riewe, F. Nonconservative Lagrangian and Hamiltonian mechanics. Phys. Rev. E 1996, 53, 1890–1899. [Google Scholar] [CrossRef] [PubMed]
  23. Agrawal, O.P. Formulation of Euler–Lagrange equations for fractional variational problems. J. Math. Anal. Appl. 2002, 272, 368–379. [Google Scholar] [CrossRef] [Green Version]
  24. Almeida, R.; Pooseh, S.; Torres, D.F.M. Computational Methods in the Fractional Calculus of Variations; Imperial College Press: London, UK, 2015. [Google Scholar]
  25. Malinowska, A.B.; Torres, D.F.M. Introduction to the Fractional Calculus of Variations; Imperial College Press: London, UK, 2012. [Google Scholar]
  26. Malinowska, A.B.; Odzijewicz, T.; Torres, D.F.M. Advanced Methods in the Fractional Calculus of Variations; Springer Briefs in Applied Sciences and Technology; Springer International Publishing: Cham, Switzerland, 2015. [Google Scholar]
  27. Muslih, S.I.; Baleanu, D. Hamiltonian formulation of systems with linear velocities within Riemann–Liouville fractional derivatives. J. Math. Anal. Appl. 2005, 304, 599–606. [Google Scholar] [CrossRef] [Green Version]
  28. Rabei, E.M.; Nawafleh, K.I.; Hijjawi, R.S.; Muslih, S.I.; Baleanu, D. The Hamilton formalism with fractional derivatives. J. Math. Anal. Appl. 2007, 327, 891–897. [Google Scholar] [CrossRef]
  29. Almeida, R.; Pooseh, S.; Torres, D.F.M. Fractional variational problems depending on indefinite integrals. Nonlinear Anal. 2012, 75, 1009–1025. [Google Scholar] [CrossRef] [Green Version]
  30. Gregory, J. Generalizing variational theory to include the indefinite integral, higher derivatives, and a variety of means as cost variables. Methods Appl. Anal. 2008, 15, 427–435. [Google Scholar] [CrossRef] [Green Version]
  31. Martins, N.; Torres, D.F.M. Generalizing the variational theory on time scales to include the delta indefinite integral. Comput. Math. Appl. 2011, 61, 2424–2435. [Google Scholar] [CrossRef] [Green Version]
  32. Almeida, R. Further properties of Osler’s generalized fractional integrals and derivatives with respect to another function. Rocky Mountain J. Math. 2019, 49, 2459–2493. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Almeida, R. Euler–Lagrange-Type Equations for Functionals Involving Fractional Operators and Antiderivatives. Mathematics 2023, 11, 3208. https://doi.org/10.3390/math11143208

AMA Style

Almeida R. Euler–Lagrange-Type Equations for Functionals Involving Fractional Operators and Antiderivatives. Mathematics. 2023; 11(14):3208. https://doi.org/10.3390/math11143208

Chicago/Turabian Style

Almeida, Ricardo. 2023. "Euler–Lagrange-Type Equations for Functionals Involving Fractional Operators and Antiderivatives" Mathematics 11, no. 14: 3208. https://doi.org/10.3390/math11143208

APA Style

Almeida, R. (2023). Euler–Lagrange-Type Equations for Functionals Involving Fractional Operators and Antiderivatives. Mathematics, 11(14), 3208. https://doi.org/10.3390/math11143208

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop