Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Analysis of an SQEIAR epidemic model with media coverage and asymptomatic infection

  • An SQEIAR model with media coverage and asymptomatic infection is proposed for populations with a certain level of immunity. Firstly, we discuss the extinction and persistence for the diseases of the model by using basic reproduction number RC. Then the parameter threshold is analyzed and the effect of parameters on the basic reproduction number is discussed. Furthermore, the optimal media coverage strategy and quarantine strategy for optimal problems under quadratic cost function are derived by applying Pontryagin's Maximum Principle.

    Citation: Xiangyun Shi, Xiwen Gao, Xueyong Zhou, Yongfeng Li. Analysis of an SQEIAR epidemic model with media coverage and asymptomatic infection[J]. AIMS Mathematics, 2021, 6(11): 12298-12320. doi: 10.3934/math.2021712

    Related Papers:

    [1] Xiaoying Pan, Longkun Tang . A new model for COVID-19 in the post-pandemic era. AIMS Mathematics, 2024, 9(8): 21255-21272. doi: 10.3934/math.20241032
    [2] Muhammad Altaf Khan, E. Bonyah, Yi-Xia Li, Taseer Muhammad, K. O. Okosun . Mathematical modeling and optimal control strategies of Buruli ulcer in possum mammals. AIMS Mathematics, 2021, 6(9): 9859-9881. doi: 10.3934/math.2021572
    [3] Chang Hou, Qiubao Wang . The influence of an appropriate reporting time and publicity intensity on the spread of infectious diseases. AIMS Mathematics, 2023, 8(10): 23578-23602. doi: 10.3934/math.20231199
    [4] Danni Wang, Hongli Yang, Liangui Yang . Research on nonlinear infectious disease models influenced by media factors and optimal control. AIMS Mathematics, 2024, 9(2): 3505-3520. doi: 10.3934/math.2024172
    [5] Yue Wu, Shenglong Chen, Ge Zhang, Zhiming Li . Dynamic analysis of a stochastic vector-borne model with direct transmission and media coverage. AIMS Mathematics, 2024, 9(4): 9128-9151. doi: 10.3934/math.2024444
    [6] Shifan Luo, Dongshu Wang, Wenxiu Li . Dynamic analysis of a SIV Filippov system with media coverage and protective measures. AIMS Mathematics, 2022, 7(7): 13469-13492. doi: 10.3934/math.2022745
    [7] Xavier Bardina, Marco Ferrante, Carles Rovira . A stochastic epidemic model of COVID-19 disease. AIMS Mathematics, 2020, 5(6): 7661-7677. doi: 10.3934/math.2020490
    [8] Houssine Zine, Abderrahim El Adraoui, Delfim F. M. Torres . Mathematical analysis, forecasting and optimal control of HIV/AIDS spatiotemporal transmission with a reaction diffusion SICA model. AIMS Mathematics, 2022, 7(9): 16519-16535. doi: 10.3934/math.2022904
    [9] Puntipa Pongsumpun, Jiraporn Lamwong, I-Ming Tang, Puntani Pongsumpun . A modified optimal control for the mathematical model of dengue virus with vaccination. AIMS Mathematics, 2023, 8(11): 27460-27487. doi: 10.3934/math.20231405
    [10] Moh. Mashum Mujur Ihsanjaya, Nanang Susyanto . A mathematical model for policy of vaccinating recovered people in controlling the spread of COVID-19 outbreak. AIMS Mathematics, 2023, 8(6): 14508-14521. doi: 10.3934/math.2023741
  • An SQEIAR model with media coverage and asymptomatic infection is proposed for populations with a certain level of immunity. Firstly, we discuss the extinction and persistence for the diseases of the model by using basic reproduction number RC. Then the parameter threshold is analyzed and the effect of parameters on the basic reproduction number is discussed. Furthermore, the optimal media coverage strategy and quarantine strategy for optimal problems under quadratic cost function are derived by applying Pontryagin's Maximum Principle.



    It is widely know that the spread of infectious diseases still poses a great threat to human life. For example, the H5N1 and H7N9 avian influenza viruses occasionally infect cells in the human upper respiratory tract and have the potential to transform into pandemic causing strains of influenza [1,2]. The outbreak of SARS (severe acute respiratory syndrom) virus in China, sent the human society into a panic [3,4]. The outbreak of COVID-19 has been announced as a Public Health Emergency of International Concern[5], which brought great losses to the society and people [6,7,8].

    As an important tool to analyze and control the spread of infectious diseases, mathematical modeling is used to predict how the disease will spread in the future, taking into account the factors that affect the spread of the disease, which is very important in the early stages of an epidemic, when no treatment or vaccination programme has been developed and therefore drug interventions are often not possible [9]. Mathematical models address the need for understanding the transmission dynamics and other significant factors of the disease that would aid policymakers to make accurate decisions and reduce the rate of transmission of the disease. Many authors have developed mathematical models to study the spread of COVID-19. In [8], Rahim ud Din et al. used an SEIR epidemic model to describe the dynamics of COVID-19 under convex incidence rate. In [10], Isys F. Mello et al. employed a spin S=1/2 Ising-like model and a (logistic) Fermi-Dirac-like function to describe the spread of COVID-19. Especially, Hafeez Aderinsayo Adekola et al. [11] examined various forms of mathematical models that are relevant for the containment, risk analysis and features of COVID-19.

    When a disease breaks out in a human population, people's response to the threat of the disease depends on their perception of the risk, and this perception is influenced by public and private information that is widely disseminated by the media [12,13,14]. Media coverage of outbreaks allows people to understand the extent of the risk and the prevention needs in the areas at risk, and encourages the public to take preventive measures such as wearing masks, avoiding public places, avoiding travel when sick and washing hands frequently. Extensive news coverage and the rapid flow of information can have a profound psychological impact on public health [15]. Media coverage can lower the rate of human exposure and the probability of transmission, Many researchers have used mathematical models to study the influence of media coverage [16,17].

    In particular, novel coronavirus, which has the highest proportion of the asymptomatic among all epidemic diseases, therefore, in addition to the symptomatic, COVID-19 is likely to excrete a large amount of virus and carry out hidden transmission during the asymptomatic phase. Individuals testing positive in serological tests or blood tests for disease without symptoms is referred to as asymptomatic [11]. A simple SEIAR model was used to explain the transmission dynamics of the swine flu outbreak in 2009 at a residential school in Maharashtra, India [18]. Therefore, the recessive transmission of asymptomatic may be one of the main factors in the spread of the epidemic disease.

    Besides, quarantine is essential measure in the spread of epidemic diseases. There are various forms of quarantine, such as home quarantine, quarantine of suspected infected persons and quarantine of infected persons. SARS is a good example, which effectively controlled the outbreak through quarantine [3,4]. The key role played by the quarantine to block the spread of COVID-19 in terms of an interacting parameter between people. In [10], Isys F. Mello et al. presented the key role played by the quarantine to block the spread of COVID-19 in terms of an interacting parameter between people. Therefore, we will establish an SQEIAR model considering media coverage, the asymptomatic and quarantine.

    The emergence of new epidemics always imposes a certain financial burden on society, which has been widespread concerned in society. Therefore, it is of great significance to prevent and control infectious diseases at the minimum cost successfully. In this paper, we further consider the role of media coverage and quarantine in the spread of infectious diseases and establish an optimal control model. We will discuss the media coverage strategy and quarantine strategy for the optimal control problems with quadratic cost functions by applying Pontrygin maximum principle [19,20] and Hamilton-Jacobi-Bellman equation.

    The paper is organized as follows. In the next section, an SQEIAR model is proposed and its basic properties are described. The the local and global stability of the disease-free equilibrium, as well as the local stability and uniform persistence of the endemic equilibrium are discussed in Section 3. In Section 4, threshold analysis is presented. In Section 5, the model is simulated numerically. In Section 6, the optimal control system is analyzed and simulated numerically. Lastly, the model is summarized and discussed.

    In epidemic disease models, the incidence function is extremely important as it determines the epidemic trend of the disease [21]. In many epidemiological models, the bilinear incidence rate βSI and standard rates βSI/N are the most commonly used to characterize the incidence of infectious diseases, where β represents the transmission rate of the disease. Considering the media coverage, Cui and Zhu [16,22] adopted function β(I)=βemI to describe the media impact factors and proposed a new transmission rate function. Liu established the EIH model proposed in [12], the transmission coefficient was assumed to be βe(α1E+α2I+α3H). In [23], Sahu established the epidemic model and took βem1I+m2HNS(I+ηH)N as the transmission term of the disease. In this paper, the total population at time t, denoted by N(t), which has been sub-divided into six mutually exclusive compartments of the susceptible S(t), the quarantined Q(t), the exposed E(t), the infectious I(t), the asymptomatic A(t) and the recovered R(t), so that N(t) = S(t)+Q(t)+E(t)+I(t)+A(t)+R(t). The susceptible transforms into the quarantined, at the same time part of the quarantined backs into the susceptible, some of the susceptible gets into the exposed, at the same time part of the exposed is divided into two major categories of the symptomatic and the asymptomatic, finally some of the two classes of the infected get into the recovered. The population recruitment rate is Λ, population in all compartments decreases at rate μ due to nature death. Inspired by the literature above, in consideration of the media coverage, we adopt βem1I+m2ANS(I+θA)N as the transmission rate function to conduct research, where parameter m1 and m2 represent vector coverage coefficients for symptomatic and asymptomatic infected persons, respectively. The schematic diagram of SQEIAR is shown as following.

    The following model is established based on the above assumptions

    {dSdt=βem1I+m2ANS(I+θA)NpS+λQ+ΛμS,dQdt=pSλQμQ,dEdt=βem1I+m2ANS(I+θA)NσEμE,dIdt=σρEγIIμIIμI,dAdt=σ(1ρ)EγAAμA,dRdt=γII+γAAμR. (2.1)

    The initial conditions of system (2.1) is: S(0)=S0>0, Q(0)=Q0>0, E(0)=E0>0, I(0)=I0>0, A(0)=A0>0, R(0)=R0>0. The description of all the parameters is summarized in Table 1. Note that

    dNdt=ΛμNμII.
    Table 1.  Description of parameters for the system (2.1).
    Parameter Description Unit
    Λ Recruitment rate day1
    μ Natural death rate day1
    β The effective contact rate day1
    θ Modification parameter for reduction in the symptomatic of
    the asymptomatic
    m1 Coefficients of media coverage corresponding to I
    m2 Coefficients of media coverage corresponding to A
    p Progression rate from susceptible to quarantined class day1
    σ Progression rate from exposed to infectious class day1
    ρ Progression rate from infectious to symptomatic infectious class day1
    λ Progression rate from quarantine to susceptible class day1
    μI Disease-induced death rate for symptomatic infectious individuals day1
    μA Disease-induced death rate for asymptomatic infectious individuals day1
    γI Recovery rate for symptomatic infectious individuals day1
    γA Recovery rate for asymptomatic infectious individuals day1

     | Show Table
    DownLoad: CSV

    We only consider the solutions with initial conditions in the feasible region

    Γ={(S,Q,E,I,A,R)6+:S,Q,E,I,A,R0,S+Q+E+I+A+RΛμ}.

    In the feasible region, the existence and uniqueness of the general solution and the continuation result are all valid. It can be proved that the feasible region Γ is a positive invariant for (2.1).

    Proposition 2.1. All the solution trajectories of system (2.1) initiating inside Γ approach enter or stay within the interior of Γ.

    Proof. Let 6+={(S,Q,E,I,A,R)6:S0,Q0,E0,I0,A0,R0} denote the non-negative cone in six-dimensional Euclidean space. From the system (2.1), we observe that

    dSdt|S=0=Λ+λQ>0,dQdt|Q=0=pS0,dEdt|E=0=βem1I+m2AN>0,
    dIdt|I=0=σρE0,dAdt|A=0=σ(1ρ)E0,dRdt|R=0=rII+rAA0.

    Since S(t),Q(t),E(t),I(t),A(t),R(t) are continuous with respect to t, and the vector field on each bounding hyperplane of 6+ is pointing inward direction of 6+, so all the solution trajectories initiating in 6+ will remain inside 6+ for all the time. Thus, the fact that 6+ is positively invariant for the system (2.1) is established. And the total population N(t) satisfies dNdt=ΛμNμII. Then, dNdtΛμN. Applying Birkhoff's and Rota's theorems on differential inequality [24,25], as t, we have 0N(t)Λμ=N0. Hence the solution of system (2.1) is bounded and any solution of the system which originated from Γ remains in Γ. Thus we can find that the solution of system (2.1) is bounded and independent of the initial condition. So the feasible region Γ is an invariant set.

    Note that, Γ is a biologically significant region for COVID-19 or other models of the same type of infectious disease. For COVID-19, as S0,Q0,E0,I0,A0,R0, i.e., the feasible region Γ is non-negative and bounded.

    Reduce the system (2.1) into non-dimensional form by using

    ˜S=SN,˜Q=QN,˜E=EN,˜I=IN,˜A=AN,˜R=RN,˜N=NN0,N0=Λμ,˜t=μt.

    Since ˜S=1˜Q˜E˜I˜A˜R, then

    d˜Sd˜t=˜βem1˜Im2˜A˜S(˜I+θ˜A)˜p˜S+1˜N˜S˜S˜Nd˜Nd˜t.

    The system (2.1) can be transformed the equivalent non-dimensional system

    {d˜Qd˜t=˜p˜S˜λ˜Q˜Q˜Q˜Nd˜Nd˜t:=f1,d˜Ed˜t=˜βem1˜Im2˜A˜S(˜I+θ˜A)˜σ˜E˜E˜E˜Nd˜Nd˜t:=f2,d˜Id˜t=˜σ˜ρ˜E˜γI˜I˜μI˜I˜I˜I˜Nd˜Nd˜t:=f3,d˜Ad˜t=˜σ(1˜ρ)˜E˜γA˜A˜A˜A˜Nd˜Nd˜t:=f4,d˜Rd˜t=˜γI˜I+˜γA˜A˜R˜R˜Nd˜Nd˜t:=f5,d˜Nd˜t=1(1+˜μI˜I)˜N:=f6 (2.2)

    with the initial conditions:

    ˜Q(0)=˜Q0>0,˜E(0)=˜E0>0,˜I(0)=˜I0>0,˜A(0)=˜A0>0,˜R(0)=˜R0>0,˜N(0)=˜N0>0, (2.3)

    where

    ˜β=βμ,˜λ=λμ,˜σ=σμ,˜ρ=ρμ,˜μI=μIμ,˜γA=rAμ,˜γI=rIμ,˜p=pμ.

    In the next sections, we will study the dynamic behavior of the system (2.2) with the initial condition (2.3).

    In this section, we calculate the basic reproduction number of the system and all feasible equilibrium states. Observe that the feasible region for the non-dimensional system (2.2) is

    Ω={(˜Q,˜E,˜I,˜A,˜R,˜N):0˜Q,˜E,˜I,˜A,˜R,˜N1}.

    For the system (2.2), Ω is positively invariant. We only consider the solution of the system (2.2) with initial conditions inside the feasible region Ω.

    The disease-free equilibrium of the system (2.2) is E0=(0,0,0,0,0,1). The local stability of the disease-free equilibrium E0 will be explored by using the basic reproduction number RC. The non-negative matrix F of the new infection terms, and the matrix V of the remaining terms are given by the following formula

    F=(˜p˜S˜βem1˜Im2˜A˜S(˜I+θ˜A)00)

    and

    V=(˜Q(1+1˜N˜μI˜I)+˜Q(˜λ+1)˜E(1+1˜N˜μI˜I)+˜E(˜σ+1)˜σ˜ρ˜E+˜I(1+1˜N˜μI˜I)+˜I(˜γI+˜μI+1)˜σ(1˜ρ)˜E+˜A(1+1˜N˜μI˜I)+˜A(˜γA+1)).

    The corresponding linearized matrices at the disease-free equilibrium E0 are

    F=(˜p˜p˜p˜p00˜β˜β˜θ00000000)

    and

    V=(˜λ+10000˜σ+1000˜σ˜ρ˜γI+˜μI+100˜σ(1˜ρ)0˜γA+1).

    It follows that

    FV1=((˜p)˜λ+1(˜p)˜σ+1(1+1˜γI+˜μI+1+˜σ(1˜ρ)˜γA+1)(˜p)˜γI+˜μI+1(˜p)˜γA+10˜σ˜β˜σ+1(˜ρ˜γI+˜μI+1+(1˜ρ)θ˜γA+1)˜β˜γI+˜μI+1˜βθ˜γA+100000000).

    Then, by applying RC=ρ(FV1), where ρ is the spectral radius, the basic reproduction number RC is given by

    RC=˜β(1+b3θ)b2(˜σ+1)+b1˜β(1+b3θ)=˜β(˜λ+1)[˜σ˜ρ(˜γA+1)+˜σ(1˜ρ)θ(˜γI+˜μI+1)](˜σ+1)(˜γA+1)(˜γI+˜μI+1)(˜p+˜λ+1),

    where

    b1=˜p˜σ˜ρ(˜σ+1)(˜rA+1)(˜γI+˜μI+1)˜σ˜ρ˜β(˜λ+1)[˜σ˜ρ(˜γA+1)+θ˜σ(1˜ρ)(˜γI+˜μI+1)],b2=˜γI+˜μI+1˜σ˜ρ,b3=˜σ(1˜ρ)b2˜rA+1.

    Using Theorem 2 in [26], we can get the following result.

    Theorem 3.1. The disease-free equilibrium E0 of the system (2.2) is locally asymptotically stable if RC<1 and E0 is unstable if RC>1.

    The basic reproduction number represents the average number of people who will catch the disease from a single infected person at the beginning of the disease [26]. If RC>1, the disease will be persistent and become endemic, otherwise, RC<1, the disease will be extinct.

    The possible positive equilibrium of the system (2.2) is obtained by solving the nonlinear equations which is obtained by equipping the derivatives of the system (2.2) to zero. The endemic equilibrium ˉE=(Q,E,I,A,R,N) of the model (2.2) is given by

    Q=b1emI,E=b2I,A=b3I,R=b4I,N=11+˜μII,

    where

    b4=γI˜σ˜ρ(˜γA+1)+˜rA˜σ(1˜ρ)(˜γI+˜μI+1),b5=1+b2+b3+b4,m=m1+m2b3.

    The value of I is given by

    1b5I=emIRC. (3.1)

    Suppose there is no media coverage, i.e., m=0, we obtain I=1b5(11RC). Therefore, I exists in a positive value if and only if RC>1, and it is the unique endemic equilibrium ˉE. Or else, the value of I is given from (3.1). Now, we establish the existence of I for RC>1. In Figure 1, we plot the curve emIRC and straight line 1b5I against I in the range [0, 1]. Note that RC<1, system (2.2) does not exist the endemic equilibrium. And when RC>1, I exists uniquely at positive level and hence the unique endemic equilibrium exists in this case (see Figure 1(b)). From the previous discussion, we can conclude the following Theorem 3.2.

    Figure 1.  (a) shows that non-existence of I when RC<1; (b) shows that the existence of I when RC>1.

    Theorem 3.2. The system (2.2) has no endemic equilibrium if RC1. And the system (2.2) has a unique endemic equilibrium ˉE if RC>1.

    In the following, we prove the local stability of the endemic equilibrium ˉE.

    Theorem 3.3. The endemic equilibrium ˉE is locally asymptotically stable if RC>1, but close to 1.

    Proof. The Jacobian matrix J0 at disease-free equilibrium is given by

    J0=(˜p˜λ1˜p˜p˜p˜p00˜σ1˜β˜βθ000˜σ˜ρ˜γI˜μI10000˜σ(1˜ρ)0˜γA10000˜γI˜γA1000˜μI001).

    Here, we establish the local stability of endemic equilibrium based on the central manifold theory and taking β as bifurcation parameter [27]. A critical value of bifurcation parameter β at RC=1 is given as

    ˜βc=(˜σ+1)(˜γA+1)(˜γI+˜μI+1)(˜p+˜λ+1)(˜λ+1)[˜σ˜ρ(˜γA+1)˜σ(1˜ρ)θ(˜γI+˜μI+1)].

    Corresponding to the zero eigenvalue, the Jacobian can be easily verified J0 at ˜β=˜βc has a right eigenvector given by W=(w1,w2,w3,w4,w5,w6)T, where

    w1=˜p(˜σ˜ρ+˜γI+˜μI+1)(˜γA+1)+˜p(˜γA+1)w4+˜p˜μ(˜γA+1)w5(˜γA+1)(˜p+˜λ+1),w2=˜γI+˜μI+1,w3=˜σ˜ρ,
    w4=˜σ(1˜ρ)(˜γI+˜μI+1)˜γA+1,w5=˜γI˜σ˜ρ(˜γA+1)+˜γA˜σ(1˜ρ)(˜γI+˜μI+1)˜μ(˜γA+1),w6=˜μI˜σ˜ρ.

    Furthermore, corresponding to the zero eigenvalue, the components of the left eigenvector, V=(v1,v2,v3,v4,v5,v6), must satisfy the equalities VJ0=0 and VW=1, so

    v1=v5=v6=0,
    v2=v4(˜λ+1)[˜σ˜ρ(˜γA+1)+˜σ(1˜ρ)θ(˜γI+˜μI+1)]θ(˜σ+1)(˜γI+˜μI+1)(˜p+˜λ+1),v3=v4(˜γA+1)θ(˜γI+˜μI+1),v4=1w4+B1+B2,

    where

    B1=˜σ˜ρ(˜γA+1)θ(˜γI+˜μI+1),B2=(˜λ+1)[˜σ˜ρ(˜γA+1)+˜σ(1˜ρ)θ(˜γI+˜μI+1)]θ(˜σ+1)(˜p+˜λ+1).

    Use the notations x1˜Q, x2˜E, x3˜I, x4˜A, x5˜R, x6˜N, we have

    a=6k,i,j=1vkwiwj2fk(0,0)xixj,b=6k,i,j=1vkwi2fk(0,0)xiβ.

    When β=βc, substituting the following values of all the second order derivatives evaluated at the disease-free equilibrium E0, we get

    a=v1(2w1w6+2w1w3˜μI)+v2[2w2w6+2w2w3(˜βc+˜μI)+2w2w4(˜βcθ)+2w1w3(˜βc)+2w1w4(˜βcθ)+w23(2˜βc2m1˜βc)+w24(2˜βcθ2m2˜βcθ)+2w3w4(˜βcm2˜βc˜βcθm1˜βcθ)+2w3w5(˜βc)+2w4w5(˜βcθ)]+v3(2w3w6+2w23˜μI)+v4(2w4w6+2w3w4˜μI)+v5(2w5w6+2w3w5˜μI)+v6(2w3w6˜μI)

    and

    b=v2(w3+θw4).

    Finally, substituting the values of V and W in a and b, we obtain

    a=2v4B2B3(˜σ+1)(˜γA+1)(˜p+˜λ+1),b=v4B22θ(˜σ+1)(˜p+˜λ+1)(˜λ+1)(˜γA+1)(˜γI+˜μI+1),

    where

    B3=(˜λ+1)(˜γI+˜μI+1)(˜γA+1)+(˜λ+1)[˜γI˜σ˜ρ(˜γA+1)+˜γA˜σ(1˜ρ)(˜γI+˜μI+1)]+(˜λ+1)[˜σ˜ρ(˜γA+1)+˜σ(1˜ρ)(˜γI+˜μI+1)]+(˜p+˜λ+1)[m1˜σ˜ρ(˜γA+1)+m2˜σ(1˜ρ)(˜γI+˜μI+1)].

    Since a<0 and b>0 at β=βc, by using Theorem 4.1 and Remark 1 stated in [27], there is a transcritical bifurcation at RC=1 and the locally asymptotically stable of the unique endemic equilibrium for RC>1 is established.

    We analyze the globally asymptotical stability of the disease-free equilibrium E0 in this subsection.

    Theorem 3.4. Suppose RC<1 and ˜μI=0. The disease-free equilibrium E0 is globally asymptotically stable.

    Proof. We present the global stability of the disease-free equilibrium by using the method in [28]. When ˜μI=0, we have d˜Nd˜t=1˜N. Then ˜N1 as ˜t. Take ˜N in the limiting case, i.e., ˜N=1, then the system (2.2) reduces to

    {d˜Qd˜t=˜p˜S(˜λ+1)˜Q,d˜Ed˜t=˜βem1˜Im2˜A˜S(˜I+θ˜A)(˜σ+1)˜E,d˜Id˜t=˜σ˜ρ˜E(˜γI+1)˜I,d˜Ad˜t=˜σ(1˜ρ)˜E(˜γA+1)˜A,d˜Rd˜t=˜γI˜I+˜γA˜A˜R. (3.2)

    Let X=(˜Q) and Z=(˜E,˜I,˜A,˜R), here U0=(X0,Z0), where X0=(˜p˜p+˜λ+1) and Z0=(0,0,0,0). We have

    dXd˜t=F(X,Z)=˜p(1˜Q˜E˜I˜A˜R)˜λ˜Q˜Q.

    At Z=Z0, G(X,0)=0. Now dXd˜t=F(X,0)=˜p(˜p+˜λ+1)X, as ˜t, XX0. Thus, X=X0 (=˜Q0=˜p˜p+˜λ+1) is globally asymptotically stable. Hence, the condition (H1) in [28] is satisfied. From (3.2), we get

    dZd˜t=G(X,Z)=BZˆG(X,Z),

    where

    B=(˜σ1˜β˜βθ0˜σ˜ρ˜γI100˜σ(1˜ρ)0˜γA100˜γI˜γA1),

    and

    ˆG(X,Z)=(˜β(˜I+θ˜A)(1em1˜Im2˜A˜S)000).

    Obviously, B is an M-matrix. We have 0<em1˜Im2˜A1, for ˜I0, ˜A0. Therefore, ˆG(X,Z)0 since 0˜S1. Thus both the conditions (H1) and (H2) in [28] are satisfied. Hence, the disease-free equilibrium E0 is globally asymptotically stable if RC<1.

    In the next, the uniform persistence for the system (2.2) will be explored. The system (2.2) is said to be uniformly persistent if there exists a constant c such that any solution (˜Q(t),˜E(t), ˜I(t),˜A(t),˜R(t),˜N(t)) satisfies

    lim inft˜Q(t)c,lim inft˜E(t)c,lim inft˜I(t)c,
    lim inft˜A(t)c,lim inft˜R(t)c,lim inft˜N(t)c

    provided that (˜Q(0),˜E(0),˜I(0),˜A(0),˜R(0),˜N(0))Ω [29]. Similar as the proof of [30], the following theorem for persistence can be stated:

    Theorem 3.5. The system (2.2) is uniformly persistent in Ω if and only if RC>1.

    Proof. From Theorem 3.1, the disease-free equilibrium of the model (3.2) is unstable when RC>1. Applying the uniform persistence result stated in [29], finally it can be proved in a similar way as Proposition 3.3 in [31]. At the beginning, the necessity of RC>1 be deduced from the result of Theorem 3.2 and the asymptotic stability of the disease-free equilibrium. The sufficiency of the condition RC>1 is given by uniformly persistent results in the Theorem 4.3 in [29]. It can be proved that the system (3.2) satisfies the hypothesis (H) of Theorem 4.3 in [29]. Eventually, the proposition is proved by observing that, in (3.2), the necessary and sufficient conditions for uniform persistence in Theorem 4.3 in [29] is equivalent to the disease-free equilibrium is unstable. The consequence shows that in limiting case, ˜Q, ˜E, ˜I and ˜A of the model (3.2) will remain a positive value and the disease will persist.

    In this section, the quarantine parameter ˜p, the progression rate ˜σ from the exposed to the infectious and the recovery rate for the symptomatic infectious ˜γI with respect to RC is measured qualitatively. The threshold analysis on the parameters ˜p is performed by calculating the partial derivatives of the effective reproduction number RC. We can observe that

    RC˜p=˜β(˜λ+1)[˜σ˜ρ(˜γA+1)+˜σ(1˜ρ)θ(˜γI+˜μI+1)](˜σ+1)(˜γA+1)(˜γI+˜μI+1)(˜p+˜λ+1)2.

    Hence

    RC˜p<0.

    From the above discussion, we have the following conclusion: the larger the quarantine parameter is, the better the quarantine effect is. The basic reproduction number RC decreases with the increase of the quarantine parameter ˜p, when RC<1 and the infectious disease will be extinct.

    Theorem 4.1. For the model (2.2), the use of quarantine of the susceptile individuals will have positive population-level impact since RC˜p<0.

    In the same way, the threshold analysis respect to the parameters (1˜ρ) is performed by calculating the partial derivatives of the effective reproduction number RC:

    RC(1˜ρ)=˜β(˜λ+1)[˜σ(˜γA+1)+˜σθ(˜γI+˜μI+1)](˜σ+1)(˜γA+1)(˜γI+˜μI+1)(˜p+˜λ+1).

    It follows that

    RC(1˜ρ)<0(or>0)ifθ<θ(1˜ρ)(orθ>θ(1˜ρ)),

    where

    θ(1˜ρ)=˜γA+1˜γI+˜μI+1.

    Therefore, it can be concluded that if θ<θ(1˜ρ), the quarantine will achieve good results.

    Theorem 4.2. For the model (2.2), the use of quarantine of the infectious individuals will have positive (negative) population-level impact since θ<θ(1˜ρ)(θ>θ(1˜ρ)).

    In this section, we provide numerical simulations to illustrate the results under the values of the biologically feasible parameters which is shown in Table 2. The system (2.1) is simulated by taking initial value S0=390,Q0=58,E0=245,I0=95,A0=45,R0=160. The initial values of the numerical simulation are arbitrary and are used only for the simulation, rather than data from a particular region. In fact, the main theoretical results in this paper are independent of the initial values, which we will also present in the subsequent numerical simulation.

    Table 2.  Parameter values used in the simulation for the system (2.1).
    Parameter Nominal values Unit Parameter Nominal values Unit
    Λ 136 day1 ρ 0.74 day1
    μ 0.0751 day1 p 1/100 day1
    β [0.09, 0.65] day1 λ 0.95 day1
    θ 0.5 μI 0.03521 day1
    m1 0.2 γI 1/6 day1
    m2 0.2 γA 1/14 day1
    σ 1/5.2 day1

     | Show Table
    DownLoad: CSV

    Based on the form of the basic regeneration number of system (2.2) calculated in Section 3, the expression form of the basic regeneration number of system (2.1) can be obtained as follows (we also note it as RC):

    RC=β(λ+μ)[σρ(γA+μ)+σ(μρ)θ(γI+μI+1)]μ(σ+μ)(γA+μ)(γI+μI+μ)(p+λ+μ).

    For β=0.09, the basic reproduction number RC=0.3438, the disease-free equilibrium is locally asymptotically stable, see Figure 2(a), the recovered is represented by the solid green line which is decreasing as t increasing but does not go to zero. When β=0.65, the basic reproduction number RC=2.4833>1, the unique endemic equilibrium is locally asymptotically stable as shown in Figure 2(b), which means that the infectious increases over time but reaches a certain level after a period of time and remains stable. The fraction of the total number of infected individuals as a function of time with different initial values for RC=0.3438<1 is presented by Figure 3(a) and RC=2.4833>1 is presented by Figure 3(b), respectively. We can conclude that when the RC is a definite value, the trend of infected individuals has no change with the different initial values of infectious. The effect of parameter m respect to infected individuals is shown separately in Figure 4, which has shown that media coverage decrease the infectious. The effect of parameter σ on the number of the infectious is shown in Figure 5(a), when σ0.03, the number of infected individuals decrease and tend to 0, while σ0.21, the number of infected individuals increase. For 1ρ, in Figure 6(b), 1ρ has both positive and negative effects on the basic reproduction number, when 1ρ>0.5, RC increases, otherwise, decrease. It can be seen from Figure 6(a) that the basic reproduction number decreases as the quarantine parameter p increases and RC<1 if p>0.7. With γI increasing, the basic regenerative number is decreasing and γI=0.1053 is the critical value.

    Figure 2.  (a) The variation of the scaled population in scaled-time, taking m1=m2=0.2 and β=0.09 with RC=0.3438<1; (b) The variation of the scaled population in scaled-time, taking m1=m2=0.2 and β=0.65 with RC=2.4833>1.
    Figure 3.  (a) The fraction of the total number of infected individuals as a function of time with different initial values for RC=0.3438<1; (b) The fraction of the total number of infected individuals as a function of time with different initial values for RC=2.4833>1.
    Figure 4.  (a) Effect of m1 and m2 on I, when β=0.09 and RC=0.3438<1; (b) Effect of m1 and m2 on I, when β=0.6 and RC=2.4833>1.
    Figure 5.  (a) Effect of σ on I, when β=0.65; (b) Effect of 1ρ on effective reproduction number RC.
    Figure 6.  (a) Effect of quarantine parameter p on effective reproduction number RC; (b) Effect of recovery rate for symptomatic infectious individuals γI on effective reproduction number RC.

    From the above sections, the influence of media coverage and quarantine measures on the number of infected individuals is discussed. To ensure the effectiveness of control while minimizing the cost of various prevention measures is the top concern of the government in the face of the epidemic. Therefore, we add media coverage control and quarantine control to the system (2.1), and the model is shown as follows:

    {dSdt=βeu1(t)I+u2(t)ANS(I+θA)Nu3(t)S+λQ+ΛμS,dQdt=u3(t)SλQμQ,dEdt=βeu1(t)I+u2(t)ANS(I+θA)NσEμE,dIdt=σρEγIIμIIμI,dAdt=σ(1ρ)EγAAμA,dRdt=γII+γAAμR. (6.1)

    The control function u1(t) represents media coverage control for the symptomatic, u2(t) represents media coverage control for the asymptomatic and the control of the quarantine is expressed by u3(t). In order to solve u1(t), u2(t), u3(t), we construct the objective functional J. Then has

    J(u1(t),u2(t),u3(t))=T0(A1I+A2A+B12u21+B22u22+B32u23)dt,

    where A1, A2 represent the positive weight. We need to find an optimal control u1(t), u2(t) and u3(t) such that

    J(u1(t),u2(t),u3(t))=minJ(u1(t),u2(t),u3(t)),(u1(t),u2(t),u3(t))U,

    where

    U={(u1(t),u2(t),u3(t))ui(t)isLebsguemeasurableon[0,1],0ui(t)1,i=1,2,3}

    is the control set.

    Optimal control existence can be demonstrated through a well known classical result: (see [20]), we must inspect that the following axioms are satisfied:

    (H1) Controls set and set of state variables are nonempty.

    (H2) The admissible control set U is convex and closed.

    (H3) R.H.S of state system is bounded by a linear function of the state variables and controls.

    (H4) The objective functional J has convex integrand on U and is bounded below by c1(3i=1|ui|2)τ2c2, where c1,c2>0 and τ>1.

    The existence of solutions for the system is established by using the result given by Lukes [32]. In this way we verify the above hypotheses. (H1) is accomplished because the coefficients are bounded. The boundedness of solutions shows that the set of controls fulfils (H2). Since the system of equations is bilinear in u1, u2 and solutions are bounded. So, R.H.S of system satisfies the criteria (H3). Because the integrand of objective functional is convex and the last condition is satisfied.

    M1I+M2A+12D1u21+12D2u22+12D3u23c1(3i=1ui2)τ2c2,

    where M1, M2, D1, D2, D3, c1, c2>0 and τ>1. Thus we have the following theorem:

    Theorem 6.1. For the objective functional J(u1(t),u2(t),u3(t))=T0(A1I+A2A+B12u21+B22u22+B32u23)dt, where U={(u1(t),u2(t),u3(t))0ui(t)1,i=1,2,3,t[0,T]} subject to system with initial conditions, there exists an optimal control u1(t),u2(t),u3(t), so that J(u1(t),u2(t),u3(t))=minJ(u1(t),u2(t),u3(t)),(u1(t),u2(t),u3(t))U.

    In order to obtain an optimal control u1(t),u2(t),u3(t) to make J(u1(t),u2(t),u3(t))=minJ(u1(t),u2(t),u3(t)),(u1(t),u2(t),u3(t))U. The optimal solution can be obtained by finding the Lagrangian as well as Hamiltonian for the system. The Lagrangian is:

    L(I,A,u1(t),u2(t),u3(t))=A1I+A2A+B12u21+B22u22+B32u23.

    We need establish the Lagrangian minimum value. To achieve this goal, we construct the Hamiltonian function H for the optimal problem as follows:

    H=A1I+A2A+B12u21+B22u22+B32u23+λ1(βeu1I+u2ANS(I+θA)Nu3S+λQ+ΛμS)+λ2(u3SλQμQ)+λ3(βeu1I+u2ANS(I+θA)NσEμE)+λ4(σρEγIIμIIμI)+λ5(σ(1ρ)EγAAμA)+λ6(γII+γAAμR).

    We can find the essential conditions for the optimal control problem by applying the Pontryagin's Maximum Principle [19]. For the system, the solution of optimal control problem, a non-zero vector function λ(t)=(λ1,λ2,λ3,λ4,λ5,λ6) exists and the subsequent conditions is established. Then the state equation, optimality condition and adjoint equation, respectively as follows:

    dxdt=λ(H(t,u1,u2,u3,λ(t))),0=u(H(t,u1,u2,u3,λ(t))),dλdt=x(H(t,u1,u2,u3,λ(t))).

    The essential conditions applied to the Hamiltonian H give the following result:

    Theorem 6.2. Given that (S,Q,E,I,A,R) are optimal state solutions and (u1,u2,u3) are associated optimal control variable for the optimal control problem (6.1), then, there exists adjoint variables λi, for i=1,2,3,4,5,6 which satisfies

    {dλ1dt=(λ1λ3)βeu1I+u2ANS(I+θA)N+(λ1λ3)u3+λ1μ,dλ2dt=(λ2λ1)λ+λ2μ,dλ3dt=(λ3λ5)σ+(λ5λ4)ρ+λ3μ,dλ4dt=A1+(λ1λ3)βeu1I+u2AN(SNu1S(I+θA)N2)+(λ4λ6)γI+λ4(μI+μ),dλ5dt=A2+(λ1λ3)βeu1I+u2AN(SθNu2S(I+θA)N2)+(λ5λ6)γA+λ5μ,dλ6dt=λ6μ, (6.2)

    with the transversality condition λ1(T)=λ2(T)==λ6(T)=0. Additionally, u1,u2,u3 are expressed as

    u1(t)=max{min{1,IB1B2NLambertW(βS(I+θA)(I2B2+A2B1)(λ3λ1)B1B2N3)(I2B2+A2B1)(λ3λ1)λ3λ1B1},0},u2(t)=max{min{1,AB1B2NLambertW(βS(I+θA)(I2B2+A2B1)(λ3λ1)B1B2N3)(I2B2+A2B1)(λ3λ1)λ3λ1B2},0},u3(t)=max{min{1,(λ1λ2)SB3},0}.

    Proof. The adjoint equations and the conditions of transversality are obtained by the Hamilton function H. By putting S=S,Q=Q,E=E,I=I,A=A,R=R and differentiating the Hamiltonian with respect to S,Q,E,I,A and R respectively, we obtain

    {dλ1dt=(λ1λ3)βeu1I+u2ANS(I+θA)N+(λ1λ3)u3+λ1μ,dλ2dt=(λ2λ1)λ+λ2μ,dλ3dt=(λ3λ5)σ+(λ5λ4)ρ+λ3μ,dλ4dt=A1+(λ1λ3)βeu1I+u2AN(SNu1S(I+θA)N2)+(λ4λ6)γI+λ4(μI+μ),dλ5dt=A2+(λ1λ3)βeu1I+u2AN(SθNu2S(I+θA)N2)+(λ5λ6)γA+λ5μ,dλ6dt=λ6μ,

    with the transversality condition λ1(T)=λ2(T)==λ6(T)=0. Further, by solving the equations Hui=0, i=1,2,3 on the interior of the control set and using the optimality condition, we obtain

    {Hu1=B1u1+(λ1λ3)βeu1Iu2ANINS(I+θA)N=0,Hu2=B2u2+(λ1λ3)βeu1Iu2ANANS(I+θA)N=0,Hu3=B3u3λ1S+λ2S=0. (6.3)

    Let's solve u1 by using the first two equations of (6.3), thus

    {B1u1=(λ3λ1)βeu1Iu2ANINS(I+θA)N,B2u2=(λ3λ1)βeu1Iu2ANANS(I+θA)N. (6.4)

    Hence,

    u2=AB1IB2u1. (6.5)

    Putting (6.5) in the first equation of (6.4), we get u1. And then we can easily get u2 and u3 as follows

    u1(t)=max{min{1,IB1B2NLambertW(βS(I+θA)(I2B2+A2B1)(λ3λ1)B1B2N3)(I2B2+A2B1)(λ3λ1)λ3λ1B1},0},u2(t)=max{min{1,AB1B2NLambertW(βS(I+θA)(I2B2+A2B1)(λ3λ1)B1B2N3)(I2B2+A2B1)(λ3λ1)λ3λ1B2},0},u3(t)=max{min{1,(λ1λ2)SB3},0}.

    Therefore, the optimal control problem of (6.1) is finished. Furthermore, some graphs is obtained by numerical simulation.

    In this section, we assess numerically by investigating the effect of control strategies on the transmission dynamics of disease. An iterative scheme of fourth order Runge-Kutta method is used for solving the optimality system (6.1). The common parameter values used in the computations are the same in Table 2 except β, and the initial value are S0=400,Q0=200,E0=300,I0=50,A0=60,R0=35 and β=0.26.

    The system is numerically simulated under the control of media coverage and quarantine. When media control and quarantine control measure are used together to optimize the objection J, it is observed from Figure 7 that the number of symptomatic and asymptomatic individuals are both reduced, while the number of susceptible individuals will increase greatly. The upper bounds for the controls u1,u2 and u3 are 140 days.

    Figure 7.  Simulation of the solutions of system (6.1) with and without control (Figures 7(a)-(f)) and control profile Figure 7(g).

    The SQEIAR model proposed in this paper is based on the media coverage and quarantine. We discussed the dynamic behavior of the proposed model. The disease-free equilibrium of the system (2.2) is locally-asymptotically stable if RC<1, and if RC>1, it is unstable. Then the existence of the endemic equilibrium, the local asymptotic stability and global asymptotic stability of the equilibria of system (2.2) related to basic reproduction number are discussed and the uniformly persistent property of the system is obtained. The threshold analysis is carried out to discuss the influence of different parameters on the effective reproduction number, such as the quarantine parameter p, the recovery rate for symptomatic infections individuals γI, the progression rate parameters σ and 1ρ, which have both positive and negative effects on the basic reproduction number. The increase of quarantine parameter p while the effective reproduction number decrease, which indicates that quarantine measures affect the extinction rate of the disease. The recovery rate parameter γI has the same effect on the basic regeneration number as the quarantine p. While the progression rate parameters 1ρ has both positive and negative effects on the basic reproduction number, when the 1ρ is greater than the threshold, RC increases, otherwise, RC decrease. The parameter σ has both positive and negative effects on the number of the infectious, when the value of σ is less than a certain small value, the number of infected individuals decrease and tend to 0, while σ is bigger than a certain big value, the number of infected individuals increase.

    Furthermore, the optimal control solutions are shown in Figure 7, showing media coverage and quarantine measures are simultaneously more conducive to controlling the spread of epidemic diseases in the early stages. By using the optimal control theory, we obtain the control strategy that minimizes infected and asymptomatic infected persons with the least cost of implementing control measures. Under this kinds of control, the number of symptomatic and asymptomatic individuals are both reduced, while the number of susceptible individuals will increase greatly. Therefore, the suggestion for the control of epidemic disease is to greatly reduce the number of infected persons through positive media influence and effective quarantine.

    The model works well for the spread of COVID-19, based on the basic regeneration number formula given by the model, we can calculate the basic regeneration number of COVID-19 transmission in a certain region, and then predict the transmission trend of COVID-19. The impact of the analyzed model parameters on the basic regeneration number or the number of infected people can provide suggestions for the prevention and treatment of COVID-19 from media reports, quarantine, treatment and recovery. In particular, our model is able to give a control strategy that minimizes the number of COVID-19 infected and asymptomatic infected people at the least cost, which is crucial for economically underdeveloped countries to control COVID-19. Compared with other conventional COVID-19 models, such as SEIR model in the literature [8,33], the proposed model not only adds the category of asymptomatic infection, but also takes into account the impact of media coverage and quarantine on the transmission of COVID-19, which is more consistent with the current actual situation of COVID-19 transmission. This makes our model more practical in revealing the transmission pattern of the COVID-19.

    This work is supported by the National Natural Science Foundation of China (No. 11701495), the Scientific and Technological Key Projects of Henan Province (No. 202102310631), the Key Scientific Research Project in Colleges and Universities of Henan Province of China (No. 19A110036) and Nanhu Scholars Program for Young Scholars of XYNU.

    The authors declare no conflict of interest.



    [1] K. Tharakaraman, A. Jayaraman, R. Raman, K. Viswanathan, N. Stebbins, D. Johnson, et al., Glycan receptor binding of the influenza a virus H7N9 hemagglutinin, Cell, 153 (2013), 1486-1493. doi: 10.1016/j.cell.2013.05.034
    [2] K. Tharakaraman, R. Raman, K. Viswanathan, N. Stebbins, A. Jayaraman, A. Krishnan, et al., Structural determinants for naturally evolving H5N1 hemagglutinin to switch its receptor specificity, Cell, 153 (2013), 1475-1485. doi: 10.1016/j.cell.2013.05.035
    [3] W. D. Wang, S. G. Ruan, Simulating the SARS outbreak in Beijingwith limited data, J. Theoret. Biol., 227 (2004), 369-379. doi: 10.1016/j.jtbi.2003.11.014
    [4] X. F. Yan, Y. Zou, Optimal and sub-optimal quarantine and isolation control in SARS epidemics, Math. Comput. Model., 47 (2008), 235-245. doi: 10.1016/j.mcm.2007.04.003
    [5] J. X. Guan, Y. Y. Wei, Y. Zhao, F. Chen, Modeling the transmission dynamics of COVID-19 epidemic: A systematic review, J. Biomed. Res., 34 (2020), 422-430. doi: 10.7555/JBR.34.20200119
    [6] X. M. Rong, L. Yang, H. D. Chu, M. Fan, Effect of delay in diagnosis on transmission of COVID-19, Math. Biosci. Eng., 17 (2020), 2725-2740. doi: 10.3934/mbe.2020149
    [7] B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Y. Tang, Y. N. Xiao, et al., Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions, J. Clin. Med., 9 (2020), 462. doi: 10.3390/jcm9020462
    [8] Rahim ud Din, Aly R. Seadawy, Kamal Shah, Aman Ullah, Dumitru Baleanu, Study of global dynamics of COVID-19 via a new mathematical model, Results Phys., 19 (2020), 103468. doi: 10.1016/j.rinp.2020.103468
    [9] C. Sun, W. Yang, J. Arino, K. Khan, Effect of media-induced social distancing on disease transmission in a two patch setting, Math. Biosci., 230 (2011), 87-95. doi: 10.1016/j.mbs.2011.01.005
    [10] I. F. Mello, L. Squillante, G. O. Gomes, A. C. Seridonio, M. de Souza, Epidemics, the Ising-model and percolation theory: A comprehensive review focused on COVID-19, Physica A, 573 (2021), 125963. doi: 10.1016/j.physa.2021.125963
    [11] H. A. Adekola, I. A. Adekunle, H. O. Egberongbe, S. A. Onitilo, Idris Nasir Abdullahi, Mathematical modeling for infectious viral disease: The COVID-19 perspective, J. Public Aff., 20 (2020), e2306.
    [12] R. Liu, J. Wu, H. Zhu, Media/psychological impact on multiple outbreaks of emerging infectious diseases, Comput. Math. Methods Med., 8 (2007), 153-164. doi: 10.1080/17486700701425870
    [13] J. H. Pang, J. A. Cui, An SIRS epidemiological model with nonlinear incidence rate incorporating media coverage, Second International Conference on Information and Computing Science, IEEE, (2009), 116-119.
    [14] J. M. Tchuenche, N. Dube, C. P. Bhunu, R. J. Smith, C. T. Bauch, The impact of media coverage on the transmission dynamics of human influenza, BMC Public Health, 11 (2011), S5.
    [15] A. Wang, Y. Xiao, A Filippov system describing media effects on the spread of infectious diseases, Nonlinear Anal. Hybrid Syst., 11 (2014), 84-97. doi: 10.1016/j.nahs.2013.06.005
    [16] J. Cui, Y. Sun, H. Zhu, The impact of media on the control of infectious diseases, J. Dyn. Differ. Equ., 20 (2007), 31-53.
    [17] Y. P. Liu, J. A. Cui, The impact of media coverage on the dynamics of infectious disease, Int. J. Biomath., 1 (2008), 65-74. doi: 10.1142/S1793524508000023
    [18] P. Shil, Mthematical modeling of viral epidemics: A review, Biomed. Res. J., 3 (2016), 195-215.
    [19] L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, E. F. Mishchenko, The mathematical theory of optimal processes, Bell Syst. Tech. J., 27 (1986), 623-656.
    [20] W. H. Fleming, R. W. Rishel, Deterministic and Stochastic Optimal Control, Springer, New York, 1975.
    [21] G. P. Sahu, J. Dhar, Analysis of an SVEIS epidemic model with partial temporary immunity and saturation incidence rate, Appl. Math. Model., 36 (2012), 908-923. doi: 10.1016/j.apm.2011.07.044
    [22] J. Cui, X. Tao, H. P. Zhu, An SIS infection on model in corporating media coverage, Rocky Mt. J. Math., 38 (2008), 1323-1334.
    [23] G. P. Sahu, J. Dhar, Dynamics of an SEQIHRS epidemic model with media coverage, quarantine and isolation in a community with pre-existing immunity, J. Math. Anal. App., 421 (2015), 1651-1672. doi: 10.1016/j.jmaa.2014.08.019
    [24] G. Birkhoff, G. C. Rota, Ordinary Differential Equations, John Wiley, New York, 1998.
    [25] H. L. Smith, P. Waltman, The Theory of the Chemostat, Cambridge University Press, Cambridge, 1995.
    [26] P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29-48.
    [27] C. Castillo-Chavez, B. Song, Dynamical models of tuberculosis and their applications, Math. Biosci. Eng., 1 (2004), 361-404. doi: 10.3934/mbe.2004.1.361
    [28] C. Castillo-Chavez, Z. L. Feng, W. Huang, On the computation of R0 and its role on global stability, Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction, Springer-Verlag, New York, (2002), 229-250.
    [29] H. I. Freedman, S. Ruan, M. Tang, Uniform persistence and flows near a closed positively invariant set, J. Dynam. Differ. Equ., 6 (1994), 583-600. doi: 10.1007/BF02218848
    [30] M. Safi, A. B. Gumel, Dynamics of a model with quarantine-adjusted incidence and quarantine of susceptible individuals, J. Math. Anal. Appl., 399 (2013), 565-575.
    [31] M. Y. Li, J. R. Graef, L. Wang, Global dynamics of a SEIR model with varying total population size, Math. Biosci., 160 (1999), 191-213. doi: 10.1016/S0025-5564(99)00030-9
    [32] D. L. Lukes, Differential Equations: Classical to Controlled, Mathematics in Science and Engineering, Academic Press, New York, 1982.
    [33] Yaqing Fang, Yiting Nie, Marshare Penny, Transmission dynamics of the COVID-19 outbreak and effectiveness of government interventions: A data-driven analysis, J. Med. Virol., 92 (2020), 645-659. doi: 10.1002/jmv.25750
  • This article has been cited by:

    1. Xueyong Zhou, Xiwen Gao, Xiangyun Shi, Analysis of an SQEIAR stochastic epidemic model with media coverage and asymptomatic infection, 2022, 15, 1793-5245, 10.1142/S1793524522500838
    2. Yaping Wang, Lin Hu, Linfei Nie, Dynamics of a Hybrid HIV/AIDS Model with Age-Structured, Self-Protection and Media Coverage, 2022, 11, 2227-7390, 82, 10.3390/math11010082
    3. Xueyong Zhou, Xiangyun Shi, Stability analysis and backward bifurcation on an SEIQR epidemic model with nonlinear innate immunity, 2022, 30, 2688-1594, 3481, 10.3934/era.2022178
    4. Dan Zhou, Xiangyun Shi, Xueyong Zhou, Dynamic Analysis of a Stochastic Delayed SEIRS Epidemic Model with Lévy Jumps and the Impact of Public Health Education, 2023, 12, 2075-1680, 560, 10.3390/axioms12060560
    5. Zohreh Abbasi, Xinzhi Liu, Switched Epidemic Models: A Contact Rate Dependent Analysis, 2024, 128, 0307904X, 168, 10.1016/j.apm.2024.01.016
    6. Lihua Dai, Xianning Liu, Yuming Chen, Global dynamics of a fractional-order SIS epidemic model with media coverage, 2023, 111, 0924-090X, 19513, 10.1007/s11071-023-08838-4
    7. Lin Hu, Linfei Nie, Stability and Hopf Bifurcation Analysis of a Multi-Delay Vector-Borne Disease Model with Presence Awareness and Media Effect, 2023, 7, 2504-3110, 831, 10.3390/fractalfract7120831
    8. Hong Qiu, Rumei Hou, Dynamics and optimal control of an SEIAQR epidemic model with media coverage, 2024, 535, 0022247X, 128192, 10.1016/j.jmaa.2024.128192
    9. Rujie Yang, Hong Qiu, A Valid Transport Related SVEIHR Stochastic Epidemic Model with Coverage and Time Delays, 2024, 23, 2224-2880, 815, 10.37394/23206.2024.23.84
    10. Mohammad Idrees, Nigar Ali, Ihtisham Ul Haq, Imtiaz Ahmad, Mohammed Daher Albalwi, Md. Haider Ali Biswas, Pramita Mishra, Nonlinear Epidemic Model With Media Coverage: Stability, Bifurcation, and Control in a Population With Pre‐Existing Immunity, 2025, 2025, 1076-2787, 10.1155/cplx/1572080
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3128) PDF downloads(132) Cited by(10)

Figures and Tables

Figures(7)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog