Incidence functions | f(I) | Source |
Bi-linear | βI | [27] |
Saturated | βI1+a1I | [28,29,30] |
Beddington-DeAngelis | βI1+a1S+a2I | [31,32] |
Specific nonlinear | βI1+a1S+a2I+a3SI | [33] |
We study the collective behaviors of the Wigner-Lohe (WL) model for quantum synchronization in phase space which corresponds to the phase description of the Schrödinger-Lohe (SL) model for quantum synchronization, and it can be formally derived from the SL model via the generalized Wigner transform. For this proposed model, we show that the WL model exhibits asymptotic aggregation estimates so that all the elements in the generalized Wigner distribution matrix tend to a common one. On the other hand, for the global unique solvability, we employ the fixed point argument together with the classical semigroup theory to derive the global unique solvability of mild and classical solutions depending on the regularity of initial data.
Citation: Seung-Yeal Ha, Gyuyoung Hwang, Dohyun Kim. On the complete aggregation of the Wigner-Lohe model for identical potentials[J]. Networks and Heterogeneous Media, 2022, 17(5): 665-686. doi: 10.3934/nhm.2022022
[1] | Ahmad A. Abubaker, Khaled Matarneh, Mohammad Faisal Khan, Suha B. Al-Shaikh, Mustafa Kamal . Study of quantum calculus for a new subclass of q-starlike bi-univalent functions connected with vertical strip domain. AIMS Mathematics, 2024, 9(5): 11789-11804. doi: 10.3934/math.2024577 |
[2] | Ekram E. Ali, Georgia Irina Oros, Rabha M. El-Ashwah, Abeer M. Albalahi . Applications of fuzzy differential subordination theory on analytic p -valent functions connected with q-calculus operator. AIMS Mathematics, 2024, 9(8): 21239-21254. doi: 10.3934/math.20241031 |
[3] | Ebrahim Amini, Mojtaba Fardi, Shrideh Al-Omari, Rania Saadeh . Certain differential subordination results for univalent functions associated with q-Salagean operators. AIMS Mathematics, 2023, 8(7): 15892-15906. doi: 10.3934/math.2023811 |
[4] | Shujaat Ali Shah, Ekram Elsayed Ali, Adriana Cătaș, Abeer M. Albalahi . On fuzzy differential subordination associated with q-difference operator. AIMS Mathematics, 2023, 8(3): 6642-6650. doi: 10.3934/math.2023336 |
[5] | Muhammad Sabil Ur Rehman, Qazi Zahoor Ahmad, H. M. Srivastava, Nazar Khan, Maslina Darus, Bilal Khan . Applications of higher-order q-derivatives to the subclass of q-starlike functions associated with the Janowski functions. AIMS Mathematics, 2021, 6(2): 1110-1125. doi: 10.3934/math.2021067 |
[6] | Syed Ghoos Ali Shah, Shahbaz Khan, Saqib Hussain, Maslina Darus . q-Noor integral operator associated with starlike functions and q-conic domains. AIMS Mathematics, 2022, 7(6): 10842-10859. doi: 10.3934/math.2022606 |
[7] | Aisha M. Alqahtani, Rashid Murtaza, Saba Akmal, Adnan, Ilyas Khan . Generalized q-convex functions characterized by q-calculus. AIMS Mathematics, 2023, 8(4): 9385-9399. doi: 10.3934/math.2023472 |
[8] | Georgia Irina Oros, Gheorghe Oros, Daniela Andrada Bardac-Vlada . Certain geometric properties of the fractional integral of the Bessel function of the first kind. AIMS Mathematics, 2024, 9(3): 7095-7110. doi: 10.3934/math.2024346 |
[9] | Ekram E. Ali, Miguel Vivas-Cortez, Rabha M. El-Ashwah . New results about fuzzy γ-convex functions connected with the q-analogue multiplier-Noor integral operator. AIMS Mathematics, 2024, 9(3): 5451-5465. doi: 10.3934/math.2024263 |
[10] | Ekram E. Ali, Nicoleta Breaz, Rabha M. El-Ashwah . Subordinations and superordinations studies using q-difference operator. AIMS Mathematics, 2024, 9(7): 18143-18162. doi: 10.3934/math.2024886 |
We study the collective behaviors of the Wigner-Lohe (WL) model for quantum synchronization in phase space which corresponds to the phase description of the Schrödinger-Lohe (SL) model for quantum synchronization, and it can be formally derived from the SL model via the generalized Wigner transform. For this proposed model, we show that the WL model exhibits asymptotic aggregation estimates so that all the elements in the generalized Wigner distribution matrix tend to a common one. On the other hand, for the global unique solvability, we employ the fixed point argument together with the classical semigroup theory to derive the global unique solvability of mild and classical solutions depending on the regularity of initial data.
Infectious diseases are massive threat for humans and also economy of the state. Proper understanding of a disease dynamics plays an important role in curtailing of the infection in a community. Implementation of a suitable strategy against the disease transmission is another challenge. Mathematical modeling approach is one of the key tools in order to handle these challenges. A number of disease models have been developed in existing literature which enables us to explore and control the spread of infectious diseases in a better way [1,2,3].
Most of these models are based on integer-order differential equations (IDEs). However, in the last few years, it is noticed that fractional-order differential equations (FDEs) can be used to model universal phenomena with greater degree of accuracy and its applications can be found in virous fields for example engineering, economics, control theory, finance and in epidemiology [4,5]. Fractional calculus is the generalization of classical integer-order calculus. The increasing interest of using FDEs in modeling of real world complex problems is due to its various properties which are not found in IDEs. In contrast of IDEs which local in nature, the FDEs are non-local and posses the memory effects which make it more superior because in many situations the future state of the model depends not only upon the current state but also on the previous history. These features enables FDEs to model the phenomena having not only the non-Gaussian but also for non-Markovian behavior. Further, the classical IDEs are unable to provide the information in between two different integer values. Various type of fractional-order operators were introduced in existing literature to over come such limitations of integer-order derivative [6,7,8]. The applications of these fractional operators can be found in various fields [9,10,11]. In epidemiology, the previous experiences and history of the past states of an epidemic have a great impact to explored its dynamics in a better way. Additionally, due to the non-Markovian behavior of any real epidemic process, models with FDEs are more realistic and therefore, play an essential role in the understanding the spread and possible control of diseases. Some recent and interesting work regarding fractional analysis to the real life problem have been studied in [12,13,14,15]. In these work, the authors studied mathematical model with delay in fractional derivatives and presented the results for stability and bifurcation.
In Mathematical models, the transmission of a disease i.e., the progression at which susceptible become infected is known as incidence rate function. This function plays a key role in epidemiological models and ensure that the model provides an acceptable qualitative analysis of the disease under consideration. A verity of incidence rate function have been implemented in the literature [16,17,18]. Most commonly, in disease models the bilinear incidence rate of the type bSI has been adopted, where b indicates the per unit contact rate. In 1927 Kermack and McKendrick [19] introduced a simple epidemic model using bilinear incidence rate and the assumptions that there are no births or deaths from the infection and each person of the population is initially susceptible to the disease. In [20] the authors introduced the saturated incidence after the cholera epidemic in Bari. Latter on, Korobeinikov and Maini[20] studied the global analysis of the SIR and SIRS models with nonlinear incidence rate. Recently, a fractional order SIR epidemic model with nonlinear incidence rate has been developed in [21].
Keeping the above discussion in view and applicability of fractional order derivatives, in the preset investigation, we develop a fractional SIR epidemic model with generalized incidence rate function of the form f(I)S. The fractional derivative is considered in the Caputo and Atangana-Baleanu sense. We simulate the proposed SIR model for different values of relevant parameters and α∈(0,1]. Further, we provide a detail stability analysis of both disease free and endemic equilibriums of the model. The reaming sections of the manuscript are organized as: The basic definitions and relevant results are provided in section 2. Model description in Caputo case and its basic properties are presented in section 3. The stability result of the Caputo model is discussed in section in 4 and its numerical simulations and brief discussion are given in section 5. Model formulation of the ABC case and its existence and uniqueness are presented in section 6. In section 7 the iterative scheme is developed for ABC model and simulation results are presented. finally, concluding remarks are given in section 8.
Here some background material for the Caputo and Atangana-Baleanu fractional derivatives are presented [7,22,23].
Definition 2.1. Consider h∈Cn be a function then the fractional Caputo derivative of order α is defined as:
CDαt(h(t))=1Γ(n−α)∫t0hn(ℑ)(t−ℑ)α−n+1dℑ, |
where α∈(n−1,n) and n∈N. Clearly, CDαt(h(t)) tends to h′(t) as α→1.
Definition 2.2. The corresponding fractional integral is described as:
Iαt(h(t))=1Γ(α)∫t0(t−ℑ)α−1h(ℑ)dℑ. |
In the above integral Γ represents the Gamma function.
Definition 2.3. A fixed point x=x∗ is said to be the equilibrium point (or solution) of the Caputo fractional system defined by:
CDαtx(t)=F(t,x(t)),α∈(0,1), | (2.1) |
if and only if F(t,x∗)=0.
To present the stability analysis of non linear fractional system in the caputo sense via Lyapunov method we first recall the following necessary results from [24,25].
Theorem 2.4. Suppose x∗ be an equilibrium point for the system (2.1) and Υ∈Rn be a domain containing x∗. Let L:[0,∞)×Υ tend to R be a continuously differentiable function such that
M1(x)≤L(t,x(t))≤M2(x), | (2.2) |
and
CDαtL(t,x(t))≤−M3(x), | (2.3) |
∀ 0<α<1 and x∈Υ. Where M1(x), M2(x) and M3(x) are continuously positive definite functions on Υ, then x∗ is uniformly asymptotically stable equilibrium point for the system (2.1).
Definition 2.5. Let F∈H1(b1,b2), b2>b1, α∈(0,1], then Atangana-Beleanu-Caputo (ABC) derivatives is defined as [7]:
ABCb1DαtF(t)=B(α)1−α∫tb1F′(ζ)Eα[−α(t−ζ)α1−α]dζ. | (2.4) |
The term B(α) represents a normalization function such that B(0)=B(1)=1.
The corresponding AB integral is given as:
ABCb1IαtF(t)=1−αB(α)F(t)+αB(α)Γ(α)∫tb1F(ζ)(t−ζ)α−1dζ. | (2.5) |
To formulate the model, the total population is represented with N(t)=S(t)+I(t)+R(t), where S(t) represents susceptible, I(t) represents number of infected and R(t) represents the number of recovered individuals. The complete model that describes a system of fractional differential equation is presented as follows
{CDαtS(t)=Λ−dS−Sf(I)+μ1I+θR,CDαtI(t)=Sf(I)−(d+μ1+μ2+δ)I,CDαtR(t)=μ2I−(d+θ)R. | (3.1) |
The initial condition is given by
S(0)=S0≥0,I(0)=I0≥0,R(0)=R0≥0. |
In the model (3.1), Λ is used for the recruitment rate, d denotes the natural mortality rate in all classes and δ is disease induced death rate. Parameter μ1 shows the rate at which the infected individuals move to the susceptible class and μ2 is the transfer rate from the infected class to the recovered class. The parameter θ is the susceptibility of recovered individuals.
To start the analysis of the system (3.1), it is necessary to present the existence, uniqueness, and non-negativity of the system solution.
To make our study more effective we have assumed that f(I) is to be non-negative and continuously differentiable in the interior of R+ and further, the same hypotheses is taken in account as in [26] that is, f(I) is a real locally Lipschitz function on R+=[0,∞) satisfying
(ⅰ) f(0)=0 and f(I)>0 for I>0;
(ⅱ) f(I)/I is continuous and monotonously non-increasing for I>0 and limI→0+f(I)I exists.
It is easy to analyze from condition (ⅱ) that
f(I)≤βIforI∈R+. |
Then ∫10+1/f(u)du=1 and hence assumption (ⅲ) in [26] is redundant. The incidence rate f(I) is dependent on the concentration of infection. Furthermore, it most famous generalized forms satisfying the above hypothesis are given in the Table 1.
Incidence functions | f(I) | Source |
Bi-linear | βI | [27] |
Saturated | βI1+a1I | [28,29,30] |
Beddington-DeAngelis | βI1+a1S+a2I | [31,32] |
Specific nonlinear | βI1+a1S+a2I+a3SI | [33] |
This section investigates the positivity of the solution of the fractional SIR model system (3.1). Let we consider
R3+={y∈R3∣y≥0}andy(t)=(S(t),I(t),R(t))T. |
To proceed further first we recall the generalized mean values theorem [34].
Lemma 3.1. Consider that h(y)∈C[a1,a2] and CDαth(y)∈(a1,a2], then
h(t)=h(a1)+1Γ(α)(h)(ℑ)(t−a1)α, |
with a1≤ℑ≤t, ∀ t∈(a1,a2].
Corollary 3.2. Consider that h(y)∈C[a1,a2] and CDαth(y)∈(a1,a2], where α∈(0,1], then if
(i).CDαth(y)≥0,∀y∈(a1,a2),thenh(y)isnon−decreasing |
(ii).CDαth(y)≤0,∀y∈(a1,a2),thenh(y)isnon−increasing |
The following main result is presented:
Theorem 3.3. A unique solution y(t) of the system (3.1) will exists and remains in R3+. Furthermore, the solution is positive.
Proof. The existence of the Caputo fractional SIR model can be shown with the help of theorem (3.1) in [35], while the uniqueness of the solution can be easily obtained by making use of the remark 3.2, [35], ∀ t>0. In order to explore the solution non-negativity, it is necessary to show that on each hyperplane bounding the positive orthant, the vector field point R3+. Form the system (3.1) we obtained
CDαtS∣S=0=Λ≥0,CDαtI∣I=0=Sf(I)=0,CDαtR∣R=0=μ2I≥0. |
Therefore, by using the above theorem, we conclude that the solution will stay in R3+
and hence biologically feasible region is constructed as:
Φ={(S,I,R)∈R3+:S,I,R≥0.} |
To obtain the equilibria of the fractional SIR model we solve the system (3.1) at the steady state we have
CDαtS=CDαtI=CDαtR=0. |
then the following is established:
Theorem 3.4. The proposed model (3.1) have at most two equilibria obtained as below:
1. The disease free equilibrium (DFE) expressed by E0=(S0,I0,R0)=(Λd,0,0),
2. and the endemic equilibrium (EE) expressed by E1=(S∗∗,I∗∗,R∗∗),
such that,
S∗∗=(d+μ1+μ2+δ)I∗∗f(I∗∗),R∗∗=μ2I∗∗(d+θ), | (3.2) |
where I∗∗ is a positive solution of the function H,
H(I)=d(d+μ1+μ2+δ)If(I)+(d+δ+dθ(d+θ))I−Λ. | (3.3) |
By condition (ⅱ) H is a strictly increasing function on R+. Furthermore,
limI→0+H(I)=d(d+μ1+μ2+δ)β−Λ and H(Λd)>0. Therefore, H has a positive zero (and only one in this case) iff d(d+μ1+μ2+δ)β−Λ<0 or equivalently
Λβd(d+μ1+μ2+δ)>1. | (3.4) |
The expression for the most biologically important threshold parameter known as the basic reproduction number R0 obtained by the next generation approach which is given as below:
R0=Λd(d+μ1+μ2+δ)∂f(0)∂I, | (3.5) |
where ∂f(0)∂I denote partial dervative of incidence functios with respect to I at disease free state.
In this section we proceed to confirm the stability results in both local and global case. The corrosoponding Jacobian of (3.1) around the DFE is given as
JE0=(−dμ1−S0∂f(0)∂Iθ0S0∂f(0)∂I−(d+μ1+μ2+δ)00μ2−(d+θ)). | (4.1) |
Theorem 4.1. For any p1,p2∈Z, such that gcd(p1,p2)=1. Let α=(p1p2) and consider M=p2, then the DFE of the system (3.1) is locally asymptotically stable (LAS) if |arg(Π)|>π2M, for all roots Π of the characteristic equation (4.2) of the matrix JE0.
det(diag[Πp1Πp1Πp1]−JE0)=0. | (4.2) |
Proof. Evaluating the determinant (4.2), we derived the below equation in term of Π.
(Πp1+d)(Π2p1+a1Πp1+a2)=0, | (4.3) |
where the coefficients are given below:
a1=(d+θ)+(d+δ+μ1+μ2)(1−R0),a2=(d+θ)(d+δ+μ1+μ2)(1−R0). |
The arguments of the roots of the equation Πp1+d=0, are as follow:
arg(Πk)=πp1+k2πp1>πM>π2M,wherek=0,1⋯,(p1−1). | (4.4) |
In similar patren, if R0<1, then the nessarry condition fulfil for all the roots of polynomial (4.3) i.e., |arg(Π)|>π2M. Thus the DFE is LAS for R0<1.
Theorem 4.2. For α∈(0,1) and R0<1 then the DFE, of the fractional SIR model (3.1) is globally asymptotically stable (GAS) and unstable for R0>1.
Proof. To confirm the GAS around the DFE, we formulate the suitable Lyapunov function as:
L(t)=I(t), |
Calculative the time fractional derivative of the above function L(t) we obtain
CDαtL(t)=CDαtI. |
Utilizing system (3.1) we get
CDαtL(t)=Sf(I)−(d+μ1+μ2+δ)I,≤{βS0−(d+μ1+μ2+δ)}I,asf(I)≤βI≤{Λβd−(d+μ1+μ2+δ)}I=(d+μ1+μ2+δ)I{R0−1}. |
Clearly once R0<1 then CDαtL(t) is −ve, and hence, by result (2.4), we ensure that the DFE, E0 is GAS in Φ.
Here, we investigate the global stability of the endemic equilibrium of the fractional SIR model (3.1). At the steady state of the model (3.1) we obtained
{Λ=dS∗∗+S∗∗f(I∗∗)−μ1I∗∗−θR∗∗,(d+μ1+μ2+δ)=S∗∗f(I∗∗)I∗∗,(d+θ)=μ2I∗∗R∗∗. | (4.5) |
Theorem 4.3. If R0>1 and with condition (4.10) satisfy, then the endemic equilibrium of the system (3.1) is globally asymptotically stable.
Proof. We define the following Lyapunov function
V(t)=(S−S∗∗−S∗∗lnSS∗∗)+(I−I∗∗−I∗∗lnII∗∗)+(R−R∗∗−R∗∗lnRR∗∗). |
The differentiation of V(t) with the solution of (3.1) is
CDαtV(t)=(1−S∗∗S)CDαtS(t)+(1−I∗∗I)CDαtI(t)+(1−R∗∗R)CDαtR(t). |
By direct calculations, we have
(1−S∗∗S)CDαtS(t)=(1−S∗∗S)[Λ−dS−Sf(I)+μ1I+θR]=(1−S∗∗S)[dS∗∗+S∗∗f(I∗∗)−μ1I∗∗−θR∗∗−dS−Sf(I)+μ1I+θR]=dS∗∗(2−SS∗∗−S∗∗S)+S∗∗f(I∗∗)(1−Sf(I)S∗∗f(I∗∗)−S∗∗S+f(I)f(I∗∗))+θR∗∗(RR∗∗−1−S∗∗RSR∗∗+S∗∗S)+μ1I∗∗(II∗∗−1−S∗∗ISI∗∗+S∗∗S), | (4.6) |
(1−I∗∗I)CDαtI(t)=(1−I∗∗I)[Sf(I)−(d+μ1+μ2+δ)I]=(1−I∗∗I)[Sf(I)−S∗∗f(I∗∗)II∗∗]=Sf(I)−S∗∗f(I∗∗)II∗∗−Sf(I)I∗∗I+S∗∗f(I∗∗), | (4.7) |
(1−R∗∗R)CDαtR(t)=(1−R∗∗R)(μ2I−(d+θ)R)=(1−R∗∗R)(μ2I−μ2I∗∗R∗∗R)=μ2I∗∗(II∗∗−RR∗∗−R∗∗IRI∗∗+1). | (4.8) |
After some arrangement we have
CDαtV(t)=dS∗∗(2−SS∗∗−S∗∗S)+μ1I∗∗(II∗∗−1−S∗∗S(II∗∗−1))+θR∗∗(RR∗∗−1−S∗∗S(RR∗∗−1))+μ2I∗∗(1−RR∗∗−II∗∗(R∗∗R−1))+S∗∗f(I∗∗)(2−S∗∗S−II∗∗−f(I)f(I∗∗)(SI∗∗S∗∗I−1)). | (4.9) |
It follows from the property arithmetic mean and we have
(2−SS∗−S∗S)≤0, |
and if
(II∗∗−1−S∗∗S(II∗∗−1))≤0,(RR∗∗−1−S∗∗S(RR∗∗−1))≤0,(1−RR∗−II∗(R∗R−1))≤0,(2−S∗∗S−II∗∗−f(I)f(I∗∗)(SI∗∗S∗∗I−1))≤0, | (4.10) |
then it ensures that the model given by (3.1) at the EE case is GAS when R0>1 and with the condition (4.10) satisfied.
In the current section, we investigate the numerical results of the Caputo SIRS model given by system (3.1) with general incidence rate. The model in each case is solved numerically using predictor corrector of molten type technique. In order to perform the numerical experiment and analyze the importance of incidence function on disease dynamics we assign different incidence rate to the generalized incidence function f and simulate the model with various fractional values of α. The parameter values considered in the numerical solution for the fractional SIRS model are depicted in each Figure caption. Initially, we consider the incidence rate f in the form as a bilinear incidence rate i.e., f(S,I)=βSI and obtain the graphical results in sub-plots (a-c) of Figure 1, while in sub-plot (d), we depicted the phase portrait analysis of the model variables. We can observed that the by decreasing the fractional order α the population in suspectable class is increasing whereas the population of infected and recovered classes is decreasing sharply. The graphical dynamics of the model (3.1) for various values of α with non-linear incidence function is shown in sub-plots (a-c) of Figure 2. The corresponding phase portrait analysis in this case is depicted in sub-plot (d) of Figure 2. The influence of parameter m (the half saturation constant) for different values of α=1,0.8,0.5, is shown in Figures 3, 4 and 5 respectively. We can see that population in the infected class is decreasing sharply with the increasing values of m. The influence of m is more significant biologically for fractional order of ABC operator rather than integer case. Finally, the dynamics of the Caputo model (3.1) with Beddington-DeAngelis incidence function with various fractional order α of the derivative are depicted in Figure 6. The corresponding phase portrait analysis in this case is presented in sub-plot (d) of Figure 6. These graphical interpretations explore the importance of the fractional derivative to the proposed model where one can see the behavior of the model for any non-integer case.
In this part of the paper we reconstruct the proposed model using the generalized fractional operator known as Atangana-Baleanu-Caputo (ABC) derivative in order to explore the dynamics of the model in more better way. The model in ABC sense is given by the following system
{ABC0DαtS(t)=Λ−dS−Sf(I)+μ1I+θR,ABC0DαtI(t)=Sf(I)−(d+μ1+μ2+δ)I,ABC0DαtR(t)=μ2I−(d+θ)R, | (6.1) |
with some appropriate nonnegative inial conditions. Before to analyze the dynamics of fractional model (6.1), first we study the existence and uniqueness od the model in the following theorem.
We apply the fixed point theory in order to explore the existence and uniqueness of the solution of fractional model (6.1). To do this we firstly re-write the system (6.1) in the following form:
{ABC0Dαtu(t)=G(t,u(t)),u(0)=u0,0<t<T<∞. | (6.2) |
The vector u(t)=(S,I,R) and G in the system (6.2) represent the state variables and a continuous vector function respectively defined as follows:
G=(G1G2G3)=(Λ−dS−Sf(I)+μ1I+θR,Sf(I)−(d+μ1+μ2+δ)I,μ2I−(d+θ)R), |
with the initial condition by u0(t)=(S(0),I(0),R(0)). Furthermore, the G satisfies the Lipschitz condition given as:
‖G(t,u1(t))−G(t,u2(t))‖≤M‖u1(t)−u2(t)‖. | (6.3) |
To confirm the existence and uniqueness for the solution of system (6.1) with fractional order, we proceed as in the following theorem.
Theorem 6.1. The exists a unique solution of the SIR epidemic model (6.1) if the below condition holds:
(1−α)ABC(α)M+αABC(α)Γ(α)MTαmax<1. | (6.4) |
Proof. On applying the ABC fractional integral to both sides of the model (6.2), we get the non-linear Volterra integral equation:
u(t)=u0+1−αABC(α)G(t,u(t))+αABC(α)Γ(α)∫t0(t−ζ)α−1G(ζ,u(ζ))dζ. | (6.5) |
Consider J=(0,T), and Ξ:C(J,R3)→C(J,R3) stated as
Ξ[u(t)]=u0+1−αABC(α)G(t,u(t))+αABC(α)Γ(α)∫t0(t−ζ)α−1G(ζ,u(ζ))dζ. | (6.6) |
Eq (6.5) becomes:
u(t)=Ξ[u(t)]. | (6.7) |
The supremum norm on u denoted by ‖.‖J and is defined as:
‖u(t)‖J=supt∈J‖u(t)‖,u(t)∈C. | (6.8) |
Clearly, C(J,R3) with the norm ‖.‖J construct a Banach space. Further, the following inequality fulfil
‖∫t0K(t,ζ)u(ζ)dζ‖≤T‖K(t,ζ)‖J‖u(t)‖J, | (6.9) |
with w(t)∈C(J,R3), H(t,ζ)∈C(J2,R) such that
‖K(t,ζ)‖J=supt,ζ∈J|K(t,ζ)|. | (6.10) |
Applying result in (6.7), we have
‖Ξ[u1(t)]−Ξ[u2(t)]‖J≤‖(1−α)ABC(α)(G(t,u1(t))−G(t,u2(t))+αABC(α)Γ(α)×∫t0(t−ζ)α−1(G(ζ,u1(ζ))−G(ζ,u2(ζ)))dζ‖J. | (6.11) |
Furthermore, taking in account equations (6.3), (6.9) and triangular inequality, we have:
‖Ξ[u1(t)]−Ξ[u2(t)]‖J≤((1−α)MABC(α)+αΓ(α)ABC(α)MTαmax)‖u1(t)−u2(t)‖J. | (6.12) |
Finally, we derived the following expression
‖Ξ[u1(t)]−Ξ[u1(t)]‖J≤ϝ‖u1(t)−u2(t)‖J, | (6.13) |
where,
ϝ=(1−α)KABC(α)+αABC(α)Γ(α)KTαmax. |
The operator Ξ is known as a contraction if it follows the condition in equation (6.4). Hence, the solution of the model (6.2) will also has a unique solution.
Here firstly, we obtain the iterative scheme for the system (6.1). The approach developed in [36] is applied in order to approximate the ABC integral. To proceed further, using the fundamental theorem on the system (6.2), we obtained
u(t)−u(0)=(1−α)ABC(α)G(t,u(t))+αABC(α)×Γ(α)∫t0G(ζ,u(ζ)(t−ζ)α−1dζ. | (7.1) |
At t=tm+1, m=0,1,2,..., we have
u(tm+1)−u(0)=1−αABC(α)G(tm,u(tm))+αΓ(α)×ABC(α)∫tm+10G(ζ,u(ζ))(tm+1−ζ)α−1dζ,=1−αABC(α)G(tm,u(tm))+αABC(α)×Γ(α)m∑j=0∫tj+1tjG(ζ,u(ζ))(tm+1−ζ)α−1dζ. | (7.2) |
Utilizing the interpolation polynomial approach, the function G(ζ,u(ζ)) can be approximated over [tj,tj+1]:
G(ζ,u(ζ))≅Pk(ζ)=G(tj,u(tj))h(ζ−tj−1)−G(tj−1,u(tj−1))h(ζ−tj). | (7.3) |
Equation (7.2) becomes:
u(tm+1)=u(0)+1−αABC(α)G(tm,u(tm))+αABC(α)×Γ(α)m∑j=0(G(tj,u(tj))h∫tj+1tj(ζ−tj−1)(tm+1−ζ)α−1dζ−G(tj−1,u(tj−1))h∫tj+1tj(ζ−tj)(tm+1−ζ)α−1dζ). | (7.4) |
On further simplifications we finally arrived the approximate solution as below:
u(tm+1)=u(t0)+1−αABC(α)G(tm,u(tm))+αABC(α)m∑j=0(hαG(tj,u(tj))Γ(α+2)((m−j+1)α(m+α−j+2)−(m−j)α(m+2α−j+2))−hαG(tj−1,u(tj−1))Γ(α+2)((m−j+1)α+1−(m−j)α(m+α−j+1))). | (7.5) |
Hence, for the presented fractional system, we have the following recursive formulae:
S(tm+1)=S(t0)+1−αABC(α)G1(tm,u(tm))+αABC(α)m∑j=0(hαG1(tj,u(tj))Γ(α+2)((m−j+1)α(m+α−j+2)−(m−j)α(m+2α−j+2))−hαG1(tj−1,u(tj−1))α(α+2)((m−j+1)α+1−(m−j)α(m+α−j+1))),I(tm+1)=I(t0)+1−αABC(α)G2(tm,u(tm))+αABC(α)m∑j=0(hαG2(tj,u(tj))α(α+2)((m+1−j)α(m+α−j+2)−(m−j)α(m+2α−j+2))−hαG2(tj−1,u(tj−1))α(α+2)((m−j+1)α+1−(m−j)α(m+α−j+1))),R(tm+1)=R(t0)+1−αABC(α)G3(tm,u(tm))+αABC(α)m∑j=0(hαG3(tj,u(tj))α(α+2)((m−j+1)α(m−j+2+α)−(m−j)α(m+2α−j+2))−hαG3(tj−1,u(tj−1))α(α+2)((m−j+1)α+1−(m−j)α(m+α−j+1))). | (7.6) |
Based on numerical scheme derived in (7.6), we present the graphical results of the proposed ABC fractional model (6.1) in Figures 7–9. To analyze the influence of fractional order and incidence function on the disease dynamics and control we simulate the ABC SIRS model (6.1) with different incidence rates as in Caputo case. In Figure 7, we have shown the dynamics of the model by taking bilinear incidence rate in account. The graphical behavior of the model with nonlinear saturated incidence rate and for different values of α∈(0,1] is analyzed in Figure 8. Finally, the graphical results of fractional SIRS model by incorporating Beddington-DeAngelis type incidence function for various values of α are depicted in Figure 9.
From the graphical analysis, we observed that the population in infected class decrease significantly by decreasing the fractional order and the disease endemic state move to disease free state in each case. Further, from the graphical results of ABC model reveal that this operator is more helpful to explore the model dynamics and provide biologically more feasible results in compression with Caputo case.
In this paper, we analyzed the dynamics of a fractional SIR model with nonlinear generalized incidence rate. Initially, we formulated a general SIR model by apply the Caputo derivative. The local and global stability for the fractional SIR model are obtained and discussed in detailed. It is shown that the model at the disease free case is locally asymptotically stable when the basic reproduction number less than 1. Further, we obtain that the fractional SIR model at the disease free case is globally asymptotically stable for the case of basic reproduction number less than 1. We used fractional Lyapunov function for the proof of the global stability of the disease free case and also we proved by defining the fractional Lyapunov function for the proof of the global stability of the endemic equilibrium when the basic reproduction number greater than 1. Numerical results with various values of the fractional order parameter are obtained and discussed in detailed by choosing some different type of incidence functions. Since the Caputo fractional operator has singular kernel and can not explore the crossover behavior of some biological phenomena. Therefore, we further reformulated the proposed model using ABC derivative having nonsingular and nonlocal kernel. The existence and uniqueness of the ABC model is proved via Picard-Landelof approach. An efficient numerical scheme is used to carried out the iterative solution of the model and perform various simulation results r to analyzed the role of different incidence function in disease dynamic and control. These results indicates that the fractional order model in ABC case is more beneficial and important than that of the ordinary integer order model.
All authors declare no conflict of interest in this paper.