Next Article in Journal
A Discrete-Time Fractional-Order Flocking Control Algorithm of Multi-Agent Systems
Next Article in Special Issue
Contributions to the Numerical Solutions of a Caputo Fractional Differential and Integro-Differential System
Previous Article in Journal
Bifurcation Analysis of Time-Delayed Non-Commensurate Caputo Fractional Bi-Directional Associative Memory Neural Networks Composed of Three Neurons
Previous Article in Special Issue
An Efficient Numerical Method Based on Bell Wavelets for Solving the Fractional Integro-Differential Equations with Weakly Singular Kernels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bang-Bang Property and Time-Optimal Control for Caputo Fractional Differential Systems

by
Shimaa H. Abel-Gaid
1,†,
Ahlam Hasan Qamlo
2,† and
Bahaa Gaber Mohamed
1,*,†
1
Department of Mathematics and Computer Science, Faculty of Sciences, Beni-Suef University, Beni-Suef 62511, Egypt
2
Department of Mathematics, Faculty of Sciences, Umm Al-Qura University, Makkah 21955, Saudi Arabia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2024, 8(2), 84; https://doi.org/10.3390/fractalfract8020084
Submission received: 2 December 2023 / Revised: 11 January 2024 / Accepted: 19 January 2024 / Published: 26 January 2024

Abstract

:
In this paper, by using the controllability method, a bang-bang property and a time optimal control problem for time fractional differential systems (FDS) are considered. First, we formulate our problem and prove the existence theorem. We then state and prove the bang-bang theorem. Finally, we state the optimality conditions that characterize the optimal control. Some application examples are given to illustrate our results.

1. Introduction

Fractional calculus is one of the most novel types of calculus used in many different scientific and technical disciplines. The order of derivatives in fractional calculus can be any real number, which distinguishes fractional calculus from ordinary calculus. Therefore, fractional calculus can be considered a generalization of ordinary calculus. Fractional order differential equations have been successfully used in modeling many different problems for various applications of fractional order derivatives and integrals in classical mechanics, quantum mechanics, image processing, earthquake engineering, biomedical engineering, physics, nuclear physics, hadron spectroscopy, viscoelasticity, bioengineering, and many other fields ([1,2,3,4,5,6,7]).
In recent decades, fractional optimal control theory for partial differential equations have found a wide range of applications in science, engineering, economics, and some other fields. For example in [8], Agrawal provided a general formulation of FOCP in the Riemann–Liouville (RL) sense and described a solution scheme for FOCP for the classical optimal control problem that was based on variational virtual work coupled with the Lagrange multiplier technique. In ([9,10]), Agrawal formulated FOCP in Caputo’s sense (CFDs) instead of the RL sense and an iterative numerical scheme was applied to find the solution of that problem numerically where the time domain was divided into small segments. CFDs permit one to include the usual initial conditions in a simple manner, and consequently they are favored choices by most researchers. The fractional derivatives of the system were approximated in [11,12] , using the Grunwald–Letnikov definition that yielded a set of algebraic equations that can be solved by numerical techniques. In [13], Gastao and Torres studied the fractional optimal control in the sense of Caputo and the fractional Noether’s Theorem. In [14], Jajarmi and Baleanu dicussed the suboptimal control of fractional-order dynamic systems with delay argument.
The time-optimal control problem for non-fractional parabolic equations with control constraints and an infinite number of variables is considered in Bahaa [15]. The fractional optimal control problem for variational inequalities with control constraints and for differential systems with control constraints for infinite order systems with control constraints are studied, respectively, in Bahaa ([16,17,18]). In Bahaa ([19,20]), fractional optimal control problems for differential systems with control constraints with delay argument, and with variable fractional order, respectively, are considered. In Bahaa [21], the optimality conditions for fractional differential inclusions with non-singular Mittag–Leffler Kernel are proved. The optimality conditions of fractional diffusion equations with weak Caputo derivatives and variational formulation are formulated in Bahaa and Tang [22]. The time-optimal control of second-order and infinite-order distributed parabolic non-fractional systems involving multiple time-varying lags are investigated in Bahaa and Kotarski [23]. In Baleanu [24], a different solution scheme was suggested where an improved Grunwald–Letnikov definition was utilized to deduce a central difference formula. In [25], Hasan et al. study the fractional optimal control of distributed systems in cylindrical and spherical coordinates. A numerical scheme and formulation for fractional optimal control of cylindrical structures subjected to general initial conditions are given in [26]. In [27], Tricaud et al. studied the time-optimal control of systems with fractional dynamics.
The bang-bang property is of great importance in the theory of optimal control, as pointed out in [28,29]. In particular, the bang-bang property can be given for certain time-optimal controls governed by parabolic equations using Pontryagin’s maximum principle [30]. In [31], Phung et al. studied the bang-bang property for time-optimal controls governed by a semilinear heat equation in a rich domain with local control in a subset. In [32], Chen et al. studied the time-varying bang-bang property of time-optimal controls for the heat equation and showed the applicability of this property.
In this paper, a Bang-Bang property and a time-optimal control problem for time fractional differential systems (FDS) are considered. First, we formulate our problem and prove the existence theorem. Using the Fre’chet differentiability, we then state and prove the bang-bang theorem. Finally, we state the optimality conditions that characterize the optimal control. Some application examples are given to illustrate our results.
This paper is organized as follows: In the second section, we introduce some fractional operators and basic definitions that we use in this paper. In the third section, we summarize the time-optimal control problem for fractional systems and then introduce the main results of this work. In the fourth section, we state and prove the existence theorem. In the fifth section, we explain and prove the Bang-Bang theorem. In the sixth section, we prove the optimality conditions. In the seventh section, we present some applications and illustrate examples of this problem. Finally, in the eighth section, we draw conclusions.

