Next Article in Journal
On Ψ-Hilfer Fractional Integro-Differential Equations with Non-Instantaneous Impulsive Conditions
Next Article in Special Issue
Generalized Beta Models and Population Growth: So Many Routes to Chaos
Previous Article in Journal
Sequential Caputo–Hadamard Fractional Differential Equations with Boundary Conditions in Banach Spaces
Previous Article in Special Issue
Boundary Value Problem for Multi-Term Nonlinear Delay Generalized Proportional Caputo Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Herglotz Variational Problems Involving Distributed-Order Fractional Derivatives with Arbitrary Smooth Kernels

Center for Research and Development in Mathematics and Applications, Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2022, 6(12), 731; https://doi.org/10.3390/fractalfract6120731
Submission received: 24 October 2022 / Revised: 4 December 2022 / Accepted: 8 December 2022 / Published: 10 December 2022
(This article belongs to the Special Issue Feature Papers in Fractal and Fractional 2022–2023)

Abstract

:
In this paper, we consider Herglotz-type variational problems dealing with fractional derivatives of distributed-order with respect to another function. We prove necessary optimality conditions for the Herglotz fractional variational problem with and without time delay, with higher-order derivatives, and with several independent variables. Since the Herglotz-type variational problem is a generalization of the classical variational problem, our main results generalize several results from the fractional calculus of variations. To illustrate the theoretical developments included in this paper, we provide some examples.

1. Introduction

Fractional calculus is as old as ordinary calculus, but it was only at the end of the 20th Century that it managed to attract the attention of many researchers, who in their studies showed that this theory is an important tool to model problems, not only in mathematics, but also in other areas, such as physics, engineering, chemistry, biology, epidemiology, and control theory, among others (see [1,2,3,4,5,6]). Many important mathematicians such as Euler, Lagrange, Fourier, Abel, Liouville, and Riemann worked in fractional calculus. Fractional derivatives and fractional integrals are generalizations of the notions of integer-order derivatives and integrals and include n-th derivatives and n-fold integrals as special cases. Several different fractional derivatives have been defined, such as Riemann–Liouville, Caputo, Riesz, Erdelyi–Kober, and Hadamard, just to mention a few [7,8]. We note that each definition has its own properties and that many of them are not equivalent to each other. In this paper, we deal with the general notions of distributed-order fractional derivatives with respect to an arbitrary kernel in the Riemann–Liouville and Caputo sense, recently introduced in [9]. One of the advantages of fractional derivatives is that these operators are non-local, thus conserving system memories, as opposed to integer-order derivatives, which are local operators.
The calculus of variations deals with the optimization of functionals involving an integral in which the Lagrangian depends on the independent variable, an unknown function, and its derivative (or derivatives). The classical problem of the calculus of variations was generalized by G. Herglotz in 1930 [10], who presented a new problem involving a first-order initial-value problem that defines a function z ( · ) in a given interval [ a , b ] and consists of finding trajectories x and z that extremize the z ( b ) value. This problem is also known as a generalized variational problem. One of the advantages of this problem is that it allows giving a variational description of non-conservative and dissipative processes, even when the Lagrangian is autonomous [11], which is not possible using the classical variational calculus. Herglotz’s problem only attracted the attention of the scientific community in the late Twentieth Century. Thereafter, many important results of the classical calculus of variations were then generalized to Herglotz problems [11,12,13,14,15,16].
The fractional calculus of variations generalizes the classical variational calculus, replacing the integer-order derivatives by fractional derivatives in the Lagrangian function given in the integral of the functional to be extremized. It was in 1996 that the fractional calculus of variations had a considerable development, being relevant to better describe non-conservative systems in mechanics. Furthermore, this theory provides a more realistic approach to physics, allowing it to consider non-conservative systems in a more natural way [6,17]. Since then, this theory has attracted much attention from a large number of researchers, with several articles published [12,18,19,20,21,22,23].
In [9], the authors introduced a new fractional operator, combining two fractional operators: fractional derivatives of distributed-order and fractional derivatives with respect to another function. The order of this new fractional derivative is not constant, and this operator is defined using a function of probability, which acts as a distribution of orders of differentiation, multiplied by a fractional derivative. Our objective in this paper is to study several Herglotz-type problems involving this new fractional derivative.
This paper is organized as follows. In Section 2, we present the classical Herglotz variational problem and some necessary background on fractional calculus. In Section 3, we study the fractional Herglotz problem, in four different cases, using distributed-order fractional derivatives with arbitrary kernels. Namely, we study fractional variational problems of the Herglotz-type for the case where the orders of differentiation belong to the interval [ 0 , 1 ] , for the higher-order case, for problems involving time delay, and with several independent variables. Finally, we present some examples to illustrate our main results.

2. Preliminaries

2.1. Herglotz’s Variational Problem

We begin this section presenting the Herglotz variational problem.
Problem 
( P H ) : Determine trajectories x C 2 ( [ a , b ] , R ) and z C 1 ( [ a , b ] , R ) that extremize
z ( b ) ,
where the pair ( x , z ) satisfies the differential equation:
z ( t ) = L ( t , x ( t ) , x ( t ) , z ( t ) ) , t [ a , b ] ,
with initial condition:
z ( a ) = γ R ,
where it is assumed that the Lagrangian L satisfies the following conditions:
(1)
L C 1 ( [ a , b ] × R 3 , R ) ;
(2)
t L x ( t , x ( t ) , x ( t ) , z ( t ) ) , t L x ( t , x ( t ) , x ( t ) , z ( t ) ) , and t L z ( t , x ( t ) , x ( t ) , z ( t ) )
are differentiable functions for any admissible trajectories ( x , z ) .
It is clear that, if the Lagrangian function L does not depend on the variable z, Problem ( P H ) reduces to the classical problem of the calculus of variations. Herglotz proved that a necessary optimality condition for a pair ( x , z ) to be a local extremizer of Problem ( P H ) is given by the following equation, known as the generalized Euler–Lagrange equation [10]:
L x ( t , x ( t ) , x ( t ) , z ( t ) ) d d t L x ( t , x ( t ) , x ( t ) , z ( t ) ) + L z ( t , x ( t ) , x ( t ) , z ( t ) ) L x ( t , x ( t ) , x ( t ) , z ( t ) ) = 0 , t [ a , b ] .
Obviously, if the Lagrangian function L does not depend on z, then we obtain the famous Euler–Lagrange equation (see [24]):
L x ( t , x ( t ) , x ( t ) ) d d t L x ( t , x ( t ) , x ( t ) ) = 0 , t [ a , b ] .

2.2. Distributed-Order Fractional Calculus with Respect to an Arbitrary Smooth Kernel

