Stability, Boundedness, and Lagrange Stability of Fractional Differential Equations with Initial Time Difference

Abstract

Differential inequalities, comparison results, and sufficient conditions on initial time difference stability, boundedness, and Lagrange stability for fractional differential systems have been evaluated.

1. Introduction

The problem of stability of solutions is one of the major problems in the theory of differential equations. Lyapunov function and the Lyapunov direct method allow us to obtain sufficient conditions for the stability of a system without explicitly solving the differential equations [1, 2]. The method generalizes the idea which shows that the system is stable if there are some Lyapunov function candidates for the system.

Only a few decades ago, it was realized that fractional calculus provides an attractive tool for modelling the real world problems. The differentiation and integration of arbitrary orders have found applications in diverse fields of science and engineering like viscoelasticity, electrochemistry, diffusion processes, control theory, heat conduction, electricity, mechanics, chaos, and fractals [3–5]. Recently, some attention has been drawn on stability analysis of fractional differential equations (FDE) [6–9].

In practical situations, it is possible to have not only a change in initial position but also in initial time because of all kinds of disturbed factors. When we do consider such a deviation in initial time, it causes measuring the difference between any two different solutions starting with different initial times. From this point of view, several studies have been made on this problem to explore the stability and boundedness, criteria for differential systems relative to initial time difference (ITD) by using variation of parameters and differential inequalities technique [10–13]. In this paper, the stability and boundedness criteria for FDE relative to initial time difference have been investigated by using comparison method. In Section 2 the differences between classical notion of stability and the notion of initial time difference (ITD) stability have been discussed and compared by giving basic definitions. In Section 3 new differential inequalities and comparison results relative to initial time difference are obtained. Then, we investigate ITD stability, boundedness and Lagrange stability by using results obtained in Section 3. Lastly the conclusions are given in Section 4.

2. Fractional Calculus

Fractional calculus generalizes the derivative and the integral of a function to a noninteger order [3, 4, 14]. Although there are several definitions of fractional derivatives and fractional integrals, only the relative ones are given below. Let 𝑓∶[𝑎,𝑏]→ℝ be a function.

Definition 1. The fractional integral (or the Riemann-Liouville (RL) type integral) of order 𝑞>0 is defined as 𝐼^𝑞_𝑎𝑓(𝑡)=1Γ𝑞^𝑡_𝑎(𝑡−𝑠)^𝑞−1𝑓(𝑠)𝑑𝑠. (1)

Definition 2. The RL fractional derivatives of order 𝑞∈[𝑛−1,𝑛) of 𝑓(𝑡) are defined as _𝑎𝐷^𝑞_𝑡𝑓(𝑡)=𝑑^𝑛𝑑𝑡^𝑛𝑎𝐼_𝑡^𝑛−𝑞𝑓(𝑡)=1Γ𝑑𝑛−𝑞𝑑𝑡^𝑛𝑛𝑡_𝑎(𝑡−𝑠)^{𝑛−𝑞−1}𝑓(𝑠)𝑑𝑠. (2)

Definition 3. The Caputo fractional derivative of order 𝑞∈(𝑛−1,𝑛] of 𝑓(𝑡) is defined as ^𝑐_𝑎𝐷^𝑞_𝑡𝑓(𝑡)=𝐼_𝑎^𝑛−𝑞𝑑^𝑛𝑑𝑡^𝑛𝑓(𝑡=1)Γ𝑛−𝑞^𝑡_𝑎(𝑡−𝑠)^{𝑛−𝑞−1}𝑓^(𝑛)(𝑠)𝑑𝑠. (3) To simplify the notations we will use 𝐼^𝑞, 𝐷^𝑞, and ^𝑐𝐷^𝑞 instead of 𝐼^𝑞_𝑎, _𝑎𝐷^𝑞_𝑡, and ^𝑐_𝑎𝐷^𝑞_𝑡, respectively. There are some relation between these fractional integral and derivatives. For details please see [3, 4, 14].

