Next Article in Journal
Applications of Caputo-Type Fractional Derivatives for Subclasses of Bi-Univalent Functions with Bounded Boundary Rotation
Previous Article in Journal
Approximate Controllability of Hilfer Fractional Stochastic Evolution Inclusions of Order 1 < q < 2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Fractional Differential Equations with Impulsive Effects

1
Department of Mathematical Analysis and Numerical Mathematics, Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, 842 48 Bratislava, Slovakia
2
Mathematical Institute of Slovak Academy of Sciences, 814 73 Bratislava, Slovakia
3
STAR-UBB Institute, Bebes-Bolyai University, 400084 Cluj-Napoca, Romania
4
Romanian Institute of Science and Technology, 400487 Cluj-Napoca, Romania
5
Department of Electrical Engineering, City University of Hong Kong, Hong Kong SAR, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(9), 500; https://doi.org/10.3390/fractalfract8090500
Submission received: 25 June 2024 / Revised: 17 August 2024 / Accepted: 22 August 2024 / Published: 26 August 2024
(This article belongs to the Section Numerical and Computational Methods)

Abstract

:
This paper discusses impulsive effects on fractional differential equations. Two approaches are taken to obtain our results: either with fixed or changing lower limits in Caputo fractional derivatives. First, we derive an existence result for periodic solutions of fractional differential equations with periodically changing lower limits. Then, the impulsive effects are modeled for fractional differential equations regarding the nonlinearities rather than the initial value conditions. The proposed impulsive model differs from common discontinuous and nonsmooth dynamical systems.
MSC:
34A08; 34A37; 34C25

1. Introduction

Fractional differential equations (FDEs) arise from engineering and scientific disciplines regarding mathematical modeling of systems in the fields of physics, aerodynamics, chemistry, electrodynamics, and so on, involving derivatives of fractional order [1,2,3,4], where the theory of impulsive differential equations (IDEs) is particularly important for its wide spectrum of applications [5,6,7]. The main reason for its applicability is due to the fact that impulsive differential problems appropriately describe many processes where, at certain moments, the states rapidly change, and these processes cannot be well modeled by classical differential equations [8,9]. Moreover, impulsive fractional differential equations (IFDEs) have found many promising applications [10,11].
Contrarily to the classical derivative, the fractional derivative is nonlocal, which leads to obstacles in studying IFDEs [12,13]. Recently, a new approach to modelling impulsive effects has been suggested and studied [14]. This approach can be used when impulses in time are non-instantaneous; that is, impulses last over certain time intervals. The present paper does not deal with this new and interesting topic of non-instantaneous IDEs; rather, it focuses on yet another way to introduce impulses into differential equations.
Now, consider a continuous mapping f : R 0 × R m R m , where R 0 = [ 0 , ) and { t k } k = 0 R 0 is an increasing sequence so that t 0 = 0 , lim k t k = and { y k } k = 1 R m .
The objective here is to solve the following ODE:
x = f ( t , x ) , t ( t k , t k + 1 )
with impulses
x ( t k + ) = x ( t k ) + y k .
A conventional approach is to start with initial value x ( t 0 ) = x 0 at t 0 = 0 to solve (1) for t [ 0 , t 1 ) and apply the impulse (2) to obtain x ( t 1 + ) = x ( t 1 ) + y 1 , and then, we use this as the initial value to solve (1) on [ t 1 , t 2 ) , and so on. During this iteration process, one can use any numerical method to solve the problem on each interval ( t k , t k + 1 ) , k = 0 , 1 , 2 , .
It follows from (1) that
x ( t ) = x ( t k + ) + t k t f ( s , x ( s ) ) d s ,
for t ( t k , t k + 1 ) , which leads to
x ( t ) = x ( t k + ) + t k t f ( s , x ( s ) ) d s = x ( t k ) + y k + t k t f ( s , x ( s ) ) d s = x ( t k 1 + ) + t k 1 t k f ( s , x ( s ) ) d s + y k + t k t f ( s , x ( s ) ) d s = x ( t 0 ) + i = 1 k y i + i = 1 k t i 1 t i f ( s , x ( s ) ) d s + t k t f ( s , x ( s ) ) d s = x ( t 0 ) + i = 1 k y i + t 0 t f ( s , x ( s ) ) d s
for t ( t 0 , t ) .
This approach cannot be extended to FDEs, however, keeping the lower limit at t 0 . Note that
x ( t ) = x ( t 0 ) + i = 1 k y i + 1 Γ ( q ) t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s ,
for t ( t 0 , ) , is a solution to the IFDE, as presented in [13]:
D t 0 q c x ( t ) = f ( t , x ) ,
for t ( t 0 , ) , with (2) and q ( 0 , 1 ) , fixing the lower limit at t 0 , where D t 0 q c x ( t ) is a generalized Caputo fractional derivative with lower limit at t 0 .
Next, if the FDE version of (1) is considered in the following form:
D t 0 q c x ( t ) = f ( t , x ) , t ( t k , t k + 1 )
where q ( 0 , 1 ) , then (3) becomes
x ( t ) = x ( t k + ) + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s ,
for t ( t k , t k + 1 ) . Thus, one can study this type of IFDE (6) with (2) when the lower limits are changing at all impulses. Consequently, D t k q c x ( t ) is a generalized Caputo fractional derivative with lower limit at t k . Similarly, it follows from (7) that
x ( t ) = x ( t k + ) + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s   = x ( t k ) + y k + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s   = x ( t k 1 + ) + 1 Γ ( q ) t k 1 t k ( t k s ) q 1 f ( s , x ( s ) ) d s     + y k + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s   = x ( t 0 ) + i = 1 k y i + 1 Γ ( q ) i = 1 k t i 1 t i ( t i s ) q 1 f ( s , x ( s ) ) d s     + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s
for t ( t k , t k + 1 ) . Because q ( 0 , 1 ) , there is a memory effect in (6), meaning that
  i = 1 k t i 1 t i ( t i s ) q 1 f ( s , x ( s ) ) d s + t k t ( t s ) q 1 f ( s , x ( s ) ) d s t 0 t ( t s ) q 1 f ( s , x ( s ) ) d s
