# Assessing the Impact of Continuous Vaccination and Voluntary Isolation on the Dynamics of COVID-19: A Mathematical Optimal Control of SEIR Epidemic Model

Authors: Yue Yu; Manman Shi; Manfeng Hu; Jingxiang Zhang
Journal: Computational Intelligence and Neuroscience
 (2022)

---

## Abstract

In order to study the impact of continuous vaccination and voluntary isolation for the COVID-19, a susceptible-exposed-infected-recovered-quarantine-vaccines (SEIR-QV) model is proposed. A basic regeneration numberR0 is defined to determine the extinction or persistence of the disease. We numerically analyze the impact of key parameters based on actual parameters of COVID-19, such as the vaccination rate, population importation rate, and natural (or causal) mortality transmission rate on the dynamics of disease transmission. Then we obtain sensitivity indices of some parameters on R0 by sensitivity analysis. Finally, the stability of the system and the effectiveness of the optimal control strategy are verified by numerical simulation.


---

## Body

## 1. Introduction

Since the outbreak of the coronavirus disease 2019 (COVID-19), it has been significantly impacting the world economy and our lives. Globally, as of 4:00 pm CET, 28 February 2022, there have been 434,154,739 confirmed cases of COVID-19, including 5,944,342 deaths, reported to the WHO. As of 26 February 2022, a total of 10,585,766,316 vaccine doses have been administered.With the emergence of multiple mutant strains of novel coronavirus, including Alpha, Beta, Gamma, Delta, and Omicron, many scholars further explore the transmission and epidemic pattern of COVID-19 from different view. These research results are helpful to learn about the transmission and mode of infectious diseases and provide reliable information for the prediction and control of infectious diseases.In the study of the dynamics of infectious diseases, the susceptibility-infection-recovery (SIR) model was mainly established by Kermack and McKendrick in 1927 with the kinetic method. In order to better study the characteristics of infectious diseases and obtain the best prevention strategies, scholars propose many improved models based on the traditional SIR model for different situations. By introducing time variation factors, Chen et al. [1] present a new SIR model, which has stronger adaptability and robustness in predicting the number of confirmed cases and inflection points of infectious diseases. Wang et al. [2] propose a SIR model with time-varying isolation protocol. In this model, the dynamic evolution of infectious diseases can be comprehensively analyzed and evaluated. In addition, Hota et al. [3] introduce a closed-loop framework combined with the SIR model and expound the significance of early detection through two feasible optimization questions.In view of the latent characteristics of some infectious diseases, SEIR models are further proposed based on the SIR model. Li et al. [4] fit the basic regeneration number curve based on the SEIR model. Trends in COVID-19 outbreaks in China, the USA, India, and Iran are predicted and analyzed. The results show that the rapid and efficient isolation measures adopted in China are significant in suppressing COVID-19. In addition, some new models take into account the influence of different factors on infectious diseases, such as the presence of temperature and humidity [5, 6]. It is concluded that COVID-19 is more infectious and lethal at low temperatures and humidity. Besides, many research studies show that the different epidemic prevention measures, such as mandatory isolation, wearing masks, vaccination, and government control policies, have significant effects on preventing the rapid spread of COVID-19 [7–9].Due to the problem of asymptomatic infected persons in infectious diseases, Yu et al. [10] propose an SEIR-AQ model by considering factors including prevention and control efforts, isolation strategies, and asymptomatic infected persons. The SEIR-AQ model can well anticipate the spread trend of COVID-19, providing technical support for the scientific assessment of the infectious disease situation. Analysis from the SEIR-AQ model shows that prevention and control isolation, medical follow-up isolation, and other measures have a significant inhibitory impact on the spread of COVID-19. In addition, some scholars use deep learning methods to study the impact of the COVID-19 outbreak on human society [11]. Some other scholars have analyzed and studied the transmission trends of the COVID-19 outbreak. The spatiotemporal evolution and transmission trends of the COVID-19 epidemic are analyzed using an ontological modeling approach by Liu et al. [12].Although control measures, vaccination, latent infectivity, and other factors are considered in these studies, the stability and existence of the model are not discussed from the mathematical view point. Therefore, we establish a new susceptible-exposed-infected-recovered-quarantine-vaccines (SEIR-QV) model based on the latent period and the influence of prevention and control measures of COVID-19. In mathematical theory, the extinction or persistence of a disease can be determined based on the basic reproduction number and the Lyapunov function. We numerically simulate and analyze the impact of key real-world parameters of COVID-19. What’s more, our sensitivity analysis indicates that the related parameters have significant effects on the stability and existence of the SEIR-QV model.


## 2. The SEIR-QV Model

In this paper, the susceptible-exposed-infectious-recovered-asymptomatic-quarantine-vaccines (SEIR-AQV) model is developed based on the traditional SEIR model. The SEIR-AQV model takes into account the effect of isolation based on the traditional model and divides the population into unisolated susceptible people (S), unisolated exposed people (E), unisolated infected people (I), isolated susceptible people (Sq), isolated exposed people (Eq), asymptomatic infected people (A), hospitalized people (H), and recovered people (R).The SEIR-AQV model has the following assumptions: first assumption: medical resources are adequate, i.e., isolated exposed people can be directly converted to hospitalized people after diagnosis; second assumption: once infected with COVID-19, one must undergo inpatient treatment to recover, i.e., infected people are unlikely to recover on their own; third assumption: both isolated exposed people and hospitalized people are isolated from the outside world, i.e., they are not infectious; fourth assumption: the infected people will be immune after recovery, i.e., he or she will not become susceptible again after recovery.The bin transformation relationship is shown in Figure1.Figure 1
Bin view of the population transformation relationship of the SEIR-AQV model.The equations of the SEIR-AQV model are as follows:(1)dSdt=Λ−cβδSI+θA+υE+λ2Sq−λ1+ηS,dSqdt=λ1S−λ2+ηSq,dEdt=c1−qβδSI+θA+υE−σ+ηE,dEqdt=cqβδSI+θA+υE−b3+ηEq,dIdt=σeE−b1+α1+ηI,dAdt=σ1−eE−b2+α2+ηA,dHdt=b1I+b2A+b3Eq−r+α3+ηH,dRdt=rH−ηR,where c means the rate of exposure; η means the natural mortality rate; q represents the proportion of isolated; δ represents the vaccine coverage rate; β represents the probability of transmission; θ stands for the ratio of A relative to the transmission capacity of I; υ stands for the ratio of E relative to the transmission capacity of I; b1,b2,b3 stand for the conversion of I, A, andEq, respectively to H; λ1 represent the rate of conversion of S into Sq; λ2 stands for the rate of conversion of Sq into S; σ stands for the ratio of E to I conversion and σ=1/Incubation period; e means the ratio of E to I conversion; α1,α2,α3 mean the rate of conversion of I, A, and the rate of cause-specific death of H; r means the rate of recovery of H; and Λ means the rate of population importation.However, given the reality of the COVID-19 epidemic, the status of isolated susceptible people is unlikely to persist over time, and isolated susceptible people have no translational relationship to hamlets other than unisolated susceptible people. Therefore, the population type of isolated susceptible people and the parameters involved in it are removed. Considering also that the proportion of asymptomatic infected people is very small and has minimal effect on the overall transmission trend of the COVID-19 outbreak, the population type of asymptomatic infected people and its involved parameters are also removed. Therefore, we improved the SEIR-AQV model in the context of the COVID-19 outbreak. We construct a new seasonal susceptible-exposed-infected-removed-quarantine-vaccines (SEIR-QV) model with the population transformation relationship bin view shown in Figure2.Figure 2
Bin view of the population transformation relationship of the SEIR-QV model.The equations of the SEIR-QV model are as follows:(2)dSdt=Λ−cβδSI+υE−ηS,dEdt=c1−qβδSI+υE−σ+ηE,dEqdt=cqβδSI+υE−b2+ηEq,dIdt=σE−b1+α1+ηI,dHdt=b1I+b2Eq−r+α2+ηH,dRdt=rH−ηR,where α1,α2 represent the rate of cause-specific death of unisolated infected I and hospitalized H, respectively, and b1,b2 represent the rate of conversion of unisolated infected I and isolated exposed Eq to hospitalized H, respectively.