2. Basic Definitions

In this section, we introduce some basic definitions related to fractional derivatives (see [33,34]).
Definition 1.
The left Riemann–Liouville fractional integral and the right Riemann–Liouville fractional integral are presented, respectively, by
I a a f ( t ) = 1 Γ ( α ) a t ( t τ ) α 1 f ( τ ) d τ ,
I b α f ( t ) = 1 Γ ( α ) t b ( τ t ) α 1 f ( τ ) d τ ,
where α > 0 , n 1 < α < n . From now on, Γ ( α ) represents the Gamma function.
The left Riemann–Liouville fractional derivative is given by
D a a f ( t ) = 1 Γ ( n α ) d d t n a t ( t τ ) n α 1 f ( τ ) d τ .
The right Riemann–Liouville fractional derivative is defined by
D b α f ( t ) = 1 Γ ( n α ) d d t n t b ( τ t ) n α 1 f ( τ ) d τ .
The fractional derivative of a constant takes the form
D a a C = C ( t a ) α Γ ( 1 α ) ,
and the fractional derivative of a power of t has the following form
D a a ( t a ) β = Γ ( α + 1 ) ( t a ) β α Γ ( β α + 1 ) ,
for β > 1 , α 0 .
The Caputo’s fractional derivatives are defined as follows:
The left Caputo fractional derivative
D a a C f ( t ) = 1 Γ ( n α ) a t ( t τ ) n α 1 d d τ n f ( τ ) d τ ,
and the right Caputo fractional derivative
D b a C f ( t ) = 1 Γ ( n α ) t b ( τ t ) n α 1 d d τ n f ( τ ) d τ ,
where α represents the order of the derivative such that n 1 < α < n . By definition, the Caputo fractional derivative of a constant is zero. For more properties and information about these derivatives, see [35,36].
Lemma 1
(see [35,36]). Let 0 < α < 1 . Then, for any ϕ C ( Q ¯ ) , we have
0 T Ω ( a C D α y ( x , t ) + A y ( x , t ) ) ϕ ( x , t ) d x d t = Ω ϕ ( x , T ) a I 1 α y ( x , T ) d x Ω ϕ ( x , 0 ) a I 1 α y ( x , 0 + ) d x
+ 0 T Ω y ϕ ν d σ d t 0 T Ω y ν ϕ d σ d t + 0 T Ω y ( x , t ) ( R D b α ϕ ( x , t ) + A ϕ ( x , t ) ) d x d t .
where A is conjugate of the operator A , which is given in the next section and
y ν A = i , j = 1 n a i j y x j c o s ( n , x j ) on Ω ,
c o s ( n , x j ) is the i-th direction cosine of n , n being the normal at Ω exterior to Ω.

3. Time-Optimal Control Problem for Caputo Fractional Differential System