For the notions of distributed-order fractional derivatives with respect to an arbitrary smooth kernel in the Riemann–Liouville and Caputo senses, in the case where α [ 0 , 1 ] , we refer the reader to [9].
Next, we recall the extensions of these two derivatives for the case of higher-order derivatives.
Let n N and ϕ : [ n 1 , n ] [ 0 , 1 ] be a continuous function such that
n 1 n ϕ ( α ) d α > 0 .
We start with some needed definitions (see e.g., [25]).
Definition 1.
Let x : [ a , b ] R be an integrable function and ψ C n ( [ a , b ] , R ) be an increasing function such that ψ ( t ) 0 , for all t [ a , b ] . The left and right Riemann–Liouville distributed-order fractional derivatives of a function x with respect to the kernel ψ are defined by:
D a + ϕ ( α ) , ψ x ( t ) : = n 1 n ϕ ( α ) D a + α , ψ x ( t ) d α a n d D b ϕ ( α ) , ψ x ( t ) : = n 1 n ϕ ( α ) D b α , ψ x ( t ) d α ,
where D a + α , ψ and D b α , ψ are the left and right ψ-Riemann–Liouville fractional derivatives of order α [ n 1 , n ] , respectively.
Definition 2.
Let x , ψ C n ( [ a , b ] , R ) be two functions such that ψ is increasing and ψ ( t ) 0 , for all t [ a , b ] . The left and right Caputo distributed-order fractional derivatives of x with respect to ψ are defined by:
C D a + ϕ ( α ) , ψ x ( t ) : = n 1 n ϕ ( α ) C D a + α , ψ x ( t ) d α a n d C D b ϕ ( α ) , ψ x ( t ) : = n 1 n ϕ ( α ) C D b α , ψ x ( t ) d α ,
where C D a + α , ψ and C D b α , ψ are the left and right ψ-Caputo fractional derivatives of order α [ n 1 , n ] , respectively.
In the following, we denote
I a + n ϕ ( α ) , ψ x ( t ) : = n 1 n ϕ ( α ) I a + n α , ψ x ( t ) d α a n d I b n ϕ ( α ) , ψ x ( t ) : = n 1 n ϕ ( α ) I b n α , ψ x ( t ) d α ,
where I a + n α , ψ and I b n α , ψ are, respectively, the left and right Riemann–Liouville fractional integrals of order n α with respect to the kernel ψ . For brevity’s sake, we will use the following notation:
y ψ [ m ] ( t ) : = 1 ψ ( t ) d d t m y ( t ) .
To finalize this section, we present the following result, which is fundamental in the proof of our results (cf. [25]).
Theorem 1
(Integration by parts formulas). Given x : [ a , b ] R a continuous function and y C n ( [ a , b ] , R ) , then
a b x ( t ) C D a + ϕ ( α ) , ψ y ( t ) d t = a b D b ϕ ( α ) , ψ x ( t ) ψ ( t ) y ( t ) ψ ( t ) d t + k = 0 n 1 1 ψ ( t ) d d t k I b n ϕ ( α ) , ψ x ( t ) ψ ( t ) y ψ [ n k 1 ] ( t ) t = a t = b
and
a b x ( t ) C D b ϕ ( α ) , ψ y ( t ) d t = a b D a + ϕ ( α ) , ψ x ( t ) ψ ( t ) y ( t ) ψ ( t ) d t + k = 0 n 1 ( 1 ) n k 1 ψ ( t ) d d t k I a + n ϕ ( α ) , ψ x ( t ) ψ ( t ) y ψ [ n k 1 ] ( t ) t = a t = b .

3. Main Results

In this section, we study four different types of variational problems of the Herglotz type involving distributed-order fractional derivatives with arbitrary smooth kernels.

3.1. Herglotz Fractional Variational Problem—Case 1