### 2.1. Proof of Stability of the SEIR-QV Model

To study the stability of equation (2), it is sufficient to study the stability of the following three formulas in (2) [13, 14]:(3)dSdt=Λ−cβδSI+υE−ηS,dEdt=c1−qβδSI+υE−σ+ηE,dIdt=σE−b1+α1+ηI.Considering the biological significance of the system, the dynamical properties of equation (3) are discussed mainly in the closed set Ω=S,E,I∈R+3S|≥0,E≥0,I≥0, where R+3 denotes the first trigonometric limit of R3 and contains the boundary.


#### 2.1.1. Existence of Equilibrium Point

The equilibrium point of (3) satisfies(4)dSdt=0,dEdt=0,dIdt=0.WhenI=0, (3) has the point of diseased equilibrium P0=S0,E0,I0=Λ/η,0,0.WhenI≠0, from dI/dt=0, we get I=σE/α1+b1+η; from dE/dt=0, we get SI+υE=δσ+ηE/cβ1−q, from dS/dt=0, we get SI+υE=δ/cβΛ−ηS, so S=Λ/η−σ+ηE/1−qη.BringI=σE/α1+b1+η and S=Λ/η−σ+ηE/1−qη into SI+υE=δ/cβΛ−ηS to obtain the endemic equilibrium point P∗=S∗,E∗,I∗, where(5)S∗=δσ+ηα1+b1+ηcβ1−qσ+υα1+b1+η,E∗=Λ1−qσ+η−δηα1+b1+ηcβσ+υα1+b1+η,I∗=σΛ1−qσ+ηα1+b1+η−σδηcβσ+υα1+b1+η.


#### 2.1.2. Basic Regeneration Number

The basic regeneration number (R0) indicates the number of people infected by a patient during the average disease period when all are susceptible at the beginning of the disease. R0=1 can be used as a threshold to decide whether the disease is extinguished or not. When R0<1, the disease will become extinct. When R0>1, the disease will persist. The basic regeneration number is closely related to the stability of the endemic equilibrium point.Next, we study the basic regeneration numberR0 of equation (3). Let x=E,I,ST, then equation (3) can be rewritten as(6)dxdt=Fx−Vx,where(7)Fx=c1−qβδSI+υE00,Vx=σ+ηEα1+b1+ηI−σEηS+cβδSI+υE−Λ.The Jacobi matrices ofFx and Vx at the disease-free equilibrium P0=Λ/η,0,0 are, respectively,(8)DFP0=F2×2000=c1−qβυΛδηc1−qβΛδη0000000,DVP0=V2×20cβυΛδηcβΛδηη=σ+η00−σα1+b1+η0cβυΛδηcβΛδηη,where(9)F=c1−qβυΛδηc1−qβΛδη00,V=σ+η0−σα1+b1+η.The basic regeneration number, denoted byR0, is thus given by the following equation:(10)R0=ρFV−1=cΛβ1−qηδη+συ+σα1+b1+η.


#### 2.1.3. Proof of Stability of the SEIR-QV Model