Property 1. (i)𝐼^𝛼(𝐼^𝛽𝑓(𝑡))=𝐼^𝛼+𝛽𝑓(𝑡).(ii)ℒ[𝐷^𝑞𝑓(𝑡)]=𝑠^𝑞𝑋(𝑠)−∑^𝑛−1_𝑘=0𝑠^𝑘[𝐷^{𝑞−𝑘−1}𝑓(𝑡)]_𝑡=0, where 𝐷^𝑞 represents ₀𝐷^𝑞_𝑡.(iii)ℒ[^𝑐𝐷^𝑞𝑓(𝑡)]=𝑠^𝑞𝑋(𝑠)−∑^𝑛−1_𝑘=0𝑠^{𝑞−𝑘−1}𝑓^(𝑘)(0), where ^𝑐𝐷^𝑞 represents ^𝑐₀𝐷^𝑞_𝑡.(iv)^𝑐𝐷^𝑞𝐼^𝑞𝑓(𝑡)=𝑓(𝑡) and 𝐷^𝑞𝐼^𝑞𝑓(𝑡)=𝑓(𝑡).(v)𝐼^𝑞𝐷^𝑞𝑓(𝑡)=𝑓(𝑡)−∑^𝑛−1_𝑘=1[𝐷^𝑞−𝑘𝑓(𝑡)]_𝑡=𝑎((𝑡−𝑎)^𝑞−𝑘/Γ(𝑘−𝑞+1)) and 𝐼^𝑞𝑐𝐷^𝑞𝑓(𝑡)=𝑓(𝑡)0𝑥02012∑^𝑛−1_𝑘=0((𝑡−𝑎)^𝑘/𝑘!)𝑓^(𝑘)(𝑎).(vi)^𝑐𝐷^𝑞𝑓(𝑡)=𝐷^𝑞(𝑓(𝑡)−∑^𝑛−1_𝑘=0((𝑡−𝑎)^𝑘/𝑘!)𝑓^(𝑘)(𝑎)).(vii)𝐷^𝑞𝐶=𝐶(𝑡−𝑎)^−𝑞/Γ(1−𝑞) and ^𝑐𝐷^𝑞𝐶=0, where 𝐶 is arbitrary constant.

After giving definition and properties of fractional integral and derivatives, we consider fractional order IVP with (RL) and Caputo derivative, respectively; 𝐷^𝑞𝑥(𝑡)(),𝑥𝑡=𝑓𝑡,𝑥₀=𝑥⁰=Γ𝑞𝑥(𝑡)𝑡−𝑡₀^1−𝑞|_𝑡=𝑡0,_𝑐𝐷^𝑞𝑥(𝑡)()𝑥𝑡=𝑓𝑡,𝑥₀=𝑥₀, (4) where 0<𝑞<1, 𝑓∈𝐶[𝐽×ℝ,ℝ] and 𝐽=[𝑡₀,𝑇]. Then IVP's are equivalent to the following Volterra fractional integral equations: 𝑥(𝑡)=𝑥⁰𝑡−𝑡₀^𝑞−1Γ𝑞+1Γ𝑞^𝑡_𝑡0(𝑡−𝑠)^𝑞−1𝑓((𝑠𝑥(𝑡)𝑠,𝑥))𝑑𝑠,=𝑥₀+1Γ𝑞_𝑡^𝑡_𝑡0(𝑡−𝑠)^𝑞−1𝑓(𝑠,𝑥(𝑠))𝑑𝑠, (5)(6) respectively. In [15] authors develop a relation between the solutions of Caputo fractional differential equations and those of (RL) fractional differential equations.

3. Stability versus Initial Time Difference Stability

Consider the IVP for the system of nonlinear fractional differential equation ^𝑐𝐷^𝑞𝑥(𝑡)(),𝑥𝑡=𝑓𝑡,𝑥₀=𝑥₀, (7) where 𝑡₀∈ℝ₊, 𝑓∈𝐶[ℝ₊×ℝ^𝑛,ℝ^𝑛] and 0<𝑞<1. Suppose that the function 𝑓 is smooth enough to guarantee existence, uniqueness, and continuous dependence of solutions of IVP (7). Denote by 𝑥(𝑡)=𝑥(𝑡,𝑡₀,𝑥₀) the solution of (7). Consider also the initial value problem at a different initial data; that is, let 𝜏₀∈ℝ₊ and 𝑦(𝑡)=𝑦(𝑡,𝜏₀,𝑦₀) be the solution of the system (7). Assume that 𝑥(𝑡)=𝑥(𝑡,𝑡₀,𝑥₀) is the solution on which we shall study stability and boundedness criteria with respect to it. Set 𝜂=𝜏₀−𝑡₀>0 and denote 𝑆(𝜌)={𝑥∈ℝ^𝑛∶‖𝑥‖<𝜌}.

Before giving our comparison theorem, stability, boundedness criteria, and Lagrange stability for FDE we need to introduce the following definitions.

Definition 4. The solution 𝑥(𝑡)=𝑥(𝑡,𝑡₀,𝑥₀) of (7) is said to be(S1)stable with ITD, if given 𝜖>0 and 𝑡₀∈ℝ₊; there exist 𝛿=𝛿(𝜖,𝑡₀)>0 and 𝛿=𝛿(𝜖,𝑡₀)>0 such that 𝑦₀−𝑥₀𝜂<𝑦<𝛿,𝛿imply𝑡+𝜂,𝜏₀,𝑦₀−𝑥𝑡,𝑡₀,𝑥₀<𝜖for𝑡≥𝑡₀; (8)(S2)uniformly stable with ITD, if (𝑆₁) holds with 𝛿 and 𝛿 independent of 𝑡₀.