Let us consider the optimization problem in the following form
D a 0 C y ( t ; v ) + A ( t ) y ( t ; v ) = f + B v , x Ω , t ( 0 , T ) ,
y ( 0 ; v ) = y 0 ,
where f , y 0 are given two elements in L 2 ( Ω ) . The second-order operator A ( t ) in the state Equation (9) takes the form
A ( t ) y ( x , t ) = i , j = 1 n x i ( a i j ( x , t ) y ( x , t ) x j ) + a 0 ( x , t ) y ( x , t ) ,
where a i j ( x , t ) , i , j = 1 , , n , be given function on Ω with the properties
a 0 ( x , t ) , a i j ( x , t ) L ( Ω ) ( with real values ) ,
a 0 ( x , t ) δ > 0 , i , j = 1 n a i j ( x , t ) ξ i ξ j δ ( ξ 1 2 + + ξ n 2 ) , ξ R n ,
a.e. on Ω , i.e., A ( t ) is a bounded second-order self-adjoint elliptic partial differential operator, which maps H 0 1 ( Ω ) onto H 1 ( Ω ) .
Let us denote by U : = L 2 ( 0 , T ; L 2 ( Ω ) ) = L 2 ( Q ) , the space of controls and by Y : = L 2 ( 0 , T ; H 0 1 ( Ω ) ) the space of states.
We suppose that U a d is a closed, convex subset of U and y 1 is a given element in L 2 ( Ω ) ,
B L ( L 2 ( Q ) , L 2 ( 0 , T , H 1 ( Ω ) ) ) .
Assume (Controllability)
there exists a v U a d , such that y ( τ ; v ) = y 1 for an appropriate τ .
The optimal time is given by
τ 0 = inf τ , τ such that ( 13 ) holds .
For the operator A ( t ) , we define the bilinear form π ( t ; y , ɸ ) as follows:
Definition 2.
On H 0 1 ( Ω ) , we define for each t ] 0 , t [ the following bilinear form
π ( t ; y , ɸ ) = ( A ( t ) y ɸ ) , L 2 ( Ω ) , y , ɸ H 0 1 ( Ω ) .
Then
π ( t ; y , ϕ ) = A y , ϕ L 2 ( Ω ) = i , j = 1 n x i a i j ( x ) y x j + a 0 ( x ) y , ϕ ( x ) L 2 ( Ω ) = Ω i , j = 1 n a i j x i y ( x ) x j ϕ ( x ) d x + Ω a 0 ( x ) y ( x ) ϕ ( x ) d x .
Lemma 2.
The bilinear form (15) is coercive on H 0 1 ( Ω ) ; that is,
π ( t ; y , y ) λ | | y | | H 0 1 { a α , 2 } ( Ω ) 2 , λ > 0 .
Proof. 
π ( t ; y , y ) = Ω i , j = 1 n a i j ( x , t ) x i y ( x , t ) x j y ( x , t ) d x + Ω a 0 ( x , t ) y ( x , t ) y ( x , t ) d x i , j = 1 n a i j ( x , t ) | | x i y ( x , t ) | | L 2 ( Ω ) 2 + | | y ( x , t ) | | L 2 ( Ω ) 2 λ | | y | | H 0 1 ( Ω ) 2 , λ > 0 .
Also, we can assume that
y , ɸ H 0 1 ( Ω ) the bilinear form t π ( t ; y , ɸ ) is continuously differentiable in [ 0 , T ] , the function t π ( t ; y , ɸ ) is measurable on [ 0 , T ] , | π ( t ; y , ɸ ) | c | | y | | . | | ɸ | | ,
the bilinear form (15) is symmetric,
π ( t ; y , ɸ ) = π ( t ; ɸ , y ) y , ɸ H 0 1 ( Ω ) ,
and there exists a λ such that
π ( t ; ɸ , ɸ ) + λ | ɸ | 2 γ | | ɸ | | 2 , γ > 0 , ɸ H 0 1 ( Ω ) , t [ 0 , T ] .
Equations (9)–(18) constitute a fractional Dirichlet problem.
Remark 1
([35,36]). The operator D t a 0 C + A ( t ) is a second-order parabolic operator that maps L 2 ( 0 , T ; H 0 1 ( Ω ) ) onto L 2 ( 0 , T ; H 1 ( Ω ) ) .
Remark 2
([35,36]). Equations (9)–(18) have the unique generalized solution
y W ( 0 , T ) : = y | y L 2 ( 0 , T ; H 0 1 ( Ω ) ) , D t a 0 C y L 2 ( 0 , T ; H 1 ( Ω ) ) ,
which continuously depends on the initial condition (10) and the right-hand side of (9). Furthermore, y W ( 0 , T ) is a continuous function [ 0 , T ] L 2 ( Ω ) (compare with Theorems 1.1 and 1.2 Chapt. 3 [29]).
We are going to study the following problems:
(i) the existence of optimal control, i.e., u U a d such that
y ( τ 0 ; u ) = y 1 ;
(ii) properties of the optimal control, if it exists. These problems will be treated in Section 4 and Section 5.

4. Existence Theorem

Theorem 1.
Let α [ 1 2 , 1 ] . We assume that (12), (13) and (17) hold and that U a d is bounded. Then, there exists an optimal control, that is u U a d , such that (19) holds.
Proof. 
Let τ n be such that
y ( τ n ; v n ) = y 1 , v n U a d ,
τ n τ 0 .
Set y n = y ( v n ) . Sin ce U a d is bounded , we may verify that y n ( resp . D a 0 C y n ) ranges in a bounded set in L 2 ( 0 , T , H 0 1 ( Ω ) ) ( resp . L 2 ( 0 , T , H 1 ( Ω ) ) ) .
We may then extract a subsequence, again denoted by { v n , y n } , such that
v n v weakly in U , v U a d y n y ( resp D a 0 C y n D a 0 C y ) weakly in L 2 ( 0 , T , H 0 1 ( Ω ) ) ( resp . L 2 ( 0 , T , H 1 ( Ω ) ) ) .
We deduce from the equality
D a 0 C y n + A y n = f + v n , x Ω , t ( 0 , T ) ,
that
D a 0 C y + A y = f + v , x Ω , t ( 0 , T ) ,
and y ( 0 ) = y 0 and hence
y = y ( v ) .
But
y ( τ n ; v n ) y ( τ 0 ; v ) = y ( τ n ; v n ) y ( τ 0 ; v n ) + y ( τ 0 ; v n ) y ( τ 0 ; v ) .
Now, from (23), y ( τ 0 ; v n ) y ( τ 0 ; v ) weakly in L 2 ( 0 , T , H 0 1 ( Ω ) ) and
| | y ( τ n ; v n ) y ( τ 0 ; v n ) | | H 1 ( Ω ) = | | I a τ 0 D a τ 0 C y ( τ n ; v n ) | | H 1 ( Ω ) 1 Γ ( α ) τ 0 τ n ( τ n τ ) α 1 | | D a τ 0 C y ( τ ; v n ) | | H 1 ( Ω ) d τ 1 Γ ( α ) ( τ 0 τ n ( τ n τ ) 2 α 2 d τ ) 1 2 ( τ 0 τ n | | D a τ 0 C y ( τ ; v n ) | | H 1 ( Ω ) 2 d τ ) 1 2 C ( τ n τ 0 ) α 1 2 ;
hence, (25) shows that
y ( τ n ; v n ) y ( τ 0 ; v ) 0 weakly in H 1 ( Ω ) .
But from (20), it follows that y ( τ 0 ; v ) = y 1 . □
Example 1
(Neumann problem with boundary control). Take U = L 2 ( Σ ) . Let the state y ( v ) be as follows
D a 0 C y ( v ) + A ( t ) y ( v ) = f , y ν A ( v ) = v , y ( x , 0 , v ) = y 0 ( x ) .
We may apply Theorem 1 to this example. Thus, the theorem deals with the case of boundary control.
Remark 3.
Theorem 1 may be modified easily to deal with the case of systems with Dirichlet boundary conditions (cf ([29]), Section 9) and where the control is exercised through the boundary.

5. Bang-Bang Theorem

We consider the same data and hypotheses as in the previous section with
A independent of t ,
U a d = { v | | v ( t ) | 1 a . e . } .
Then A is the infinitesimal generator of a semi-group G ( t ) in L 2 ( Q ) .
Theorem 2
(Bang-Bang Theorem). We assume that (12), (17), (26), (27) and (13) hold. Let u be an optimal control, that is, an element of U a d satisfying (19) (we know from Theorem 1 that such elements exist). Then
| u ( t ) | = 1 a . e . on [ 0 , τ 0 ] .
We shall establish some lemmas before giving a proof of the Theorem 2. Let us at once isolate.
Corollary 1.
Under the hypotheses of Theorem 2, there exists a unique optimal control.
Proof. 
Let u 1 , u 2 U a d ; y ( τ 0 ; u i ) = y i , i = 1 , 2 Then
y ( τ 0 ; ( 1 θ ) u 1 + θ u 2 ) = y 1 , 0 < θ < 1 and u = ( 1 θ ) u 1 + θ u 2
does not satisfy (19) unless u 1 ( t ) = u 2 ( t ) a . e .
Notation
K s = { h | h = G v ( s ) = 0 s G ( s t ) v ( t ) d t , v L ( 0 , s , L 2 ( Ω ) ) } K s ( e ) = { h | h = G v ( s ) = 0 s G ( s t ) v ( t ) d t , v L ( 0 , s , L 2 ( Ω ) ) with support in e [ 0 , s ] , e = measurable subset of [ 0 , T ] } .
Lemma 3.
With the notation of (29), we have:
K s = K s , s , s .
We then set
K = K s , s .
We have
G ( s ) H K .
Proof. 
  • Let τ > s . Then we may easily verify
    G v ( s ) = G v ¯ ( τ ) , where v ¯ ( τ ) = 0 if 0 < t < τ s , v ( t ( τ s ) ) if t > τ s .
    Consequently
    K s K τ .
  • Let τ < s . We may again verify that
    G v ( s ) = G w ( s ) , where w ( s ) = v ( t + τ s ) ) + 1 s G ( t ) 0 τ s G ( τ s t ) v ( t ) d t ,
    which proves the reverse inclusion of (34), whence (30) and (31).
  • We check that
    G ( s ) h = G × ( 1 s G h ) ( s ) , whence ( 32 ) .