Theorem 1.
For equation (3), the disease-free equilibrium P0=Λ/η,0,0 is locally asymptotically stable if R0<1.Proof.
Linearizing equation (3) at the disease-free equilibrium point P0=Λ/η,0,0, we obtain the linearization matrix at point P0=Λ/η,0,0 as the following equation:(11)JP0=−η−cβυΛδη−cβΛδη0c1−qβυΛδη−σ+ηc1−qβΛδη0σ−α1+b1+η.
The characteristic equation of this matrix isdetλL−JP0, where L is a 3×3 unit matrix. Expanding it gives(12)λ+ηλ2+α1+b1+η+σ+η−c1−qβυΛδηλ+α1+b1+ησ+η1−R0=0.
Obviously, this characteristic equation has a negative characteristic rootλ=−η. The other characteristic roots satisfy the following equation:(13)λ2+α1+b1+η+σ+η−c1−qβυΛδηλ+α1+b1+ησ+η1−R0=0.
It is known thatR0=ρFV−1=cΛβ1−q/ηδη+συ+σ/α1+b1+η≤1.
LetR1=R0−cΛβ1−q/ηδη+σ⋅σ/α1+b1+η=cΛβ1−qυ/ηδη+σ<1, then(14)α1+b1+η+σ+η−c1−qβυΛδη=α1+b1+η+σ+η1−c1−qβυΛδησ+η=α1+b1+η+σ+η1−R1.
So we can obtain(15)λ2+α1+b1+η+σ+η1−R1λ+α1+b1+ησ+η1−R0=0.
It is easy to verify forR0<1 and R1<1(16)α1+b1+η+σ+η1−R1α1+b1+ησ+η1−R0>0,that the roots of the quadratic equation are as follows:(17)λ2+α1+b1+η+σ+η1−R1λ+α1+b1+ησ+η1−R0=0.
All have negative real parts, i.e., the disease-free equilibrium pointP0=Λ/η,0,0 of equation (3) is proved to be locally asymptotically stable.Theorem 2.
For equation (3), the disease-free equilibrium P0=Λ/η,0,0 is globally asymptotically stable if R0≤1.Proof.
By considering the Lyapunov function,(18)Vt=x1Et+x2It.
Clearly, the solution ofVt along equation (3) has Vt≥0 and Vt=0 when and only when Et=0 and It=0. The full derivative of the solution of the function Vt along equation (3) is as follows:(19)dVtdt=x1dEtdt+x2dItdt=x1c1−qβδSI+υE−σ+ηE+x2σE−b1+α1+ηI=x1c1−qβδSI+x1c1−qβυδS−x1σ+η+x2σE−x2b1+α1+ηI≤c1−qβΛδηx1I−b1+α1+ηx2I+c1−qβυΛδηx1−σ+ηx1+σx2E.
Takex1=σ, x2=σ+ηc1−qβυΛ/δη.
Then,(20)dVtdt≤c1−qβΛσδηI−b1+α1+ησ+η−c1−qβυΛδηI=c1−qβΛσδη+b1+α1+ηc1−qβυΛδη−b1+α1+ησ+ηI=c1−qβΛδησ+b1+α1+ηυ−b1+α1+ησ+ηI=b1+α1+ησ+ηc1−qβΛδησ+ηυ+σb1+α1+η−b1+α1+ησ+ηI,=b1+α1+ησ+ηR0−1I,where the inequality sign is obtained based on S≤S0. When R0≤1, there is dVt/dt≤0. Thus, dVt/dt≤0 when and only when It=0. Therefore, if R0≤1, the maximum tight invariant set S,E,I∈Ω|dVt/dt=0 is the single point set P0. Therefore, if R0≤1, the disease-free equilibrium point P0 of equation (3) is globally asymptotically stable in Ω.Theorem 3.
S∗,E∗,I∗ are positive when R0>1.Proof.
(21)S∗=δσ+ηα1+b1+ηcβ1−qσ+υα1+b1+η.
It is known that the parameters are positive andq<1, so S∗>0 obviously holds the following equation:(22)E∗=Λ1−qσ+η−δηα1+b1+ηcβσ+υα1+b1+η=cβΛ1−qσ+υα1+b1+η−δηα1+b1+ησ+ηcβσ+ησ+υα1+b1+η=cβΛ1−q/σ+ησ+υα1+b1+η−δηα1+b1+ηcβσ+υα1+b1+η=cβΛ1−q/δησ+ησ+υα1+b1+η−α1+b1+ηcβ/δησ+υα1+b1+η=α1+b1+ηcβΛ1−q/δησ+ηυ+σ/α1+b1+η−1cβ/δησ+υα1+b1+η=α1+b1+ηR0−1cβ/δησ+υα1+b1+η.
It is known thatR0>1, q<1 and the parameters are positive, so E∗>0.(23)I∗=σΛ1−qσ+ηα1+b1+η−σδηcβσ+υα1+b1+η=σ⋅cβΛ1−qσ+υα1+b1+η−δησ+ηα1+b1+ηcβσ+ηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+υα1+b1+η−σ+ηα1+b1+ηcβ/δησ+ηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+ηα1+b1+ηυ+σ/α1+b1+η−α1+b1+ηcβ/δηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+ηυ+σ/α1+b1+η−1cβ/δησ+υα1+b1+η=σ⋅R0−1cβ/δησ+υα1+b1+η.
It is known thatR0>1, q<1 and the parameters are positive, so I∗>0.Theorem 4.
For equation (3), the endemic equilibrium P∗=S∗,E∗,I∗ is locally asymptotically stable if R0>1.Proof.
Linearize equation (3) at the endemic equilibrium point P∗=S∗,E∗,I∗ and obtain the linearization matrix at point P∗=S∗,E∗,I∗ as follows:(24)JP∗=−cβδI∗+υE∗−η−cβυδS∗−cβδS∗c1−qβδI∗+υE∗c1−qβυδS∗−σ+ηc1−qβδS∗0σ−α1+b1+η.
The characteristic equation of this matrix isdetλL−JP∗, where L is a 3×3 unit matrix. Expanding it gives λ3+Aλ2+Bλ+C=0, where(25)A=−c1−qβυS∗δ+σ+η+α1+b1+η+η+cβI∗+υE∗δ,B=−cη1−qβυS∗δ+ησ+η+ηα1+b1+η−cσ1−qβS∗δ+α1+b1+η−c1−qβυS∗δ+σ+η+cβα1+b1+η+σ+ηI∗+υE∗δ,C=cη+cβα1+b1+ησ+ηI∗+υE∗δ.
Then,(26)AB−C≥c2β21−q2υυη+σ+υα1+b1+ηδ2S∗2+f1+f2+ηf1f2+ηf1+f2−c2β21−q2υα1+b1+η+2υη+υσ+σδ2−cβf1f2δ−cη−cβ1−qΛυf1f2+υηf1+f2+υη+σ+υf1f1+f2+ηδη,where f1=α1+b1+η and f2=σ+η.
Due toI∗+υE∗≥I0+υE0 and S∗≤S0, AB−C>0.
Therefore, the endemic equilibrium pointP∗=S∗,E∗,I∗ of equation (3) is locally asymptotically stable on Ω0.Theorem 5.
For equation (3), the endemic equilibrium P∗=S∗,E∗,I∗ is globally asymptotically stable if R0>1.Proof.
Considering the Lyapunov function, we can obtain(27)Vt=S−S∗−S∗lnSS∗+x1E−E∗−E∗lnEE∗+x2I−I∗−I∗lnII∗.
The full derivative of the solution of the functionVt along equation (3) is as follows:(28)dVtdt=1−S∗SdSdt+x11−E∗EdEdt+x21−I∗IdIdt=1−S∗SΛ−cβδSI+υE−ηS+x21−I∗IσE−b1+α1+ηI+x11−E∗Ec1−qβδSI+υE−σ+ηE=I+υEcβS∗δ+x1c1−qβSδ−cβSδ−x1c1−qβE∗SδE+Λ−ηS−ΛS∗S+ηS∗−x1σ+ηE+x1σ+ηE∗+x2σE−x2b1+α1+ηI−x2σEI∗I+x2b1+α1+ηI∗.
Bringing inI∗=σE∗/α1+b1+η, S∗=Λ/η−σ+ηE∗/1−qη yields(29)dVtdt=I+υEcβδΛη−σ+ηE∗1−qη+x1c1−qβSδ−cβSδ−x1c1−qβE∗SδE+Λ−ηS+Λη−σ+ηE∗1−qηη−ΛS−x1σ+ηE+x1σ+ηE∗+x2σE−x2b1+α1+ηI−x2σ2EE∗Ib1+α1+η+x2b1+α1+ησE∗b1+α1+η=cβδI+υES∗+x11−q−1S−x11−qE∗SE−Λ−ηS2ηS+Λ−ηSσ+ηE∗ηS1−q+x2σ−x1σ+ηE+x2σ+x1σ+ηE∗−x2b1+α1+ηI−x2σ2EE∗Ib1+α1+η.
Takex1=2/1−q and x2=2σ+η/σ1−q.
Then,(30)dVtdt=cβδI+υES∗+S−2E∗SE−Λ−ηS2ηS+Λ−ηSσ+ηE∗ηS1−q+4σ+η1−qE∗−2σ+ηb1+α1+ηIσ1−q−2σσ+ηEE∗Ib1+α1+η1−q=cβδI+υES∗+1−2E∗ES−Λ−ηSSΛ−ηSη−σ+ηE∗η1−q+22σb1+α1+ησ+ηE∗I−σ+ηb1+α1+ηI2−σ2σ+ηEE∗σ1−qb1+α1+ηI.
SinceΛ−ηS<Λ and when E>α1+b1+η, σ2E2<σ2EE∗/α1+b1+η.
At this time,(31)dVtdt≥cβδI+υES∗+1−2E∗ES+ηS−ΛSS∗+2σ+ηI−σE2σ1−qI,and when E>2E∗ and S>Λ/η, 1−2E∗/E>0, and ηS−Λ>0 holds, so dVt/dt≥0. Therefore, if R0>1, the endemic equilibrium point P∗=S∗,E∗,I∗ of equation (3) is globally asymptotically stable on Ω0.


## 2.1. Proof of Stability of the SEIR-QV Model

To study the stability of equation (2), it is sufficient to study the stability of the following three formulas in (2) [13, 14]:(3)dSdt=Λ−cβδSI+υE−ηS,dEdt=c1−qβδSI+υE−σ+ηE,dIdt=σE−b1+α1+ηI.Considering the biological significance of the system, the dynamical properties of equation (3) are discussed mainly in the closed set Ω=S,E,I∈R+3S|≥0,E≥0,I≥0, where R+3 denotes the first trigonometric limit of R3 and contains the boundary.