For this problem, we restrict ourselves to the case where α [ 0 , 1 ] , that is considering the definitions introduced in [9].
Consider two continuous functions ϕ , φ : [ 0 , 1 ] [ 0 , 1 ] satisfying the following conditions:
0 1 ϕ ( α ) d α > 0 a n d 0 1 φ ( α ) d α > 0 .
In what follows, we use the notation:
[ x , z ] ( t ) : = t , x ( t ) , C D a + ϕ ( α ) , ψ x ( t ) , C D b φ ( α ) , ψ x ( t ) , z ( t )
and we denote the partial derivative of L with respect to its ith-coordinate by i L .
We can formulate the problem as follows:
Problem 
( P H 1 ) : Determine trajectories x C 1 ( [ a , b ] , R ) and z C 1 ( [ a , b ] , R ) that extremize (minimize or maximize)
z ( b )
where the pair ( x , z ) satisfies the differential equation:
z ( t ) = L t , x ( t ) , C D a + ϕ ( α ) , ψ x ( t ) , C D b φ ( α ) , ψ x ( t ) , z ( t ) , t [ a , b ] ,
and
z ( a ) = γ R .
It is assumed that C D a + ϕ ( α ) , ψ x and C D b φ ( α ) , ψ x are of class C 1 , L : [ a , b ] × R 4 R is of class C 1 , and the maps exist and are continuous on [ a , b ] :
t D b ϕ ( α ) , ψ λ ( t ) · 3 L [ x , z ] ( t ) ψ ( t ) and t D a + φ ( α ) , ψ λ ( t ) · 4 L [ x , z ] ( t ) ψ ( t ) ,
for all admissible pairs ( x , z ) , where
λ ( t ) : = e a t 5 L [ x , z ] ( s ) d s , t [ a , b ] .
The following result gives a necessary condition of the Euler–Lagrange type and natural boundary conditions, for an admissible pair ( x , z ) to be a solution of the problem ( P H 1 ) .
Theorem 2.
(Necessary optimality conditions for Problem ( P H 1 ) ) If the pair ( x , z ) is a solution of Problem ( P H 1 ) , then ( x , z ) satisfies the generalized fractional Euler–Lagrange equation:
λ ( t ) 2 L [ x , z ] ( t ) + D b ϕ ( α ) , ψ λ ( t ) 3 L [ x , z ] ( t ) ψ ( t ) ψ ( t ) + D a + φ ( α ) , ψ λ ( t ) 4 L [ x , z ] ( t ) ψ ( t ) ψ ( t ) = 0 ,
for all t [ a , b ] . Furthermore, if x ( a ) is free, then ( x , z ) satisfies the following condition:
I b 1 ϕ ( α ) , ψ λ ( t ) 3 L [ x , z ] ( t ) ψ ( t ) = I a + 1 φ ( α ) , ψ λ ( a ) 4 L [ x , z ] ( t ) ψ ( t ) , a t   t = a ,
and if x ( b ) is free, then ( x , z ) satisfies the following condition:
I b 1 ϕ ( α ) , ψ λ ( t ) 3 L [ x , z ] ( t ) ψ ( t ) = I a + 1 φ ( α ) , ψ λ ( t ) 4 L [ x , z ] ( t ) ψ ( t ) , a t   t = b .
Proof. 
Suppose that the pair ( x , z ) is a solution of Problem ( P H 1 ) and h C 1 ( [ a , b ] , R ) is an arbitrary function, such that its Caputo distributed-order fractional derivatives, C D a + ϕ ( α ) , ψ h and C D b φ ( α ) , ψ h , are continuously differentiable. Define the function β : [ a , b ] R by
β ( t ) : = d d ϵ z [ x + ϵ h ] ( t ) ϵ = 0 .
Since z ( a ) is fixed, we have that β ( a ) = 0 . Now, we define g : ] r , r [ R , where r > 0 , by
g ( ϵ ) = z [ x + ϵ h ] ( b ) .
We have that zero is a local extremizer of g, since z ( b ) is a local extremum, and therefore,
β ( b ) : = d d ϵ z [ x + ϵ h ] ( b ) ϵ = 0 = g ( 0 ) = 0 .
Since
β ( t ) = d d t d d ϵ z [ x + ϵ h ] ( t ) ϵ = 0 = d d ϵ d d t z [ x + ϵ h ] ( t ) ϵ = 0 = d d ϵ L [ x + ϵ h , z ] ( t ) ϵ = 0 = 2 L [ x , z ] ( t ) · h ( t ) + 3 L [ x , z ] ( t ) · C D a + ϕ ( α ) , ψ h ( t ) + 4 L [ x , z ] ( t ) · C D b φ ( α ) , ψ h ( t ) + 5 L [ x , z ] ( t ) · d d ϵ z [ x + ϵ h ] ( t ) ϵ = 0 = 2 L [ x , z ] ( t ) · h ( t ) + 3 L [ x , z ] ( t ) · C D a + ϕ ( α ) , ψ h ( t ) + 4 L [ x , z ] ( t ) · C D b φ ( α ) , ψ h ( t ) + 5 L [ x , z ] ( t ) · β ( t ) ,
then
β ( t ) 5 L [ x , z ] ( t ) · β ( t ) = 2 L [ x , z ] ( t ) · h ( t ) + 3 L [ x , z ] ( t ) · C D a + ϕ ( α ) , ψ h ( t ) + 4 L [ x , z ] ( t ) · C D b φ ( α ) , ψ h ( t ) .
Solving this equation, we obtain
e a t 5 L [ x , z ] ( s ) d s β ( t ) β ( a ) = a t e a s 5 L [ x , z ] ( p ) d p · ( 2 L [ x , z ] ( s ) · h ( s ) + 3 L [ x , z ] ( s ) · C D a + ϕ ( α ) , ψ h ( s ) + 4 L [ x , z ] ( s ) · C D b φ ( α ) , ψ h ( s ) ) d s .
Considering t = b , we obtain
a b λ ( s ) · ( 2 L [ x , z ] ( s ) · h ( s ) + 3 L [ x , z ] ( s ) · C D a + ϕ ( α ) , ψ h ( s ) + 4 L [ x , z ] ( s ) · C D b φ ( α ) , ψ h ( s ) ) d s = 0 .
Using fractional integration by parts (Theorem 1) in Equation (5), we obtain
a b ( λ ( s ) 2 L [ x , z ] ( s ) + D b ϕ ( α ) , ψ λ ( s ) 3 L [ x , z ] ( s ) ψ ( s ) ψ ( s ) + D a + φ ( α ) , ψ λ ( s ) 4 L [ x , z ] ( s ) ψ ( s ) ψ ( s ) ) h ( s ) d s + [ h ( s ) I b 1 ϕ ( α ) , ψ λ ( s ) 3 L [ x , z ] ( s ) ψ ( s ) h ( s ) I a + 1 φ ( α ) , ψ λ ( s ) 4 L [ x , z ] ( s ) ψ ( s ) ] s = a s = b = 0 .
Considering h ( a ) = h ( b ) = 0 in Equation (6), we have
a b ( λ ( s ) 2 L [ x , z ] ( s ) + D b ϕ ( α ) , ψ λ ( s ) 3 L [ x , z ] ( s ) ψ ( s ) ψ ( s ) + D a + φ ( α ) , ψ λ ( s ) 4 L [ x , z ] ( s ) ψ ( s ) ψ ( s ) ) h ( s ) d s = 0 .
From the fundamental lemma of the calculus of variations (see [24]), we obtain
λ ( s ) 2 L [ x , z ] ( s ) + D b ϕ ( α ) , ψ λ ( s ) 3 L [ x , z ] ( s ) ψ ( s ) ψ ( s ) + D a + φ ( α ) , ψ λ ( s ) 4 L [ x , z ] ( s ) ψ ( s ) ψ ( s ) = 0 ,
for all s [ a , b ] , proving the generalized fractional Euler–Lagrange Equation (2). Since h ( a ) is arbitrary if x ( a ) is free, using (2) and considering h ( a ) 0 and h ( b ) = 0 in (6), we obtain
I a + 1 φ ( α ) , ψ λ ( s ) 4 L [ x , z ] ( s ) ψ ( s ) = I b 1 ϕ ( α ) , ψ λ ( s ) 3 L [ x , z ] ( s ) ψ ( s ) , a t s = a ,
proving the natural boundary condition (3). Similarly, since h ( b ) is arbitrary if x ( b ) is free, considering h ( a ) = 0 and h ( b ) 0 in (6) and using (2), we obtain the natural boundary condition (4). □
Remark 1.
We note that if the Lagrangian L does not depend on z, then we obtain as a corollary Theorem 3.2 of [9].

3.2. Herglotz Fractional Variational Problem—Case 2

For this problem, let us consider the case where the Lagrangian depends on higher-order distributed-order fractional derivatives (see Definition 2).
Consider the distribution functions ϕ i , φ i with domains [ i 1 , i ] , i = 1 , , n , where n N is fixed, satisfying the following conditions
i 1 i ϕ i ( α ) d α > 0 a n d i 1 i φ i ( α ) d α > 0 , for all i .
For the simplicity of notation, we consider the following:
[ x , z ] n ( t ) : = t , x ( t ) , C D a + ϕ 1 ( α ) , ψ x ( t ) , C D b φ 1 ( α ) , ψ x ( t ) , , C D a + ϕ n ( α ) , ψ x ( t ) , C D b φ n ( α ) , ψ x ( t ) , z ( t ) .
Problem 
( P H n ) : Determine trajectories x C n ( [ a , b ] , R ) and z C 1 ( [ a , b ] , R ) that extremize
z ( b )
where ( x , z ) satisfies the differential equation:
z ( t ) = L [ x , z ] n ( t ) , t [ a , b ] ,
subject to the boundary condition:
z ( a ) = γ R .
We assume that, for each i = 1 , , n , C D a + ϕ i ( α ) , ψ x and C D b φ i ( α ) , ψ x are all of class C 1 , the Lagrangian function L : [ a , b ] × R 2 n + 2 R is of class C 1 , and the maps exist and are continuous on [ a , b ] :
t D b ϕ i ( α ) , ψ λ ( t ) · 2 i + 1 L [ x , z ] n ( t ) ψ ( t ) and t D a + φ i ( α ) , ψ λ ( t ) · 2 i + 2 L [ x , z ] n ( t ) ψ ( t ) ,
for all admissible pairs ( x , z ) , where
λ ( t ) : = e a t 2 n + 3 L [ x , z ] ( s ) d s , t [ a , b ] .
We are now in a position to present our second result.
Theorem 3.
(Necessary optimality conditions for Problem ( P H n ) ) If the pair ( x , z ) is a solution of Problem ( P H n ) , then ( x , z ) satisfies the generalized fractional Euler–Lagrange equation:
λ ( t ) 2 L [ x , z ] n ( t ) + i = 1 n [ D b ϕ i ( α ) , ψ λ ( t ) 2 i + 1 L [ x , z ] n ( t ) ψ ( t ) ψ ( t ) + D a + φ i ( α ) , ψ λ ( t ) 2 i + 2 L [ x , z ] n ( t ) ψ ( t ) ψ ( t ) ] = 0 ,
for all t [ a , b ] . Furthermore,
(i) 
For a given i = 0 , , n 1 , if x ψ [ i ] ( a ) is free, then
k = i + 1 n [ 1 ψ ( t ) 1 d t k i 1 I b k ϕ k ( α ) , ψ λ ( t ) 2 k + 1 L [ x , z ] n ( t ) ψ ( t ) + ( 1 ) i + 1 1 ψ ( t ) 1 d t k i 1 I a + k φ k ( α ) , ψ λ ( t ) 2 k + 2 L [ x , z ] n ( t ) ψ ( t ) ] = 0 , a t   t = a ;
(ii) 
For a given i = 0 , , n 1 , if x ψ [ i ] ( b ) is free, then
k = i + 1 n [ 1 ψ ( t ) 1 d t k i 1 I b k ϕ k ( α ) , ψ λ ( t ) 2 k + 1 L [ x , z ] n ( t ) ψ ( t ) + ( 1 ) i + 1 1 ψ ( t ) 1 d t k i 1 I a + k φ k ( α ) , ψ λ ( t ) 2 k + 2 L [ x , z ] n ( t ) ψ ( t ) ] = 0 , a t   t = b .
Proof. 
Let h C n ( [ a , b ] , R ) be an arbitrary function such that h ψ [ i ] ( a ) = 0 or h ψ [ i ] ( b ) = 0 , if x ψ [ i ] ( a ) or x ψ [ i ] ( b ) are fixed, respectively, for each i = 0 , , n 1 . Defining β : [ a , b ] R by
β ( t ) : = d d ϵ z [ x + ϵ h ] ( t ) ϵ = 0 ,
then β ( a ) = 0 , β ( b ) = 0 , and
β ( t ) = 2 L [ x , z ] n ( t ) · h ( t ) + i = 1 n ( 2 i + 1 L [ x , z ] n ( t ) · C D a + ϕ i ( α ) , ψ h ( t ) + 2 i + 2 L [ x , z ] n ( t ) · C D b φ i ( α ) , ψ h ( t ) ) + 2 n + 3 L [ x , z ] n ( t ) · β ( t ) .
The solution of Equation (12) is defined by
e a t 2 n + 3 L [ x , z ] ( s ) d s · β ( t ) β ( a ) = a t e a s 2 n + 3 L [ x , z ] ( p ) d p · ( 2 L [ x , z ] n ( s ) · h ( s ) + i = 1 n 2 i + 1 L [ x , z ] n ( s ) · C D a + ϕ i ( α ) , ψ h ( s ) + 2 i + 2 L [ x , z ] n ( s ) · C D b φ i ( α ) , ψ h ( s ) ) d s .
Considering λ ( t ) = e a t 2 n + 3 L [ x , z ] ( s ) d s and taking t = b in (13), Theorem 1 allows us to prove that
a b ( λ ( s ) 2 L [ x , z ] n ( s ) + i = 1 n [ D b ϕ i ( α ) , ψ λ ( s ) 2 i + 1 L [ x , z ] n ( s ) ψ ( s ) ψ ( s ) + D a + φ i ( α ) , ψ λ ( s ) 2 i + 2 L [ x , z ] n ( s ) ψ ( s ) ψ ( s ) ] ) h ( s ) d s + i = 1 n k = 0 i 1 [ ( 1 ψ ( s ) 1 d s k I b i ϕ i ( α ) , ψ λ ( s ) 2 i + 1 L [ x , z ] n ( s ) ψ ( s ) + ( 1 ) i k 1 ψ ( s ) 1 d s k I a + i φ i ( α ) , ψ λ ( s ) 2 i + 2 L [ x , z ] n ( s ) ψ ( s ) ) h ψ [ i k 1 ] ( s ) ] s = a s = b = 0 .
Since
i = 1 n k = 0 i 1 [ ( 1 ψ ( s ) 1 d s k I b i ϕ i ( α ) , ψ λ ( s ) 2 i + 1 L [ x , z ] n ( s ) ψ ( s ) + ( 1 ) i k 1 ψ ( s ) 1 d s k I a + i φ i ( α ) , ψ λ ( s ) 2 i + 2 L [ x , z ] n ( s ) ψ ( s ) ) h ψ [ i k 1 ] ( s ) ] s = a s = b = i = 0 n 1 h ψ [ i ] ( s ) k = i + 1 n [ ( 1 ψ ( s ) 1 d s k i 1 I b k ϕ k ( α ) , ψ λ ( s ) 2 k + 1 L [ x , z ] n ( s ) ψ ( s ) + ( 1 ) i + 1 1 ψ ( s ) 1 d s k i 1 I a + k φ k ( α ) , ψ λ ( s ) 2 k + 2 L [ x , z ] n ( s ) ψ ( s ) ) ] s = a s = b ;
from the arbitrariness of h and using the fundamental lemma of calculus of variations, we have proven the generalized fractional Euler–Lagrange Equation (9) and the necessary conditions (10) and (11). □
Remark 2.
It is easy to see that Theorem 9 of [25] is a corollary of Theorem 3.

3.3. Herglotz Fractional Variational Problem—Case 3

It is well known that time delay is a common phenomenon that occurs in many engineering and nature problems. Therefore, it is extremely important to consider when formulating such problems the explicit dependence of a time delay, in order to better understand the evolution of the dynamical systems under observation. Motivated by the importance of considering a time delay in the formulation of variational problems, we now study the Herglotz variational problem with time delay. For the simplicity of presentation, we restrict ourselves to the case where α [ 0 , 1 ] .
In what follows, τ is a fixed real number such that 0 τ < b a , and in order to simplify the notation, we write:
[ x , z ] τ ( t ) : = t , x ( t ) , x ( t τ ) , C D a + ϕ ( α ) , ψ x ( t ) , C D b φ ( α ) , ψ x ( t ) , z ( t ) .
Problem 
( P H τ ) : Determine x C 1 ( [ a τ , b ] , R ) and z C 1 ( [ a , b ] , R ) that extremize
z ( b )
where
z ( t ) = L [ x , z ] τ ( t ) , t [ a , b ] ;
z ( a ) = γ R and x ( t ) = μ ( t ) on [ a τ , a ] , where μ C 1 ( [ a τ , a ] , R ) is a given initial function.
It is assumed that C D a + ϕ ( α ) , ψ x and C D b φ ( α ) , ψ x are of class C 1 and L satisfies the following conditions:
  • L : [ a , b ] × R 5 R is of class C 1 ;
  • The functions exist and are continuous:
    t D ( b τ ) ϕ ( α ) , ψ λ ( t ) · 4 L [ x , z ] τ ( t ) ψ ( t ) and t D a + φ ( α ) , ψ λ ( t ) · 5 L [ x , z ] τ ( t ) ψ ( t )
    on [ a , b τ ] , and
    t D b ϕ ( α ) , ψ λ ( t ) · 4 L [ x , z ] τ ( t ) ψ ( t ) and t D ( b τ ) + φ ( α ) , ψ λ ( t ) · 5 L [ x , z ] τ ( t ) ψ ( t )
    on [ b τ , b ] , where
    λ ( t ) : = e a t 6 L [ x , z ] τ ( s ) d s , t [ a , b ] .
The next result presents the necessary optimality conditions for the fractional variational problem of the Herglotz type with time delay.
Theorem 4
(Necessary optimality conditions for Problem ( P H τ ) ). If the pair ( x , z ) is a solution of Problem ( P H τ ) , then ( x , z ) satisfies the generalized fractional Euler–Lagrange equations:
λ ( t ) 2 L [ x , z ] τ ( t ) + λ ( t + τ ) 3 L [ x , z ] τ ( t + τ ) + D ( b τ ) ϕ ( α ) , ψ λ ( t ) 4 L [ x , z ] τ ( t ) ψ ( t ) ψ ( t ) + D a + φ ( α ) , ψ λ ( t ) 5 L [ x , z ] τ ( t ) ψ ( t ) ψ ( t ) 0 1 ϕ ( α ) Γ ( 1 α ) d d t b τ b ( ψ ( s ) ψ ( t ) ) α λ ( s ) 4 L [ x , z ] τ ( s ) d s d α = 0 , t [ a , b τ ]
and
λ ( t ) 2 L [ x , z ] τ ( t ) + D b ϕ ( α ) , ψ λ ( t ) 4 L [ x , z ] τ ( t ) ψ ( t ) ψ ( t ) + D ( b τ ) + φ ( α ) , ψ λ ( t ) 5 L [ x , z ] τ ( t ) ψ ( t ) ψ ( t ) + 0 1 φ ( α ) Γ ( 1 α ) d d t a b τ ( ψ ( t ) ψ ( s ) ) α λ ( s ) 5 L [ x , z ] τ ( s ) d s d α = 0 , t [ b τ , b ] .
Furthermore, if x ( b ) is free, then ( x , z ) satisfies the natural boundary condition:
I b 1 ϕ ( α ) , ψ λ ( t ) 4 L [ x , z ] τ ( t ) ψ ( t ) = I a + 1 φ ( α ) , ψ λ ( t ) 5 L [ x , z ] τ ( t ) ψ ( t ) , a t   t = b .
Proof. 
Let h C 1 ( [ a τ , b ] , R ) be an arbitrary function such that h ( t ) = 0 , a τ t a , C D a + ϕ ( α ) , ψ h , and C D b φ ( α ) , ψ h are of class C 1 . Defining function β : [ a , b ] R by
β ( t ) : = d d ϵ z [ x + ϵ h ] τ ( t ) ϵ = 0
we have that β ( a ) = β ( b ) = 0 . Hence, we obtain
β ( t ) = d d t d d ϵ z [ x + ϵ h ] τ ( t ) ϵ = 0 = d d ϵ d d t z [ x + ϵ h ] τ ( t ) ϵ = 0 = d d ϵ L [ x + ϵ h , z ] τ ( t ) ϵ = 0 = 2 L [ x , z ] τ ( t ) · h ( t ) + 3 L [ x , z ] τ ( t ) · h ( t τ ) + 4 L [ x , z ] τ ( t ) · C D a + ϕ ( α ) , ψ h ( t ) + 5 L [ x , z ] τ ( t ) · C D b φ ( α ) , ψ h ( t ) + 6 L [ x , z ] τ ( t ) · β ( t ) .
Solving the differential Equation (17), we obtain
λ ( t ) β ( t ) β ( a ) = a t λ ( s ) ( 2 L [ x , z ] τ ( s ) · h ( s ) + 3 L [ x , z ] τ ( s ) · h ( s τ ) + 4 L [ x , z ] τ ( s ) · C D a + ϕ ( α ) , ψ h ( s ) + 5 L [ x , z ] τ ( s ) · C D b φ ( α ) , ψ h ( s ) ) d s .
Considering t = b and replacing β ( a ) = β ( b ) = 0 in (18), we have
a b ( λ ( s ) 2 L [ x , z ] τ ( s ) · h ( s ) + λ ( s ) 3 L [ x , z ] τ ( s ) · h ( s τ ) + λ ( s ) 4 L [ x , z ] τ ( s ) · C D a + ϕ ( α ) , ψ h ( s ) + λ ( s ) 5 L [ x , z ] τ ( s ) · C D b φ ( α ) , ψ h ( s ) ) d s = 0 .
Since
a b λ ( s ) 3 L [ x , z ] τ ( s ) · h ( s τ ) d s = a b τ λ ( s + τ ) 3 L [ x , z ] τ ( s + τ ) · h ( s ) d s ,
then, from (19), we obtain
a b τ λ ( s ) 2 L [ x , z ] τ ( s ) + λ ( s + τ ) 3 L [ x , z ] τ ( s + τ ) · h ( s ) d s + b τ b λ ( s ) 2 L [ x , z ] τ ( s ) · h ( s ) d s + a b λ ( s ) 4 L [ x , z ] τ ( s ) · C D a + ϕ ( α ) , ψ h ( s ) + λ ( s ) 5 L [ x , z ] τ ( s ) · C D b φ ( α ) , ψ h ( s ) d s = 0 .
Note also that
D b ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) = D ( b τ ) ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) 0 1 ϕ ( α ) Γ ( 1 α ) 1 ψ ( s ) d d s b τ b ( ψ ( p ) ψ ( s ) ) α λ ( p ) 4 L [ x , z ] τ ( p ) d p d α ,
for all s [ a , b τ ] , and
D a + φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) = D ( b τ ) + φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) + 0 1 φ ( α ) Γ ( 1 α ) 1 ψ ( s ) d d s a b τ ( ψ ( s ) ψ ( p ) ) α λ ( p ) 5 L [ x , z ] τ ( p ) d p d α = 0 ,
for all s [ b τ , b ] . Using Equation (21) and Theorem 1, we conclude that
a b λ ( s ) 4 L [ x , z ] τ ( s ) · C D a + ϕ ( α ) , ψ h ( s ) d s = a b τ ( D ( b τ ) ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) 0 1 ϕ ( α ) Γ ( 1 α ) d d s b τ b ( ψ ( p ) ψ ( s ) ) α λ ( p ) 4 L [ x , z ] τ ( p ) d p d α ) h ( s ) d s + b τ b D b ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) h ( s ) d s + I b 1 ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) h ( s ) s = a s = b .
Similarly, using (22) and Theorem 1, we obtain
a b λ ( s ) 5 L [ x , z ] τ ( s ) · C D b φ ( α ) , ψ h ( s ) d s = b τ b ( D ( b τ ) + φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) + 0 1 φ ( α ) Γ ( 1 α ) d d s a b τ ( ψ ( s ) ψ ( p ) ) α λ ( p ) 5 L [ x , z ] τ ( p ) d p d α ) h ( s ) d s + a b τ D a + φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) h ( s ) d s I a + 1 φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) h ( s ) s = a s = b .
Therefore, we obtain
a b τ ( λ ( s ) 2 L [ x , z ] τ ( s ) + λ ( s + τ ) 3 L [ x , z ] τ ( s + τ ) + D ( b τ ) ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) 0 1 ϕ ( α ) Γ ( 1 α ) d d s b τ b ( ψ ( p ) ψ ( s ) ) α λ ( p ) 4 L [ x , z ] τ ( p ) d p d α + D a + φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) ) h ( s ) d s + b τ b ( λ ( s ) 2 L [ x , z ] τ ( s ) + D b ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) + D ( b τ ) + φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) ψ ( s ) + 0 1 φ ( α ) Γ ( 1 α ) d d s a b τ ( ψ ( s ) ψ ( p ) ) α λ ( p ) 5 L [ x , z ] τ ( p ) d p d α ) h ( s ) d s + I b 1 ϕ ( α ) , ψ λ ( s ) 4 L [ x , z ] τ ( s ) ψ ( s ) h ( s ) I a + 1 φ ( α ) , ψ λ ( s ) 5 L [ x , z ] τ ( s ) ψ ( s ) h ( s ) s = a s = b = 0 ,
introducing (23) and (24) into (20). Therefore, choosing the appropriate h and using the fundamental lemma of calculus of variations, we obtain the generalized fractional Euler–Lagrange Equations (14) and (15) and the natural boundary condition (16). □
Remark 3.
It is clear that Theorem 2 of [25] can be obtained from Theorem 4 in the particular case where the Lagrangian is independent of z.