The following lemma is fundamental.
Lemma 4.
For almost all t e , we have
K t ( e ) = K .
Proof. 
1. Clearly, K t ( e ) K t and, hence, it suffices to prove that for almost all t e , K K t ( e ) . Let h K t . Then h K s with s arbitrary small and therefore for any t 1 < t , we may represent h in the form (cf. (33))
h = t 1 t G ( t σ ) v ( σ ) d σ , v L ( 0 , T ; H ) .
Now, suppose that we can find a sequence t n such that
t n < t n + 1 < < t , t n t , measure ( [ t n , t n + 1 ] e ) ρ ( t n + 1 t n ) , ρ > 0 , t n + 1 t n t n + 2 t n + 1 c .
We shall see in part 2 of the proof that
for almost all t e , we may find a sequence t n such that ( 39 ) holds .
In (38), let us choose as t 1 the first element of the sequence { t n } . From (38), we deduce that we may write
h = n = 1 G ( t t n + 1 ) t n t n + 1 G ( t n + 1 σ ) v ( σ ) d σ = n = 1 e [ t n + 1 , t n + 2 ] 1 measure ( e [ t n + 1 , t n + 2 ] ) G ( t σ 1 ) × G ( σ 1 t n + 1 ) d σ 1 t n t n + 1 G ( t n + 1 σ 2 ) v ( σ 2 ) d σ 2 = 0 t G ( t σ ) w ( σ ) d σ ,
where
w ( s ) = 1 measure ( e [ t n + 1 , t n + 2 ] ) G ( σ t n + 1 ) t n t n + 1 G ( t n + 1 σ 2 ) v ( σ 2 ) d σ 2 in the set e [ t n + 1 , t n + 2 ] , n = 1 , 2 , ; 0 otherwise .
From (39), we have
| w ( σ ) | 1 ρ ( t n + 2 t n + 1 ) C ( t n + 1 t n ) constant .
Hence,
h = G w ( t ) , w L ( 0 , T , H ) , w has support in e , and h K t ( e ) .
Then the lemma is proved except (40).
2. Proof of (40). This is a result of the measure theory. We first define
e m = σ | σ e , m e a s u r e ( e [ σ 1 k , σ ] ) 1 2 k , k m
and d m = set of points of density of e m . It is known that m e a s u r e e m 1 d m = 0 and, hence, it suffices to prove (40) for t d m .
However, then we may construct { t n } , t n d m n , t n + 1 = t n + s n < t , s n > 0 , s n s n + 1 e (in a manner such that the last property of (39) holds) and since t n d m e m n , there exists a ρ > 0 such that
measure ( [ t n , t n + 1 ] e ) ρ ( t n + 1 t n ) .
Lemma 5.
Let u U a d be optimal with respect to { y 0 , y 1 } , that is, y ( τ ; u ) = y 1 with τ minimum. Then for any s < τ , u is optimal with respect to { y 0 , y ( s ; u ) } .
In other words, if w U a d satisfies
y ( σ ; w ) = y ( s ; u ) for some appropriate σ ,
we necessarily have σ s .
Proof. 
Assume that (43) holds with σ < s . Then define v as.
v ( t ) = w ( t ) i n [ 0 , σ ] , u ( t + s σ ) i n [ σ , T ( s σ ) ] .
Then, from (43),
y ( σ ; v ) = y ( σ ; w ) = y ( s ; u )
and
D a 0 C y ( t ; v ) + A y ( t ; v ) = u ( t + s σ ) , for t σ .
Hence,
y ( t ( s σ ) ; v ) = y ( t , u ) , t s ,
and hence,
y ( τ ( s σ ) ; v ) = y ( τ , u ) = y .
But since τ ( s σ ) < τ , this contradicts the hypothesis that u is optimal with respect to { y 0 , y 1 } . □
Lemma 6.
Let v U a d be such that
| v ( t ) | 1 ε almost everywhere , ε > 0 ,
y ( τ ; v ) = y 1 for some appropriate τ .
Then τ > τ 0 .
Proof. 
To prove this, we will verify that there exists a s < τ and w U a d such that
y ( s ; w ) = y ( τ ; v )
(which proves that τ is not optimal). For this, we note that (46) may be written as
0 τ G ( τ σ ) v ( σ ) d ( σ ) + G ( τ ) y 0 G ( s ) y 0 = 0 s G ( s σ ) w ( σ ) d σ .
But it may be easily verified that the left-hand side of (47) may be written as
0 s G ( s σ ) v ¯ ( σ ) d σ ,
with
v ¯ ( σ ) = v ( σ + τ s ) + 1 s G ( σ ) 0 τ s G ( τ s σ 1 ) v ( σ 1 ) d σ 1 + 1 s [ G ( σ + τ s ) y 0 G ( σ ) y 0 ]
and for τ s sufficiently small, we have, by virtue of (44), | v ¯ ( σ ) | 1 . Therefore, we may take w = v ¯ . □
Proof of Bang-Bang Theorem (Theorem 2). 
Assume that (28) did not hold. Then there would exist e [ 0 , τ 0 ] , measure ( e ) > 0 such that
| u ( t ) | 1 ε , ε > 0 , t e .
It follows from the fundamental Lemma (4) that we can find an s such that K s ( e ) = K = K s . In other words, there exists a g ¯ L ( 0 , T , H ) , with support in e, such that
0 s G ( s σ ) g ¯ ( σ ) d σ = 0 s G ( s σ ) u ( σ ) d σ .
Let us introduce the control
v = ( 1 δ ) u + δ g ¯ i n ( 0 , τ 0 ) , δ > 0 almost everywhere , 0 for t > τ 0 .
We shall choose δ in a manner such that
| v ( t ) | 1 ε 1 almost everywhere , ε 1 > 0 appropriate .
This is possible. To see this, we note that on e, we have from (48)
| v ( t ) | ( 1 δ ) ( 1 ε ) + δ | g ¯ ( t ) | 1 ε 2
for δ sufficiently small, and outside e,
| v ( t ) | = ( 1 δ ) | u ( t ) | ( 1 δ ) ( s i n c e u U a d ) .
But
y ( s ; v ) = ( 1 δ ) 0 s G ( s σ ) u ( σ ) d σ + δ 0 s G ( s σ ) g ¯ ( σ ) d σ + G ( s ) y 0
and using (49), we have
y ( s ; v ) = y ( s ; u ) .
However, then, from Lemma (6), (since we have (51) u is not optimal with respect to { y 0 , y ( s ; u ) } and, hence, from Lemma (5), u is not optimal with respect to { y 0 , y 1 } , contradicting the hypotheses. □
Remark 4.
It is possible to proceed further using much simpler arguments when the semi-group G is a group—in other words when we can reverse time. Indeed we have the following result.
Theorem 3.
We assume that A is the infinitesimal generator of a group G ( t ) . We assume that
U a d = { v | v ( t ) H a d H , H a d = neighborhood of 0 i n H } .
We further suppose that there exists a v U a d such that (13) holds and that there exists an optimal control u (that is, y ( τ 0 ; u ) = y 1 , τ 0 = optimal time defined in (12).
Then
u ( t ) H a d ( boundary of H a d ) , almost everywhere .
Proof. 
If (54) does not hold, there exists an e [ 0 , τ 0 ] such that
measure ( e ) > 0 , distance ( u ( t ) , H a d ) c 1 > 0 , almost everywhere for t e .
Let χ e be the characteristic function of e. Then
y ( τ 0 ; u ) = y 1 = G ( τ 0 ) y 0 + 0 τ 0 G ( τ 0 σ ) u ( σ ) d σ
may be written as—by virtue of the fact that G is a group—
y ( τ 0 ; u ) = G ( τ 1 ) y 0 + 0 τ 1 G ( τ 1 σ ) u ¯ ( σ ) d σ ,
with
u ¯ ( σ ) = u ( σ ) + 1 m e a s u r e ( e ) χ e ( σ ) G ( σ τ 1 ) [ y ( τ 0 ; u ) y ( τ 1 ; u ) ] , τ 1 < τ 2 .
We may choose τ 1 sufficiently near to τ 0 such that
| u ( σ ) u ¯ ( σ ) | c 1 2 for σ e ,
and, therefore, from (55), u U a d . Then, from (56),
y ( τ 1 ; u ¯ ) = y ( τ 0 ; u ) = y 1 ,
which contradicts the fact that τ 0 is optimal. □

6. Optimality Conditions

Theorem 4.
Assume that the hypotheses of Theorem 3 hold. We further assume that H a d is convex. Then there exists an h 0 H such that
( h 0 , y ( τ 0 ; v ) y ( τ 0 ; u ) ) 0 v U a d .
Proof. 
1. Define the set - which is clearly convex:
K = { h | h H , h = y ( τ 0 ; v ) y ( τ 0 ; 0 ) ) , v U a d } .
Let us verify that
0 is in the i n t e r i o r o f K ,
y 1 y ( τ 0 ; 0 ) ( which belongs to K ) is not in the interior of K .
To prove (60), let h H be arbitrary. We have to prove that for sufficiently small values of | ε | , ε h K . Now,
h = 0 τ 0 G ( τ 0 σ ) G ( σ τ 0 ) h d σ ,
and therefore,
ε h = y ( τ 0 ; v G ( σ τ 0 ) h ) ) y ( τ 0 ; 0 ) K f o r | ε | sufficiently small .
We prove (61) by contradiction. If y 1 y ( τ 0 ; 0 ) were in the interior of K, for an appropriate ε > 0 , we would have
y 1 y ( τ 0 ; 0 ) + ε ( y 1 G ( τ 0 ) y 0 ) K
and there would exist a v U a d such that
y 1 y ( τ 0 ; 0 ) + ε ( y 1 G ( τ 0 ) y 0 ) = y ( τ 0 ; v ) y ( τ 0 ; 0 ) ,
whence
y 1 = y ( τ 0 ; 0 ) + 0 τ 0 G ( τ 0 σ ) 1 1 + ε v ( σ ) d ( σ ) = y ( τ 0 ; w )
and therefore, w ( σ ) interior of U a d , contradicting (54).
2. The result now follows as a consequence of (60), (61) and of Corollary 5, Section 9 of Dunford–Schwartz [37]. □
Remark 5
(58). We may interpret this as follows: let us introduce the adjoint state by
D τ 0 a R p + A p = 0 , in [ 0 , τ 0 ] ,
I τ 0 1 α p ( τ 0 ) = h 0 .
Then, by using the Green formula given in Lemma 1, (58) is equivalent to
0 τ 0 ( p ( t ) , v ( t ) u ( t ) ) d t 0 v U a d ,
whence, we may pass to local conditions in t:
( p ( t ) , h u ( t ) ) 0 almost everywhere in [ 0 , τ 0 ] h H a d .