Definition 5. The system (7) is said to be(B1)equibounded with ITD, if given 𝛼>0 and 𝑡₀∈ℝ₊; there exist 𝛿=𝛿(𝛼,𝑡₀)>0 and 𝛽=𝛽(𝛼,𝑡₀)>0 such that 𝑦₀−𝑥₀𝜂<𝑦≤𝛼,𝛿imply𝑡+𝜂,𝜏₀,𝑦₀−𝑥𝑡,𝑡₀,𝑥₀<𝛽𝑡≥𝑡₀; (9)(B2)uniformly bounded with ITD, if (𝐵₁) holds with 𝛿 and 𝛽 independent of 𝑡₀;(A1)attractive in the large with ITD, if for every 𝜖>0 and each 𝛼>0 there exists 𝛿=𝛿(𝜖,𝛼,𝑡₀) and 𝑇=𝑇(𝜖,𝛼,𝑡₀) such that 𝑦₀−𝑥₀𝜂<𝑦≤𝛼,𝛿imply𝑡+𝜂,𝜏₀,𝑦₀−𝑥𝑡,𝑡₀,𝑥₀<𝜖for𝑡≥𝑡₀+𝑇; (10)(L1)Lagrange stable if (B1) and (A1) hold.

Definition 6. A function 𝑎 is said to belong to the class 𝒦 such that 𝒦=𝑎∈𝐶ℝ₊,ℝ₊∶𝑎(0)=0and𝑎isstrictlyincreasing. (11)

Definition 7. We define the generalized derivative with respect to the system (7) as follows: ^𝑐𝐷^𝑞₊𝑉𝑡,𝑦−𝑥=lim_ℎ→0+1supℎ^𝑞𝑉𝑡,𝑦−𝑥−𝑉𝑡−ℎ,𝑦−𝑥−ℎ^𝑞×𝑓𝑡+𝜂,𝑦−𝑓(𝑡,𝑥) (12) for (𝑡,𝑦−𝑥)∈ℝ₊×ℝ^𝑛, where 𝑥=𝑥(𝑡)=𝑥(𝑡,𝑡₀,𝑥₀) is the solution of (7) and 𝑦=𝑦(𝑡,𝑡₀,𝑦₀)=𝑦(𝑡+𝜂,𝜏₀,𝑦₀), where 𝑦(𝑡,𝜏₀,𝑦₀) is the solution of (7) and 𝜂=𝜏₀−𝑡₀>0.
The stability with ITD gives us an opportunity to compare solutions of FDE where both initial time and position are different. In the case of differential equation, stability with ITD is studied in [10–13, 16]. We will give a brief overview of both concepts of stability.

Case 1 (classical notion of stability). Let 𝑥(𝑡)=𝑥(𝑡,𝑡₀,𝑥₀) be a solution of (7). Study the stability of 𝑥(𝑡). Consider the solution 𝑌(𝑡)=𝑌(𝑡,𝑡₀,𝑦₀) of (7). Set the IVP as ^𝑐𝐷^𝑞𝑓(),𝑧𝑡𝑧=𝑡,𝑧₀=𝑦₀−𝑥₀, (13) where 𝑓(𝑡,𝑧)=𝑓(𝑡,𝑧+𝑥(𝑡))−𝑓(𝑡,𝑥(𝑡)). The function 𝑧(𝑡;𝑡₀,𝑦₀,𝑥₀)=𝑌(𝑡)−𝑥(𝑡) is a solution of the IVP (13).
The IVP (13) has a zero solution and the study of stability properties of the nonzero solution 𝑥(𝑡) of (7) is reduced to the stability of the zero solution of transformed system (13).

Case 2 (stability with ITD). Study the stability with initial time difference of 𝑥(𝑡). Consider the solution of (7) with different initial data as 𝑦(𝑡)=𝑦(𝑡,𝜏₀,𝑦₀). Set the IVP as ^𝑐𝐷^𝑞𝑓(),𝑧𝑡𝑧=𝑡,𝑧₀=𝑦₀−𝑥₀, (14) where 𝑓(𝑡,𝑧;𝜂)=𝑓(𝑡+𝜂,𝑧+𝑥(𝑡))−𝑓(𝑡,𝑥(𝑡)). Then 𝑧(𝑡;𝑡₀,𝜏₀,𝑦₀,𝑥₀)=𝑦(𝑡+𝜂)−𝑥(𝑡) is a solution of (14). The system (14) has no zero solution since 𝑓(𝑡,0)=𝑓(𝑡+𝜂,𝑥(𝑡))−𝑓(𝑡,𝑥(𝑡))≠0. Therefore, in this case the study of stability with ITD of 𝑥(𝑡) could not be reduced to the study of stability of the zero solution of an appropriate fractional differential system.

4. Main Results

4.1. Differential Inequalities and Comparison Results

4.1.1. Differential Inequalities

Definition 8 (see [15]). 𝑚 is said to be 𝐶^𝑞 continuous; that is, 𝑚∈𝐶^𝑞([𝑡₀,𝑇],ℝ), if and only if the Caputo derivative of ^𝑐𝐷^𝑞𝑚(𝑡) exists and satisfies ^𝑐𝐷^𝑞𝑚(𝑡)=1Γ1−𝑞^𝑡_𝑡0(𝑡−𝑠)^−𝑞𝑚^′(𝑠)𝑑𝑠. (15)