### 2.1.1. Existence of Equilibrium Point

The equilibrium point of (3) satisfies(4)dSdt=0,dEdt=0,dIdt=0.WhenI=0, (3) has the point of diseased equilibrium P0=S0,E0,I0=Λ/η,0,0.WhenI≠0, from dI/dt=0, we get I=σE/α1+b1+η; from dE/dt=0, we get SI+υE=δσ+ηE/cβ1−q, from dS/dt=0, we get SI+υE=δ/cβΛ−ηS, so S=Λ/η−σ+ηE/1−qη.BringI=σE/α1+b1+η and S=Λ/η−σ+ηE/1−qη into SI+υE=δ/cβΛ−ηS to obtain the endemic equilibrium point P∗=S∗,E∗,I∗, where(5)S∗=δσ+ηα1+b1+ηcβ1−qσ+υα1+b1+η,E∗=Λ1−qσ+η−δηα1+b1+ηcβσ+υα1+b1+η,I∗=σΛ1−qσ+ηα1+b1+η−σδηcβσ+υα1+b1+η.


### 2.1.2. Basic Regeneration Number

The basic regeneration number (R0) indicates the number of people infected by a patient during the average disease period when all are susceptible at the beginning of the disease. R0=1 can be used as a threshold to decide whether the disease is extinguished or not. When R0<1, the disease will become extinct. When R0>1, the disease will persist. The basic regeneration number is closely related to the stability of the endemic equilibrium point.Next, we study the basic regeneration numberR0 of equation (3). Let x=E,I,ST, then equation (3) can be rewritten as(6)dxdt=Fx−Vx,where(7)Fx=c1−qβδSI+υE00,Vx=σ+ηEα1+b1+ηI−σEηS+cβδSI+υE−Λ.The Jacobi matrices ofFx and Vx at the disease-free equilibrium P0=Λ/η,0,0 are, respectively,(8)DFP0=F2×2000=c1−qβυΛδηc1−qβΛδη0000000,DVP0=V2×20cβυΛδηcβΛδηη=σ+η00−σα1+b1+η0cβυΛδηcβΛδηη,where(9)F=c1−qβυΛδηc1−qβΛδη00,V=σ+η0−σα1+b1+η.The basic regeneration number, denoted byR0, is thus given by the following equation:(10)R0=ρFV−1=cΛβ1−qηδη+συ+σα1+b1+η.


### 2.1.3. Proof of Stability of the SEIR-QV Model