3.4. Herglotz Fractional Variational Problem—Case 4

For this the last problem, we consider the case where the state function depends on several independent variables. Here, we have the case where the fractional orders belong to the interval [ 0 , 1 ] .
We consider U = [ a , b ] × Λ and Λ = i = 1 n [ a i , b i ] . We denote by t [ a , b ] the time variable and s = ( s 1 , , s n ) Λ the spacial coordinates. We use the notation:
[ x , z ] ( t , s ) : = t , s , x ( t , s ) , C D + ϕ ( α ) , ψ x ( t , s ) , C D φ ( α ) , ψ x ( t , s ) , z ( t ) ,
where
C D + ϕ ( α ) , ψ x ( t , s ) = C D a + ϕ ( α ) , ψ x ( t , s ) , C D a 1 + ϕ ( α ) , ψ x ( t , s ) , , C D a n + ϕ ( α ) , ψ x ( t , s ) R n + 1
and
C D φ ( α ) , ψ x ( t , s ) = C D b φ ( α ) , ψ x ( t , s ) , C D b 1 φ ( α ) , ψ x ( t , s ) , , C D b n φ ( α ) , ψ x ( t , s ) R n + 1 ,
and C D a + ϕ ( α ) , ψ x and C D b φ ( α ) , ψ x denote the left and right partial distributed-order fractional derivatives of x with respect to the variable t; C D a i + ϕ ( α ) , ψ x and C D b i φ ( α ) , ψ x , for any i { 1 , , n } , denote the left and right partial distributed-order fractional derivatives of x with respect to the variable s i . Furthermore, it is assumed that the domain of function ψ contains the intervals [ a , b ] and [ a i , b i ] , for i = 1 , , n .
Problem 
( P H * ) : Determine trajectories x C 1 ( U , R ) and z C 1 ( [ a , b ] , R ) that extremize
z ( b )
such that the pair ( x , z ) satisfies the differential equation:
z ( t ) = Λ L [ x , z ] ( t , s ) d n s , t [ a , b ] , d n s = d s 1 d s n ,
and
z ( a ) = γ R .
Furthermore, we assume that x ( t , s ) is fixed when s i = a i and s i = b i , for all t [ a , b ] . We suppose that C D a + ϕ ( α ) , ψ x , C D b φ ( α ) , ψ x , C D a i + ϕ ( α ) , ψ x , and C D b i φ ( α ) , ψ x are of class C 1 for i = 1 , , n , L : [ a , b ] × R 3 n + 4 R is a continuously differentiable function, and the maps exist and are continuous on U:
( t , s ) D b ϕ ( α ) , ψ λ ( t ) · n + 3 L [ x , z ] ( t , s ) ψ ( t ) , ( t , s ) D a + φ ( α ) , ψ λ ( t ) · 2 n + 4 L [ x , z ] ( t , s ) ψ ( t ) ,
and
( t , s ) D b i ϕ ( α ) , ψ λ ( t ) · n + 3 + i L [ x , z ] ( t , s ) ψ ( s i ) , ( t , s ) D a i + φ ( α ) , ψ λ ( t ) · 2 n + 4 + i L [ x , z ] ( t , s ) ψ ( s i ) ,
for all admissible pairs ( x , z ) and for all i = 1 , , n , where
λ ( t ) : = e a t Λ 3 n + 5 L [ x , z ] ( p , s ) d n s d p , t [ a , b ] .
Under these assumptions, we can prove our last result.
Theorem 5
(Necessary optimality conditions for Problem ( P H * ) ). If the pair ( x , z ) is a solution of Problem ( P H * ) , then ( x , z ) satisfies the generalized fractional Euler–Lagrange equation:
λ ( t ) n + 2 L [ x , z ] ( t , s ) + D b ϕ ( α ) , ψ λ ( t ) n + 3 L [ x , z ] ( t , s ) ψ ( t ) ψ ( t ) + D a + φ ( α ) , ψ λ ( t ) 2 n + 4 L [ x , z ] ( t , s ) ψ ( t ) ψ ( t ) + i = 1 n [ D b i ϕ ( α ) , ψ λ ( t ) n + 3 + i L [ x , z ] ( t , s ) ψ ( s i ) ψ ( s i ) + D a i + φ ( α ) , ψ λ ( t ) 2 n + 4 + i L [ x , z ] ( t , s ) ψ ( s i ) ψ ( s i ) ] = 0 ,
for all ( t , s ) U . Furthermore, if x ( a , · ) is free, then
Λ I b 1 ϕ ( α ) , ψ λ ( t ) n + 3 L [ x , z ] ( t , s ) ψ ( t ) d n s = Λ I a + 1 φ ( α ) , ψ λ ( t ) 2 n + 4 L [ x , z ] ( t , s ) ψ ( t ) d n s , a t   t = a ,
and if x ( b , · ) is free, then
Λ I a + 1 φ ( α ) , ψ λ ( t ) 2 n + 4 L [ x , z ] ( t , s ) ψ ( t ) d n s = Λ I b 1 ϕ ( α ) , ψ λ ( t ) n + 3 L [ x , z ] ( t , s ) ψ ( t ) d n s , a t   t = b .
Proof. 
Let h C 1 ( U , R ) be an arbitrary function such that its Caputo distributed-order fractional derivatives are continuously differentiable. Because the state function is fixed when s i = a i and s i = b i , we suppose that, for any i { 1 , , n } , if s i = a i or s i = b i , then h ( t , s ) = 0 for all t [ a , b ] . Defining function β : [ a , b ] R by
β ( t ) : = d d ϵ z [ x + ϵ h ] ( t ) ϵ = 0 ,
then β ( a ) = β ( b ) = 0 and
β ( t ) = d d t d d ϵ z [ x + ϵ h ] ( t ) ϵ = 0 = d d ϵ d d t z [ x + ϵ h ] ( t ) ϵ = 0 = d d ϵ Λ L [ x + ϵ h , z ] ( t , s ) d n s ϵ = 0 = Λ d d ϵ L [ x + ϵ h , z ] ( t , s ) d n s ϵ = 0 = Λ ( n + 2 L [ x , z ] ( t , s ) · h ( t , s ) + n + 3 L [ x , z ] ( t , s ) · C D a + ϕ ( α ) , ψ h ( t , s ) + 2 n + 4 L [ x , z ] ( t , s ) · C D b φ ( α ) , ψ h ( t , s ) + i = 1 n [ n + 3 + i L [ x , z ] ( t , s ) · C D a i + ϕ ( α ) , ψ h ( t , s ) + 2 n + 4 + i L [ x , z ] ( t , s ) · C D b i φ ( α ) , ψ h ( t , s ) ] ) d n s + β ( t ) · Λ 3 n + 5 L [ x , z ] ( t ) d n s ;
we obtain
a b λ ( t ) · Λ ( n + 2 L [ x , z ] ( t , s ) · h ( t , s ) + n + 3 L [ x , z ] ( t , s ) · C D a + ϕ ( α ) , ψ h ( t , s ) + 2 n + 4 L [ x , z ] ( t , s ) · C D b φ ( α ) , ψ h ( t , s ) + i = 1 n [ n + 3 + i L [ x , z ] ( t , s ) · C D a i + ϕ ( α ) , ψ h ( t , s ) + 2 n + 4 + i L [ x , z ] ( t , s ) · C D b i φ ( α ) , ψ h ( t , s ) ] ) d n s d t = 0 .
Using Theorem 1, we obtain
a b Λ ( λ ( t ) n + 2 L [ x , z ] ( t , s ) + D b ϕ ( α ) , ψ λ ( t ) n + 3 L [ x , z ] ( t , s ) ψ ( t ) ψ ( t ) + D a + φ ( α ) , ψ λ ( t ) 2 n + 4 L [ x , z ] ( t , s ) ψ ( t ) ψ ( t ) + i = 1 n [ D b i ϕ ( α ) , ψ λ ( t ) n + 3 + i L [ x , z ] ( t , s ) ψ ( s i ) ψ ( s i ) + D a i + φ ( α ) , ψ λ ( t ) 2 n + 4 + i L [ x , z ] ( t , s ) ψ ( s i ) ψ ( s i ) ] ) h ( t , s ) d n s d t + [ Λ I b 1 ϕ ( α ) , ψ λ ( t ) n + 3 L [ x , z ] ( t , s ) ψ ( t ) h ( t , s ) d n s Λ I a + 1 φ ( α ) , ψ λ ( t ) 2 n + 4 L [ x , z ] ( t , s ) ψ ( t ) h ( t , s ) d n s ] t = a t = b = 0 .
From the arbitrariness of h and using the fundamental lemma of the calculus of variations, we obtain the generalized fractional Euler–Lagrange Equation (25) and the natural boundary conditions (26) and (27). □