Lemma 9. Let 𝑚∈𝐶^𝑞([𝑡₀,𝑇],ℝ). Suppose that for any 𝑡₁∈(𝑡₀,𝑇], one has 𝑚(𝑡₁)=0 and 𝑚(𝑡)<0 for 𝑡₀≤𝑡<𝑡₁; then it follows that ^𝑐𝐷^𝑞𝑚𝑡₁>0. (16)

Proof. From the relation between Riemann-Liouville and Caputo fractional derivatives we write ^𝑐𝐷^𝑞𝑚(𝑡)=𝐷^𝑞𝑚(𝑡)−𝑚𝑡₀=𝐷^𝑞𝑚(𝑡)−𝑚𝑡₀𝑡−𝑡₀^−𝑞Γ.1−𝑞 (17) Since 𝑚(𝑡₀)<0 we have 𝐷^𝑞𝑚(𝑡₀)<0. Therefore, we obtain ^𝑐𝐷^𝑞𝑚(𝑡)>𝐷^𝑞𝑚(𝑡). (18) Finally using the lemma for R-L derivative from [6] and (18) implies that ^𝑐𝐷^𝑞𝑚𝑡₁>𝐷^𝑞𝑚𝑡₁≥0. (19)

Theorem 10. Let 𝑣,𝑤∈𝐶^𝑞([𝑡₀,𝑇],ℝ), 𝑓∈𝐶([𝑡₀,𝑇]×ℝ,ℝ) and (𝑖)^𝑐𝐷^𝑞𝑣(𝑡)((𝑡,()≤𝑓𝑡,𝑣))𝑖𝑖^𝑐𝐷^𝑞𝑤(𝑡)>𝑓(𝑡,𝑤(𝑡)),𝑡₀≤𝑡≤𝑇. (20) Then 𝑣(𝑡₀)<𝑤(𝑡₀) implies 𝑣(𝑡)<𝑤(𝑡),𝑡₀≤𝑡≤𝑇. (21)

Proof. Suppose that relation (21) is false. Then since 𝑣(𝑡₀)<𝑤(𝑡₀) and 𝑣(𝑡), 𝑤(𝑡) are continuous, there exists a 𝑡₁∈(𝑡₀,𝑇] with 𝑣(𝑡₁)=𝑤(𝑡₁) and 𝑣(𝑡)<𝑤(𝑡) for 𝑡₀≤𝑡<𝑡₁. Set 𝑚(𝑡)=𝑣(𝑡)−𝑤(𝑡). Then 𝑚(𝑡₁)=0 and 𝑚(𝑡)<0 for 𝑡₁∈[𝑡₀,𝑡₁). Hence, the hypothesis of Lemma 9 holds and we conclude that ^𝑐𝐷^𝑞𝑚(𝑡₁)>0, which means that ^𝑐𝐷^𝑞𝑣(𝑡₁)>^𝑐𝐷^𝑞𝑤(𝑡₁). Consider 𝑓𝑡₁,𝑣𝑡₁≥^𝑐𝐷^𝑞𝑣𝑡₁>^𝑐𝐷^𝑞𝑤𝑡₁>𝑓𝑡₁,𝑤𝑡₁ (22) which is a contradiction. Thus the conclusion of the theorem holds and the proof is complete.

Theorem 11. Let 𝑣,𝑤∈𝐶^𝑞([𝑡₀,𝑇],ℝ), 𝑓∈𝐶([𝑡₀,𝑇]×ℝ,ℝ) and (𝑖)^𝑐𝐷^𝑞𝑣(𝑡)((𝑡,()≤𝑓𝑡,𝑣))𝑖𝑖^𝑐𝐷^𝑞𝑤(𝑡)≥𝑓(𝑡,𝑤(𝑡)),𝑡₀≤𝑡≤𝑇. (23) Assume 𝑓 satisfies the Lipschitz condition 𝑓(𝑡,𝑥)−𝑓𝑡,𝑦≤𝐿𝑥−𝑦,𝑥≥𝑦,𝐿>0. (24) Then, 𝑣(𝑡₀)≤𝑤(𝑡₀) implies 𝑣(𝑡)≤𝑤(𝑡),𝑡₀≤𝑡≤𝑇. (25)