Theorem 1.
For equation (3), the disease-free equilibrium P0=Λ/η,0,0 is locally asymptotically stable if R0<1.Proof.
Linearizing equation (3) at the disease-free equilibrium point P0=Λ/η,0,0, we obtain the linearization matrix at point P0=Λ/η,0,0 as the following equation:(11)JP0=−η−cβυΛδη−cβΛδη0c1−qβυΛδη−σ+ηc1−qβΛδη0σ−α1+b1+η.
The characteristic equation of this matrix isdetλL−JP0, where L is a 3×3 unit matrix. Expanding it gives(12)λ+ηλ2+α1+b1+η+σ+η−c1−qβυΛδηλ+α1+b1+ησ+η1−R0=0.
Obviously, this characteristic equation has a negative characteristic rootλ=−η. The other characteristic roots satisfy the following equation:(13)λ2+α1+b1+η+σ+η−c1−qβυΛδηλ+α1+b1+ησ+η1−R0=0.
It is known thatR0=ρFV−1=cΛβ1−q/ηδη+συ+σ/α1+b1+η≤1.
LetR1=R0−cΛβ1−q/ηδη+σ⋅σ/α1+b1+η=cΛβ1−qυ/ηδη+σ<1, then(14)α1+b1+η+σ+η−c1−qβυΛδη=α1+b1+η+σ+η1−c1−qβυΛδησ+η=α1+b1+η+σ+η1−R1.
So we can obtain(15)λ2+α1+b1+η+σ+η1−R1λ+α1+b1+ησ+η1−R0=0.
It is easy to verify forR0<1 and R1<1(16)α1+b1+η+σ+η1−R1α1+b1+ησ+η1−R0>0,that the roots of the quadratic equation are as follows:(17)λ2+α1+b1+η+σ+η1−R1λ+α1+b1+ησ+η1−R0=0.
All have negative real parts, i.e., the disease-free equilibrium pointP0=Λ/η,0,0 of equation (3) is proved to be locally asymptotically stable.Theorem 2.
For equation (3), the disease-free equilibrium P0=Λ/η,0,0 is globally asymptotically stable if R0≤1.Proof.
By considering the Lyapunov function,(18)Vt=x1Et+x2It.
Clearly, the solution ofVt along equation (3) has Vt≥0 and Vt=0 when and only when Et=0 and It=0. The full derivative of the solution of the function Vt along equation (3) is as follows:(19)dVtdt=x1dEtdt+x2dItdt=x1c1−qβδSI+υE−σ+ηE+x2σE−b1+α1+ηI=x1c1−qβδSI+x1c1−qβυδS−x1σ+η+x2σE−x2b1+α1+ηI≤c1−qβΛδηx1I−b1+α1+ηx2I+c1−qβυΛδηx1−σ+ηx1+σx2E.
Takex1=σ, x2=σ+ηc1−qβυΛ/δη.
Then,(20)dVtdt≤c1−qβΛσδηI−b1+α1+ησ+η−c1−qβυΛδηI=c1−qβΛσδη+b1+α1+ηc1−qβυΛδη−b1+α1+ησ+ηI=c1−qβΛδησ+b1+α1+ηυ−b1+α1+ησ+ηI=b1+α1+ησ+ηc1−qβΛδησ+ηυ+σb1+α1+η−b1+α1+ησ+ηI,=b1+α1+ησ+ηR0−1I,where the inequality sign is obtained based on S≤S0. When R0≤1, there is dVt/dt≤0. Thus, dVt/dt≤0 when and only when It=0. Therefore, if R0≤1, the maximum tight invariant set S,E,I∈Ω|dVt/dt=0 is the single point set P0. Therefore, if R0≤1, the disease-free equilibrium point P0 of equation (3) is globally asymptotically stable in Ω.Theorem 3.
S∗,E∗,I∗ are positive when R0>1.Proof.
(21)S∗=δσ+ηα1+b1+ηcβ1−qσ+υα1+b1+η.
It is known that the parameters are positive andq<1, so S∗>0 obviously holds the following equation:(22)E∗=Λ1−qσ+η−δηα1+b1+ηcβσ+υα1+b1+η=cβΛ1−qσ+υα1+b1+η−δηα1+b1+ησ+ηcβσ+ησ+υα1+b1+η=cβΛ1−q/σ+ησ+υα1+b1+η−δηα1+b1+ηcβσ+υα1+b1+η=cβΛ1−q/δησ+ησ+υα1+b1+η−α1+b1+ηcβ/δησ+υα1+b1+η=α1+b1+ηcβΛ1−q/δησ+ηυ+σ/α1+b1+η−1cβ/δησ+υα1+b1+η=α1+b1+ηR0−1cβ/δησ+υα1+b1+η.
It is known thatR0>1, q<1 and the parameters are positive, so E∗>0.(23)I∗=σΛ1−qσ+ηα1+b1+η−σδηcβσ+υα1+b1+η=σ⋅cβΛ1−qσ+υα1+b1+η−δησ+ηα1+b1+ηcβσ+ηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+υα1+b1+η−σ+ηα1+b1+ηcβ/δησ+ηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+ηα1+b1+ηυ+σ/α1+b1+η−α1+b1+ηcβ/δηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+ηυ+σ/α1+b1+η−1cβ/δησ+υα1+b1+η=σ⋅R0−1cβ/δησ+υα1+b1+η.
It is known thatR0>1, q<1 and the parameters are positive, so I∗>0.Theorem 4.
For equation (3), the endemic equilibrium P∗=S∗,E∗,I∗ is locally asymptotically stable if R0>1.Proof.
Linearize equation (3) at the endemic equilibrium point P∗=S∗,E∗,I∗ and obtain the linearization matrix at point P∗=S∗,E∗,I∗ as follows:(24)JP∗=−cβδI∗+υE∗−η−cβυδS∗−cβδS∗c1−qβδI∗+υE∗c1−qβυδS∗−σ+ηc1−qβδS∗0σ−α1+b1+η.
The characteristic equation of this matrix isdetλL−JP∗, where L is a 3×3 unit matrix. Expanding it gives λ3+Aλ2+Bλ+C=0, where(25)A=−c1−qβυS∗δ+σ+η+α1+b1+η+η+cβI∗+υE∗δ,B=−cη1−qβυS∗δ+ησ+η+ηα1+b1+η−cσ1−qβS∗δ+α1+b1+η−c1−qβυS∗δ+σ+η+cβα1+b1+η+σ+ηI∗+υE∗δ,C=cη+cβα1+b1+ησ+ηI∗+υE∗δ.
Then,(26)AB−C≥c2β21−q2υυη+σ+υα1+b1+ηδ2S∗2+f1+f2+ηf1f2+ηf1+f2−c2β21−q2υα1+b1+η+2υη+υσ+σδ2−cβf1f2δ−cη−cβ1−qΛυf1f2+υηf1+f2+υη+σ+υf1f1+f2+ηδη,where f1=α1+b1+η and f2=σ+η.
Due toI∗+υE∗≥I0+υE0 and S∗≤S0, AB−C>0.
Therefore, the endemic equilibrium pointP∗=S∗,E∗,I∗ of equation (3) is locally asymptotically stable on Ω0.Theorem 5.
For equation (3), the endemic equilibrium P∗=S∗,E∗,I∗ is globally asymptotically stable if R0>1.Proof.
Considering the Lyapunov function, we can obtain(27)Vt=S−S∗−S∗lnSS∗+x1E−E∗−E∗lnEE∗+x2I−I∗−I∗lnII∗.
The full derivative of the solution of the functionVt along equation (3) is as follows:(28)dVtdt=1−S∗SdSdt+x11−E∗EdEdt+x21−I∗IdIdt=1−S∗SΛ−cβδSI+υE−ηS+x21−I∗IσE−b1+α1+ηI+x11−E∗Ec1−qβδSI+υE−σ+ηE=I+υEcβS∗δ+x1c1−qβSδ−cβSδ−x1c1−qβE∗SδE+Λ−ηS−ΛS∗S+ηS∗−x1σ+ηE+x1σ+ηE∗+x2σE−x2b1+α1+ηI−x2σEI∗I+x2b1+α1+ηI∗.
Bringing inI∗=σE∗/α1+b1+η, S∗=Λ/η−σ+ηE∗/1−qη yields(29)dVtdt=I+υEcβδΛη−σ+ηE∗1−qη+x1c1−qβSδ−cβSδ−x1c1−qβE∗SδE+Λ−ηS+Λη−σ+ηE∗1−qηη−ΛS−x1σ+ηE+x1σ+ηE∗+x2σE−x2b1+α1+ηI−x2σ2EE∗Ib1+α1+η+x2b1+α1+ησE∗b1+α1+η=cβδI+υES∗+x11−q−1S−x11−qE∗SE−Λ−ηS2ηS+Λ−ηSσ+ηE∗ηS1−q+x2σ−x1σ+ηE+x2σ+x1σ+ηE∗−x2b1+α1+ηI−x2σ2EE∗Ib1+α1+η.
Takex1=2/1−q and x2=2σ+η/σ1−q.
Then,(30)dVtdt=cβδI+υES∗+S−2E∗SE−Λ−ηS2ηS+Λ−ηSσ+ηE∗ηS1−q+4σ+η1−qE∗−2σ+ηb1+α1+ηIσ1−q−2σσ+ηEE∗Ib1+α1+η1−q=cβδI+υES∗+1−2E∗ES−Λ−ηSSΛ−ηSη−σ+ηE∗η1−q+22σb1+α1+ησ+ηE∗I−σ+ηb1+α1+ηI2−σ2σ+ηEE∗σ1−qb1+α1+ηI.
SinceΛ−ηS<Λ and when E>α1+b1+η, σ2E2<σ2EE∗/α1+b1+η.
At this time,(31)dVtdt≥cβδI+υES∗+1−2E∗ES+ηS−ΛSS∗+2σ+ηI−σE2σ1−qI,and when E>2E∗ and S>Λ/η, 1−2E∗/E>0, and ηS−Λ>0 holds, so dVt/dt≥0. Therefore, if R0>1, the endemic equilibrium point P∗=S∗,E∗,I∗ of equation (3) is globally asymptotically stable on Ω0.


## 2.1.1. Existence of Equilibrium Point

The equilibrium point of (3) satisfies(4)dSdt=0,dEdt=0,dIdt=0.WhenI=0, (3) has the point of diseased equilibrium P0=S0,E0,I0=Λ/η,0,0.WhenI≠0, from dI/dt=0, we get I=σE/α1+b1+η; from dE/dt=0, we get SI+υE=δσ+ηE/cβ1−q, from dS/dt=0, we get SI+υE=δ/cβΛ−ηS, so S=Λ/η−σ+ηE/1−qη.BringI=σE/α1+b1+η and S=Λ/η−σ+ηE/1−qη into SI+υE=δ/cβΛ−ηS to obtain the endemic equilibrium point P∗=S∗,E∗,I∗, where(5)S∗=δσ+ηα1+b1+ηcβ1−qσ+υα1+b1+η,E∗=Λ1−qσ+η−δηα1+b1+ηcβσ+υα1+b1+η,I∗=σΛ1−qσ+ηα1+b1+η−σδηcβσ+υα1+b1+η.


## 2.1.2. Basic Regeneration Number

The basic regeneration number (R0) indicates the number of people infected by a patient during the average disease period when all are susceptible at the beginning of the disease. R0=1 can be used as a threshold to decide whether the disease is extinguished or not. When R0<1, the disease will become extinct. When R0>1, the disease will persist. The basic regeneration number is closely related to the stability of the endemic equilibrium point.Next, we study the basic regeneration numberR0 of equation (3). Let x=E,I,ST, then equation (3) can be rewritten as(6)dxdt=Fx−Vx,where(7)Fx=c1−qβδSI+υE00,Vx=σ+ηEα1+b1+ηI−σEηS+cβδSI+υE−Λ.The Jacobi matrices ofFx and Vx at the disease-free equilibrium P0=Λ/η,0,0 are, respectively,(8)DFP0=F2×2000=c1−qβυΛδηc1−qβΛδη0000000,DVP0=V2×20cβυΛδηcβΛδηη=σ+η00−σα1+b1+η0cβυΛδηcβΛδηη,where(9)F=c1−qβυΛδηc1−qβΛδη00,V=σ+η0−σα1+b1+η.The basic regeneration number, denoted byR0, is thus given by the following equation:(10)R0=ρFV−1=cΛβ1−qηδη+συ+σα1+b1+η.