where it is important to note that the integral kernels of the first term on the left-hand side are ( t i s ) q 1 but not ( t s ) q 1 and that [ t 0 , t ) = i = 0 k 1 [ t i , t i + 1 ) [ t k , t ) holds only for t ( t k , t k + 1 ) . Thus, (8) shows that the memory effect of (5) is not a sum of the memory effects of (6) on each of these subintervals.
A comparative numerical study using examples of the above two IFDEs, with fixed lower limit (5)–(2) and with changing lower limit (6), can be found in [15].
One can study two types of IFDEs: the first is (5)–(2) with a solution formula (4); the second is (6)–(2) with a solution formula (7). One can see that gluing (7) from 1 to k together does not give a solution formula (4). This is an essential difference between integer-order and fractional-order differential equations. This also implies that IFDE has no nonconstant periodic solutions when the lower limit is fixed at t 0 = 0 [16]. However, (6)–(2) may have periodic solutions when the impulses are periodic. As a matter of fact, periodic IFDEs create discrete dynamical systems, to which a general theory of difference equations can be applied [17,18].
The rest of this paper is organized as follows: Periodic IFDEs are discussed in Section 2, while other possibilities of impulsive effects on FDEs are discussed in Section 3 and Section 4, inspired by [19], wherein impulsive effects were not on the initial values but on the right-hand sides of FDEs. Finally, our conclusions and a discussion are presented in the last section.

2. Periodic IFDEs