Proof. We set 𝑤_𝜖(𝑡)=𝑤(𝑡)+𝜖𝜆(𝑡), where 𝜆(𝑡)=𝜆(𝑡₀)𝐸_𝑞(2𝐿(𝑡−𝑡₀)^𝑞) is the solution of linear fractional differential equation ^𝑐𝐷^𝑞𝜆(𝑡)=2𝐿𝜆(𝑡), 𝜆(𝑡₀)=𝜆₀. Then from 𝑤_𝜖(𝑡₀)=𝑤(𝑡₀)+𝜖𝜆(𝑡₀), we get 𝑤_𝜖(𝑡₀)>𝑤(𝑡₀)≥𝑣(𝑡₀). In order to get (25), we need to apply Theorem 10 to 𝑣(𝑡) and 𝑤_𝜖(𝑡). Using (𝑖𝑖) and Lipschitz condition we have ^𝑐𝐷^𝑞𝑤_𝜖(𝑡)=^𝑐𝐷^𝑞𝑤(𝑡)+^𝑐𝐷^𝑞𝜆(𝑡)≥𝑓(𝑡,𝑤(𝑡))−𝑓𝑡,𝑤_𝜖(𝑡)+𝑓𝑡,𝑤_𝜖(𝑡)+𝜖2𝐿𝜆_𝑐(𝑡),𝐷^𝑞𝑤_𝜖(𝑡)>𝑓𝑡,𝑤_𝜖(𝑡). (26)
Applying now Theorem 10 to 𝑣(𝑡) and 𝑤_𝜖(𝑡), we get 𝑣(𝑡)<𝑤_𝜖(𝑡) for every 𝜖>0 and consequently making 𝜖→0, we get the desired estimate (25).

4.1.2. Comparison Results

The most commonly used technique in the theory of differential equations is related to the estimation of a function satisfying a differential inequality by the extremal solutions of the related differential equation. The following theorems give such estimate with initial time difference.

Theorem 12. Assume that 𝑚∈𝐶^𝑞([𝑡₀,𝑇],ℝ) and ^𝑐𝐷^𝑞𝑚(𝑡)≤𝑔(𝑡,𝑚(𝑡)),𝑡₀≤𝑡≤𝑇, (27) where 𝑔∈𝐶([𝑡₀,𝑇]×ℝ,ℝ). Let 𝑟(𝑡) be the maximal solution of the IVP ^𝑐𝐷^𝑞𝑢=𝑔(𝑡,𝑢),𝑢𝑡₀=𝑢₀, (28) existing on [𝑡₀,𝑇] such that 𝑚(𝑡₀)≤𝑢₀. Then one has 𝑚(𝑡)≤𝑟(𝑡),𝑡₀≤𝑡≤𝑇. (29)

Proof. In view of the definition of the maximal solution 𝑟(𝑡), it is enough to prove that 𝑚(𝑡)<𝑢(𝑡,𝜖), 𝑡₀≤𝑡≤𝑇, where 𝑢(𝑡,𝜖) is any solution of the IVP ^𝑐𝐷^𝑞𝑢=𝑔(𝑡,𝑢)+𝜖,withtheinitialvalue𝑢𝑡₀+𝜖,𝜖>0. (30) Using Lemma 9, we get 𝑚(𝑡)<𝑢(𝑡,𝜖). And from lim_𝜖→0𝑢(𝑡,𝜖)=𝑟(𝑡) uniformly on each compact set 𝑡₀≤𝑡≤𝑇₀<𝑇, we get the desired estimate (29).

Theorem 13. Assume that(i)𝑚∈𝐶^𝑞([𝑡₀,𝑇],ℝ₊), 𝑔∈𝐶([𝑡₀,𝑇]×ℝ₊,ℝ) and ^𝑐𝐷^𝑞𝑚(𝑡)≤𝑔(𝑡,𝑚(𝑡)),𝑚𝑡₀≤𝑤₀; (31)(ii)the maximal solution 𝑟(𝑡)=𝑟(𝑡,𝜏₀,𝑤₀) of the IVP ^𝑐𝐷^𝑞𝑤=𝑔(𝑡,𝑤),𝑤𝜏₀=𝑤₀≥0 (32) exists for 𝑡≥𝜏₀;(iii)𝑔(𝑡,𝑤) is nondecreasing in 𝑡 for each 𝑤 and 𝜏₀>𝑡₀.Then (a) 𝑚(𝑡)≤𝑟(𝑡+𝜂), 𝑡≥𝑡₀ and (b) 𝑚(𝑡−𝜂)≤𝑟(𝑡), 𝑡≥𝜏₀.