## 2.1.3. Proof of Stability of the SEIR-QV Model

Theorem 1.
For equation (3), the disease-free equilibrium P0=Λ/η,0,0 is locally asymptotically stable if R0<1.Proof.
Linearizing equation (3) at the disease-free equilibrium point P0=Λ/η,0,0, we obtain the linearization matrix at point P0=Λ/η,0,0 as the following equation:(11)JP0=−η−cβυΛδη−cβΛδη0c1−qβυΛδη−σ+ηc1−qβΛδη0σ−α1+b1+η.
The characteristic equation of this matrix isdetλL−JP0, where L is a 3×3 unit matrix. Expanding it gives(12)λ+ηλ2+α1+b1+η+σ+η−c1−qβυΛδηλ+α1+b1+ησ+η1−R0=0.
Obviously, this characteristic equation has a negative characteristic rootλ=−η. The other characteristic roots satisfy the following equation:(13)λ2+α1+b1+η+σ+η−c1−qβυΛδηλ+α1+b1+ησ+η1−R0=0.
It is known thatR0=ρFV−1=cΛβ1−q/ηδη+συ+σ/α1+b1+η≤1.
LetR1=R0−cΛβ1−q/ηδη+σ⋅σ/α1+b1+η=cΛβ1−qυ/ηδη+σ<1, then(14)α1+b1+η+σ+η−c1−qβυΛδη=α1+b1+η+σ+η1−c1−qβυΛδησ+η=α1+b1+η+σ+η1−R1.
So we can obtain(15)λ2+α1+b1+η+σ+η1−R1λ+α1+b1+ησ+η1−R0=0.
It is easy to verify forR0<1 and R1<1(16)α1+b1+η+σ+η1−R1α1+b1+ησ+η1−R0>0,that the roots of the quadratic equation are as follows:(17)λ2+α1+b1+η+σ+η1−R1λ+α1+b1+ησ+η1−R0=0.
All have negative real parts, i.e., the disease-free equilibrium pointP0=Λ/η,0,0 of equation (3) is proved to be locally asymptotically stable.Theorem 2.
For equation (3), the disease-free equilibrium P0=Λ/η,0,0 is globally asymptotically stable if R0≤1.Proof.
By considering the Lyapunov function,(18)Vt=x1Et+x2It.
Clearly, the solution ofVt along equation (3) has Vt≥0 and Vt=0 when and only when Et=0 and It=0. The full derivative of the solution of the function Vt along equation (3) is as follows:(19)dVtdt=x1dEtdt+x2dItdt=x1c1−qβδSI+υE−σ+ηE+x2σE−b1+α1+ηI=x1c1−qβδSI+x1c1−qβυδS−x1σ+η+x2σE−x2b1+α1+ηI≤c1−qβΛδηx1I−b1+α1+ηx2I+c1−qβυΛδηx1−σ+ηx1+σx2E.
Takex1=σ, x2=σ+ηc1−qβυΛ/δη.
Then,(20)dVtdt≤c1−qβΛσδηI−b1+α1+ησ+η−c1−qβυΛδηI=c1−qβΛσδη+b1+α1+ηc1−qβυΛδη−b1+α1+ησ+ηI=c1−qβΛδησ+b1+α1+ηυ−b1+α1+ησ+ηI=b1+α1+ησ+ηc1−qβΛδησ+ηυ+σb1+α1+η−b1+α1+ησ+ηI,=b1+α1+ησ+ηR0−1I,where the inequality sign is obtained based on S≤S0. When R0≤1, there is dVt/dt≤0. Thus, dVt/dt≤0 when and only when It=0. Therefore, if R0≤1, the maximum tight invariant set S,E,I∈Ω|dVt/dt=0 is the single point set P0. Therefore, if R0≤1, the disease-free equilibrium point P0 of equation (3) is globally asymptotically stable in Ω.Theorem 3.
S∗,E∗,I∗ are positive when R0>1.Proof.
(21)S∗=δσ+ηα1+b1+ηcβ1−qσ+υα1+b1+η.
It is known that the parameters are positive andq<1, so S∗>0 obviously holds the following equation:(22)E∗=Λ1−qσ+η−δηα1+b1+ηcβσ+υα1+b1+η=cβΛ1−qσ+υα1+b1+η−δηα1+b1+ησ+ηcβσ+ησ+υα1+b1+η=cβΛ1−q/σ+ησ+υα1+b1+η−δηα1+b1+ηcβσ+υα1+b1+η=cβΛ1−q/δησ+ησ+υα1+b1+η−α1+b1+ηcβ/δησ+υα1+b1+η=α1+b1+ηcβΛ1−q/δησ+ηυ+σ/α1+b1+η−1cβ/δησ+υα1+b1+η=α1+b1+ηR0−1cβ/δησ+υα1+b1+η.
It is known thatR0>1, q<1 and the parameters are positive, so E∗>0.(23)I∗=σΛ1−qσ+ηα1+b1+η−σδηcβσ+υα1+b1+η=σ⋅cβΛ1−qσ+υα1+b1+η−δησ+ηα1+b1+ηcβσ+ηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+υα1+b1+η−σ+ηα1+b1+ηcβ/δησ+ηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+ηα1+b1+ηυ+σ/α1+b1+η−α1+b1+ηcβ/δηα1+b1+ησ+υα1+b1+η=σ⋅cβΛ1−q/δησ+ηυ+σ/α1+b1+η−1cβ/δησ+υα1+b1+η=σ⋅R0−1cβ/δησ+υα1+b1+η.
It is known thatR0>1, q<1 and the parameters are positive, so I∗>0.Theorem 4.
For equation (3), the endemic equilibrium P∗=S∗,E∗,I∗ is locally asymptotically stable if R0>1.Proof.
Linearize equation (3) at the endemic equilibrium point P∗=S∗,E∗,I∗ and obtain the linearization matrix at point P∗=S∗,E∗,I∗ as follows:(24)JP∗=−cβδI∗+υE∗−η−cβυδS∗−cβδS∗c1−qβδI∗+υE∗c1−qβυδS∗−σ+ηc1−qβδS∗0σ−α1+b1+η.
The characteristic equation of this matrix isdetλL−JP∗, where L is a 3×3 unit matrix. Expanding it gives λ3+Aλ2+Bλ+C=0, where(25)A=−c1−qβυS∗δ+σ+η+α1+b1+η+η+cβI∗+υE∗δ,B=−cη1−qβυS∗δ+ησ+η+ηα1+b1+η−cσ1−qβS∗δ+α1+b1+η−c1−qβυS∗δ+σ+η+cβα1+b1+η+σ+ηI∗+υE∗δ,C=cη+cβα1+b1+ησ+ηI∗+υE∗δ.
Then,(26)AB−C≥c2β21−q2υυη+σ+υα1+b1+ηδ2S∗2+f1+f2+ηf1f2+ηf1+f2−c2β21−q2υα1+b1+η+2υη+υσ+σδ2−cβf1f2δ−cη−cβ1−qΛυf1f2+υηf1+f2+υη+σ+υf1f1+f2+ηδη,where f1=α1+b1+η and f2=σ+η.
Due toI∗+υE∗≥I0+υE0 and S∗≤S0, AB−C>0.
Therefore, the endemic equilibrium pointP∗=S∗,E∗,I∗ of equation (3) is locally asymptotically stable on Ω0.Theorem 5.
For equation (3), the endemic equilibrium P∗=S∗,E∗,I∗ is globally asymptotically stable if R0>1.Proof.
Considering the Lyapunov function, we can obtain(27)Vt=S−S∗−S∗lnSS∗+x1E−E∗−E∗lnEE∗+x2I−I∗−I∗lnII∗.
The full derivative of the solution of the functionVt along equation (3) is as follows:(28)dVtdt=1−S∗SdSdt+x11−E∗EdEdt+x21−I∗IdIdt=1−S∗SΛ−cβδSI+υE−ηS+x21−I∗IσE−b1+α1+ηI+x11−E∗Ec1−qβδSI+υE−σ+ηE=I+υEcβS∗δ+x1c1−qβSδ−cβSδ−x1c1−qβE∗SδE+Λ−ηS−ΛS∗S+ηS∗−x1σ+ηE+x1σ+ηE∗+x2σE−x2b1+α1+ηI−x2σEI∗I+x2b1+α1+ηI∗.
Bringing inI∗=σE∗/α1+b1+η, S∗=Λ/η−σ+ηE∗/1−qη yields(29)dVtdt=I+υEcβδΛη−σ+ηE∗1−qη+x1c1−qβSδ−cβSδ−x1c1−qβE∗SδE+Λ−ηS+Λη−σ+ηE∗1−qηη−ΛS−x1σ+ηE+x1σ+ηE∗+x2σE−x2b1+α1+ηI−x2σ2EE∗Ib1+α1+η+x2b1+α1+ησE∗b1+α1+η=cβδI+υES∗+x11−q−1S−x11−qE∗SE−Λ−ηS2ηS+Λ−ηSσ+ηE∗ηS1−q+x2σ−x1σ+ηE+x2σ+x1σ+ηE∗−x2b1+α1+ηI−x2σ2EE∗Ib1+α1+η.
Takex1=2/1−q and x2=2σ+η/σ1−q.
Then,(30)dVtdt=cβδI+υES∗+S−2E∗SE−Λ−ηS2ηS+Λ−ηSσ+ηE∗ηS1−q+4σ+η1−qE∗−2σ+ηb1+α1+ηIσ1−q−2σσ+ηEE∗Ib1+α1+η1−q=cβδI+υES∗+1−2E∗ES−Λ−ηSSΛ−ηSη−σ+ηE∗η1−q+22σb1+α1+ησ+ηE∗I−σ+ηb1+α1+ηI2−σ2σ+ηEE∗σ1−qb1+α1+ηI.
SinceΛ−ηS<Λ and when E>α1+b1+η, σ2E2<σ2EE∗/α1+b1+η.
At this time,(31)dVtdt≥cβδI+υES∗+1−2E∗ES+ηS−ΛSS∗+2σ+ηI−σE2σ1−qI,and when E>2E∗ and S>Λ/η, 1−2E∗/E>0, and ηS−Λ>0 holds, so dVt/dt≥0. Therefore, if R0>1, the endemic equilibrium point P∗=S∗,E∗,I∗ of equation (3) is globally asymptotically stable on Ω0.