Let N 0 = { 0 , 1 , 2 , } and consider
D t k q c x ( t ) = f ( t , x ( t ) ) , t ( t k , t k + 1 ) , k N 0 , x ( t k + ) = x ( t k ) + y k , k N , x ( 0 ) = x 0 ,
where q ( 0 , 1 ) . Assume that
(i)
f : R 0 × R m R m is continuous and T-periodic in t;
(ii)
there is a constant K 0 , such that f ( t , x ) f ( t , y )   K x y for all t R 0 and x , y R m ;
(iii)
{ t k } k = 1 R 0 is increasing, and { y k } k = 1 R m , such that there is an N N with T = t N + 1 , t k + N + 1 = t k + T and y k + N + 1 = y k , for all k N .
As covered in [2], under assumptions (i) and (ii), (9) has a unique solution on R 0 .
Now, consider the Poincaré mapping
P ( x 0 ) = x ( T + ) .
Clearly, the fixed points of P determine the T-periodic solutions of (9) ([20], Lemma 2.2). The following existence and uniqueness result follows ([21], Lemma 2.1 and Theorem 2.2).
Theorem 1.
Under assumptions (i) and (ii), it holds that
P ( x ) P ( y ) Θ x y , x , y R m
with
Θ = k = 0 N E q ( K ( t k + 1 t k ) q ) ,
where t 0 = 0 , and E q is the Mittag–Leffler function (see ([2], p. 40)). If (iii) holds as well and Θ < 1 , then (9) has a unique T-periodic solution x ¯ 0 R m , which is asymptotically stable, namely
x ¯ 0 P n ( y 0 )   Θ n x ¯ 0 y 0
for all n N and y 0 R m .
To this end, the following existence can be established.
Theorem 2.
Suppose that assumptions (i) and (iii) hold, and moreover,
(iv) 
there are constants K 0 and L 0 , such that f ( t , x )   K x + L for all t R 0 and x R m .
If Θ < 1 in (10), then (9) has a T-periodic solution x ¯ 0 R m , satisfying
x ¯ 0   r 0 = 1 1 Θ ( k = 0 N j = N k N E q ( K ( t j + 1 t j ) q ) × y N k + 1   + ( t N k + 1 t N k ) q Γ ( q + 1 ) ) .
Proof. 
On each interval ( t k , t k + 1 ) , k = 0 , 1 , , N , Equation (9) is equivalent to
x ( t ) = x ( t k + ) + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s .
This implies that
x ( t ) x ( t k + )   + 1 Γ ( q ) t k t ( t s ) q 1 f ( s , x ( s ) ) d s x ( t k + )   + 1 Γ ( q ) t k t ( t s ) q 1 ( K x ( s )   +   L ) d s = x ( t k + )   + L ( t t k ) q Γ ( q + 1 ) + K Γ ( q ) t k t ( t s ) q 1 x ( s ) d s .
Applying the Gronwall fractional inequality ([3], Theorem 1) to (12) leads to
x ( t ) x ( t k + )   + L ( t t k ) q Γ ( q + 1 ) E q ( K ( t t k ) q ) , t ( t k , t k + 1 ) ,
which gives
x ( t k + 1 + ) y k + 1   + L ( t k + 1 t k ) q Γ ( q + 1 ) E q ( K ( t k + 1 t k ) q ) + E q ( K ( t k + 1 t k ) q ) x ( t k + )
for k = 0 , 1 , , N . Consequently,
P ( x 0 ) = x ( t N + 1 + ) k = 0 N j = N k N E q ( K ( t j + 1 t j ) q ) × y N k + 1   + ( t N k + 1 t N k ) q Γ ( q + 1 ) + Θ x 0 .
If x 0   r 0 , then (11) implies P ( x 0 )   r 0 . Hence, P has a fixed point in the ball B r 0 ( 0 ) by the Brouwer fixed point theorem. The proof is complete. □
Remark 1.
Clearly, ( i i ) ( i v ) , but ( i v ) ( i i ) . Hence, (iv) does not guarantee the uniqueness of solutions.
From the above discussion, it follows that
D 0 q c x ( t ) = f ( t , x ( t ) ) , t ( t k , t k + 1 ) , k N 0 , x ( t k + ) = x ( t k ) + y k , k N , x ( 0 ) = x 0
has no periodic solutions. So, it suffices to consider only the periodic condition
x ( T ) = x ( 0 ) .
More related results can be found in [21].
Example 1.
Consider the simplest case of a linear FDE with a period-one impulse sequence, described by
D k q c x ( t ) = A x ( t ) , t ( k , k + 1 ) , k N 0 , x ( k + ) = x ( k ) + y , k N , x ( 0 ) = x 0 ,
where q ( 0 , 1 ) and A is an m × m matrix, with
P ( x 0 ) = x ( 1 + ) = E q ( A ) x 0 + y .
A fixed point x 0 of (14) is determined by
( I E q ( A ) ) x 0 = y ,
which is uniquely solvable, and
x 0 = ( I E q ( A ) ) 1 = y
if and only if 1 σ ( E q ( A ) ) .
If 1 σ ( E q ( A ) ) , then the Lyapunov–Schmidt method [22] can be applied. Specifically, for the case of m = 3 and A x = x × v with a vector product × and v = ( 0 , 0 , η ) , where η R , one has
D k q c x ( t ) = x ( t ) × v , t ( k , k + 1 ) , k N 0 , x ( k + ) = x ( k ) + y , k N , x ( 0 ) = x 0 .
Note that
( x 1 , x 2 , x 3 ) × v = η ( x 2 , x 1 , 0 ) .
Thus,
A = 0 η 0 η 0 0 0 0 0 .
Setting z = x 1 + x 2 ı , where ı = 1 , and identifying R 3 with C × R , i.e., ( x 1 , x 2 , x 3 ) = ( z , u ) , one obtains
( z , u ) × v = ( η ı z , 0 ) .
Then, (15) is transformed to
D k q c z ( t ) = η ı z ,   D k q c u ( t ) = 0 , t ( k , k + 1 ) , k N 0 , z ( k + ) = z ( k ) + y 1 , u ( k + ) = u ( k ) + y 2 , k N , z ( 0 ) = z 0 , u ( 0 ) = u 0
for y = ( y 1 , y 2 ) C × R . Thus, (14) becomes
P ( z 0 , u 0 ) = ( E q ( η ı ) z 0 + y 1 , u 0 + y 2 ) .
Since E q ( η ı ) is a complex number, the only fixed point of (18) is
z 0 = y 1 1 E q ( η ı ) , u 0
which exists only if y 2 = 0 .
Note that, for the matrix A defined in (16), σ ( A ) = { ± η ı , 0 } .
Now, consider a small perturbation ε R to (17), yielding
D k q c z ( t ) = η ı z ( t ) + ε u ( t ) z ( t ) , D k q c u ( t ) = ε u ( t ) z ( t ) , t ( k , k + 1 ) , k N 0 , z ( k + ) = z ( k ) + y 1 , u ( k + ) = u ( k ) + ε y 2 , k N , z ( 0 ) = z 0 , u ( 0 ) = u 0
Then, in (19) expanding
z ( t ) = z 0 ( t ) + ε z 1 ( t ) + O ( ε 2 ) , u ( t ) = u 0 ( t ) + ε u 1 ( t ) + O ( ε 2 ) ,
yields
D k q c z 0 ( t ) = η ı z 0 ( t ) , D k q c u 0 ( t ) = 0 , t ( k , k + 1 ) , k N 0 , z 0 ( k + ) = z 0 ( k ) + y 1 , u 0 ( k + ) = u 0 ( k ) , k N , z 0 ( 0 ) = z 0 , u 0 ( 0 ) = u 0
and
D k q c z 1 ( t ) = η ı z 1 ( t ) + u 0 ( t ) z 0 ( t ) , D k q c u 1 ( t ) = u 0 ( t ) ( z 0 ( t ) ) , t ( k , k + 1 ) , k N 0 , z 1 ( k + ) = z 1 ( k ) , u 1 ( k + ) = u 1 ( k ) + y 2 , k N , z ( 0 ) = 0 , u ( 0 ) = 0 .
Thus, (20) and (21) together give
z 0 ( t ) = E q ( η t q ı ) z 0 , u 0 ( t ) = u 0 , z 1 ( t ) = u 0 z 0 0 t ( t s ) q 1 E q , q ( η ( t s ) q ) E q ( η s q ı ) d s , u 1 ( t ) = u 0 Γ ( q ) 0 t ( t s ) q 1 ( E q ( η s q ı ) z 0 ) d s
for t [ 0 , 1 ) . Consequently, (14) becomes
P ( z 0 , u 0 ) = ( z 0 ( 1 ) + ε z 1 ( 1 ) + y 1 + O ( ε 2 ) , u 0 ( 1 ) + ε u 1 ( 1 ) + ε y 2 + O ( ε 2 ) ) = ( E q ( η ı ) + ε u 0 0 1 ( 1 s ) q 1 E q , q ( η ( 1 s ) q ) E q ( η s q ı ) d s z 0 + y 1 , 1 + ε 1 Γ ( q ) 0 1 ( 1 s ) q 1 ( E q ( η s q ı ) z 0 ) d s u 0 + ε y 2 ) + O ( ε 2 ) .
A fixed point of (22) is given by
z 0 = E q ( η ı ) + ε u 0 0 1 ( 1 s ) q 1 E q , q ( η ( 1 s ) q ) E q ( η s q ı ) d s z 0 + y 1 + O ( ε 2 ) , u 0 = 1 + ε 1 Γ ( q ) 0 1 ( 1 s ) q 1 ( E q ( η s q ı ) z 0 ) d s u 0 + ε y 2 + O ( ε 2 ) .
Solving the first equation of (23) gives
z 0 = 1 1 E q ( η ı ) ε u 0 0 1 ( 1 s ) q 1 E q , q ( η ( 1 s ) q ) E q ( η s q ı ) d s z 0 + y 1 1 E q ( η ı ) + O ( ε 2 ) = y 1 1 E q ( η ı ) + ε y 1 ( 1 E q ( η ı ) ) 2 u 0 0 1 ( 1 s ) q 1 E q , q ( η ( 1 s ) q ) E q ( η s q ı ) d s + O ( ε 2 ) .
It follows from the second equation of (23) with (24) that
0 = 1 Γ ( q ) 0 1 ( 1 s ) q 1 ( E q ( η s q ı ) z 0 ) d s u 0 + y 2 + O ( ε ) = u 0 Γ ( q ) 0 1 ( 1 s ) q 1 y 1 E q ( η s q ı ) 1 E q ( η ı ) d s + y 2 + O ( ε ) .
If
0 1 ( 1 s ) q 1 y 1 E q ( η s q ı ) 1 E q ( η ı ) d s 0 ,
then (25) gives
u 0 = y 2 Γ ( q ) 0 1 ( 1 s ) q 1 y 1 E q ( η s q ı ) 1 E q ( η ı ) d s + O ( ε ) .
Finally, inserting (27) into (24) leads to
z 0 = y 1 1 E q ( η ı ) ε y 1 y 2 Γ ( q ) 0 1 ( 1 s ) q 1 E q , q ( η ( 1 s ) q ) E q ( η s q ı ) d s ( 1 E q ( η ı ) ) 2 0 1 ( 1 s ) q 1 y 1 E q ( η s q ı ) 1 E q ( η ı ) d s + O ( ε 2 ) .
In summary, the following result is obtained.
Theorem 3.
Assuming (26), for any small ε R , IFDE (19) has a fixed point given by (27) and (28).
To this end, one can take specific values of q, η, ε, y 2 and y 1 in Theorem 3 to numerically study the dynamics of (18), which will not be further discussed here.