7. Application

In this section, we state some examples to explain our abstract conclusions.
Example 2.
If H a d = unit ball in H, we would have
u ( t ) = p ( t ) | p ( t ) | .
Example 3.
Let us consider the system with the following state
D a 0 R y ( t ; v ) + A ( t ) y ( t ; v ) = f + B v
I 1 a 0 y ( 0 , v ) = y 0
and the adjoint state is given by
D τ 0 a C p + A p = 0 , in [ 0 , τ 0 ] ,
p ( τ 0 ) = h 0 .
Then (58) is equivalent to
0 τ 0 ( p ( t ) , v ( t ) u ( t ) ) d t 0 v U a d .
Remark 6.
If we take α = 1 in the previous sections, we obtain the classical results in the optimal control with integer derivatives.

8. Open Problems

1—In a similar manner, we can also study the time-fractional optimal control of the above systems, where the time derivative is considered as the left Atangana–Baleanu fractional derivative in the Caputo sense:
D a 0 A B C y ( t ; v ) + A ( t ) y ( t ; v ) = f + B v , x Ω , t ( 0 , T ) ,
y ( 0 ; v ) = y 0 ,
where D a 0 A B C is the left Atangana–Baleanu fractional derivative in the sense of Caputo (see [1]).
2—The problem can be extended to the time-space fractional derivative as the following:
D a 0 A B C y ( t ; v ) + ( Δ ) β ( t ) y ( t ; v ) = f + B v , x Ω , t ( 0 , T ) ,
y ( 0 ; v ) = y 0 ,
where ( Δ ) β is is the fractional Laplacian operator for β ( 0 , 1 ) (see [35,36,38]).