## 3. Numerical Experiments

### 3.1. Numerical Simulation

The stability of disease-free and endemic equilibrium points is discussed in this section, and this section uses the software Matlab to numerically simulate equation (3) to verify the above conclusions.The initial values of the system are (0.7,0.2,0.1), and the parameters are assigned asc=2, η=0.00714, β=0.1, υ=1, b1=0.8, σ=1/7, α1=0.04, and Λ=0.00334, where q and δ are two variables. In this paper, we will take some random values to simulate the change of S, E, and I ratio when R0<1 and R0>1, respectively. When q=0.9, δ takes a random value between 0 and 1 so that R0<1 holds. When δ=1, q takes a random value between 0 and 1 so that R0<1 holds. Its value of R0 is less than 1 (R0<1) when q and delta take the values in Table 1. Its value of R0 is greater than 1 (R0<1) when q and delta take the values in Table 2.Table 1
Values of disease-free equilibrium pointsq and δ.
qδR00.90.10.72890.90.50.14580.910.07290.110.65600.510.36450.910.5729Table 2
Values of endemic equilibrium pointsq and δ.
qδR00.10.0513.120.10.16.56020.10.51.31200.10.0513.120.50.057.28910.90.051.4578The global asymptotic stability of both disease-free and endemic equilibrium points can be seen from Theorems1–5. Figures 3–6 illustrates the correctness of the results obtained from the above theorem.Figure 3
Stability of disease-free equilibrium point (q=0.9).
(a)(b)(c)Figure 4
Stability of disease-free equilibrium point (δ=1).
(a)(b)(c)Figure 5
Stability of disease-free equilibrium point (q=0.1).
(a)(b)(c)Figure 6
Stability of disease-free equilibrium point (δ=0.05).
(a)(b)(c)IfR0<1, by numerically simulating the stability of the disease-free equilibrium point, we can analytically conclude the following points. Figure 3 shows the variation curves of S, E, and I as the parameter δ varies when q is constant. The variation of S increases, and the rate of decrease and the rate of increase both increase with the increase of δ. The rate of decrease of E and I decreases with the increase of δ. Finally, S, E, and I all stabilized (see Figures 3(a)–3(c)). The change of S decreases, and the rate of decrease of E decreases with the increase of q. But the change of I has no obvious pattern with the increase of q. Finally, S, E, and I all tend to be stable (see Figures 4(a)–4(c)). Figures 3 and 4 illustrate that when R0<1, the changes of q and δ have significant effects on the trends of S, E, and I. And when R0<1, the disease-free equilibrium point is tending to be stable.In Figures5 and 6, the curve of S, E, and I variation with q and δ is simulated. By comparing the trend of S, E, and I with q in Figure 5 and the trend of S, E, and I with δ in Figure 6, we can draw the following conclusions. Increases in both q and δ lead to larger fluctuations in S and delayed stabilization times. It will lead to an increase in the value of S as it stabilizes. What is more, as δ increases, the rate of increase and decrease of E decreases. The monotonically increasing interval of I decreases and the rate of decrease of the monotonically decreasing interval decreases with the increase of δ. The maximum peak value of both E and I decreases with the increase of δ. And the time of reaching the maximum peak value of E is delayed and the time of reaching the maximum peak value of I is advanced. Furthermore, as q increases, the rate of increase and decrease of E decreases and E reaches its maximum peak earlier and at a lower value. As q changes, there is no obvious pattern of change in I, but it can be concluded that q has a significant effect on the change in the value of I. Figures 5 and 6 illustrate that when R0>1, the changes of q and δ have significant effects on the trends of S, E, and I. And when R0>1, the endemic equilibrium point is tending to be stable.


### 3.2. Model Simulation