Proof. (a) It is well known that if 𝑤(𝑡,𝜖) is any solution of ^𝑐𝐷^𝑞𝑤=𝑔(𝑡,𝑤)+𝜖,𝑤𝜏₀=𝑤₀+𝜖for𝜖>0 (33) sufficiently small then lim_𝜖→0𝑤(𝑡,𝜖)=𝑟(𝑡,𝜏₀,𝑤₀) on every compact set [𝜏₀,𝜏₀+𝑇].
Setting 𝑤₀(𝑡,𝜖)=𝑤(𝑡+𝜂,𝜖), we have 𝑤₀(𝑡₀,𝜖)=𝑤(𝑡₀+𝜂,𝜖)=𝑤(𝜏₀,𝜖)=𝑤₀+𝜖>𝑤₀≥𝑚(𝑡₀). Thus we get 𝑚(𝑡₀)<𝑤₀(𝑡₀,𝜖). On the other hand, using (iii) ^𝑐𝐷^𝑞𝑤₀(𝑡,𝜖)=𝑔𝑡+𝜂,𝑤₀(𝑡,𝜖)+𝜖>𝑔𝑡+𝜂,𝑤₀(𝑡,𝜖)>𝑔𝑡,𝑤₀(𝑡,𝜖),𝑡≥𝑡₀. (34) Then we get 𝑚(𝑡)<𝑤₀(𝑡,𝜖),  𝑡≥𝑡₀ and hence it follows that 𝑚(𝑡)≤𝑟(𝑡+𝜂),  𝑡≥𝑡₀.
(b) We set 𝑚₀(𝑡)=𝑚(𝑡−𝜂) so that 𝑚₀(𝜏₀)=𝑚(𝑡₀)≤𝑤₀<𝑤₀+𝜖 and ^𝑐𝐷^𝑞𝑚₀(𝑡)≤𝑔𝑡−𝜂,𝑚₀(𝑡)<𝑔𝑡,𝑚₀(𝑡),𝑡≥𝜏₀. (35) Then we get 𝑚₀(𝑡)<𝑤(𝑡,𝜖), 𝑡≥𝜏₀. The conclusion follows taking the limit as 𝜖→0. The proof of theorem is complete.

In the following theorem, we obtain a comparison result in terms of Lyapunov-like functions with ITD.

Theorem 14. Assume that(i)𝑉∈𝐶[ℝ₊×ℝ^𝑛,ℝ₊], 𝑉(𝑡,𝑥) is locally Lipschitzian in 𝑥∈ℝ^𝑛, 𝑔∈𝐶[ℝ²₊,ℝ] and ^𝑐𝐷^𝑞₊𝑉𝑡,𝑦𝑡+𝜂−𝑥(𝑡)≤𝑔𝑡,𝑉𝑡,𝑦𝑡+𝜂−𝑥(𝑡); (36)(ii)the maximal solution 𝑟(𝑡)=𝑟(𝑡,𝜏₀,𝑢₀) of the fractional scalar differential equation ^𝑐𝐷^𝑞𝑢=𝑔(𝑡,𝑢),𝑢𝜏₀=𝑢₀≥0, (37) exists for 𝑡≥𝜏₀;(iii)𝑔(𝑡,𝑢) is nondecreasing in 𝑡 for each 𝑢.Then 𝑉(𝑡₀,𝑦₀−𝑥₀)≤𝑢₀ implies 𝑉𝑡,𝑦𝑡+𝜂,𝜏₀,𝑦₀−𝑥𝑡,𝑡₀,𝑥₀≤𝑟𝑡+𝜂,𝜏₀,𝑢₀𝑓𝑜𝑟𝑡≥𝑡₀. (38)

Proof. Define 𝑚(𝑡)=𝑉(𝑡,𝑦(𝑡+𝜂,𝜏₀,𝑦₀)−𝑥(𝑡,𝑡₀,𝑥₀)) so that 𝑚𝑡₀=𝑉𝑡₀,𝑦₀−𝑥₀≤𝑢₀. (39) Let 𝑧(𝑡,𝑡₀,𝑦₀−𝑥₀)=𝑦(𝑡+𝜂,𝜏₀,𝑦₀)−𝑥(𝑡,𝑡₀,𝑥₀) so that ^𝑐𝐷^𝑞𝑓()=𝑓((𝑡,𝑧𝑡𝑧=𝑡,𝑧𝑡+𝜂,𝑦𝑡+𝜂−𝑓𝑡,𝑥))₀=𝑦₀−𝑥₀,𝑚(𝑡)−𝑚(𝑡−ℎ)ℎ^𝑞=𝑉(𝑡,𝑧(𝑡))−𝑉𝑡−ℎ,𝑧(𝑡)−ℎ^𝑞𝑓(𝑡,𝑧)ℎ^𝑞+𝑉𝑡−ℎ,𝑧(𝑡)−ℎ^𝑞𝑓(𝑡,𝑧)−𝑉𝑡−ℎ,𝑆𝑧,ℎ,𝑞ℎ^𝑞, (40) where 𝑆(𝑧,ℎ,𝑞)=𝑧(𝑡)−ℎ^𝑞(𝑓(𝑡,𝑧))−𝜖(ℎ^𝑞). Since 𝑉 is locally Lipschitzian in 𝑥 and 𝐿>0 is the Lipschitz constant and 𝜖(ℎ^𝑞)/ℎ^𝑞→0 as ℎ→0, we have ^𝑐𝐷^𝑞₊𝑚(𝑡)≤lim_ℎ→0+𝐿𝜖ℎ^𝑞ℎ^𝑞+limsup_ℎ→0+𝑉𝑡,𝑦𝑡+𝜂−𝑥(𝑡)ℎ^𝑞−𝑉𝑡−ℎ,𝑦𝑡+𝜂−𝑥(𝑡)ℎ^𝑞−ℎ^𝑞𝑓𝑡+𝜂,𝑦𝑡+𝜂−𝑓(𝑡,𝑥(𝑡))ℎ^𝑞≤𝑔(𝑡,𝑚(𝑡)). (41) By using Theorem 13 we obtain that 𝑚(𝑡)=𝑉𝑡,𝑦𝑡+𝜂,𝜏₀,𝑦₀−𝑥𝑡,𝑡₀,𝑥₀≤𝑟𝑡+𝜂,𝜏₀,𝑢₀,for𝑡≥𝑡₀. (42)

