Identifying influential spreaders in complex networks is a crucial issue that can help control the propagation process in complex networks. An aviation network is a typical complex network, and accurately identifying the key city nodes in the aviation network can help us better prevent network attacks and control the spread of diseases. In this paper, a method for identifying key nodes in undirected weighted networks, called weighted Laplacian energy centrality, was proposed and applied to an aviation network constructed from real flight data. Based on the analysis of the topological structure of the network, the paper recognized critical cities in this network, then simulation experiments were conducted on key city nodes from the perspectives of network dynamics and robustness. The results indicated that, compared with other methods, weighted Laplacian energy centrality can identify the city nodes with the most spreading influence in the network. From the perspective of network robustness, the identified key nodes also have the characteristics of accurately and quickly destroying network robustness.
Citation: Shuying Zhao, Shaowei Sun. A study on centrality measures in weighted networks: A case of the aviation network[J]. AIMS Mathematics, 2024, 9(2): 3630-3645. doi: 10.3934/math.2024178
[1] | 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 |
[2] | Kirti Kaushik, Anoop Kumar, Aziz Khan, Thabet Abdeljawad . Existence of solutions by fixed point theorem of general delay fractional differential equation with p-Laplacian operator. AIMS Mathematics, 2023, 8(5): 10160-10176. doi: 10.3934/math.2023514 |
[3] | Kishor D. Kucche, Sagar T. Sutar, Kottakkaran Sooppy Nisar . Analysis of nonlinear implicit fractional differential equations with the Atangana-Baleanu derivative via measure of non-compactness. AIMS Mathematics, 2024, 9(10): 27058-27079. doi: 10.3934/math.20241316 |
[4] | Luchao Zhang, Xiping Liu, Zhensheng Yu, Mei Jia . The existence of positive solutions for high order fractional differential equations with sign changing nonlinearity and parameters. AIMS Mathematics, 2023, 8(11): 25990-26006. doi: 10.3934/math.20231324 |
[5] | Yige Zhao, Yibing Sun, Zhi Liu, Yilin Wang . Solvability for boundary value problems of nonlinear fractional differential equations with mixed perturbations of the second type. AIMS Mathematics, 2020, 5(1): 557-567. doi: 10.3934/math.2020037 |
[6] | Sunisa Theswan, Sotiris K. Ntouyas, Jessada Tariboon . Coupled systems of ψ-Hilfer generalized proportional fractional nonlocal mixed boundary value problems. AIMS Mathematics, 2023, 8(9): 22009-22036. doi: 10.3934/math.20231122 |
[7] | 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 |
[8] | 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 |
[9] | Iyad Suwan, Mohammed S. Abdo, Thabet Abdeljawad, Mohammed M. Matar, Abdellatif Boutiara, Mohammed A. Almalahi . Existence theorems for Ψ-fractional hybrid systems with periodic boundary conditions. AIMS Mathematics, 2022, 7(1): 171-186. doi: 10.3934/math.2022010 |
[10] | Naimi Abdellouahab, Keltum Bouhali, Loay Alkhalifa, Khaled Zennir . Existence and stability analysis of a problem of the Caputo fractional derivative with mixed conditions. AIMS Mathematics, 2025, 10(3): 6805-6826. doi: 10.3934/math.2025312 |
Identifying influential spreaders in complex networks is a crucial issue that can help control the propagation process in complex networks. An aviation network is a typical complex network, and accurately identifying the key city nodes in the aviation network can help us better prevent network attacks and control the spread of diseases. In this paper, a method for identifying key nodes in undirected weighted networks, called weighted Laplacian energy centrality, was proposed and applied to an aviation network constructed from real flight data. Based on the analysis of the topological structure of the network, the paper recognized critical cities in this network, then simulation experiments were conducted on key city nodes from the perspectives of network dynamics and robustness. The results indicated that, compared with other methods, weighted Laplacian energy centrality can identify the city nodes with the most spreading influence in the network. From the perspective of network robustness, the identified key nodes also have the characteristics of accurately and quickly destroying network robustness.
This article is devoted to the following nonlinear fractional differential equation with periodic boundary condition
{cDα0+x(t)−λx(t)=f(t,x(t)), 0<t≤ω,x(0)=x(ω), | (1.1) |
where λ≤0, 0<α≤1 and cDα0+ is Caputo fractional derivative
cDα0+x(t)=1Γ(1−α)∫t0(t−s)−αx′(s)ds. |
Differential equations of fractional order occur more frequently on different research areas and engineering, such as physics, economics, chemistry, control theory, etc. In recent years, boundary value problems for fractional differential equation have become a hot research topic, see [2,3,4,5,6,7,9,10,11,12,13,15,16,17,18,19,20,21,24,25]. In [27], Zhang studied the boundary value problem for nonlinear fractional differential equation
{cDα0+u(t)=f(t,u(t)), 0<t<1,u(0)+u′(0)=0,u(1)+u′(1)=0, | (1.2) |
where 1<α≤2, f:[0,1]×[0,+∞)→[0,+∞) is continuous and cDα0+ is Caputo fractional derivative
cDα0+u(t)=1Γ(2−α)∫t0(t−s)1−αu″(s)ds. |
The author obtained the existence of the positive solutions by using the properties of the Green function, Guo-Krasnosel'skill fixed point theorem and Leggett-Williams fixed point theorem.
Ahmad and Nieto [1] studied the anti-periodic boundary value problem of fractional differential equation
{cDq0+u(t)=f(t,u(t)), 0≤t≤T, 1<q≤2,u(0)=−u(T), cDp0+u(0)=−cDp0+u(T), 0<p<1, | (1.3) |
where f:[0,T]×R→R is continuous. The authors obtained some existence and uniqueness results by applying fixed point principles. The anti-periodic boundary value condition in this article corresponds to the anti-periodic condition u(0)=−u(T),u′(0)=−u′(T) in ordinary differential equation.
In [26], Zhang studied the following fractional differential equation
{Dδ0+u(t)=f(t,u), 0<t≤T,limt→0+t1−αu(t)=u0, | (1.4) |
where 0<δ<1, T>0, u0∈R and Dδ0+ is Riemann-Liouville fractional derivative
Dδ0+u(t)=1Γ(1−δ)ddt∫t0(t−s)−δu(s)ds. |
The author obtained the existence and uniqueness of the solutions by the method of upper and lower solutions and monotone iterative method.
In [7], Belmekki, Nieto and Rodriguez-Lopez studied the following equation
{Dδ0+u(t)−λu(t)=f(t,u(t)),0<t≤1,limt→0+t1−δu(t)=u(1), | (1.5) |
where 0<δ<1, λ∈R, f is continuous. The authors obtained the existence and uniqueness of the solutions by using the fixed point theorem. Cabada and Kisela [8] studied the following equation
{Dδ0+u(t)−λu(t)=f(t,t1−αu(t)),0<t≤1,limt→0+t1−δu(t)=u(1), | (1.6) |
where 0<δ<1, λ≠0(λ∈R), f is continuous. The authors studied the existence and uniqueness of periodic solutions by using Krasnosel'skii fixed point theorem and monotone iterative method. In [7,8], the boundary condition limt→0+t1−δu(t)=u(1) was called as periodic boundary value condition of Riemann-Liouville fractional differential equation, which is different from the periodic condition for ordinary differential equation. The boundary value condition u(0)=u(1) is not suitable for Riemann-Liouville fractional differential equation.
For the ordinary differential equation, the periodic boundary value problem is closely related to the periodic solution. For the Caputo fractional differential equation, the periodic boundary value condition u(0)=u(w) is meaningful. As far as we know, few work involves the periodic boundary value problem for Caputo fractional. The aim of this paper is to show the existence of positive solutions of (1.1) by using Krasnosel'skii fixed point theorem. Meanwhile, we also use the monotone iterative method to study the extremal solutions problem
{cDα0+u(t)=f(t,u(t)),0<t≤ω,u(0)=u(ω). | (1.7) |
The paper is organized as follows. In Section 2, we recall and derive some results on Mittag-Leffler functions. In Section 3, we use the Laplace transform to obtain the solution of a linear problem and discuss some properties of Green's function. In Section 4, the existence of positive solution is studied by using the Krasnosel'skii fixed point theorem. In Section 5, the existence of extremal solutions is proved by utilizing the monotone iterative technique. Section 6 is conclusion of the paper.
A key role in the theory of linear fractional differential equation is played by the well-known two-parameter Mittag-Leffler function
Eα,β(z)=Σ∞k=0zkΓ(αk+β), z∈R, α, β>0. | (2.1) |
We recall and derive some of their properties and relationships summarized in the following.
Proposition 2.1. Let α∈(0,1],β>0,λ∈R and ξ>0. Then it holds
(C1) limt→0+Eα,β(λtα)=1Γ(β),limt→0+Eα,1(λtα)=1.
(C2) Eα,α+1(λtα)=λ−1t−α(Eα,1(λtα)−1).
(C3) Eα,α(λtα)>0, Eα,1(λtα)>0 for all t≥0.
(C4) Eα,α(λtα) is decreasing in t for λ<0and increasing for λ>0 for all t>0.
(C5) Eα,1(λtα) is decreasing in t for λ<0 and increasing forλ>0 for all t>0.
(C6) ∫ξ0tβ−1Eα,β(λtα)dt=ξβEα,β+1(λξα).
Proof. (C1) It is obtained by an immediate calculation from (2.1).
(C2) By (2.1), we get
Eα,1(λtα)=Σ∞k=0(λtα)kΓ(αk+1)=1+λtαΓ(α+1)+(λtα)2Γ(2α+1)+(λtα)3Γ(3α+1)+⋯, |
Eα,α+1(λtα)=Σ∞k=0(λtα)kΓ(αk+α+1)=1Γ(α+1)+λtαΓ(2α+1)+(λtα)2Γ(3α+1)+⋯. |
Hence,
Eα,α+1(λtα)=λ−1t−α(Eα,1(λtα)−1). |
(C3) It follows from [23,Lemma 2.2].
(C4) It follows from [8,Proposition 1].
(C5) By a direct calculation, we get
ddtEα,1(λtα)=λtα−1Eα,α(λtα), |
since α∈(0,1], t>0 and Eα,α(λtα) is positive by Proposition (C3), the assertion is proved.
(C6) It follows from (1.99) of [20].
In this section, we deal with the linear case that f(t,x)=f(t) is a continuous function by mean of the Laplace transform for caputo fractional derivative
(LcDα0+x)(s)=sαX(s)−sα−1x(0), 0<α≤1, | (3.1) |
where L denotes the Laplace transform operator, X(s) denotes the Laplace transform of x(t).
From Lemma 3.2 of [14], we get
(LEα(λtα))(s)=sα−1sα−λ,Re(s)>0,λ∈C,|λs−α|<1, | (3.2) |
and
(Ltβ−1Eα,β(λtα))(s)=sα−βsα−λ,Re(s)>0,λ∈C,|λs−α|<1. | (3.3) |
We do Laplace transform to the equation
cDα0+x(t)−λx(t)=f(t), x(0)=x(ω). | (3.4) |
By (3.1), we obtain
sαX(s)−λX(s)=F(s)+x(0)⋅sα−1, |
X(s)=F(s)sα−λ+sα−1sα−λ⋅x(0), |
where F denotes the Laplace transform of f. By (3.2) and (3.3), we obtain that
x(t)=∫t0(t−s)α−1Eα,α(λ(t−s)α)f(s)ds+x(0)⋅Eα,1(λtα). | (3.5) |
Hence,
x(ω)=x(0)Eα,1(λwα)+∫ω0(ω−s)α−1Eα,α(λ(ω−s)α)f(s)ds=x(0), |
which implies that
x(0)=∫ω0(ω−s)α−1Eα,α(λ(ω−s)α)f(s)ds1−Eα,1(λωα) |
if Eα,1(λωα)≠1. Therefore, if Eα,1(λωα)≠1, the solution of the problem (3.4) is
x(t)=∫ω0(ω−s)α−1Eα,α(λ(ω−s)α)f(s)ds1−Eα,1(λωα)⋅Eα,1(λtα)+∫t0Eα,α(λ(t−s)α)(t−s)1−αf(s)ds=∫t0Eα,1(λtα)Eα,α(λ(ω−s)α)(1−Eα,1(λωα))(ω−s)1−αf(s)ds+∫t0Eα,α(λ(t−s)α)(t−s)1−αf(s)ds+∫ωtEα,1(λtα)Eα,α(λ(ω−s)α)(1−Eα,1(λωα))(ω−s)1−αf(s)ds. |
Theorem 3.1. Let Eα,1(λωα)≠1, the periodic boundary value problem (3.4) has a unique solutiongiven by
x(t)=∫ω0Gα,λ(t,s)f(s)ds, |
where
Gα,λ(t,s)={Eα,1(λtα)Eα,α(λ(ω−s)α)(1−Eα,1(λωα))(ω−s)1−α+Eα,α(λ(t−s)α)(t−s)1−α, 0≤s<t≤ω,Eα,1(λtα)Eα,α(λ(ω−s)α)(1−Eα,1(λωα))(ω−s)1−α, 0≤t≤s<ω. | (3.6) |
Remark 3.2. The unique solution x of (3.4) is continuous on [0,ω].
Lemma 3.3. Let 0<α≤1,λ≠0 and sign(η) denotes the signumfunction. Then
(F1) limt→0+Gα,λ(t,s)=Eα,α(λ(ω−s)α)(1−Eα,1(λωα))(ω−s)1−α for anyfixed s∈[0,ω),
(F2)lims→ω−Gα,λ(t,s)=sign(−λ)⋅∞ for any fixed t∈[0,ω],
(F3)limt→s+Gα,λ(t,s)=∞ for any fixeds∈[0,ω),
(F4)Gα,λ(t,s)>0 for λ<0 and for all t∈[0,ω] and s∈[0,ω),
(F5)Gα,λ(t,s) changes its sign for λ>0 for t∈[0,ω] and s∈[0,ω).
Proof. (F1) When 0≤t≤s<ω, by Proposition 2.1 (C1) we can get (F1).
(F2) When 0≤t≤s<ω, it follows by Proposition 2.1 (C5) that 1−Eα,1(λωα) is positive for λ<0 and negative for λ>0. The unboundedness is implied by continuity of Mittag-Leffler function, Proposition 2.1 (C1) and the relation limt→0+t−r=∞ for r>0.
(F3) When 0≤s<t≤ω, the first term of (3.6) is finite due to the continuity of the involved functions. And by a similar argument as in the previous point of this proof we have the second term tends to infinity.
(F4) It is obtained by the positivity of all involved functions (Proposition 2.1 (C3)) and the inequation 1−Eα,1(λωα)>0 for λ<0.
(F5) When 0≤s<t≤ω, the second term of (3.6) is positive due to
lims→t−(t−s)α−1Eα,α(λ(t−s)α)=+∞, |
and by the positivity of all involved functions (Proposition 2.1 (C3)) we get the proof. When 0≤t≤s<ω, it is obtained by 1−Eα,1(λωα)<0 for λ>0 and the positivity of all involved functions (Proposition 2.1 (C3)).
Proposition 3.4. Let α∈(0,1] and λ<0. Then the Green's function (3.6) satisfies
(K1) Gα,λ(t,s)≥m=:Eα,1(λωα)Eα,α(λωα)|λ|ωEα,α+1(λωα)>0,
(K2)∫ω0Gα,λ(t,s)ds=M=:1|λ| for all t∈[0,ω].
Proof. (K1) For 0≤t≤s<ω, we deduce from Proposition 2.1 (C4),(C5) that Gα,λ has the minimum on the line t=s. Hence,
Gα,λ(t,s)≥Gα,λ(t,t)=Eα,1(λtα)Eα,α(λ(ω−t)α)(1−Eα,1(λωα))(ω−t)1−α≥Eα,1(λωα)Eα,α(λωα)(1−Eα,1(λωα))ω1−α=Eα,1(λωα)Eα,α(λωα)[1−(Eα,α+1(λωα)λωα+1)]ω1−α=Eα,1(λωα)Eα,α(λωα)|λ|ωEα,α+1(λωα). |
For 0≤s<t≤ω, we have
Gα,λ(t,s)≥Eα,1(λωα)Eα,α(λωα)|λ|ωEα,α+1(λωα)+Eα,α(λωα)ω1−α≥Eα,1(λωα)Eα,α(λωα)|λ|ωEα,α+1(λωα). |
(K2) Employing Proposition 2.1, we get
∫ω0Gα,λ(t,s)ds=∫ω0Eα,1(λtα)Eα,α(λ(ω−s)α)(ω−s)1−α(1−Eα,1(λωα))ds+∫t0Eα,α(λ(t−s)α)(t−s)1−αds=Eα,1(λtα)1−Eα,1(λωα)∫ω0(ω−s)α−1Eα,α(λ(ω−s)α)ds+∫t0(t−s)α−1Eα,α(λ(t−s)α)ds=Eα,1(λtα)1−Eα,1(λωα)⋅ωαEα,α+1(λωα)+tαEα,α+1(λtα)=Eα,1(λtα)1−Eα,1(λωα)⋅ωα⋅λ−1ω−α(Eα,1(λωα)−1)+tα⋅λ−1t−α(Eα,1(λtα)−1)=1|λ|, |
which completes the proof.
Let C[0,ω] be the space continuous function on [0,ω] with the norm ‖x‖=sup{|x(t)|:t∈[0,ω]}.
In this section, we always assume that λ<0. Clearly, x is a solution of (1.1) if and only if
x(t)=∫ω0Gα,λ(t,s)f(s,x(s))ds, | (4.1) |
where Gα,λ is Green's function defined in Theorem 3.1.
The following famous Krasnosel'skii fixed point theorem, which is main tool of this section.
Theorem 4.1. [22] Let B be a Banach space, and let P⊂B be a cone. AssumeΩ1,Ω2 two open and bounded subsets of B with0∈Ω1,Ω1⊂Ω2 and letA:P∩(¯Ω2∖Ω1)→P bea completely continuous operator such that one of the followingconditions is satisfied:
(L1) ‖Ax‖≤‖x‖, if x∈P∩∂Ω1, and‖Ax‖≥‖x‖, if x∈P∩∂Ω2,
(L2) ‖Ax‖≥‖x‖, if x∈P∩∂Ω1, and‖Ax‖≤‖x‖, if x∈P∩∂Ω2.
Then, A has at least one fixed point inP∩(¯Ω2∖Ω1).
Proposition 4.2. Assume that there exist 0<r<R,0<c1<c2 such that
f:[0,ω]×[mc1ωMc2r,R]→R is continuous, | (4.2) |
c1≤f(t,u)≤c2, ∀(t,u)∈[0,ω]×[mc1ωMc2r,R]. | (4.3) |
Let P⊂C[0,ω] be the cone
P={x∈C[0,ω]:mint∈[0,ω]x(t)≥mc1ωMc2‖x‖}. |
Then the operator A:¯PR∖Pr→P givenby
Ax(t)=∫ω0Gα,λ(t,s)f(s,x(s))ds | (4.4) |
is completely continuous, where Pl={u∈P:‖u‖<l}.
Proof. Let x∈¯PR∖Pr, then
mc1ωMc2r≤x(t)≤R for all t∈[0,ω]. | (4.5) |
We first show that A is well-defined, i.e. that A:¯PR∖Pr→P. Note that
Ax(t)=∫ω0Gα,λ(t,s)f(s,x(s))ds=∫ω0Eα,1(λtα)Eα,α(λ(ω−s)α)(ω−s)1−α(1−Eα,1(λωα))f(s,x(s))ds +∫t0Eα,α(λ(t−s)α)(t−s)1−αf(s,x(s))ds=kq(ω)Eα,1(λtα)+q(t), | (4.6) |
where k=11−Eα,1(λωα) and
q(t)={∫t0Eα,α(λ(t−s)α)(t−s)1−αf(s,x(s))ds,0<t≤ω,0,t=0. | (4.7) |
Clearly, for t∈(0,ω]
0<q(t)≤c2∫t0Eα,α(λ(t−s)α)(t−s)1−αds≤c2Γ(α)∫t0(t−s)α−1ds=c2Γ(α+1)tα, | (4.8) |
which implies that q is continuous at t=0. On the other hand,
q(t)=∑k<1α−1λk∫t0(t−s)αk+α−1Γ(αk+α)f(s,x(s))ds+∫t0∑k≥1α−1λkΓ(αk+α)(t−s)αk+α−1f(s,x(s))ds=∑k<1α−1λkΓ(αk+α)∫t0uαk+α−1f(t−u,x(t−u))du+∫t0∑k≥1α−1λkΓ(αk+α)(t−s)αk+α−1f(s,x(s))ds=:H1(t)+H2(t), |
where
H1(t)=∑k<1α−1λkΓ(αk+α)∫t0uαk+α−1f(t−u,x(t−u))du,H2(t)=∫t0∑k≥1α−1λkΓ(αk+α)(t−s)αk+α−1f(s,x(s))ds. |
Since
|uαk+α−1f(t−u,x(t−u))|≤c2uαk+α−1,u>0, t∈(0,ω], x∈¯PR/Pr,|λk(t−s)αk+α−1Γ(αk+α)f(s,x(s))|≤λkΓ(αk+α)tαk+α−1c2, 0≤s≤t, x∈¯PR/Pr,∫t0uαk+α−1du<+∞, t∈(0,ω],∑k≥1α−1λktαk+α−1Γ(αk+α)<+∞, t∈(0,ω], |
we obtain that H1∈C[0,ω] and
H2(t)=∑k≥1α−1λkΓ(αk+α)∫t0(t−s)αk+α−1f(s,x(s))ds=:∑k≥1α−1λkΓ(αk+α)uk(t). |
Noting that uk∈C(0,ω]
|uk(t)|≤c2ωαk+ααk+α, t∈(0,ω], ∑λkΓ(αk+α)ωαk+ααk+α<+∞, |
we have H2∈C(0,ω]. Hence, q∈C[0,ω].
Moreover,
mc1ωMc2‖Ax‖=mc1ωMc2sup∫ω0Gα,λ(t,s)f(s,x(s))ds≤mc1ωMsup∫ω0Gα,λ(t,s)ds=mc1ω≤min∫ω0Gα,λ(t,s)f(s,x(s))ds=mint∈[0,ω]Ax(t), | (4.9) |
which means that A:¯PR∖Pr→P.
Next, we show that A is continuous on ¯PR∖Pr. Let xn,x∈¯PR∖Pr and ‖xn−x‖→0. From (4.2), we have ‖f(t,xn(t))−f(t,x(t))‖→0,
‖Axn−Ax‖=supt∈[0,ω]|∫ω0Gα,λ(t,s)(f(s,x(s))−f(s,y(s)))ds|≤supt∈[0,ω]∫ω0Gα,λ(t,s)ds‖f(t,xn(t))−f(t,x(t))‖≤M‖f(t,xn(t))−f(t,x(t))‖→0, |
which implies that A is continuous. From (4.6), we get that Ax(t) is uniformly bounded. Finally, we show that {Ax|x∈¯PR/Pr} is an equicontinuity in C[0,ω]. By (4.6), we have
|Ax(t1)−Ax(t2)|=|∫ω0(Gα,λ(t1,s)−Gα,λ(t2,s))f(s,x(s))ds|≤kq(ω)|Eα,1(λtα1)−Eα,1(λtα2)|+|q(t1)−q(t2)|. |
Since Eα,1(λtα)∈C[0,ω], q(t)∈C[0,ω] are uniformly continuous, |Eα,1(λtα1)−Eα,1(λtα2)| and |q(t1)−q(t2)| tend to zero as |t1−t2|→0. Hence, {Ax(t)|x∈¯PR∖Pr} is equicontinuous in C[0,ω].
Finally, by Arzela-Ascoli theorem, we can obtain that A is compact. Hence, it is completely continuous.
Theorem 4.3. Assume that there exist 0<r<R,0<c1<c2 such that (4.2)and (4.3) hold. Further suppose one of the followingconditions is satisfied
(i)f(t,u)≥Mc2m2ω2c1u, ∀(t,u)∈[0,ω]×[mc1ωMc2r,r],
f(t,u)≤|λ|u, ∀(t,u)∈[0,ω]×[mc1ωMc2R,R],
(ii)f(t,u)≤|λ|u, ∀(t,u)∈[0,ω]×[mc1ωMc2r,r],
f(t,u)≥Mc2m2ω2c1u, ∀(t,u)∈[0,ω]×[mc1ωMc2R,R].
Then (1.1) has at least a positive solution x with r≤‖x‖≤R.
Proof. Here we only consider the case (i). By Proposition 4.2, A:¯PR∖Pr→P is completely continuous. For x∈∂Pr, we have
‖x‖=r, mc1ωMc2r≤x(t)≤r, ∀t∈[0,ω] |
and
Ax(t)≥m∫ω0f(s,x(s))ds≥Mc2mω2c1∫ω0x(s)ds≥r=‖x‖. |
Similarly, if x∈∂PR,
mc1ωMc2R≤x(t)≤R, t∈[0,ω], |
0≤Ax(t)≤∫ω0Gα,λ(t,s)|λ|x(s)ds≤|λ|R∫ω0Gα,λ(t,s)ds=R=‖x‖. |
By Theorem 4.1, there exists x∈¯PR∖Pr such that Ax=x and x is a solution of (1.1). Moreover,
mc1ωMc2r≤x(t)≤R. |
Corollary 4.4. Let c1<c2 be positive reals and f(t,x) satisfy the conditions
(i) c1≤f(t,x)≤c2 for all x≥0,
(ii) f:[0,w]×(0,+∞)→R is a continuous function.
Then problem (1.1) has a positive solution.
Proof. Let 0<r<c21m2ω2Mc2, R>c22M|λ|m1c1ω, then (4.2) and (4.3) are satisfied. Clearly, for (t,u)∈[0,ω]×[mc1ωMc2r,r],
f(t,u)≥c1≥Mc2m2ω2c1r≥Mc2m2ω2c1u |
and for (t,u)∈[0,ω]×[mc1ωMc2R,R],
f(t,u)≤c2≤|λ|mc1ωMc2R≤|λ|u. |
Hence, by Theorem 4.3 (1.1) has at least a positive solution.
Example 4.5. Consider the equation
{cDα0+x(t)−λx(t)=1+x1β(t), 0<x≤ω,x(0)=x(ω), | (4.10) |
where 0<α≤1,β>1 and
Λ={λ<0,|λ|Eα,1(λωα)Eα,α(λωα)≥4Eα,α+1(λωα)}≠∅. |
Choosing c1=1,c2=2,r=110min{1,m2ω22M},R=1. It is easy to check that (4.2) and (4.3) hold. For λ∈Λ,
f(t,u)=1+u1β≥1≥Mc2m2ωc1r≥Mc2m2ωc1u,∀(t,u)∈[0,ω]×[mc1ωMc2r,r], |
f(t,u)≤2≤|λ|⋅mc1ωMc2=|λ|⋅mω2M≤|λ|u,∀(t,u)∈[0,ω]×[mc1ωMc2R,R]. |
Hence, (4.10) has at least one positive solution for λ∈Λ.
In this section, by using the monotone iterative method, we discuss the existence of solutions when λ=0 in (1.1). Firstly, we give the definition of the upper and lower solutions and get monotone iterative sequences with the help of the corresponding linear equation. Finally, we prove the limits of the monotone iterative sequences are solutions of (1.7).
Definition 5.1. Let h,k∈C1[0,ω]. h and k are called lower solution and upper solution of problem (1.7), respectively if h and k satisfy
cDα0+h(t)≤f(t,h(t)), 0<t≤ω, h(0)≤h(ω), | (5.1) |
cDα0+k(t)≥f(t,k(t)), 0<t≤ω, k(0)≥k(ω), | (5.2) |
Clearly, if g the lower solution or upper solution of (1.7), then cDα0+g is continuous on [0,ω].
Lemma 5.2. Let δ∈C[0,ω] with δ≥0 and p∈R with p≤0. Then
{cDα0+z(t)−λz(t)=δ(t), 0<t≤ω,z(ω)−z(0)=p, | (5.3) |
has a unique solution z(t)≥0 for t∈[0,ω], where 0<α≤1,λ<0.
Proof. Let z1, z2 are two solutions of (5.3) and v=z1−z2, then
{cDα0+v(t)−λv(t)=0, 0<t≤ω,v(ω)=v(0). | (5.4) |
Using Theorem 3.1, (5.4) has trivial solution v=0.
By (3.5), we can verify that problem (5.3) has a unique solution
z=∫ω0Gα,λ(t,s)δ(t)ds+p⋅Eα,1(λωα)Eα,1(λωα)−1. |
As consequence, by Proposition 2.1 (C3), (C5) and Lemma 3.3 (F4), we conclude that z(t)≥0. This completes the proof.
Theorem 5.3. Assume that h,k are the lower and upper solutions ofproblem (1.7) and h≤k. Moreover, supposethat f satisfies the following properties:
(M) there is λ<0 such that for all fixed t∈[0,ω], f(t,x)−λx is nondecreasing in h(t)≤x≤k(t),
(J) f:[0,ω]×[h(t),k(t)]→R is a continuous function.
Then there are two monotone sequences {hn} and {kn} are nonincreasingand nondecreasing, respectively with h0=h and k0=ksuch that limn→∞hn=¯h(t), limn→∞kn=¯k(t)uniformly on [0,ω], and ¯h,¯k are the minimal and the maximalsolutions of (1.7) respectively, such that
h0≤h1≤h2≤...≤hn≤¯h≤x≤¯k≤kn≤...≤k2≤k1≤k0 |
on [0,ω], where x is any solution of (1.7) such that h(t)≤x(t)≤k(t)on [0,ω].
Proof. Let [h,k]={u∈C[0,ω]:h(t)≤u(t)≤k(t),t∈[0,ω]}. For any η∈[h,k], we consider the equation
{cDα0+x(t)−λx(t)=f(t,η(t))−λη(t), 0<t≤ω,x(0)=x(ω), |
Theorem 3.1 implies the above problem has a unique solution
x(t)=∫ω0Gα,λ(t,s)(f(s,η(s))−λη(s))ds. | (5.5) |
Define an operator B by x=Bη, we shall show that
(a) k≥Bk,Bh≥h,
(b) B is nondecreasing on [h,k].
To prove (a). Denote θ=k−Bk, we have
cDα0+θ(t)−λθ(t)= cDα0+k(t)− cDα0+Bk(t)−λ(k(t)−Bk(t))≥f(t,k(t))−((f(t,k(t))−λk(t))−λk(t)=0, |
and θ(w)−θ(0)≤0. Since k∈C1[0,ω],
cDα0+k∈C[0,ω], cDα0+Bk∈C[0,ω]. |
By Lemma 5.2, θ≥0, i.e. k≤Bk. In an analogous way, we can show that Bh≥h.
To prove (b). We show that Bη1≤Bη2 if h≤η1≤η2≤k. Let z1=Bη1, z2=Bη2 and z=z2−z1, then by (M), we have
cDα0+z(t)−λz(t)= cDα0+z2(t)− cDα0+z1(t)−λ(z2(t)−z1(t))=f(t,η2(t))−λη2(t)−(f(t,η1(t))−λη1(t))≥0, |
and v(ω)=v(0). By Lemma 5.2, z(t)≥0, which implies Bη1≤Bη2.
Define the sequence {hn}, {kn} with h0=h,k0=k such that hn+1=Bhn,kn+1=Bkn for n=0,1,2,.... From (a) and (b), we have
h0≤h1≤h2≤...≤hn≤kn≤...≤k2≤k1≤k0 |
on t∈[0,ω], and
hn(t)=∫ω0Gα,λ(t,s)(f(s,hn−1(s))−λhn−1(s))ds, |
kn(t)=∫ω0Gα,λ(t,s)(f(s,kn−1(s))−λkn−1(s))ds. |
Therefore, there exist ¯h,¯k such that limn→∞hn=¯h, limn→∞kn=¯k.
Similar to the proof of Proposition 4.2, we can show that B:[h,k]→[h,k] is a completely continuous operator. Therefore, ¯h,¯k are solutions of (1.7).
Finally, we prove that if x∈[h0,k0] is one solution of (1.7), then ¯h(t)≤x(t)≤¯k(t) on [0,ω]. To this end, we assume, without loss of generality, that hn(t)≤x(t)≤kn(t) for some n. From property (b), we can get that hn+1(t)≤x(t)≤kn+1(t),t∈[0,ω]. Since h0(t)≤x(t)≤k0(t), we can conclude that
hn(t)≤x(t)≤kn(t), forall n. |
Passing the limit as n→∞, we obtain ¯h(t)≤x(t)≤¯k(t), t∈[0,ω]. This completes the proof.
Example 5.4. Consider the equation
{cDα0+x(t)=t+1−x2(t), 0<x≤1,x(0)=x(1). | (5.6) |
It easy to check that h=1,k=2 are the low solution and upper solution of (5.6), respectively. Let λ=−10. For all t∈[0,ω],
f(t,u)−λu=t2+1−u2+10u |
is nondecreasing on u∈[1,2] and
f(t,u)=t+1−u2 |
is continuous on [0,ω]×[1,2].
Hence, there exist two monotone sequences {hn} and {kn}, nonincreasing and nondecreasing respectively, that converge uniformly to the extremal solutions of (5.6) on [h,k].
This paper focuses on the existence of solutions for the Caputo fractional differential equation with periodic boundary value condition. We use Green's function to transform the problem into the existence of the fixed points of some operator, and we prove the existence of positive solutions by using the Krasnosel'skii fixed point theorem. On the other hand, the existence of the extremal solutions for the special case of the problem is obtained from monotone iterative technique and lower and upper solutions method. Since the fractional differential equation is nonlocal equation, the process of verifying the compactness of operator is very tedious, and we will search for some better conditions to prove the compactness of the operator A in the follow-up research. Meanwhile, since the existence result for 0<α≤1 is obtained in present paper, we will discuss the existence of solutions for the Caputo fractional differential equation when n−1<α≤n in follow-up research.
The authors declare that there are no conflicts of interest regarding the publication of this paper.
[1] |
M. Ouyang, Z. Pan, L. Hong, L. Zhao, Correlation analysis of different vulnerability metrics on power grids, Physica A, 396 (2014), 204–211. https://doi.org/10.1016/j.physa.2013.10.041 doi: 10.1016/j.physa.2013.10.041
![]() |
[2] |
B. S. Kerner, Criticism of generally accepted fundamentals and methodologies of traffic and transportation theory: A brief review, Physica A, 391 (2013), 5261–5282. https://doi.org/10.1016/j.physa.2013.06.004 doi: 10.1016/j.physa.2013.06.004
![]() |
[3] |
D. Chen, H. Gao, L. Lü, T. Zhou, Identifying influential nodes in large-scale directed networks: The role of clustering, PLoS One, 8 (2013), e77455. https://doi.org/10.1371/journal.pone.0077455 doi: 10.1371/journal.pone.0077455
![]() |
[4] |
D. Chen, L. Lü, M. Shang, Y. Zhang, T. Zhou, Identifying influential nodes in complex networks, Physica A, 391 (2012), 1777–1887. https://doi.org/10.1016/j.physa.2011.09.017 doi: 10.1016/j.physa.2011.09.017
![]() |
[5] |
B. Michele, C. Davide, V. Simone, Efficiency of attack strategies on complex model and real-world networks, Physica A, 414 (2014), 174–180. https://doi.org/10.1016/j.physa.2014.06.079 doi: 10.1016/j.physa.2014.06.079
![]() |
[6] |
S. P. Borgatti, Identifying sets of key players in a social network, Comput. Math. Organ. Theory, 12 (2006), 21–34. https://doi.org/10.1007/s10588-006-7084-x doi: 10.1007/s10588-006-7084-x
![]() |
[7] |
T. Wen, D. Pelusi, Y. Deng, Vital spreaders identification in complex networks with multi-local dimension, Knowl. Based Syst., 195 (2020), 105717. https://doi.org/10.1016/j.knosys.2020.105717 doi: 10.1016/j.knosys.2020.105717
![]() |
[8] |
J. Zhao, Y. Song, F. Liu, Y. Deng, The identification of influential nodes based on structure similarity, Connect. Sci., 33 (2021), 201–218. https://doi.org/10.1080/09540091.2020.1806203 doi: 10.1080/09540091.2020.1806203
![]() |
[9] |
L. Freeman, Centrality in social networks conceptual clarification, Soc. Networks, 1 (1978), 215–239. https://doi.org/10.1016/0378-8733(78)90021-7 doi: 10.1016/0378-8733(78)90021-7
![]() |
[10] |
A. Barrat, M. Barthelemyt, R. Pastor-Satorrast, A. Vespignani, The architecture of complex weighted networks, PNAS, 101 (2004), 3747–3752. https://doi.org/10.1073/pnas.0400087101 doi: 10.1073/pnas.0400087101
![]() |
[11] | G. Sabidussi, The centrality index of a graph, Psychometrika, 31 (1966), 581–603. |
[12] |
L. Freeman, A set of measures of centrality based upon betweenness, Sociometry, 40 (1977), 35–41. https://doi.org/10.2307/3033543 doi: 10.2307/3033543
![]() |
[13] |
T. Opsahl, F. Agneessens, J. Skvoretzc, Node centrality in weighted networks: Generalizing degree and shortest paths, Soc. Networks, 32 (2010), 245–251. https://doi.org/10.1016/j.socnet.2010.03.006 doi: 10.1016/j.socnet.2010.03.006
![]() |
[14] |
Y. Ma, Z. Cao, X. Qi, Quasi-Laplacian centrality: A new vertex centrality measurement based on Quasi-Laplacian energy of networks, Physica A, 527 (2019), 121130. https://doi.org/10.1016/j.physa.2019.121130 doi: 10.1016/j.physa.2019.121130
![]() |
[15] |
S. Zhao, S. Sun, Identification of node centrality based on Laplacian energy of networks, Physica A, 609 (2023), 128353. https://doi.org/10.1016/j.physa.2022.128353 doi: 10.1016/j.physa.2022.128353
![]() |
[16] |
S. Bansal, J. Sen, Network assessment of Tier-Ⅱ Indian cities airports in terms of type, accessibility, and connectivity, Transp. Policy, 124 (2022), 221–232. https://doi.org/10.1016/j.tranpol.2021.05.009 doi: 10.1016/j.tranpol.2021.05.009
![]() |
[17] |
R. W. Daniel, R. Soumen, M. D. S. Raissa, Resilience and rewiring of the passenger airline networks in the United States, Phys. Rev. E, 82 (2010), 056101. https://doi.org/10.1103/PhysRevE.82.056101 doi: 10.1103/PhysRevE.82.056101
![]() |
[18] | A. Reggiani, S. Signoretti, P. Nijkamp, A. Cento, Network measures in civil air transport: A case study of lufthansa, 1 Eds., Berlin: Springer Press, 2009. https://doi.org/10.1007/978-3-540-68409-1-14 |
[19] |
J. Shen, H. Zong, Identification of critical transportation cities in the multimodal transportation network of China, Physica A, 628 (2023), 129174. https://doi.org/10.1016/j.physa.2023.129174 doi: 10.1016/j.physa.2023.129174
![]() |
[20] | F. Gao, Y. Dang, Analysis on distribution property of an international air transport network, Sci. Sci. Manage. S. T., 7 (2009), 75–79. |
[21] |
R. Guimera, L. Amaral, Modeling the world-wide airport network, Eur. Phys. J. B, 38 (2004), 381–385. https://doi.org/10.1140/epjb/e2004-00131-0 doi: 10.1140/epjb/e2004-00131-0
![]() |
[22] |
J. Hu, Y. Wang, X. He, Analysis and application of global aviation network based on complex network, Comput. Sci., 48 (2021), 321–325. https://doi.org/10.11896/jsjkx.200900112 doi: 10.11896/jsjkx.200900112
![]() |
[23] | W. Liu, M. Han, Z. Xie, Connectivity characteristics and community identification of worldcity network based on global airline, Econ. Geogr., 40 (2020), 34–40. |
[24] |
O. Lordan, J. Sallan, P Simo, Robustness of the air transport network, Transp. Res. Part E, 68 (2014), 155–163. https://doi.org/10.1016/j.tre.2014.05.011 doi: 10.1016/j.tre.2014.05.011
![]() |
[25] |
H. Mo, F. Jin, Y. Liu, J. Wang, Network analysis on centrality of airport system, Sci. Geol. Sin., 30 (2010), 204–212. https://doi.org/10.13249/j.cnki.sgs.2010.02.204 doi: 10.13249/j.cnki.sgs.2010.02.204
![]() |
[26] |
J. Li, X. Wen, M. Wu, F. Liu, S. Li, Identification of key nodes and vital edges in aviation network based on minimum connected dominating set, Physica A, 87 (2020), 123340. https://doi.org/10.1016/j.physa.2019.123340 doi: 10.1016/j.physa.2019.123340
![]() |
[27] | X. Luo, J. Wen, J. Zhong, Structural characteristics and robustness analysis of state-owned airline networks, Aeronaut. Comput. Tech., 51 (2021), 55–59. |
[28] |
D. J. Watts, S. H. Strogatz, Collective dynamics of 'small-world' networks, Nature, 393 (1998), 440–442. https://doi.org/10.1038/30918 doi: 10.1038/30918
![]() |
[29] | X. Feng, H. Jia, Aviation network robustness considering node failure and edge failure, J. Beijing Jiaotong Univ., 45 (2021), 84–92. |
[30] |
M. E. J. Newman, Y. Liu, Scientific collaboration networks. Ⅱ. Shortest paths, weighted networks, and centrality, Phys. Rev. E, 64 (2001), 016132. https://doi.org/10.1103/PhysRevE.64.016132 doi: 10.1103/PhysRevE.64.016132
![]() |
[31] |
P. Bonacich, Power and centrality: A family of measures, Am. J. Sociol., 92 (1987), 1170–1182. https://doi.org/10.1086/228631 doi: 10.1086/228631
![]() |
[32] |
K. Dietz, Infectious diseases of humans: Dynamics and control, Ann. Inter. Med., 117 (1992), 179. https://doi.org/10.1016/0169-4758(92)90018-W doi: 10.1016/0169-4758(92)90018-W
![]() |
[33] |
V. Latora, M. Marchiori, Efficient behavior of small-world networks, Phys. Rev. Lett., 87 (2001), 198701. https://doi.org/10.1103/PhysRevLett.87.198701 doi: 10.1103/PhysRevLett.87.198701
![]() |
[34] |
C. Castellano, R. Pastor-Satorras, Thresholds for epidemic spreading in networks, Phys. Rev. Lett., 105 (2010), 218701. https://doi.org/10.1103/PhysRevLett.105.218701 doi: 10.1103/PhysRevLett.105.218701
![]() |
1. | Zhilin Li, Guoping Chen, Weiwei Long, Xinyuan Pan, Variational approach to p-Laplacian fractional differential equations with instantaneous and non-instantaneous impulses, 2022, 7, 2473-6988, 16986, 10.3934/math.2022933 | |
2. | Saleh Fahad Aljurbua, Extended existence results for FDEs with nonlocal conditions, 2024, 9, 2473-6988, 9049, 10.3934/math.2024440 | |
3. | Saleh Fahad Aljurbua, Extended existence results of solutions for FDEs of order 1<γ≤2, 2024, 9, 2473-6988, 13077, 10.3934/math.2024638 | |
4. | Limin Guo, Weihua Wang, Cheng Li, Jingbo Zhao, Dandan Min, Existence results for a class of nonlinear singular p-Laplacian Hadamard fractional differential equations, 2024, 32, 2688-1594, 928, 10.3934/era.2024045 |