4. Illustrative Examples

In order to illustrate our results, we present three examples.
Example 1.
Consider the following fractional differential equation:
z ( t ) = C D 0 + ϕ ( α ) , ψ x ( t ) · ( ψ ( t ) ψ ( 0 ) ) α + 1 ( ψ ( t ) ψ ( 0 ) ) α + 3 4 + ( t cos t ) z ( t ) ,
for t [ 0 , 1 ] , where z ( 0 ) = 4 , x ( 0 ) = 0 and x ( 1 ) = ( ψ ( 1 ) ψ ( 0 ) ) α + 2 . Define ϕ : [ 0 , 1 ] [ 0 , 1 ] by
ϕ ( α ) = 2 Γ ( α + 3 ) .
If we consider x ¯ ( t ) = ( ψ ( t ) ψ ( 0 ) ) α + 2 , t [ 0 , 1 ] , then by Lemma 1 of [26], we obtain
C D 0 + α , ψ x ¯ ( t ) = Γ ( α + 3 ) 2 ( ψ ( t ) ψ ( 0 ) ) 2 .
Thus,
C D 0 + ϕ ( α ) , ψ x ¯ ( t ) = 0 1 ϕ ( α ) C D 0 + α , ψ x ¯ ( t ) d α = ( ψ ( t ) ψ ( 0 ) ) 2 .
Note that x ¯ satisfies the necessary optimality condition (2):
D 1 ϕ ( α ) , ψ 4 ( ψ ( t ) ψ ( 0 ) ) α + 1 λ ( t ) C D 0 + ϕ ( α ) , ψ x ( t ) · ( ψ ( t ) ψ ( 0 ) ) α + 1 ( ψ ( t ) ψ ( 0 ) ) α + 3 3 ψ ( t ) · ψ ( t ) = 0 , t [ 0 , 1 ] ,
where λ ( t ) = e t 2 2 + sin t . Therefore,
x ¯ ( t ) = ( ψ ( t ) ψ ( 0 ) ) α + 2 and z ( t ) = 4 e t 2 2 sin t
is a candidate to be a local extremizer of the value z ( 1 ) .
Example 2.
Let x C 4 ( [ 1 , 5 ] , R ) and z C 1 ( [ 1 , 5 ] , R ) such that the pair ( x , z ) is a solution of the fractional differential equation:
z ( t ) = 4 t + C D 1 + ϕ 4 ( α ) , ψ x ( t ) ( ψ ( t ) ψ ( 1 ) ) 2 ψ ( t ) + ψ ( 1 ) ln ( ψ ( t ) ψ ( 1 ) ) 5 2 t z ( t ) ,
and for i = 0 , 1 , 2 , 3 , we have x ψ [ i ] ( 1 ) = 0 , x ψ [ i ] ( 5 ) is free, and z ( 1 ) = 2 , where ϕ 4 : [ 3 , 4 ] [ 0 , 1 ] is defined by
ϕ 4 ( α ) = Γ ( 6 α ) 120 .
Note that, if x ¯ ( t ) = ( ψ ( t ) ψ ( 1 ) ) 5 , t [ 1 , 5 ] , then, by Lemma 1 of [26], we obtain
C D 1 + α , ψ x ¯ ( t ) = 120 Γ ( 6 α ) ( ψ ( t ) ψ ( 1 ) ) 5 α .
Hence,
C D 1 + ϕ 4 ( α ) , ψ x ¯ ( t ) = 3 4 ϕ 4 ( α ) C D 1 + α , ψ x ¯ ( t ) d α = ( ψ ( t ) ψ ( 1 ) ) 2 ψ ( t ) + ψ ( 1 ) ln ( ψ ( t ) ψ ( 1 ) ) .
Note that x ¯ satisfies the generalized fractional Euler–Lagrange Equation (9):
D 5 ϕ 4 ( α ) , ψ 5 ( t 2 1 ) C D 1 + ϕ 4 ( α ) , ψ x ( t ) ( ψ ( t ) ψ ( 1 ) ) 2 ψ ( t ) + ψ ( 1 ) ln ( ψ ( t ) ψ ( 1 ) ) 4 ψ ( t ) ψ ( t ) = 0 , t [ 1 , 5 ] ,
and the natural boundary conditions (11), for each i = 0 , 1 , 2 , 3 ,
k = i + 1 4 1 ψ ( t ) d d t k i 1 I 5 k ϕ 4 ( α ) , ψ 5 ( t 2 1 ) C D 1 + ϕ 4 ( α ) , ψ x ( t ) ( ψ ( t ) ψ ( 1 ) ) 2 ψ ( t ) + ψ ( 1 ) ln ( ψ ( t ) ψ ( 1 ) ) 4 ψ ( t ) = 0 , a t t = 5 .
Therefore, by Theorem 3, the pair ( x ¯ , z ) where x ¯ ( t ) = ( ψ ( t ) ψ ( 1 ) ) 5 and z ( t ) = t 2 + 1 t 2 , t [ 1 , 5 ] , is a candidate to be a local extremizer of the value z ( 5 ) .
Example 3.
Consider x C 1 ( [ 2 , 4 ] , R ) and z C 1 ( [ 0 , 4 ] , R ) such that the pair ( x , z ) is a solution of the following fractional differential equation:
z ( t ) = x ( t 2 ) ( ψ ( 4 ) ψ ( t 2 ) ) 3 2 + C D 4 φ ( α ) , ψ x ( t ) · ( ψ ( 4 ) ψ ( t ) ) α ( ψ ( 4 ) ψ ( t ) ) α + 2 2 ln ( ψ ( 4 ) ψ ( t ) ) 2 + 1 cos ( t / 4 ) z ( t ) ,
where z ( 0 ) = 3 , μ ( t ) = ( ψ ( 4 ) ψ ( t ) ) 3 , t [ 2 , 0 ] , and x ( 4 ) = 0 . Let φ : [ 0 , 1 ] [ 0 , 1 ] defined by
φ ( α ) = Γ ( 4 α ) 6 .
If x ¯ ( t ) = ( ψ ( 4 ) ψ ( t ) ) 3 , t [ 2 , 4 ] , then, by Lemma 1 of [26], we have
C D 4 α , ψ x ¯ ( t ) = 6 Γ ( 4 α ) ( ψ ( 4 ) ψ ( t ) ) 3 α ,
and the distributed-order fractional derivative of x ¯ is given by
C D 4 φ ( α ) , ψ x ¯ ( t ) = ( ψ ( 4 ) ψ ( t ) ) 3 ( ψ ( 4 ) ψ ( t ) ) 2 ln ( ψ ( 4 ) ψ ( t ) ) .
Note that x ¯ satisfies the generalized fractional Euler–Lagrange Equations (14) and (15):
2 λ ( t + 2 ) x ( t + 2 ) ( ψ ( 4 ) ψ ( t + 2 ) ) 3 + D 0 + φ ( α ) , ψ 2 λ ( t ) C D 4 φ ( α ) , ψ x ( t ) ( ψ ( 4 ) ψ ( t ) ) 3 ( ψ ( 4 ) ψ ( t ) ) 2 ln ( ψ ( 4 ) ψ ( t ) ) ( ψ ( 4 ) ψ ( t ) ) α ψ ( t ) ψ ( t ) = 0 ,
for all t [ 0 , 2 ] , and
D 2 + φ ( α ) , ψ 2 λ ( t ) C D 4 φ ( α ) , ψ x ( t ) ( ψ ( 4 ) ψ ( t ) ) 3 ( ψ ( 4 ) ψ ( t ) ) 2 ln ( ψ ( 4 ) ψ ( t ) ) ( ψ ( 4 ) ψ ( t ) ) α ψ ( t ) ψ ( t ) + 0 1 φ ( α ) Γ ( 1 α ) ( d d t 0 2 ( ψ ( t ) ψ ( s ) ) α 2 λ ( s ) C D 4 φ ( α ) , ψ x ( s ) ( ψ ( 4 ) ψ ( s ) ) 3 ( ψ ( 4 ) ψ ( s ) ) 2 ln ( ψ ( 4 ) ψ ( s ) ) ( ψ ( 4 ) ψ ( s ) ) α d s ) d α = 0 ,
for all t [ 2 , 4 ] , where λ ( t ) = e 0 t sec ( 4 s ) d s . Therefore, by Theorem 4, ( x ¯ , z ) , where
x ¯ ( t ) = ( ψ ( 4 ) ψ ( t ) ) 3 and z ( t ) = 3 ( sec ( t / 4 ) + tan ( t / 4 ) ) 4
is a candidate to be a local extremizer of the value z ( 4 ) .