9. Conclusions

In this work, we stated and proved the bang-bang theorem (Theorem 2) for (FDS) with Dirichlet and Neumann boundary conditions. The fractional derivatives were defined in the weak Caputo and Riemann–Liouville sense. The analytical results were obtained in the form of Euler–Lagrange equations for the (TOCFP). The formulation presented and the resulting equations are similar to those for classical optimal control problems. The optimization problem presented in this paper represents a generalization of the time-optimal control problems of parabolic systems with Dirichlet and Neumann boundary conditions considered in Lion’s [29] to fractional time-optimal control problems. Moreover, the main result of the paper contains necessary optimality conditions for fractional systems of non-integer order, which allows for the characterization of optimality problems.

Author Contributions

These authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Abdeljawad, T.; Baleanu, D. Integration by parts and its applications of a new nonlocal fractional derivative with Mittag-Leffler nonsingular kernel. J. Nonlinear Sci. Appl. 2017, 10, 1098–1107. [Google Scholar] [CrossRef]
  2. Agrawal, O.P. Formulation of Euler-Lagrange equations for fractional variational problems. J. Math. Anal. Appl. 2002, 272, 368–379. [Google Scholar] [CrossRef]
  3. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; A Wiley- Interscience Pub.; John Wiley & Sons, Inc.: New York, NY, USA, 1993. [Google Scholar]
  4. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  6. Riewe, F. Nonconservative Lagrangian and Hamiltonian mechanics. Phys. Rev. E 1996, 53, 1890–1899. [Google Scholar] [CrossRef] [PubMed]
  7. Riewe, F. Mechanics with fractional derivatives. Phys. Rev. E 1997, 55, 81–92. [Google Scholar] [CrossRef]
  8. Agrawal, O.P. A General formulation and solution scheme for fractional optimal control problems. Nonlinear Dyn. 2004, 38, 323–337. [Google Scholar] [CrossRef]
  9. Agrawal, O.P. Fractional optimal control of a distributed system using eigenfunctions. J. Comput. Nonlinear Dyn. 2008, 3, 1–6. [Google Scholar] [CrossRef]
  10. Baleanu, D.; Agrawal, O.P. Fractional Hamilton formalism within Caputo’s derivative. Czechoslov. J. Phys. 2006, 56, 1087–1092. [Google Scholar] [CrossRef]
  11. Agrawal, O.P.; Baleanu, D. A Hamiltonian formulation and direct numerical scheme for fractional optimal control problems. J. Vib. Control 2007, 13, 1269–1281. [Google Scholar] [CrossRef]
  12. Agrawal, O.P.; Defterli, O.; Baleanu, D. Fractional optimal control problems with several state and control variables. J. Vib. Control 2010, 16, 1967–1976. [Google Scholar] [CrossRef]
  13. Gastao, F.F.; Torres Delfim, F.M. Fractional optimal control in the sense of Caputo and the fractional Noether’s Theorem. Int. Math. 2008, 3, 479–493. [Google Scholar]
  14. Jajarmi, A.; Baleanu, D. Suboptimal control of fractional-order dynamic systems with delay argument. J. Vib. Control 2017, 23, 1–17. [Google Scholar] [CrossRef]
  15. Bahaa, G.M. Time-optimal control problem for parabolic equations with control constraints and infinite number of variables. IMA J. Math. Control Inform. 2005, 22, 364–375. [Google Scholar] [CrossRef]
  16. Bahaa, G.M. Fractional optimal control problem for variational inequalities with control constraints. IMA J. Math. Control Inform. 2018, 35, 107–122. [Google Scholar] [CrossRef]
  17. Bahaa, G.M. Fractional optimal control problem for differential system with control constraints. Filomat 2016, 30, 2177–2189. [Google Scholar] [CrossRef]
  18. Bahaa, G.M. Fractional optimal control problem for infinite order system with control constraints. Adv. Differ. Eq. 2016, 250, 1–16. [Google Scholar] [CrossRef]
  19. Bahaa, G.M. Fractional optimal control problem for differential system with control constraints with delay argument. Adv. Differ. Eq. 2017, 69, 1–19. [Google Scholar]
  20. Bahaa, G.M. Fractional optimal control problem for variable-order differential systems. Fract. Calc. Appl. Anal. 2017, 20, 1447–1470. [Google Scholar] [CrossRef]
  21. Bahaaa, G.M.; Hamiaz, A. Optimality conditions for fractional differential inclusions with nonsingular Mittag-Leffler Kernel. Adv. Differ. 2018, 2018, 257. [Google Scholar] [CrossRef]
  22. Bahaa, G.M. Generalized variational calculus in terms of multi-parameters involving Atangana-Baleanu’s derivatives and application. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 485–501. [Google Scholar]
  23. Bahaa, G.M.; Kotarski, W. Time-optimal control of infinite-order distributed parabolic systems involving multiple time-varying lags. J. Numer. Funct. Opt. 2016, 37, 1066–1088. [Google Scholar] [CrossRef]
  24. Baleanu, D.; Defterli, O.; Agrawal, O.P. Central difference numerical scheme for fractional optimal control problems. J. Vib. Control. 2009, 15, 583–597. [Google Scholar] [CrossRef]
  25. Hasan, M.M.; Tangpong, X.W.; Agrawal, O.P. Fractional optimal control of distributed systems in spherical and cylindrical coordinates. J. Vib. Control. 2011, 18, 1506–1525. [Google Scholar] [CrossRef]
  26. Hasan, M.M.; Tangpong, X.W.; Agrawal, O.P. A Formulation and Numerical Scheme for Fractional Optimal Control of Cylindrical Structures Subjected to General Initial Conditions. In Fractional Dynamics and Control; Baleanu, D., Machado, J., Luo, A., Eds.; Springer: New York, NY, USA, 2012; pp. 3–17. [Google Scholar]
  27. Tricaud, C.; Chen, Y.Q. Time-optimal control of systems with fractional dynamics. Int. J. Differ. Eq. 2010, 2010, 461048. [Google Scholar] [CrossRef]
  28. Fattorini, H.O. Infinite Dimensional Linear Control Sytems: The Time Optimal and Norm Optimal Problems; North-Holland Math. Stud; Elsevier: Amsterdam, The Netherlands, 2005; Volume 201. [Google Scholar]
  29. Lions, J.L. Optimal Control of Systems Governed by Partial Differential Equations; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1971; Volume 170. [Google Scholar]
  30. Wang, G.; Wang, L. The bang-bang principle of time optimal controls for the heat equation with internal controls. Syst. Control Lett. 2007, 56, 709–713. [Google Scholar] [CrossRef]
  31. Phung, K.D.; Zhang, C. Bang-bang property for time-optimal control of similinear heat equation. In Annales de l’Institut Henri Poincare (C) Non Linear Analysis; Elsevier: Amsterdam, The Netherlands, 2014; Volume 31, pp. 477–499. [Google Scholar]
  32. Chen, N.; Yang, D.H. Time-varying bang-bang property of time optimal controls for the heat equation and its applications. Syst. Control Lett. 2018, 112, 18–23. [Google Scholar] [CrossRef]
  33. Jarad, F.; Maraba, T.; Baleanu, D. Fractional variational optimal control problems with delayed arguments. Nonlinear Dyn. 2010, 62, 609–614. [Google Scholar] [CrossRef]
  34. Jarad, F.; Maraba, T.; Baleanu, D. Higher order fractional variational optimal control problems with delayed arguments. Appl. Math. Comput. 2012, 218, 9234–9240. [Google Scholar] [CrossRef]
  35. Mophou, G.M. Optimal control of fractional diffusion equation. Comput. Math. Appl. 2011, 61, 68–78. [Google Scholar] [CrossRef]
  36. Mophou, G.M. Optimal control of fractional diffusion equation with state constraints. Comput. Math. Appl. 2011, 62, 1413–1426. [Google Scholar] [CrossRef]
  37. Dunford, N.; Schwartz, J.T. Linear Operators. Part I. Interscience; John Wiley & Sons: Hoboken, NJ, USA, 1958. [Google Scholar]
  38. Mophou, G.M.; Fotsing, J.M. Optimal control of a fractional diffusion equation with delay. J. Adv. Math. 2014, 6, 1017–1037. [Google Scholar]
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

Abel-Gaid, S.H.; Qamlo, A.H.; Mohamed, B.G. Bang-Bang Property and Time-Optimal Control for Caputo Fractional Differential Systems. Fractal Fract. 2024, 8, 84. https://doi.org/10.3390/fractalfract8020084

AMA Style

Abel-Gaid SH, Qamlo AH, Mohamed BG. Bang-Bang Property and Time-Optimal Control for Caputo Fractional Differential Systems. Fractal and Fractional. 2024; 8(2):84. https://doi.org/10.3390/fractalfract8020084

Chicago/Turabian Style

Abel-Gaid, Shimaa H., Ahlam Hasan Qamlo, and Bahaa Gaber Mohamed. 2024. "Bang-Bang Property and Time-Optimal Control for Caputo Fractional Differential Systems" Fractal and Fractional 8, no. 2: 84. https://doi.org/10.3390/fractalfract8020084

APA Style

Abel-Gaid, S. H., Qamlo, A. H., & Mohamed, B. G. (2024). Bang-Bang Property and Time-Optimal Control for Caputo Fractional Differential Systems. Fractal and Fractional, 8(2), 84. https://doi.org/10.3390/fractalfract8020084

Article Metrics

Back to TopTop