The real data of China in this paper are from the official website of the WHO. We preprocess the data, remove some bad points, and make the time continuous so that the data can be better adapted to the SEIR-QV model. The Chinese COVID-19 outbreak is fitted based on the SEIR-QV model. The parameters used in the model are listed as follows:c=0.219, η=0.00714, q=0.5, δ=0.2, β=1×10−8, υ=0.275, b1=0.05, σ=1/7, α1=0.035, Λ=0.00334, b2=0.05, α2=0.035, and r=0.9493. The simulation results and the actual data are shown in Figure 7.Figure 7
Model of the amount of confirmed cases in China per day.The Chinese government has decisively taken strict preventive and control measures following the emergence of the COVID-19 outbreak. These prevention and control measures, along with vigorous advocacy by the Chinese government, have reduced exposure ratesc and increased isolation rates q and minimized the development of aggregated outbreaks. The COVID-19 outbreak has become manageable and stable. As can be seen in Figure 7, the simulated data fit better with the real data from January 20, 2020 to May 20, 2021. However, there is a small outbreak of the epidemic in China in June 2021. The main reason is that the epidemic is out of control in many places abroad, and people at home are taking it lightly and weakening the efforts of prevention and control, thus giving the opportunity for the epidemic to spread.As can be seen in Figure8, both the isolation ratio and vaccine coverage play an important role in the control of the outbreak. The highest number of confirmed cases per day in China decreases as both q and δ increase.Figure 8
Impact of control strategies on the maximum number of daily confirmed cases of the epidemic in China: (a)δ=0.2 and (b) q=0.5.
(a)(b)The correlation coefficient is used to measure the accuracy of the fit, and the value ofR is found to be 0.9394. Therefore, the results show that the SEIR-QV model can be used for COVID-19 epidemic development status assessment and has important implications.


## 3.1. Numerical Simulation

The stability of disease-free and endemic equilibrium points is discussed in this section, and this section uses the software Matlab to numerically simulate equation (3) to verify the above conclusions.The initial values of the system are (0.7,0.2,0.1), and the parameters are assigned asc=2, η=0.00714, β=0.1, υ=1, b1=0.8, σ=1/7, α1=0.04, and Λ=0.00334, where q and δ are two variables. In this paper, we will take some random values to simulate the change of S, E, and I ratio when R0<1 and R0>1, respectively. When q=0.9, δ takes a random value between 0 and 1 so that R0<1 holds. When δ=1, q takes a random value between 0 and 1 so that R0<1 holds. Its value of R0 is less than 1 (R0<1) when q and delta take the values in Table 1. Its value of R0 is greater than 1 (R0<1) when q and delta take the values in Table 2.Table 1
Values of disease-free equilibrium pointsq and δ.
qδR00.90.10.72890.90.50.14580.910.07290.110.65600.510.36450.910.5729Table 2
Values of endemic equilibrium pointsq and δ.
qδR00.10.0513.120.10.16.56020.10.51.31200.10.0513.120.50.057.28910.90.051.4578The global asymptotic stability of both disease-free and endemic equilibrium points can be seen from Theorems1–5. Figures 3–6 illustrates the correctness of the results obtained from the above theorem.Figure 3
Stability of disease-free equilibrium point (q=0.9).
(a)(b)(c)Figure 4
Stability of disease-free equilibrium point (δ=1).
(a)(b)(c)Figure 5
Stability of disease-free equilibrium point (q=0.1).
(a)(b)(c)Figure 6
Stability of disease-free equilibrium point (δ=0.05).
(a)(b)(c)IfR0<1, by numerically simulating the stability of the disease-free equilibrium point, we can analytically conclude the following points. Figure 3 shows the variation curves of S, E, and I as the parameter δ varies when q is constant. The variation of S increases, and the rate of decrease and the rate of increase both increase with the increase of δ. The rate of decrease of E and I decreases with the increase of δ. Finally, S, E, and I all stabilized (see Figures 3(a)–3(c)). The change of S decreases, and the rate of decrease of E decreases with the increase of q. But the change of I has no obvious pattern with the increase of q. Finally, S, E, and I all tend to be stable (see Figures 4(a)–4(c)). Figures 3 and 4 illustrate that when R0<1, the changes of q and δ have significant effects on the trends of S, E, and I. And when R0<1, the disease-free equilibrium point is tending to be stable.In Figures5 and 6, the curve of S, E, and I variation with q and δ is simulated. By comparing the trend of S, E, and I with q in Figure 5 and the trend of S, E, and I with δ in Figure 6, we can draw the following conclusions. Increases in both q and δ lead to larger fluctuations in S and delayed stabilization times. It will lead to an increase in the value of S as it stabilizes. What is more, as δ increases, the rate of increase and decrease of E decreases. The monotonically increasing interval of I decreases and the rate of decrease of the monotonically decreasing interval decreases with the increase of δ. The maximum peak value of both E and I decreases with the increase of δ. And the time of reaching the maximum peak value of E is delayed and the time of reaching the maximum peak value of I is advanced. Furthermore, as q increases, the rate of increase and decrease of E decreases and E reaches its maximum peak earlier and at a lower value. As q changes, there is no obvious pattern of change in I, but it can be concluded that q has a significant effect on the change in the value of I. Figures 5 and 6 illustrate that when R0>1, the changes of q and δ have significant effects on the trends of S, E, and I. And when R0>1, the endemic equilibrium point is tending to be stable.


## 3.2. Model Simulation

The real data of China in this paper are from the official website of the WHO. We preprocess the data, remove some bad points, and make the time continuous so that the data can be better adapted to the SEIR-QV model. The Chinese COVID-19 outbreak is fitted based on the SEIR-QV model. The parameters used in the model are listed as follows:c=0.219, η=0.00714, q=0.5, δ=0.2, β=1×10−8, υ=0.275, b1=0.05, σ=1/7, α1=0.035, Λ=0.00334, b2=0.05, α2=0.035, and r=0.9493. The simulation results and the actual data are shown in Figure 7.Figure 7
Model of the amount of confirmed cases in China per day.The Chinese government has decisively taken strict preventive and control measures following the emergence of the COVID-19 outbreak. These prevention and control measures, along with vigorous advocacy by the Chinese government, have reduced exposure ratesc and increased isolation rates q and minimized the development of aggregated outbreaks. The COVID-19 outbreak has become manageable and stable. As can be seen in Figure 7, the simulated data fit better with the real data from January 20, 2020 to May 20, 2021. However, there is a small outbreak of the epidemic in China in June 2021. The main reason is that the epidemic is out of control in many places abroad, and people at home are taking it lightly and weakening the efforts of prevention and control, thus giving the opportunity for the epidemic to spread.As can be seen in Figure8, both the isolation ratio and vaccine coverage play an important role in the control of the outbreak. The highest number of confirmed cases per day in China decreases as both q and δ increase.Figure 8
Impact of control strategies on the maximum number of daily confirmed cases of the epidemic in China: (a)δ=0.2 and (b) q=0.5.
(a)(b)The correlation coefficient is used to measure the accuracy of the fit, and the value ofR is found to be 0.9394. Therefore, the results show that the SEIR-QV model can be used for COVID-19 epidemic development status assessment and has important implications.


## 4. Conclusion

An SEIR-QV model based on vaccination, isolation strategies, and the impact of different parameters on the development of infectious diseases is developed. The equilibrium point and stability of the new model are proved by using the basic regeneration number and the Lyapunov function theory. Simulation experiments show that the new method has certain theoretical value for analyzing and predicting the development of the COVID-19 epidemic. We can get the conclusion that there is a significant impact on the development of infectious diseases by different vaccination and isolation strategies.



---

*Source: 3309420-2022-06-02.xml*