5. Concluding Remarks

In this paper, we studied four cases of the fractional-Herglotz-variational-type problems, where the Lagrangian depends on distributed-order fractional derivatives with arbitrary smooth kernels. In the first case, the distributed order belongs to [ 0 , 1 ] , by considering the definitions introduced in [9]. In the second case, we considered the higher-order case, that is when α [ n 1 , n ] for a given n N , considering the definitions recently introduced in [25]. In the third case, we studied the Herglotz variational problem with time delay, and in the last case, we considered the Herglotz variational problem with several independent variables. We proved the necessary optimality conditions for all of these Herglotz-type problems, and three examples were presented to illustrate our results. To finalize this paper, we point out that our theoretical contributions generalize several results recently proven in the context of the fractional calculus of variations.

Author Contributions

Conceptualization, F.C., R.A. and N.M.; methodology, F.C., R.A. and N.M.; formal analysis, F.C., R.A. and N.M.; writing—original draft preparation, F.C.; writing—review and editing, F.C., R.A. and N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work is 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.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agrawal, O.P. A general formulation and solution scheme for fractional optimal control problems. Nonlinear Dyn. 2004, 38, 323–337. [Google Scholar] [CrossRef]
  2. Bergounioux, M.; Bourdin, L. Pontryagin maximum principle for general Caputo fractional optimal control problems with Bolza cost and terminal constraints. ESAIM Control Optim. Calc. Var. 2020, 26, 35. [Google Scholar] [CrossRef] [Green Version]
  3. Magin, R.L. Fractional calculus models of complex dynamics in biological tissues. Comput. Math. Appl. 2010, 59, 1586–1593. [Google Scholar] [CrossRef] [Green Version]
  4. Odzijewicz, T.; Malinowska, A.B.; Torres, D.F.M. Fractional variational calculus with classical and combined Caputo derivatives. Nonlinear Anal. 2012, 75, 1507–1515. [Google Scholar] [CrossRef] [Green Version]
  5. Oldham, K.B. Fractional differential equations in electrochemistry. Adv. Eng. Softw. 2010, 41, 9–12. [Google Scholar] [CrossRef]
  6. Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E 1997, 55, 3581–3592. [Google Scholar] [CrossRef]
  7. Podlubny, I. Fractional Differential Equations, Mathematics in Science and Engineering; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  8. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives, Translated from the 1987 Russian Original; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  9. Cruz, F.; Almeida, R.; Martins, N. Optimality conditions for variational problems involving distributed-order fractional derivatives with arbitrary kernels. Aims Math. 2021, 6, 5351–5369. [Google Scholar] [CrossRef]
  10. Herglotz, G. Berührungstransformationen; Lectures at the University of Göttingen: Göttingen, Germany, 1930. [Google Scholar]
  11. Georgieva, B.; Guenther, R. First Noether-type theorem for the generalized variational principle of Herglotz. Topol. Methods Nonlinear Anal. 2002, 20, 261–273. [Google Scholar] [CrossRef] [Green Version]
  12. Almeida, R.; Malinowska, A.B. Fractional variational principle of Herglotz. Discrete Contin. Dyn. Syst. Ser. B 2014, 19, 2367–2381. [Google Scholar] [CrossRef]
  13. Almeida, R.; Martins, N. Fractional variational principle of Herglotz for a new class of problems with dependence on the boundaries and a real parameter. J. Math. Phys. 2020, 61, 102701. [Google Scholar] [CrossRef]
  14. Georgieva, B.; Guenther, R. Second Noether-type theorem for the generalized variational principle of Herglotz. Topol. Methods Nonlinear Anal. 2005, 26, 307–314. [Google Scholar] [CrossRef] [Green Version]
  15. Santos, S.P.S.; Martins, N.; Torres, D.F.M. Variational problems of Herglotz type with time delay: DuBois-Reymond condition and Noether’s first theorem. Discrete Contin. Dyn. Syst. 2015, 35, 4593–4610. [Google Scholar] [CrossRef]
  16. Zhang, J.; Yin, L.; Zhou, C. Fractional Herglotz variational problems with Atangana-Baleanu fractional derivatives. J. Inequalities Appl. 2018, 44, 1–16. [Google Scholar] [CrossRef] [Green Version]
  17. Riewe, F. Nonconservative Lagrangian and Hamiltonian mechanics. Phys. Rev. E 1996, 53, 1890–1899. [Google Scholar] [CrossRef]
  18. 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]
  19. Almeida, R.; Morgado, M.L. The Euler–Lagrange and Legendre equations for functionals involving distributed-order fractional derivatives. Appl. Math. Comput. 2018, 331, 394–403. [Google Scholar] [CrossRef]
  20. Malinowska, A.B.; Torres, D.F.M. Introduction to the Fractional Calculus of Variations; World Scientific Publishing Company: London, UK, 2012. [Google Scholar]
  21. Malinowska, A.B.; Torres, D.F.M. Generalized natural boundary conditions for fractional variational problems in terms of the Caputo derivative. Comput. Math. Appl. 2010, 59, 3110–3116. [Google Scholar] [CrossRef] [Green Version]
  22. 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]
  23. Muslih, S.I.; Baleanu, D. Quantization of classical fields with fractional derivatives. Nuovo Cimento Soc. Ital. Fis. B 2005, 120, 507–512. [Google Scholar]
  24. Van Brunt, B. The Calculus of Variations, Universitext; Springer: New York, NY, USA, 2004. [Google Scholar]
  25. Cruz, F.; Almeida, R.; Martins, N. Variational Problems with Time Delay and Higher-order Distributed-order Fractional Derivatives with Arbitrary Kernels. Mathematics 2021, 9, 1665. [Google Scholar] [CrossRef]
  26. 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]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cruz, F.; Almeida, R.; Martins, N. Herglotz Variational Problems Involving Distributed-Order Fractional Derivatives with Arbitrary Smooth Kernels. Fractal Fract. 2022, 6, 731. https://doi.org/10.3390/fractalfract6120731

AMA Style

Cruz F, Almeida R, Martins N. Herglotz Variational Problems Involving Distributed-Order Fractional Derivatives with Arbitrary Smooth Kernels. Fractal and Fractional. 2022; 6(12):731. https://doi.org/10.3390/fractalfract6120731

Chicago/Turabian Style

Cruz, Fátima, Ricardo Almeida, and Natália Martins. 2022. "Herglotz Variational Problems Involving Distributed-Order Fractional Derivatives with Arbitrary Smooth Kernels" Fractal and Fractional 6, no. 12: 731. https://doi.org/10.3390/fractalfract6120731

APA Style

Cruz, F., Almeida, R., & Martins, N. (2022). Herglotz Variational Problems Involving Distributed-Order Fractional Derivatives with Arbitrary Smooth Kernels. Fractal and Fractional, 6(12), 731. https://doi.org/10.3390/fractalfract6120731

Article Metrics

Back to TopTop