3. Weighted FDEs

Now, consider the case where the impulsive effects are not on the initial values but on the right-hand side of the FDEs.
Specifically, consider t k = k and take { η k } k N 0 R 0 . Motivated by [19], introduce the following FDE:
D 0 q c x ( t ) = f ( t , x ( t ) , ( 1 { t } ) η [ x ( t ) ] x ( t ) ) , t R 0 , x ( 0 ) = x 0 ,
where { t } = t [ t ] is the fractional part of t.
Example 2.
Consider
D 0 q c x = a x + b ( 1 { t } ) η [ x ] x + c x 3 + d cos t ,
which gives
x ( t ) = x ( 0 ) + 1 Γ ( q ) j = 0 k 1 j j + 1 ( t s ) q 1 ( a x ( s ) + b ( 1 + j s ) η j x ( s ) + c x 3 ( s ) + d cos s ) d s + 1 Γ ( q ) k t ( t s ) q 1 ( a x ( s ) + b ( 1 + k s ) η + c x 3 ( s ) + d cos s ) d s .
The solvability of (29) can be studied either by applying fixed point theorems or by numerical methods, which will not be further discussed here.

4. More General FDEs with Impulsive Effects

Consider another situation where the impulsive effects are not on the initial values but on the right-hand sides of the FDEs.
Let q ( 0 , 1 ] and take an increasing sequence P = { p k } k = 0 with p 0 = 0 and lim k p k = . Consider
D 0 q c x ( t ) = f ( t , x ( t ) , x ( p t ) ) , t R 0 , x ( 0 ) = x 0 ,
where · is an induced floor function given by
t = k , t [ p k , p k + 1 ) .
Thus, (30) can be rewritten as
D 0 q c x ( t ) = f ( t , x ( t ) , x ( p k ) ) , t [ p k , p k + 1 ) .
Consequently, (30) has “impulses” at the nonlinearity but not at the argument t. This is a kind of impulsive delay differential equation. The solution of (30) is given by
x ( t ) = x ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 f ( s , x ( s ) , x ( p s ) ) d s .
Now, add control parameters { y k } k = 0 R , { z k } k = 0 R and modify (30) to obtain the following:
D 0 q c x ( t ) = f ( t , x ( t ) , y t x ( p t ) + z t ) , x ( 0 ) = x 0 .
Thus, (32) becomes
D 0 q c x ( t ) = f ( t , x ( t ) , y k x ( p k ) + z k ) , t [ p k , p k + 1 ) ,
which gives
x ( t ) = x ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 f ( s , x ( s ) , y s x ( p s ) + z s ) d s .
A simple form of (32) is
D 0 q c x ( t ) = f ( t , x ( t ) ) + y t x ( p t ) + z t , x ( 0 ) = x 0 ,
where the unperturbed part is given by
D 0 q c x ( t ) = f ( t , x ( t ) ) , x ( 0 ) = x 0 .
One may consider controlling the solutions of (35) with (34), which can be written as
D 0 q c x ( t ) = f ( t , x ( t ) ) + y k x ( p k ) + z k , t [ p k , p k + 1 ) , x ( 0 ) = x 0 .
Furthermore, one may consider a weighted control problem:
D 0 q c x ( t ) = f ( t , x ( t ) ) + y t x ( p t ) e α ( t p t ) + z t e β ( t p t ) , x ( 0 ) = x 0 ,
namely,
D 0 q c x ( t ) = f ( t , x ( t ) ) + y k x ( p k ) e α ( t p k ) + z k e β ( t p k ) , t [ p k , p k + 1 ) , x ( 0 ) = x 0
for α , β R 0 .
Finally, one may consider a more general form of (30), such as
D 0 q c x ( t ) = f ( t , x ( t ) , x ( p t ) , x ( p t 1 ) ) , x ( 0 ) = x 0 ,
and extend the ideas described above, from (30) to (37).
Remark 2.
If P = N 0 , i.e., p k = k with k N 0 , then p t = t for the standard floor function · . This kind of equation is studied in [19] for a scalar case and for integer-order derivatives.
The simplest case of (30) is its linear version
D 0 q c x ( t ) = A x ( t ) + B x ( p t ) + g ( t ) , t R 0 , x ( 0 ) = x 0
with n × n matrices A , B M ( n ) and g C ( R 0 , R n ) . The homogeneous case
D 0 q c x ( t ) = A x ( t ) + B x ( p t ) , t R 0 , x ( 0 ) = x 0
defines a linear problem. This means that (39) has a fundamental matrix solution X ( t ) M ( n ) that solves
D 0 q c X ( t ) = A X ( t ) + B X ( p t ) , t R 0 , X ( 0 ) = I .
Therefore, the solution of (39) is x ( t ) = X ( t ) x 0 , t R 0 . A semilinear case is
D 0 q c x ( t ) = A x ( t ) + B x ( p t ) + g ( t , x ( t ) ) , t R 0 , x ( 0 ) = x 0
with g C ( R 0 × R n , R n ) .
Example 3.
Consider P = N 0 and a scalar case of (41):
D 0 q c x ( t ) = a x ( t ) + b x ( t ) + c x 3 ( t ) + d cos t .
By (31), one has
x ( t ) = x ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 ( a x ( s ) + b x ( s ) + c x 3 ( s ) + d cos s ) d s = x ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 ( a x ( s ) + c x 3 ( s ) + d cos s ) d s + b Γ ( q ) j = 0 k 1 x ( j ) j j + 1 ( t s ) q 1 d s + b Γ ( q ) x ( k ) k t ( t s ) q 1 d s = x ( 0 ) + 1 Γ ( q ) 0 t ( t s ) q 1 ( a x ( s ) + c x 3 ( s ) + d cos s ) d s + b j = 0 k 1 x ( j ) ( t j ) q ( t j 1 ) q Γ ( q + 1 ) + b x ( k ) ( t k ) q Γ ( q + 1 ) .
for t [ k , k + 1 ) .
Note that (43) implies that
x ( k + 1 ) = x ( 0 ) + 1 Γ ( q ) 0 k + 1 ( k + 1 s ) q 1 ( a x ( s ) + c x 3 ( s ) + d cos s ) d s + b j = 0 k 1 x ( j ) ( k + 1 j ) q ( k j ) q Γ ( q + 1 ) + b x ( k ) Γ ( q + 1 ) ,
which gives a recurrent formula for computing the x ( j ) in (43). One may perform a numerical simulation of (42) for certain values of q , a , b , c , d , which will not be further discussed here.
The influence of the impulsive delay in (42) can be studied by introducing its weight dependence, which is as follows:
D 0 q c x ( t ) = a x ( t ) + b e η { t } x ( t ) + c x 3 ( t ) + d cos t ,
where η R 0 and { t } = t ( t ) is the fractional part of t. Thus, (43) becomes
x ( t ) = x ( 0 ) + 1 Γ ( q ) j = 0 k 1 j j + 1 ( t s ) q 1 ( a x ( s ) + b e η ( s j ) x ( j ) + c x 3 ( s ) + d cos s ) d s + 1 Γ ( q ) k t ( t s ) q 1 ( a x ( s ) + b e η ( s k ) x ( k ) + c x 3 ( s ) + d cos s ) d s .
Consequently, the larger the η, the more x ( j ) is concentrated near j.
Example 4.
Motivated by the above example, but with non-fixed lower limits in Caputo derivatives, consider
D k q c x ( t ) = a x ( t ) + b e η ( k + 1 t ) x ( k ) p , t [ k , k + 1 ) , k N 0 , p { 1 , 2 } .
This gives
x ( t ) = E q ( a ( t k ) q ) x ( k ) + b x ( k ) p k t E q , q ( a ( t s ) q ) e η ( k + 1 s ) d s ,
where E q and E q , q are Mittag–Leffler functions, which implies that
x ( k + 1 ) = E q ( a ) x ( k ) + b x ( k ) p k k + 1 E q , q ( a ( k + 1 s ) q ) e η ( k + 1 s ) d s = E q ( a ) x ( k ) + b x ( k ) p 0 1 E q , q ( a s q ) e η s d s .
The dynamics of { x ( k ) } k N 0 are presented by a function
F ( x ) = E q ( a ) x ( k ) ( 1 b ^ x ( k ) p 1 ) , b ^ = b 0 1 E q , q ( a s q ) e η s d s E q ( a ) .
When p = 1 , (46) is a linear map, and it is stable for
E q ( a ) | 1 b ^ | 1 .
When p = 2 , (46) is the well-known logistic map [17,23] with complex dynamics. The equation
E q ( a ) = 3
characterizes a border between simple and non-simple dynamics—the period-doubling bifurcation—as shown by Figure 1.
Note that the equation
E q ( a ) = 1 + 2 2
characterizes a border between non-simple and chaotic dynamics. A period-three cycle emerges, and system (45) is chaotic in the sense of Li–Yorke chaos for every a E q 1 ( 1 + 2 2 ) , as shown by Figure 1.