4.2. Stability and Boundedness Criteria

A comparison principle is obtained by employing the notion of Lyapunov function together with the theory of differential inequalities. In this part, one can see Lyapunov-like function as a transformation which reduces the study of stability, boundedness, and Lagrange stability properties relative to ITD of a given complicated system to a relatively simpler scalar equation.

4.2.1. Stability Criteria

Theorem 15. Assume that(i)𝑉∈𝐶[ℝ₊×ℝ^𝑛,ℝ₊], 𝑉(𝑡,𝑥) is locally Lipschitzian in 𝑥∈ℝ^𝑛, 𝑔∈𝐶[ℝ²₊,ℝ] and ^𝑐𝐷^𝑞₊𝑉𝑡,𝑦𝑡+𝜂−𝑥(𝑡)≤𝑔𝑡,𝑉𝑡,𝑦𝑡+𝜂−𝑥(𝑡) (43)(ii)the maximal solution 𝑟(𝑡)=𝑟(𝑡,𝜏₀,𝑢₀) of (37) exists for 𝑡≥𝜏₀;(iii)there exists 𝑎,𝑏∈𝒦 such that 𝑏(‖𝑥‖)≤𝑉(𝑡,𝑥)≤𝑎(‖𝑥‖) for (𝑡,𝑥)∈ℝ₊×𝑆(𝜌);(iv)𝑔(𝑡,𝑢) is nondecreasing in 𝑡 for each 𝑢 and 𝑔(𝑡,0)=0;Then the stability properties of the null solution of (37) imply the corresponding initial time difference stability properties of the solution 𝑥(𝑡,𝑡₀,𝑥₀).

Proof. Assume that the null solution of (37) is equistable. Let 0<𝜖<𝜌 be given. Then by definition of equistability given 𝑏(𝜖)>0, 𝜏₀∈ℝ₊, there exist a 𝛿₁=𝛿₁(𝜖,𝜏₀) such that 𝑢(𝑡)<𝑏(𝜖)providedthat𝑢₀<𝛿₁,𝑡≥𝜏₀, (44) where 𝑢(𝑡,𝜏₀,𝑢₀) is any solution of the (37). Choose 𝛿₂=𝛿₂(𝜖,𝜏₀)>0 as 0<𝑎(𝛿₂)<𝛿₁. Obviously lim_{(𝜏0,𝑦0)→(𝑡0,𝑥0)}‖𝑦(𝑡+𝜂,𝜏₀,𝑦₀)−𝑥(𝑡,𝑡₀,𝑥₀)‖=0. Then given 𝜖>0 and 𝜏₀∈ℝ₊, 𝛿=𝛿(𝜖,𝜏₀)>0 and 𝛿₃=𝛿₃(𝜖,𝜏₀)>0 such that 𝑦₀−𝑥₀<𝛿₃,𝜂<𝛿implies𝑦𝑡+𝜂−𝑥(𝑡)<𝜖for𝑡₀≤𝑡≤𝜏₀ (45) Let 𝑢₀=𝑎(‖𝑦₀−𝑥₀‖) and choose 𝛿=min(𝛿₂,𝛿₃). Then we claim that 𝑦𝑡+𝜂−𝑥(𝑡)<𝜖providedthat𝑦₀−𝑥₀𝜂<<𝛿,𝛿for𝑡≥𝑡₀. (46) If it is not true, from (45), there exist a 𝑡₁>𝜏₀ and a solution 𝑦(𝑡,𝜏₀,𝑦₀) with ‖𝑦₀−𝑥₀‖<𝛿 and |𝜂|<𝛿 such that 𝑦𝑡₁+𝜂−𝑥𝑡₁𝑦(𝑡≥𝜖,𝑡+𝜂−𝑥)<𝜖,for𝑡₀≤𝑡<𝑡₁. (47) Moreover since ‖𝑦₀−𝑥₀‖<𝛿, by (iii) we have 𝑉𝑡₀,𝑦₀−𝑥₀≤𝑎𝑦₀−𝑥₀<𝑎(𝛿)<𝛿₁. (48) Hence, by (i), (ii), (47) and, Theorem 14, we obtain the following estimate: 𝑉𝑡,𝑦𝑡+𝜂−𝑥(𝑡)≤𝑟𝑡+𝜂,𝜏₀,𝑢₀,𝑡₀≤𝑡<𝑡₁. (49) Consequently, the relations (44), (47), (49), and (iii) lead to the contradiction 𝑏(𝜖)≤𝑏𝑦𝑡₁+𝜂−𝑥𝑡₁≤𝑉𝑡₁,𝑦𝑡₁+𝜂−𝑥𝑡₁≤𝑟𝑡₁+𝜂,𝜏₀,𝑢₀<𝑏(𝜖). (50) which proves that the solution 𝑥(𝑡,𝑡₀,𝑥₀) of (7) is equistable with ITD.

