Citation: Xinwei Su, Shuqin Zhang, Lixin Zhang. Periodic boundary value problem involving sequential fractional derivatives in Banach space[J]. AIMS Mathematics, 2020, 5(6): 7510-7530. doi: 10.3934/math.2020481
[1] | Yujun Cui, Chunyu Liang, Yumei Zou . Existence and uniqueness of solutions for a class of fractional differential equation with lower-order derivative dependence. AIMS Mathematics, 2025, 10(2): 3797-3818. doi: 10.3934/math.2025176 |
[2] | Muneerah Al Nuwairan, Ahmed Gamal Ibrahim . Solutions and anti-periodic solutions for impulsive differential equations and inclusions containing Atangana-Baleanu fractional derivative of order $ \zeta \in (1, 2) $ in infinite dimensional Banach spaces. AIMS Mathematics, 2024, 9(4): 10386-10415. doi: 10.3934/math.2024508 |
[3] | Cuiying Li, Rui Wu, Ranzhuo Ma . Existence of solutions for Caputo fractional iterative equations under several boundary value conditions. AIMS Mathematics, 2023, 8(1): 317-339. doi: 10.3934/math.2023015 |
[4] | Chengbo Zhai, Yuanyuan Ma, Hongyu Li . Unique positive solution for a p-Laplacian fractional differential boundary value problem involving Riemann-Stieltjes integral. AIMS Mathematics, 2020, 5(5): 4754-4769. doi: 10.3934/math.2020304 |
[5] | Yitao Yang, Dehong Ji . Properties of positive solutions for a fractional boundary value problem involving fractional derivative with respect to another function. AIMS Mathematics, 2020, 5(6): 7359-7371. doi: 10.3934/math.2020471 |
[6] | Choukri Derbazi, Hadda Hammouche . Caputo-Hadamard fractional differential equations with nonlocal fractional integro-differential boundary conditions via topological degree theory. AIMS Mathematics, 2020, 5(3): 2694-2709. doi: 10.3934/math.2020174 |
[7] | Karim Guida, Lahcen Ibnelazyz, Khalid Hilal, Said Melliani . Existence and uniqueness results for sequential $ \psi $-Hilfer fractional pantograph differential equations with mixed nonlocal boundary conditions. AIMS Mathematics, 2021, 6(8): 8239-8255. doi: 10.3934/math.2021477 |
[8] | Choukri Derbazi, Zidane Baitiche, Mohammed S. Abdo, Thabet Abdeljawad . Qualitative analysis of fractional relaxation equation and coupled system with Ψ-Caputo fractional derivative in Banach spaces. AIMS Mathematics, 2021, 6(3): 2486-2509. doi: 10.3934/math.2021151 |
[9] | Bashir Ahmad, Ahmed Alsaedi, Mona Alsulami, Sotiris K. Ntouyas . Existence theory for coupled nonlinear third-order ordinary differential equations with nonlocal multi-point anti-periodic type boundary conditions on an arbitrary domain. AIMS Mathematics, 2019, 4(6): 1634-1663. doi: 10.3934/math.2019.6.1634 |
[10] | Jiqiang Zhang, Siraj Ul Haq, Akbar Zada, Ioan-Lucian Popa . Stieltjes integral boundary value problem involving a nonlinear multi-term Caputo-type sequential fractional integro-differential equation. AIMS Mathematics, 2023, 8(12): 28413-28434. doi: 10.3934/math.20231454 |
In this work, we consider the following periodic boundary value problem (PBVP for short) in a Banach space E
{D2αx(t)=f(t,x(t),Dαx(t)), t∈(0,1], 0<α≤1,x(0)=x(1), Dαx(0)=Dαx(1), | (1.1) |
where f(t,x,y) is a continuous E-value function on [0,1]×E×E, Dα is the conformable fractional derivative of order α, D2α=DαDα is the conformable sequential fractional derivative.
Sequential fractional derivative for a sufficiently smooth function g(t) due to Miller and Ross [1] is defined as Dδg(t)=Dδ1Dδ2…Dδkg(t), here δ=(δ1,δ2,…,δk) is a multi-index. In general, the operator Dδ can be Riemann–Liouville or Caputo or any other kind of differential operators. There is a close connection between the sequential fractional derivatives and the non-sequential derivatives [2]. Many research papers have appeared recently concerning the existence of solutions for fractional differential equations involving Riemann-Liouville or Caputo sequential fractional derivatives by techniques of nonlinear analysis such as fixed point theorems, coincidence degree continuation theorems and nonlinear alternatives, see, for example, the papers [3,4,5,6,7,8,9,10,11] and the references therein.
In recent years, there has been a significant development in ordinary and partial differential equations involving fractional derivatives due to their wide range of applications in varied fields of science and engineering. Generally, for most of the fractional differential equations, it is difficult to find exact solutions in closed forms. In most cases, only approximate solutions or numerical solutions can be expected. Therefore, many iterative methods have been designed to be one of the suitable and successful classes of numerical techniques for obtaining the solutions of numerous types of fractional differential equations, see, for instance, [12,13,14,15,16] and the references therein. The monotone iterative technique, combined with the method of upper and lower solutions, is an effective technique for proving the existence of solutions for initial and boundary value problems of nonlinear differential equations. The basic idea of this method is that by choosing upper and lower solutions as two initial iterations, one can construct the monotone sequences for a corresponding linear equation and that converge monotonically to the extremal solutions of the nonlinear equation. Not only does this method give constructive proof for existence theorems but also the monotone behavior of iterative sequences is useful in the treatment of numerical solutions of various initial and boundary value problems. So many authors developed the upper and lower solutions method to investigate fractional differential equations, see, for example, [17,18,19,20,21,22,23] and the references therein.
In [2], using the method of upper and lower solutions and its associated monotone iterative technique, the authors considered the existence of minimal and maximal solutions and uniqueness of solution of the following initial value problem for fractional differential equation involving Riemann-Liouville sequential fractional derivative
{D2αy(t)=f(t,y(t),Dαy(t)), t∈(0,T], 0<α≤1,t1−αy(t)|t=0=y0, t1−αDαx(t)|t=0=y1, |
where f(t,x,y) is continuous on [0,1]×R×R.
The nonlinear impulsive fractional differential equation with periodic boundary conditions
{D2αu(t)=f(t,u,Dαu), t∈(0,T]∖{t1,…,tm}, 0<α≤1,limt→0+t1−αu(t)=u(1), limt→0+t1−αDαu(t)=Dαu(1),limt→t+j(t−tj)1−α(u(t)−u(tj))=Ij(u(tj)),limt→t+j(t−tj)1−α(Dαu(t)−Dαu(tj))=¯Ij(u(tj)), |
is studied in [24], where Dα is the standard Riemann-Liouville fractional derivative and D2α=DαDα is the Riemann-Liouville sequential fractional derivative. Ij,¯Ij∈C(R,R),j=1,2,…,m. f is continuous at every point (t,u,v)∈[0,1]×R×R. The existence and uniqueness results of solutions are obtained by the method of upper and lower solutions and its associated monotone iterative technique.
A new simple well-behaved definition of the fractional derivative called conformable fractional derivative has been presented very recently in [25]. This new definition is a natural extension of the usual derivative, and it satisfies some similar properties to the integer order calculus such as derivative of the product of two functions, derivative of the quotient of two functions and the chain rule. In [26] the author developed further the definitions and properties of conformable fractional derivative and integral.
conformable fractional derivative: Given a function f:[0,∞)→R, the conformable fractional derivative of f of order α∈(0,1] is defined by
Dαf(t)=limε→0f(t+εt1−α)−f(t)ε, t>0. |
If there exists conformable fractional derivative of f of order α∈(0,1] in some (0,a),a>0, and limt→0+Dαf(t) exists, then define Dαf(0)=limt→0+Dαf(t). Evidently, Dαf(t)=f′(t) for α=1.
conformable fractional integral: Given a function f:[0,∞)→R, the conformable fractional integral of f of order α∈(0,1] is defined by
Iαf(t)=∫t0τα−1f(τ)dτ. |
The conformable sequential fractional derivative is proposed in [26]. Given a function f:[0,∞)→R and α∈(0,1], the conformable sequential fractional derivative of f of order n is expressed by
Dnαf(t)=DαDα…Dα⏟n−timesf(t), t>0. |
The physical and geometrical meaning of the conformable fractional derivative is interpreted in [27]. Several applications of the definition have shown the significance of conformable fractional derivative. For example, [28] discussed the potential conformable quantum mechanics, [29] discussed the conformable Maxwell equations, and [30,31] showed that the conformable fractional derivative models present good agreements with experimental data. For recent results on the existence, stability and oscillation of solutions for conformable fractional differential equations, we refer the reader to [32,33,34,35,36,37,38,39,40,41,42,43,44].
The existence of solutions for periodic boundary value problem of impulsive conformable fractional integro-differential equation
{tkDαx(t)=f(t,x,(Fx)(t),(Sx)(t)), t∈J=[0,T]∖{t1,…,tm}, 0<α≤1,x(0)=x(T), x(t+k)−x(t−k)=Ik(x(tk)), k=,2,…,m, |
is studied in [33], where tkDα denotes the conformable fractional derivative of order α starting from tk, f∈C(J×R3,R), (Fx)(t)=∫t0l(t,s)x(s)ds, (Sx)(t)=∫T0h(t,s)x(s)ds, l∈C(D,R+),D={(t,s)∈J2:t≥s}, h∈C(J2,R+), Ik∈C(R,R). By the method of upper and lower solutions in reversed order coupled with the monotone iterative technique, some sufficient conditions for the existence of solutions are established.
In [42], applying the upper and lower solutions method and the monotone iterative technique, the authors investigated the existence of solutions to antiperiodic boundary value problem for impulsive conformable fractional functional differential equation
{tkDαx(t)=f(t,x,x(ω(t))), t∈J=[0,T]∖{t1,…,tm}, 0<α≤1,x(0)=−x(T), x(t+k)−x(t−k)=Ik(x(tk)), k=,2,…,m,x(t)=x(0), t∈[−r,0], |
where f∈C(J×R2,R), tkDα denotes the conformable fractional derivative of order α starting from tk, ω∈C(J,J+), J+=[−r,T],r>0, t−r≤ω(t),t∈J and tk<ω(t)≤t,t∈(tk,tk+1], Ik∈C(R,R).
However, to the best of our knowledge, the existence of minimal and maximal solutions and uniqueness of solution for fractional PBVP (1.1) involving conformable sequential fractional derivatives in ordered Banach spaces have not been considered up to now. Inspired by above works, we apply the theory of noncompactness measure and the method of upper and lower solutions coupled with the monotone iterative technique to construct two monotone iterative sequences, and then prove that the sequences converge to the extremal periodic solutions of PBVP (1.1), respectively, under some monotonicity conditions and noncompactness measure conditions of f. Also, we prove minimal and maximal solutions are equal and thus we obtain the uniqueness of solution. Moreover, we give the existence and uniqueness results of periodic solutions of the non-sequential fractional differential equation.
In this section, we introduce some notations, definitions and preliminary facts which are used throughout this paper. Let J=[0,1] and E be an ordered Banach space with the norm ‖⋅‖ and the partial order "≤", whose positive cone P={x∈E,x≥θ} is normal with normal constant L, where θ denotes the zero element of E. Generally, C(J,E) denotes the ordered Banach space of all continuous E-value functions on the interval J with the norm ‖x‖c=maxt∈J‖x(t)‖ and the partial order "≤" deduced by the positive cone Pc={x:x∈C(J,E),x(t)≥θ}. Pc is also normal with the same normal constant L. Let Cα(J,E)={x:x∈C(J,E),Dαx∈C(J,E)}, evidently, Cα(J,E) also is a Banach space with the norm ‖x‖α=max{‖x‖c,‖Dαx‖c}. A function x∈Cα(J,E) is called a solution of PBVP (1.1) if it satisfies the equation and the boundary conditions in (1.1).
Now let us recall some fundamental facts of the notion of Kuratowski measure of noncompactness.
Definition 2.1. ([45]) Let E be a Banach space and let ΩE be the family of bounded subsets of E. The Kuratowski measure of noncompactness is the map μ:ΩE→[0,∞] defined by
μ(B)=inf{ε>0:B⊆∪ni=1Bi,diam(Bi)≤ε},here B∈ΩE. |
Property 2.2. The Kuratowski measure of noncompactness satisfies some properties (for more details see [45]).
(1) μ(B)=0⇔¯B is compact (B is relatively compact), where ¯B denotes the closure of B.
(2) μ(A∪B)=max{μ(A),μ(B)}.
(3) μ(cB)=|c|μ(B),c∈R.
Denote the Kuratovski noncompactness measures of bounded sets in C(J,E) and Cα(J,E) by μc and μcα, respectively. For any H⊂C(J,E) and t∈J, set H(t)={x(t):x∈H}⊂E. If H is bounded in C(J,E), then H(t) is bounded in E, and μ(H(t))≤μc(H). Furthermore, we have the following well-known result.
Lemma 2.3. ([45]) Let H⊂C(J,E) be bounded and equicontinuous. Then μ(H(t)) is continuous on J and μc(H)=maxt∈Jμ(H(t)).
We can deduce the following useful result by Lemma 2.3.
Lemma 2.4. Let H⊂Cα(J,E) be bounded and equicontinuous. Then
μcα(H)=max{maxt∈Jμ(H(t)),maxt∈Jμ(DαH(t))}, |
where DαH(t)={Dαx(t)|x∈H}⊂E, t∈J.
Proof. Firstly, we prove that μcα(H)≤d=:max{maxt∈Jμ(H(t)),maxt∈Jμ(DαH(t))}. Noting that H⊂C(J,E) and DαH⊂C(J,E) are bounded and equicontinuous, by Lemma 2.3, we know
μc(H)=maxt∈Jμ(H(t))≤d, μc(DαH)=maxt∈Jμ(DαH(t))≤d. |
Therefore, there exist V1,V2,…,Vn⊂H and W1,W2,…,Wm⊂H such that H=⋃ni=1Vi=⋃mj=1Wj and
diamcVi<d+ε, diamcDαWj<d+ε, i=1,2,…,n,j=1,2,…,m, | (2.1) |
where diamc(⋅) denotes the diameter of the bounded subset of C(J,E). At the same time, for any x1,x2∈Vi, by (2.1) we obtain
‖x1(t)−x2(t)‖≤d+ε. | (2.2) |
Similarly, for y1,y2∈Wj, we have
‖Dαy1(t)−Dαy2(t)‖≤d+ε. | (2.3) |
Let Hij={x:x∈Vi,Dαx∈DαWj}, i=1,2,…,n,j=1,2,…,m. According to (2.1), (2.2) and (2.3), we can get
diamcHij≤d+ε, diamcDαHij≤d+ε, |
this means diamcαHij≤d+ε,i=1,2,…,n,j=1,2,…,m. Then it follows from H=⋃i=1,2,⋯,nj=1,2,⋯,mHij that μcα(H)≤d.
On the other hand, for any ε>0, there exist Hi⊂H,i=1,2,…,k such that H=⋃ki=1Hi and diamcα(Hi)≤μcα(H)+ε. Hence, for any t∈J and any x1,x2∈Hi, i=1,2,…,k, we have
‖x1(t)−x2(t)‖≤‖x1−x2‖c≤‖x1−x2‖cα≤μcα(H)+ε. | (2.4) |
From the fact H(t)=⋃ki=1Hi(t) and (2.4), we have μ(H(t))≤μcα(H)+ε for any t∈J, so maxt∈Jμ(H(t))≤μcα(H)+ε. Thus, maxt∈Jμ(H(t))≤μcα(H) since ε>0 is arbitrary. Similarly, it follows that maxt∈Jμ(DαH(t))≤μcα(H). Consequently, the proof of Lemma 2.4 is completed.
The following lemma also will be used in the proof of our main results.
Lemma 2.5. ([46,47,48]) Let E be a Banach space, H={xn}⊂L(J,E) be a countable set with ‖xn(t)‖≤ρ(t) for a.a. t∈J and every xn∈H, where ρ(t)∈L(J). Then μ(H(t)) is Lebesgue integrable on J, and μ({∫Jxn(t)dt})≤2∫Jμ(H(t))dt.
Conformable calculus satisfies the following properties.
Property 2.6. ([25]) Let α∈(0,1] and f,g be conformable differentiable of order α. Then
(1) DαC=0 for all constant functions f(t)=C.
(2) Dα(af+bg)=aDα(f)+bDα(g) for all a,b∈R.
(3) Dα(fg)=fDα(g)+gDα(f).
(4) DαIαf(t)=f(t) for t>0, where f is any continuous function in the domain of Iα.
Property 2.7. By the definition of conformable fractional integral, it is easy to know Iαf:C(J)→C(J) for α∈(0,1].
Remark 2.8 ([26,39,41]) The function exp(λtαα)=eλtαα, t≥0, is called fractional conformable exponential function, where α∈(0,1] and λ∈R. The conformable derivative of fractional conformable exponential function is Dαeλtαα=λeλtαα.
Lemma 2.9. ([35,39])For 0<α≤1, the general solution of the fractional nonhomogeneous equation
Dαx(t)−λx(t)=σ(t),t>0, |
is expressed by
x(t)=eλtαα[C+Iα(e−λtαασ(t))], |
where C is a constant. Further, the unique solution of the linear initial value problem
{Dαx(t)−λx(t)=σ(t), t∈(0,1],x(0)=x0, |
has the following form
x(t)=eλtαα[x0+Iα(e−λtαασ(t))]. |
From Lemma 2.9 we can obtain the solution of linear boundary value problem.
Lemma 2.10. For σ(t)∈C(J) and 0<α≤1, the unique solution x∈C(J) of the linear boundary value problem
{Dαx(t)−λx(t)=σ(t), t∈(0,1], λ≠0,x(0)−x(1)=d, |
has the following form
x(t)=deλtαα1−eλα+∫10Gλ(t,s)σ(s)ds, | (2.5) |
where
Gλ(t,s)={eλαeλtααe−λsααsα−11−eλα+eλtααe−λsααsα−1,s≤t,eλαeλtααe−λsααsα−11−eλα,s≥t. |
The continuity of the solution x(t) in Lemma 2.10 is guaranteed by Property 2.7.
Remark 2.11. It is easy to know from the expression of Gλ(t,s) that for λ<0,
0≤Gλ(t,s)≤{(eλα1−eλα+1)sα−1=11−eλαsα−1,s≤t,11−eλαsα−1,s≥t. |
Thus ∫10Gλ(t,s)ds≤1α(1−eλα).
Furthermore, from Lemma 2.10 we can deduce the solution of linear PBVP with sequential derivative.
Lemma 2.12. For σ(t)∈C(J), 0<α≤1 and M,N>0 such that M2≥4N, the solution of the linear PBVP
{D2αx(t)+MDαx(t)+Nx(t)=σ(t), t∈(0,1],x(0)=x(1),Dαx(0)=Dαx(1), | (2.6) |
is
x(t)=∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(τ)dτ)ds, | (2.7) |
where
λ1=−M−√M2−4N2≤λ2=−M+√M2−4N2<0. |
Proof. Let (Dα−λ2)x(t)=y(t),t∈(0,1]. Then the problem (2.6) is equivalent to
{(Dα−λ2)x(t)=y(t),t∈(0,1],x(0)=x(1), | (2.8) |
and
{(Dα−λ1)y(t)=σ(t),t∈(0,1],y(0)=y(1). | (2.9) |
By Lemma 2.10, we obtain that the problems (2.8) and (2.9) have the following representation of solutions
x(t)=∫10Gλ2(t,s)y(s)ds, | (2.10) |
and
y(t)=∫10Gλ1(t,s)σ(s)ds, | (2.11) |
respectively. Substituting (2.11) into (2.10), we get (2.7).
The following comparison result plays an important role in the proofs of our main results.
Lemma 2.13. Let M2≥4N, M,N>0 and
{D2αx(t)+MDαx(t)+Nx(t)≥0, t∈(0,1],x(0)≥x(1),Dαx(0)≥Dαx(1). | (2.12) |
Then x(t)≥0 on J.
Proof. Rewrite (2.12) as follows
{D2αx(t)+MDαx(t)+Nx(t)=σ(t)≥0, t∈(0,1],x(0)−x(1)=k≥0,Dαx(0)−Dαx(1)=l≥0. | (2.13) |
Then the problem (2.13) is equivalent to
{(Dα−λ2)x(t)=y(t),t∈(0,1],x(0)−x(1)=k, |
and
{(Dα−λ1)y(t)=σ(t),t∈(0,1],y(0)−y(1)=l−λ2k. |
In view of λ2<0,l−λ2k≥0, σ(t)≥0, λ1<0 and (2.5), we can obtain y(t)≥0. Furthermore, λ2<0,k≥0, y(t)≥0 and (2.5) ensure x(t)≥0 on J.
Definition 3.1. Let v,w∈Cα(J,E). v is called a lower solution of PBVP (1.1) if it satisfies
{D2αv(t)≤f(t,v(t),Dαv(t)), t∈(0,1],v(0)≤v(1), Dαv(0)≤Dαv(1), | (3.1) |
and w is called an upper solution of PBVP (1.1) if it satisfies
{D2αw(t)≥f(t,w(t),Dαw(t)), t∈(0,1],w(0)≥w(1), Dαw(0)≥Dαw(1). | (3.2) |
In the following, we assume that v(t)≤w(t),t∈J. Define the ordered interval in space C(J,E)
[v,w]={x(t)∈C(J,E),v(t)≤x(t)≤w(t), t∈J}. |
Denote D1(t)=Dαv(t)+λ2(w(t)−v(t)) and D2(t)=Dαw(t)−λ2(w(t)−v(t)). We work with the following conditions on the function f in (1.1).
(H1) There exist constants M,N>0 with M2≥4N such that
f(t,w,Dαw)−f(t,v,Dαv)≥−M(Dαw−Dαv)−N(w−v), t∈J, |
where v,w∈Cα(J,E) are lower and upper solutions of (1.1).
(H2) There exist constants M,N>0 with M2≥4N such that
f(t,x2,y2)−f(t,x1,y1)≥−M(y2−y1)−N(x2−x1), t∈J, |
where v≤x1≤x2≤w and D1(t)≤yi(t)≤D2(t), i=1,2.
(H3) There exists a constant K≥0 such that
μ({f(t,xn,yn)+Nxn(t)+Myn(t)})≤K(μ({xn(t)})+μ({yn(t)})), t∈J, |
for any monotonic sequence {xn}⊂[v,w] and any sequence {yn} such that D1(t)≤yn(t)≤D2(t). Moreover,
8K(1−λ2)α21(1−eλ1α)(1−eλ2α)<1. |
(H4) There exist constants ¯M,¯N>0 such that
f(t,x2,y2)−f(t,x1,y1)≤¯M(y2−y1)+¯N(x2−x1), t∈J, |
where v≤x1≤x2≤w and D1(t)≤yi(t)≤D2(t), i=1,2. Moreover,
8¯K(1−λ2)α21(1−eλ1α)(1−eλ2α)<1, |
here ¯K=max{L(¯M+M),L(¯N+N)} and L is the normal constant of cone P.
Remark 3.2. The condition (H1) and Lemma 2.10 guarantee that Dα(w(t)−v(t))−λ2(w(t)−v(t))≥θ,t∈J. In fact, let z(t)=Dα(w(t)−v(t))−λ2(w(t)−v(t)), then by (3.1), (3.2) and (H1), for t∈(0,1],
Dαz(t)−λ1z(t)=D2α(w(t)−v(t))+MDα(w(t)−v(t))+N(w(t)−v(t))≥(f(t,w,Dαw)−f(t,v,Dαv))+MDα(w(t)−v(t))+N(w(t)−v(t))≥θ, |
and z(0)−z(1)≥θ. By (2.5) of Lemma 2.10, x(t)≥θ if d≥θ and σ(t)≥θ. Thus we obtain that z(t)≥θ, this implies D1(t)≤Dαv(t)≤D2(t) and D1(t)≤Dαw(t)≤D2(t).
Set Ω={η∈[v,w]∩Cα(J,E):D1(t)≤Dαη(t)≤D2(t)}. Obviously, Ω is well defined by Remark 3.2 when the condition (H1) holds. Now we are in the position to state our main results.
Theorem 3.3. Assume that f∈C([0,1]×E×E), v and w are lower and upper solutions of BPVP (1.1) and the conditions (H1), (H2) and (H3) are valid. Then there exist p(t),q(t)∈Cα(J,E) such that p(t),q(t) are minimal and maximal solutions on the ordered interval [v,w] for BPVP (1.1), respectively, that is, for any solution x(t) of BPVP (1.1) such that x∈[v,w], we have v(t)≤p(t)≤x(t)≤q(t)≤w(t),t∈J.
Proof. Let σ(η)(t)=f(t,η(t),Dαη(t))+MDαη(t)+Nη(t). For any η∈Ω, consider the linear PBVP
{D2αx(t)+MDαx(t)+Nx(t)=σ(η)(t), t∈(0,1],x(0)=x(1), Dαx(0)=Dαx(1). | (3.3) |
By Lemma 2.12, (3.3) has exactly one solution x(t) given by
x(t)=∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(η)(τ)dτ)ds. | (3.4) |
For any η∈Ω, define the operator T as Tη(t)=x(t), then the fixed points of T are exactly the solutions of PBVP (1.1). For clarity, we divide the proof into several steps.
Step 1: Firstly, T:Ω→Cα(J,E) is well defined. Indeed, for any η∈Ω, by (3.4),
Tη(t)=∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(η)(τ)dτ)ds=Aηeλ2tαα+Bηeλ2tααIαe(λ1−λ2)tαα+eλ2tααIα(e(λ1−λ2)tααIαe−λ1tαασ(η)(t)), | (3.5) |
where
Aη=eλ1+λ2α(1−eλ1α)(1−eλ2α)(Iαe−λ1tαασ(η)(t))|t=1(Iαe(λ1−λ2)tαα)|t=1+eλ2α1−eλ2αIα(e(λ1−λ2)tααIαe−λ1tαασ(η)(t))|t=1, |
and
Bη=eλ1α1−eλ1α(Iαe−λ1tαασ(η)(t))|t=1. |
Furthermore, Property 2.6, Property 2.7 and Remark 2.8 indicate for t>0,
DαTη(t)=Aηλ2eλ2tαα+Bηλ2eλ2tααIαe(λ1−λ2)tαα+λ2eλ2tααIα(e(λ1−λ2)tααIαe−λ1tαασ(η)(t))+Bηeλ1tαα+eλ1tααIαe−λ1tαασ(η)(t)=λ2∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(η)(τ)dτ)ds+∫10Gλ1(t,s)σ(η)(s)ds. | (3.6) |
The continuity of f together with Property 2.7 ensures that the right side of (3.6) is continuous on J. Hence, (3.6) is also valid for t=0. By Property 2.7, (3.5) and (3.6) we have Tη(t)∈Cα(J,E).
Step 2: For any η∈Ω, {Tη}⊂Cα(J,E) is bounded. In fact, from condition (H2) we can derive
f(t,v,Dαv)+MDαv+Nv≤f(t,η,Dαη)+MDαη+Nη≤f(t,w,Dαw)+MDαw+Nw. |
In view of the normality of the cone P, there exists ¯L>0 such that ‖f(t,η,Dαη)+MDαη+Nη‖≤¯L. Since λ1<λ2<0, by Remark 2.11, (3.5) and (3.6) we can find
‖Tη(t)‖≤¯Lα21(1−eλ1α)(1−eλ2α), |
and
‖DαTη(t)‖≤|λ2|¯Lα21(1−eλ1α)(1−eλ2α)+¯Lα11−eλ1α. |
Hence, {Tη}⊂Cα(J,E) is bounded.
Step 3: For any η∈Ω, {Tη}⊂Cα(J,E) is equicontinuous. Noticing that
|eλ2tα1α∫t10sα−1e(λ1−λ2)sααds−eλ2tα2α∫t20sα−1e(λ1−λ2)sααds|≤|eλ2tα1α∫t10sα−1e(λ1−λ2)sααds−eλ2tα2α∫t10sα−1e(λ1−λ2)sααds|+|eλ2tα2α∫t10sα−1e(λ1−λ2)sααds−eλ2tα2α∫t20sα−1e(λ1−λ2)sααds|≤1α|eλ2tα1α−eλ2tα2α|+|∫t2t1sα−1e(λ1−λ2)sααds|≤1α|eλ2tα1α−eλ2tα2α|+1α|tα2−tα1|, | (3.7) |
‖eλ2tα1α∫t10sα−1e(λ1−λ2)sαα(∫s0τα−1e−λ1ταασ(η)(τ)dτ)ds−eλ2tα2α∫t20sα−1e(λ1−λ2)sαα(∫s0τα−1e−λ1ταασ(η)(τ)dτ)ds‖≤|eλ2tα1α−eλ2tα2α|‖∫t10sα−1e(λ1−λ2)sαα(∫s0τα−1e−λ1ταασ(η)(τ)dτ)ds‖+eλ2tα2α‖∫t2t1sα−1e(λ1−λ2)sαα(∫s0τα−1e−λ1ταασ(η)(τ)dτ)ds‖≤¯Lα2e−λ1α|eλ2tα1α−eλ2tα2α|+¯Lα2e−λ1α|tα1−tα2|, | (3.8) |
and
‖eλ1tα1α∫t10sα−1e−λ1sαασ(η)(s)ds−eλ1tα2α∫t20sα−1e−λ1sαασ(η)(s)ds‖≤|eλ1tα1α−eλ1tα2α|‖∫t10sα−1e−λ1sαασ(η)(s)ds‖+eλ1tα2α‖∫t2t1sα−1e−λ1sαασ(η)(s)ds‖≤¯Lαe−λ1α|eλ1tα1α−eλ1tα2α|+¯Lαe−λ1α|tα1−tα2|, | (3.9) |
we can arrive at from (3.5)–(3.9) that {Tη}⊂Cα(J,E) is equicontinuous.
Step 4: Now we prove v(t)≤Tv(t), w(t)≥Tw(t). Set v1=Tv, then we have
{D2αv1(t)+MDαv1(t)+Nv1(t)=f(t,v(t),Dαv(t)) +MDαv(t)+Nv(t), t∈(0,1],v1(0)=v1(1), Dαv1(0)=Dαv1(1). | (3.10) |
Set z(t)=v1(t)−v(t), then by (3.1) and (3.10) we can know
{D2αz(t)+MDαz(t)+Nz(t)≥θ, t∈(0,1],z(0)≥z(1), Dαz(0)≥Dαz(1). |
Thus, z(t)≥θ by Lemma 2.13 and this means v(t)≤Tv(t). Similarly, w(t)≥Tw(t).
Next, let η1,η2∈Ω be such that η1(t)≤η2(t). (H2) and (3.5) show that σ(η1)(t)≤σ(η2)(t), Tη1(t)≤Tη2(t). For each η∈Ω, we get by the similar method to Remark 3.2 that Dα(Tη(t)−v(t))−λ2(Tη(t)−v(t))≥θ,t∈J. This implies
D1(t)≤Dαv(t)+λ2(Tη(t)−v(t))≤DαTη(t). |
Similarly, Dα(w(t)−Tη(t))−λ2(w(t)−Tη(t))≥θ,t∈J, and hence
DαTη(t)≤Dαw(t)−λ2(w(t)−Tη(t))≤D2(t). |
Therefore, T(Ω)⊂Ω and T is a monotone operator on [v,w]. Consequently, let vn=Tvn−1,wn=Twn−1(n=1,2,…), we obtain
v=v0≤v1≤…≤vn≤…≤wn…≤w1≤w0=w, |
and
D1(t)≤DαTvn(t),DαTwn(t)≤D2(t), t∈J. |
Step 5: Let V={vn:n=1,2,…}. In the following, we will show that V is a relatively compact set in Cα(J,E). Note that μc(V)=μc(V∪{v0}) and μc(DαV)=μc(DαV∪{Dαv0}) by Property 2.2. In view of Lemma 2.5, Property 2.2, condition (H3), Remark 2.11 and Lemma 2.3, it follows that
μ({∫10Gλ1(t,s)σ(vn−1)(s)ds})≤2∫10Gλ1(t,s)μ({σ(vn−1)(s)})ds≤2K∫10Gλ1(t,s)(μ({vn−1(s)})+μ({Dαvn−1(s)}))ds≤2Kα11−eλ1α(μc(V)+μc(DαV)). | (3.11) |
Furthermore,
μ(V(t))=μ({∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(vn−1)(τ)dτ)ds})≤2∫10Gλ2(t,s)μ({∫10Gλ1(s,τ)σ(vn−1)(τ)dτ})ds≤4Kα21(1−eλ1α)(1−eλ2α)(μc(V)+μc(DαV)). | (3.12) |
Hence, we obtain from Lemma 2.4
μ(V(t))≤8Kα21(1−eλ1α)(1−eλ2α)μcα(V). | (3.13) |
At the same time, by (3.6), (3.11), (3.12) and (3.13), we get
μ(DαV(t))≤|λ2|8Kα21(1−eλ1α)(1−eλ2α)μcα(V)+4Kα11−eλ1αμcα(V)≤8K(1−λ2)α21(1−eλ1α)(1−eλ2α)μcα(V). | (3.14) |
Combining Lemma 2.4, (3.13) and (3.14), we have
μcα(V)≤8K(1−λ2)α21(1−eλ1α)(1−eλ2α)μcα(V). |
Therefore, μcα(V)=0 by (H3), then Property 2.2 shows {vn} is a relatively compact set of Cα(J,E), and thus there exists subsequence converging uniformly to p∈Cα(J,E). By the monotone property of {vn} and the assumption (H2) of function f which implies {σ(vn)} is monotone, we obtain that {vn} converges uniformly to p∈Cα(J,E) and {σ(vn)} converges to σ(p). Similarly, we can show that {wn} converges uniformly to q∈Cα(J,E) and {σ(wn)} converges to σ(q). Moreover, the limits p,q satisfy
v=v0≤v1≤…≤vn≤p≤q≤wn…≤w1≤w0=w, |
and
D1(t)≤Dαp(t),Dαq(t)≤D2(t), t∈J. |
Step 6: Let n→∞ in the relation
vn(t)=Tvn−1(t)=∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(vn−1)(τ)dτ)ds. |
Applying the dominated convergence theorem, we have p satisfies the equation
p(t)=∫10Gλ2(t,s)(∫10Gλ1(s,τ)σ(p)(τ)dτ)ds, t∈J, |
which implies that p(t) is an integral representation of the solution to the problem (3.3), and thus by the definition of function σ, p is a solution of PBVP (1.1). Similarly, we can prove that q is a solution of PBVP (1.1).
Finally, we prove p and q are extremal solutions of PBVP (1.1). Assume that ¯x is a fixed point of T in Ω, then by the monotonicity of T proved in Step 4, it is easy to see that Tv≤T¯x≤Tw, that is, v1≤¯x≤w1. Furthermore, we have vn≤¯x≤wn for n=1,2,…. Let n→∞ we get p≤¯x≤q. Therefore, p and q are the minimal and maximal solutions of PBVP (1.1) in [v,w], respectively.
This completes the proof of Theorem 3.3.
Theorem 3.4. Assume that f∈C([0,1]×E×E), v and w are lower and upper solutions of BPVP (1.1) and the conditions (H1), (H2) and (H4) are valid. Then the BPVP (1.1) has a unique solution x∈[v,w].
Proof. First of all, we prove that (H2) and (H4) imply (H3). Let {xn}⊂[v,w] be increasing sequence and {yn} be such that {yn}⊂[D1(t),D2(t)], t∈J. For m,n∈N with m>n, in view of (H2) and (H4), we have
θ≤f(t,xm,ym)−f(t,xn,yn)+M(ym(t)−yn(t))+N(xm(t)−xn(t))≤(¯M+M)(ym(t)−yn(t))+(¯N+N)(xm(t)−xn(t)). | (3.15) |
By (3.15) and the normality of cone P, we can derive that
‖f(t,xm,ym)−f(t,xn,yn)+M(ym(t)−yn(t))+N(xm(t)−xn(t))‖≤L‖(¯M+M)(ym(t)−yn(t))+(¯N+N)(xm(t)−xn(t))‖≤L(¯M+M)‖ym(t)−yn(t)‖+L(¯N+N)‖xm(t)−xn(t)‖≤¯K(‖ym(t)−yn(t)‖+‖xm(t)−xn(t)‖). | (3.16) |
From (3.16) and the definition of Kuratowski measure of noncompactness, it follows that
μ({f(t,xn,yn)+Nxn(t)+Myn(t)})≤¯K(μ({xn(t)})+μ({yn(t)})), t∈J. |
If {xn}⊂[v,w] is a decreasing sequence and {yn}⊂[D1(t),D2(t)], the above inequality is also valid. Hence (H3) is satisfied. Therefore, Theorem 3.3 asserts that the PBVP (1.1) has minimal and maximal solutions p and q between v and w. In the following we show that p=q.
Since by (H2) and (H4),
θ≤q(t)−p(t)=∫10Gλ2(t,s)(∫10Gλ1(s,τ)[f(τ,q(τ),Dαq(τ))−f(τ,p(τ),Dαp(τ))+M(Dαq(τ)−Dαp(τ))+N(q(τ)−p(τ))]dτ)ds≤∫10Gλ2(t,s)(∫10Gλ1(s,τ)[(¯M+M)(Dαq(τ)−Dαp(τ))+(¯N+N)(q(τ)−p(τ))]dτ)ds, | (3.17) |
(3.17) together with the normality of cone P ensures
‖q(t)−p(t)‖≤∫10Gλ2(t,s)(∫10Gλ1(s,τ)[L(¯M+M)‖Dαq(τ)−Dαp(τ)‖+L(¯N+N)‖q(τ)−p(τ)‖]dτ)ds≤[L(¯M+M)+L(¯N+N)]1α21(1−eλ1α)(1−eλ2α)‖q−p‖Cα. | (3.18) |
Furthermore, if Dαq(t)≥Dαp(t), again by (H2) and (H4), we have
θ≤Dαq(t)−Dαp(t)=λ2∫10Gλ2(t,s)(∫10Gλ1(s,τ)[f(τ,q(τ),Dαq(τ))−f(τ,p(τ),Dαp(τ))+M(Dαq(τ)−Dαp(τ))+N(q(τ)−p(τ))]dτ)ds+∫10Gλ1(t,s)[f(s,q(s),Dαq(s))−f(s,p(s),Dαp(s))+M(Dαq(s)−Dαp(s))+N(q(s)−p(s))]ds≤∫10Gλ1(t,s)[(¯M+M)(Dαq(s)−Dαp(s))+(¯N+N)(q(s)−p(s))]ds. | (3.19) |
From (3.19) and the normality of cone P we get
‖Dαq(t)−Dαp(t)‖≤∫10Gλ1(t,s)[L(¯M+M)‖Dαq(s)−Dαp(s)‖+L(¯N+N)‖q(s)−p(s)‖]ds≤[L(¯M+M)+L(¯N+N)]1α11−eλ1α‖q−p‖Cα. | (3.20) |
On the other hand, if Dαq(t)≤Dαp(t), also by (H2) and (H4) we get
θ≤Dαp(t)−Dαq(t)=λ2∫10Gλ2(t,s)(∫10Gλ1(s,τ)[f(τ,p(τ),Dαp(τ))−f(τ,q(τ),Dαq(τ))+M(Dαp(τ)−Dαq(τ))+N(p(τ)−q(τ))]dτ)ds+∫10Gλ1(t,s)[f(s,p(s),Dαp(s))−f(s,q(s),Dαq(s))+M(Dαp(s)−Dαq(s))+N(p(s)−q(s))]ds≤λ2∫10Gλ2(t,s)(∫10Gλ1(s,τ)[(¯M+M)(Dαp(τ)−Dαq(τ))+(¯N+N)(p(τ)−q(τ))]dτ)ds. | (3.21) |
From (3.21) and the normality of cone P we know
‖Dαp(t)−Dαq(t)‖≤|λ2|∫10Gλ2(t,s)(∫10Gλ1(s,τ)[L(¯M+M)‖Dαp(τ)−Dαq(τ)‖+L(¯N+N)‖p(τ)−q(τ)‖]dτ)ds≤|λ2|[L(¯M+M)+L(¯N+N)]1α21(1−eλ1α)(1−eλ2α)‖p−q‖Cα. | (3.22) |
Consequently, from (3.18), (3.20) and (3.22) we can conclude
‖p−q‖Cα≤2¯K(1−λ2)α21(1−eλ1α)(1−eλ2α)‖p−q‖Cα. |
Thus, p=q by (H4), which means that there exists a unique solution of PBVP (1.1) in [v,w].
This completes the proof of Theorem 3.4.
Remark 3.5. Using the methods of our main results Theorem 3.3 and Theorem 3.4, we can easy to obtain the existence of solutions of the following PBVP for fractional differential equation
{Dαx(t)=f(t,x(t)), t∈(0,1], 0<α≤1,x(0)=x(1). | (3.23) |
Let v,w∈C(J,E). We say that the function v is a lower solution of problem (3.23) if
{Dαv(t)≤f(t,v(t)), t∈(0,1],v(0)≤v(1). |
Analogously, w is an upper solution for problem (3.23) if it verifies similar conditions for the inequalities reversed.
(G1) There exists a constant N>0 such that f(t,x2)−f(t,x1)≥−N(x2−x1),t∈J, where v≤x1≤x2≤w.
(G2) There exists a constant K≥0 such that μ({f(t,xn)+Nxn(t)})≤Kμ({xn(t)}),t∈J for any monotonic sequence {xn}⊂[v,w]. Moreover, 2Kα11−e−Nα<1.
(G3) There exists a constant ¯N>0 such that f(t,x2)−f(t,x1)≤¯N(x2−x1),t∈J, where v≤x1≤x2≤w. Moreover, 2¯Kα11−e−Nα<1, where ¯K=L(¯N+N) and L is the normal constant of cone P.
Theorem 3.6. Assume that f∈C(J×E), v,w are lower and upper solutions of BPVP (3.23) and v≤w. The conditions (G1) and (G2) are valid. Then there exist p(t),q(t)∈C(J,E) such that p(t),q(t) are minimal and maximal solutions on the ordered interval [v,w] for BPVP (3.23), respectively, that is, for any solution x(t) of BPVP (3.23) such that x∈[v,w], we have v(t)≤p(t)≤x(t)≤q(t)≤w(t),t∈J.
Theorem 3.7. Assume that f∈C(J×E), v,w are lower and upper solutions of BPVP (3.23), v≤w, and the conditions (G1) and (G3) hold. Then the BPVP (3.23) has a unique solution x∈[v,w].
Let E={x:x=(x1,x2,…,xn,…),xn→0} with the norm ‖x‖=supn|xn| and P={x∈E:xn≥0,n=1,2,3,…}. Then P is a normal cone in E. Consider the PBVP of infinite system for differential equations in E
{D2αxn(t)=18ent[(1−xn(t))2−2n2−12n2] +(n−1)(2n+1)34n2entsin3x2n+1(t)−Dαxn(t), t∈(0,1],xn(0)=xn(1), Dαxn(0)=Dαxn(1). | (4.1) |
Evidently, (4.1) can be regarded as a PBVP of the form (1.1) in E. In this situation, x=(x1,x2,…,xn,…), y=(y1,y2,…,yn,…) and f=(f1,f2,…,fn,…), in which
fn(t,x,y)=18ent[(1−xn(t))2−2n2−12n2]+(n−1)(2n+1)34n2entsin3x2n+1(t)−yn(t). |
It is clear f∈C(J×E×E,E). Let v=(0,0,…,0,…) and w=(1,12,…,1n,…). Then v,w∈Cα(J,E),v(t)≤w(t),t∈J, and
vn(0)=vn(1)=Dαvn(0)=Dαvn(1)=0,wn(0)=wn(1)=1n,Dαwn(0)=Dαwn(1)=0. |
Moreover,
fn(t,v(t),Dαv(t))=18ent(1−2n2−12n2)=18ent12n2>0, |
and
fn(t,w(t),Dαw(t))=18ent[(1−1n)2−2n2−12n2]+(n−1)(2n+1)34n2entsin312n+1≤18ent[(1−1n)2−2n2−12n2]+18ent2n−2n2=−18ent12n2<0. |
Hence, v and w are lower and upper solutions of (4.1).
The conditions (H1) and (H2) are satisfied with M=1,N=14. Let us first verify the condition (H2). Noticing that λ2=−12,D1(t)=−w(t)2 and D2(t)=w(t)2, for v(t)≤x(1)(t)≤x(2)(t)≤w(t) and D1(t)≤y(i)(t)≤D2(t), i=1,2, we have
0≤x(1)n(t)≤x(2)n(t)≤1n, −12n≤y(i)n(t)≤12n,i=1,2, n=1,2,3,…. |
Therefore,
fn(t,x(2),y(2))−fn(t,x(1),y(1))=18ent[(1−x(2)n(t))2−(1−x(1)n(t))2]−(y(2)n(t)−y(1)n(t))+(n−1)(2n+1)34n2ent[sin3x(2)2n+1(t)−sin3x(1)2n+1(t)]≥18ent[(1−x(2)n(t))2−(1−x(1)n(t))2]−(y(2)n(t)−y(1)n(t))=−18ent(2−x(1)n(t)−x(2)n(t))(x(2)n(t)−x(1)n(t))−(y(2)n(t)−y(1)n(t))≥−14ent(x(2)n(t)−x(1)n(t))−(y(2)n(t)−y(1)n(t))≥−14(x(2)n(t)−x(1)n(t))−(y(2)n(t)−y(1)n(t)). |
This implies that (H2) is satisfied. Obviously, by the same method we can verify (H1).
Finally, we check condition (H3). Let the sequences {x(m):x(m)=(x(m)1,x(m)2,…,x(m)n,…)} and {y(m):y(m)=(y(m)1,y(m)2,…,y(m)n,…)} be given such that {x(m)}⊂[v,w] is monotonous and −12n≤y(m)n(t)≤12n,n=1,2,3,…. Let z(m)n(t)=fn(t,x(m),y(m))+Nx(m)n(t)+My(m)n(t). In view of
|z(m)n(t)|≤184n−32n2+n−14n2+12n+Nn+M2n, t∈J, n,m=1,2,3,…, |
it follows that {z(m)n(t)} is bounded, so we can choose a subsequence {mi}⊂{m} such that z(mi)n(t)→zn(t) as i→∞, n=1,2,3,… and
|zn(t)|≤184n−32n2+n−14n2+12n+Nn+M2n, t∈J, n=1,2,3,…. |
Hence, z(t)=(z1(t),z2(t),z3(t),…)∈E for any t∈J, and it is easy to see that
‖f(t,x(mi),y(mi))+Nx(mi)(t)+My(mi)(t)−z(t)‖=supn|z(mi)n(t)−zn(t)|→0, i→∞. |
Consequently, we conclude that condition (H3) is satisfied for K=0. Therefore, Theorem 3.3 ensures that PBVP (4.1) has extremal solutions in [v,w], which can be obtained by taking limits from the iterative sequences {v(m):v(m)=(v(m)1,v(m)2,…,v(m)n,…)} and {w(m):w(m)=(w(m)1,w(m)2,…,w(m)n,…)}, here
v(0)=v,w(0)=w, |
v(m)n(t)=∫10Gλ2(t,s){∫10Gλ1(s,τ)[fn(τ,v(m−1)(τ),Dαv(m−1)(τ))+MDαv(m−1)n(τ)+Nv(m−1)n(τ)]dτ}ds=∫10G−12(t,s){∫10G−12(s,τ) [18enτ((1−v(m−1)n(τ))2−2n2−12n2)+(n−1)(2n+1)34n2enτsin3v(m−1)2n+1(τ)+14v(m−1)n(τ)]dτ}ds, |
and
w(m)n(t)=∫10Gλ2(t,s){∫10Gλ1(s,τ)[fn(τ,w(m−1)(τ),Dαw(m−1)(τ))+MDαw(m−1)n(τ)+Nw(m−1)n(τ)]dτ}ds=∫10G−12(t,s){∫10G−12(s,τ) [18enτ((1−w(m−1)n(τ))2−2n2−12n2)+(n−1)(2n+1)34n2enτsin3w(m−1)2n+1(τ)+14w(m−1)n(τ)]dτ}ds. |
This paper explores periodic solutions of some nonlinear fractional differential equations. The problem discussed involves sequential conformable fractional derivative. Under suitable monotonicity conditions and noncompactness measure conditions, the existence and uniqueness of solutions are derived from monotone iterative technique and upper and lower solutions method. Further, it is analyzed that the similar methods are well suited for investigating the existence and uniqueness of periodic solutions of the non-sequential fractional differential equation. In particular, for α=1, the classical results corresponding to ordinary differential equations of integer order are yielded. An example is given to illustrate an application of our theoretical work. To the best of our knowledge, the results obtained throughout this article are not recorded in any published literature. It is worth to be pointed out that the techniques applied in the main results of this paper can be used to investigate initial value problems or the differential equations with Riemann-Liouville sequential fractional derivatives.
The authors are grateful to the editor and referees for their valuable suggestions. This work is supported by the National Natural Science Foundation of China (11671181) and the Fundamental Research Funds for the Central Universities (2020YJSLX01).
The authors declare that they have no competing interests.
[1] | K. S. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley and Sons, New York, 1993. |
[2] | Z. L. Wei, Q. D. Li, J. L. Che, Initial value problems for fractional differential equations involving Riemann-Liouville sequential fractional derivative, J. Math. Anal. Appl., 367 (2010), 260-272. |
[3] | B. Ahmad, R. Luca, Existence of solutions for a sequential fractional integro-differential system with coupled integral boundary conditions, Chaos Soliton. Fract., 104 (2017), 378-388. |
[4] | B. Ahmad, J. J. Nieto, Sequential fractional differential equations with three-point boundary conditions, Comput. Math. Appl., 64 (2012), 3046-3052. |
[5] | B. Ahmad, J. J. Nieto, Boundary value problems for a class of sequential integrodifferential equations of fractional order, J. Funct. Space. Appl., 2013 (2013), 1-8. |
[6] | B. Ahmad, S. K. Ntouyas, Existence results for a coupled system of Caputo type sequential fractional differential equations with nonlocal integral boundary conditions, Appl. Math. Comput., 266 (2015), 615-622. |
[7] | Y. Chen, X. H. Tang, Solvability of sequential fractional order multi-point boundary value problems at resonance, Appl. Math. Comput., 218 (2012), 7638-7648. |
[8] | Q. D. Li, H. Su, Z. L. Wei, Existence and uniqueness result for a class of sequential fractional differential equations, J. Appl. Math. Comput., 38 (2012), 641-652. |
[9] | C. Promsakon, N. Phuangthong, S. K. Ntouyas, et al. Nonlinear sequential Riemann-Liouville and Caputo fractional differential equations with generalized fractional integral conditions, Adv. Differ. Equ., 2018 (2018), 1-17. |
[10] | H. Ye, R. Huang, Initial value problem for nonlinear fractional differential equations with sequential fractional derivative, Adv. Differ. Equ., 2015 (2015), 1-13. |
[11] | H. Ye, R. Huang, On the nonlinear fractional differential equations with Caputo sequential fractional derivative, Adv. Math. Phys., 2015 (2015), 1-9. |
[12] | O. Abu Arqub, Application of residual power series method for the solution of time-fractional Schrödinger equations in one-dimensional space, Fundam. Inf., 166 (2019), 87-110. |
[13] | O. Abu Arqub, Numerical algorithm for the solutions of fractional order systems of Dirichlet function types with comparative analysis, Fundam. Inf., 166 (2019), 111-137. |
[14] | O. Abu Arqub, N. Shawagfeh, Application of reproducing kernel algorithm for solving Dirichlet time-fractional diffusion-Gordon types equations in porous media, J. Porous Med., 22 (2019), 411-434. |
[15] | A. M. A. El-Sayed, H. H. G. Hashem, E. A. A. Ziada, Picard and Adomian decomposition methods for a quadratic integral equation of fractional order, Comput. Appl. Math., 33 (2014), 95-109. |
[16] | A. Lotfi, S. A. Yousefi, A numerical technique for solving a class of fractional variational problems, J. Comput. Appl. Math., 237 (2013), 633-643. |
[17] | P. Chen, Y. Li, Monotone iterative technique for a class of semilinear evolution equations with nonlocal conditions, Results Math., 63 (2013), 731-744. |
[18] | H. D. Gou, Y. X. Li, Upper and lower solution method for Hilfer fractional evolution equations with nonlocal conditions, Bound. Value Probl., 2019 (2019), 1-25. |
[19] | C. Hu, B. Liu, S. Xie, Monotone iterative solutions for nonlinear boundary value problems of fractional differential equations with deviating arguments, Appl. Math. Comput., 222 (2013), 72-81. |
[20] | F. A. Mcrae, Monotone iterative technique and existence results for fractional differential equations, Nonlinear Anal., 71 (2009), 6093-6096. |
[21] | J. Mu, Y. Li, Monotone interactive technique for impulsive fractional evolution equations, J. Inequal. Appl., 2011 (2011), 1-12. |
[22] | G. Wang, Monotone iterative technique for boundary value problems of a nonlinear fractional differential equation with deviating arguments, J. Comput. Appl. Math., 236 (2012), 2425-2430. |
[23] | L. Zhang, B. Ahmad, G. Wang, Explicit iterations and extremal solutions for fractional differential equation with nonlinear integral boundary conditions, Appl. Math. Comput., 268 (2015), 338-392. |
[24] | C. Z. Bai, Impulsive periodic boundary value problems for fractional differential equation involving Riemann-Liouville sequential fractional derivative, J. Math. Anal. Appl., 384 (2011), 211-231. |
[25] | R. Khalil, M. Al-Horani, A. Yousef, et al. A new definition of fractional derivative, J. Comput. Appl. Math., 264 (2014), 65-70. |
[26] | T. Abdeljawad, On conformable fractional calculus, J. Comput. Appl. Math., 279 (2015), 57-66. |
[27] | D. Zhao, M. Luo, General conformable fractional derivative and its physical interpretation, Calcolo, 54 (2017), 903-917. |
[28] | D. R. Anderson, D. J. Ulness, Properties of the Katugampola fractional derivative with potential application in quantum mechanics, J. Math. Phys., 56 (2015), 1-15. |
[29] | D. Zhao, X. Pan, M. Luo, A new framework for multivariate general conformable fractional calculus and potential applications, Phys. A, 15 (2018), 271-280. |
[30] | S. Yang, L. Wang, S. Q. Zhang, Conformable derivative: application to non-Darcian flow in lowpermeability porous media, Appl. Math. Lett., 79 (2018), 105-110. |
[31] | H. W. Zhou, S. Yang, S. Q. Zhang, Conformable derivative approach to anomalous diffusion, Phys. A, 491 (2018), 1001-1013. |
[32] | M. Al-Refai, T. Abdeljawad, Fundamental results of conformable Sturm-Liouville eigenvalue problems, Complexity, 2017 (2017), 1-7. |
[33] | S. Asawasamrit, S. K. Ntouyas, P. Thiramanus, et al. Periodic boundary value problems for impulsive conformable fractional integro-differential equations, Bound. Value Probl., 2016 (2016), 1-18. |
[34] | H. Batarfi, J. Losada, J. J. Nieto, et al. Three-point boundary value problems for conformable fractional differential equations, J. Funct. Space., 2015 (2015), 1-6. |
[35] | B. Bayour, D. F. M. Torres, Existence of solution to a local fractional nonlinear differential equation, J. Comput. Appl. Math., 312 (2017), 127-133. |
[36] | M. Bouaouid, M. Hannabou, K. Hilal, Nonlocal conformable fractional differential equations with a measure of noncompactness in Banach spaces, J. Math., 2020 (2020), 1-6. |
[37] | L. M. Feng, S. R. Sun, Oscillation theorems for three classes of conformable fractional differential equations, Adv. Differ. Equ., 2019 (2019), 1-30. |
[38] | A. Jaiswal, D. Bahuguna, Semilinear conformable fractional differential equations in Banach spaces, Differ. Equ. Dyn. Syst., 27 (2019), 313-325. |
[39] | M. M. Li, J. R. Wang, D. O'Regan, Existence and Ulam's stability for conformable fractional differential equations with constant coefficients, Bull. Malays. Math. Sci. Soc., 42 (2019), 1791-1812. |
[40] | S. M. Meng, Y. J. Cui, Multiplicity results to a conformable fractional differential equations involving integral boundary condition, Complexity, 2019 (2019), 1-8. |
[41] | E. Ünal, A. Gökdoğan, Solution of conformable fractional ordinary differential equations via differential transform method, Optik, 128 (2017), 264-273. |
[42] | J. F. Wang, C. Z. Bai, Antiperiodic boundary value problems for impulsive fractional functional differential equations via conformable derivative, J. Funct. Space., 2018 (2018), 1-11. |
[43] | W. Y. Zhong, L. F. Wang, Positive solutions of conformable fractional differential equations with integral boundary conditions, Bound. Value Probl., 2018 (2018), 1-12. |
[44] | W. Y. Zhong, L. F. Wang, Basic theory of initial value problems of conformable fractional differential equations, Adv. Differ. Equ., 2018 (2018), 1-14. |
[45] | K. Deimling, Nonlinear Functional Analysis, Springer-Verlag, New York, 1985. |
[46] | T. Cardinali, R. Precup, P. Rubbioni, A unified existence theory for evolution equations and systems under nonlocal conditions, J. Math. Anal. Appl., 432(2015), 1039-1057. |
[47] | H. P. Heinz, On the behaviour of measure of noncompactness with respect to differentiation and integration of vector-valued functions, Nonlinear Anal., 7 (1983), 1351-1371. |
[48] | D. O'Regan, R. Precup, Existence criteria for integral equations in Banach spaces, J. Inequal. Appl., 6 (2001), 77-97. |
1. | K. Rajendra Prasad, Mahammad Khuddush, D. Leela, Existence, uniqueness and Hyers–Ulam stability of a fractional order iterative two-point boundary value Problems, 2021, 1012-9405, 10.1007/s13370-021-00895-5 | |
2. | Alexandru Tudorache, Rodica Luca, Positive Solutions of a Fractional Boundary Value Problem with Sequential Derivatives, 2021, 13, 2073-8994, 1489, 10.3390/sym13081489 | |
3. | Ymnah Alruwaily, Bashir Ahmad, Sotiris K. Ntouyas, Ahmed S. M. Alzaidi, Existence Results for Coupled Nonlinear Sequential Fractional Differential Equations with Coupled Riemann–Stieltjes Integro-Multipoint Boundary Conditions, 2022, 6, 2504-3110, 123, 10.3390/fractalfract6020123 | |
4. | Ymnah Alruwaily, Lamya Almaghamsi, Kulandhaivel Karthikeyan, El-sayed El-hady, Existence and uniqueness for a coupled system of fractional equations involving Riemann-Liouville and Caputo derivatives with coupled Riemann-Stieltjes integro-multipoint boundary conditions, 2023, 8, 2473-6988, 10067, 10.3934/math.2023510 | |
5. | Bashir Ahmad, Ymnah Alruwaily, Ahmed Alsaedi, Sotiris K. Ntouyas, Sequential fractional differential equations with nonlocal integro-multipoint boundary conditions , 2021, Accepted, 1450-5444, 10.30755/NSJOM.12668 | |
6. | Ymnah Alruwaily, Shorog Aljoudi, Lamya Almaghamsi, Abdellatif Ben Makhlouf, Najla Alghamdi, Existence and Uniqueness Results for Different Orders Coupled System of Fractional Integro-Differential Equations with Anti-Periodic Nonlocal Integral Boundary Conditions, 2023, 15, 2073-8994, 182, 10.3390/sym15010182 | |
7. | Mohamed Bouaouid, Mild Solutions of a Class of Conformable Fractional Differential Equations with Nonlocal Conditions, 2023, 39, 0168-9673, 249, 10.1007/s10255-022-1015-6 | |
8. | Mohamed Houas, Mohamed Abbas, Francisco Martínez, Existence and Mittag-Leffler-Ulam-Stability results of sequential fractional hybrid pantograph equations, 2023, 37, 0354-5180, 6891, 10.2298/FIL2320891H | |
9. | Suzana Aleksic, Alberto Cabada, Sladjana Dimitrijevic, Tatjana Tomovic-Mladenovic, The existence of a solution for nonlinear fractional differential equations where nonlinear term depends on the fractional and first order derivative of an unknown function, 2023, 37, 0354-5180, 3871, 10.2298/FIL2312871A |