5. Conclusions and Discussion

Ordinary differential equations (ODEs) play a crucial role in modelling evolutional processes with orbits depending continuously on time. When orbits have discontinuities in time for certain values, ODEs with impulses are used, forming the impulsive ODEs (IODEs). Another extension of ODEs is to generalize their integer derivatives to noninteger ones, which leads to fractional calculus and fractional differential equations (FDEs), which have been well developed [2,5,6,8,9]. The next natural development is the combination of the theories of IODEs and FDEs to create a new research direction of impulsive FDEs (IFDEs) [13]. This paper proposes some new forms of modelling impulsive effects in FDEs. Usually, impulses are presented in the time evolution of solutions. Here, impulsive effects are not imposed on initial values but on the right-hand sides of FDEs. The proposed impulsive models differ from the common discontinuous and nonsmooth dynamical systems in many aspects, as discussed in [24,25]. Clearly, a combination of different impulses can be an interesting topic to study. Hopefully, the ideas presented here could contribute to the further development of impulsive evolution equations. Note that other types of impulsive models exist. This is another good topic for future research. For example, changing the lower limit at t k instead of t 0 opens up the opportunity to use this kind of IFDE in various applications (see, e.g., [15]).

Author Contributions

Investigation, M.F., M.-F.D. and G.C., methodology, M.F., M.-F.D. and G.C.; writing—original draft preparation, M.F., M.-F.D. and G.C.; writing—review and editing, M.F., M.-F.D. and G.C.; visualization, M.F., M.-F.D. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was partially supported by the Slovak Research and Development Agency under contract No. APVV-23-0039 and by the Slovak Grant Agency VEGA No. 1/0084/23 and No. 2/0062/24.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Atkinson, C.; Osseiran, A. Rational solutions for the time-fractional diffusion equation. SIAM J. Appl. Math. 2011, 71, 92–106. [Google Scholar] [CrossRef]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  3. Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef]
  4. Benchohra, M.; Bouriah, S.; Salim, A.; Zhou, Y. Fractional Differential Equations: A Coincidence Degree Approach; Walter de Gruyter GmbH & Co KG: Berlin, Germany; Boston, MA, USA, 2023. [Google Scholar]
  5. Bainov, D.D.; Simeonov, P.S. Impulsive Differential Equations: Periodic Solutions and Applications; CRC Press: Boca Raton, FL, USA, 1993. [Google Scholar]
  6. Samoilenko, A.M.; Perestyuk, N.A. Impulsive Differential Equations; World Scientific: Singapore, 1995. [Google Scholar]
  7. Liu, X.; Ballinger, G. Boundedness for impulsive delay differential equations and applications to population growth models. Nonlinear Anal. Theory Methods Appl. 2003, 53, 1041–1062. [Google Scholar] [CrossRef]
  8. Chicone, C. Ordinary Differential Equations with Applications, 2nd ed.; Springer Science+Bussiness Media, Inc.: New York, NY, USA, 2006. [Google Scholar]
  9. Teschl, G. Ordinary Differential Equations and Dynamical Systems; American Mathematical Society: Providence, RI, USA, 2010. [Google Scholar]
  10. Stamov, T.; Stamov, G.; Stamova, I. Fractional-order impulsive delayed reaction-diffusion gene regulatory networks: Almost periodic solutions. Fractal Fract. 2023, 7, 384. [Google Scholar] [CrossRef]
  11. Wang, Q.; Lu, D.; Fang, Y. Stability analysis of impulsive fractional differential systems with delay. Appl. Math. Lett. 2015, 40, 1–6. [Google Scholar] [CrossRef]
  12. Abbas, S.; Benchohra, M.; Nieto, J.J. Caputo-Fabrizio fractional differential equations with instantaneous impulses. AIMS Math. 2021, 6, 2932–2946. [Google Scholar] [CrossRef]
  13. Wang, J.R.; Fečkan, M.; Zhou, Y. On the concept and existence of solution for impulsive fractional differential equations. Commun. Nonlinear Sci. Numer. Simulations 2012, 17, 3050–3060. [Google Scholar] [CrossRef]
  14. Wang, J.R.; Fečkan, M. Non-Instantaneous Impulsive Differential Equations, Basic Theory and Computation; IOP Publishing Ltd.: Bristol, UK, 2018. [Google Scholar]
  15. Danca, M.-F.; Fečkan, M. Memory principle of the Matlab code for Lyapunov exponents of fractional order. Int. J. Bifurc. Chaos, 2024, accepted. [CrossRef]
  16. Kaslik, E.; Sivasundaram, S. Non-existence of periodic solutions in fractional-order dynamical systems and a remarkable difference between integer and fractional-order derivatives of periodic functions. Nonlinear Anal. Real World Appl. 2012, 13, 1489–1497. [Google Scholar] [CrossRef]
  17. Elaydi, S.N. An Introduction to Difference Equations, 3rd ed.; Springer Science+Business Media, Inc.: New York, NY, USA, 2005. [Google Scholar]
  18. Irwin, M.C. Smooth Dynamical Systems; Academic Press: London, UK, 1980. [Google Scholar]
  19. Cooke, K.L.; Wiener, J. Retarded differential equations with piecewise constant delays. J. Math. Anal. Appl. 1984, 99, 265–297. [Google Scholar] [CrossRef]
  20. Fečkan, M.; Wang, J.R.; Zhou, Y. Periodic solutions for nonlinear evolution equations with non-instantaneous impulses. Nonautonomous Dyn. Sys. 2014, 1, 93–101. [Google Scholar]
  21. Fečkan, M.; Wang, J.R. Periodic impulsive fractional differential equations. Adv. Nonlinear Anal. 2019, 8, 482–496. [Google Scholar] [CrossRef]
  22. Deimling, K. Nonlinear Functional Analysis; Springer: Berlin/Heidelberg, Germany, 1985. [Google Scholar]
  23. Smítal, J. On Functions and Functional Equations; Adam Hilger: Bristol, UK, 1988. [Google Scholar]
  24. di Bernardo, M. Piecewise-Smooth Dynamical Systems: Theory and Applications; Springer: London, UK, 2008. [Google Scholar]
  25. Filippov, A.F. Differential Equations with Discontinuous Right-Hand Sides; Kluwer Academic: Dordrecht, The Netherlands, 1988. [Google Scholar]
Figure 1. The curves of (47) and (48).
Figure 1. The curves of (47) and (48).
Fractalfract 08 00500 g001
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

Fečkan, M.; Danca, M.-F.; Chen, G. Fractional Differential Equations with Impulsive Effects. Fractal Fract. 2024, 8, 500. https://doi.org/10.3390/fractalfract8090500

AMA Style

Fečkan M, Danca M-F, Chen G. Fractional Differential Equations with Impulsive Effects. Fractal and Fractional. 2024; 8(9):500. https://doi.org/10.3390/fractalfract8090500

Chicago/Turabian Style

Fečkan, Michal, Marius-F. Danca, and Guanrong Chen. 2024. "Fractional Differential Equations with Impulsive Effects" Fractal and Fractional 8, no. 9: 500. https://doi.org/10.3390/fractalfract8090500

APA Style

Fečkan, M., Danca, M. -F., & Chen, G. (2024). Fractional Differential Equations with Impulsive Effects. Fractal and Fractional, 8(9), 500. https://doi.org/10.3390/fractalfract8090500

Article Metrics

Back to TopTop