4.2.2. Boundedness Criteria

In this section, Lagrange stability, which includes boundedness criteria, of fractional dynamic systems are discussed by employing comparison method.

Theorem 16. Let the assumption of Theorem 15 hold. Then the Lagrange stability properties of (37) imply the corresponding initial time difference Lagrange stability properties of (7).

Proof. We need to prove (B1) and (S3) for (7). Let 𝛼≥0 be given, and let ‖𝑦₀−𝑥₀‖≤𝛼. In view of (iii) 𝑉(𝑡₀,𝑦₀−𝑥₀)≤𝑎(𝛼)=𝛼₁. Assume that (37) is Lagrange stable. It follows that (37) is bounded. Then, given 𝛼₁≥0 and 𝜏₀∈ℝ₊ there exists a 𝛽₁=𝛽₁(𝜏₀,𝛼) such that 𝑢₀≤𝛼₁imply𝑟𝑡,𝜏₀,𝑢₀<𝛽₁. (51) Moreover, 𝑏(𝑢)→∞ as 𝑢→∞; we can choose a 𝛽=𝛽(𝜏₀,𝛼) verifying the relation 𝑏𝛽≥𝛽₁. (52) Then from (45), given 𝛽>0 and 𝜏₀∈ℝ₊, 𝛿=𝛿(𝜖,𝜏₀)>0 and 𝛿₁=𝛿₁(𝜖,𝜏₀)>0 such that 𝑦₀−𝑥₀<𝛿₁,𝜂<𝛿implies𝑦𝑡+𝜂−𝑥(𝑡)<𝛽for𝑡₀≤𝑡≤𝜏₀. (53) Let 𝑢₀=𝑉(𝜏₀,𝑦₀−𝑥₀). Then we claim that 𝑦𝑡+𝜂−𝑥(𝑡)<𝛽providedthat𝑦₀−𝑥₀𝜂<<𝛼,𝛿for𝑡≥𝑡₀. (54) If this were false, there would exist a 𝑡^∗>𝜏₀ and a solution 𝑦(𝑡,𝜏₀,𝑦₀) of (7) such that ‖𝑥(𝑡^∗+𝜂)−𝑥(𝑡^∗)‖=𝛽. Consequently, the relations (49), (51), (52), and (iii) lead to the contradiction 𝑏𝛽≤𝑏𝑦𝑡^∗+𝜂−𝑥𝑡^∗≤𝑉𝑡^∗,𝑦𝑡^∗+𝜂−𝑥𝑡^∗≤𝑟𝑡^∗+𝜂,𝜏₀,𝑢₀<𝛽₁≤𝑏𝛽 (55) which proves that (B1) holds for (7).
To prove attractivity, we let 𝜖>0, 𝛼≥0, 𝜏₀∈ℝ₊ be given and ‖𝑦₀−𝑥₀‖≤𝛼. In view of (iii), 𝑉(𝑡₀,𝑦₀−𝑥₀)≤𝑎(𝛼)=𝛼₁. Since (37) is attractive in the large with ITD, given 𝛼₁≥0, 𝑏(𝜖) and 𝜏₀∈ℝ₊ there exist a 𝑇=𝑇(𝜏₀,𝛼,𝜖) such that 𝑢₀≤𝛼₁implies𝑢𝑡,𝜏₀,𝑢₀<𝑏(𝜖),for𝑡≥𝜏₀+𝑇. (56) We have (49) from the proof of Theorem (45). Now, suppose that there exists a sequence {𝑡_𝑘}∈ℝ₊, 𝑡_𝑘→∞ as 𝑘→∞, 𝑡_𝑘>𝜏₀+𝑇 and a solution 𝑦(𝑡,𝜏₀,𝑦₀) of (7) such that 𝑦𝑡_𝑘+𝜂−𝑥𝑡_𝑘≥𝜖. (57) The relations (iii), (49), (56), and (57) lead to the contradiction, 𝑏(𝜖)≤𝑏𝑦𝑡_𝑘+𝜂−𝑥𝑡_𝑘≤𝑉𝑡_𝑘,𝑦𝑡_𝑘+𝜂−𝑥𝑡_𝑘≤𝑟𝑡_𝑘+𝜂,𝜏₀,𝑢₀<𝑏(𝜖) (58) which proves that the system is attractive in the large with ITD.

5. Conclusion

Firstly, differential inequalities and a new comparison principle for fractional differential equations relative to initial time difference have been developed and then stability, boundedness criteria, and Lagrange stability relative to initial time difference have been proved by employing comparison method.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.