Control of a multigroup COVID-19 model with immunity: treatment and test elimination.

This paper introduces a multigroup COVID-19 model with immunity, in which the total population of each group is partitioned into five compartments, that is, susceptible, exposed, infective, infective in treatment and recovered compartment. If the basic reproduction number is less than or equal to one, and the infection graph is strongly connected, then the disease-free equilibrium is globally asymptotically stable and the disease dies out. However, the COVID-19 is already in a pandemic state, and the basic reproduction number is large than one. Hence, in order to make the COVID-19 die out in some groups in an area, we design some appropriate control strategies which reduce the number of exposed people and increase the number of people treated. These two methods have been proved to be the most effective methods at present. An effective algorithm is proposed to identify the groups that need to be controlled. Finally, we use the actual limited data of Hubei, Guangdong and Zhejiang provinces in China to illustrate the effectiveness of the obtained results. © Springer Nature B.V. 2020.

Introduction

In December 31, 2019, the Health Commission of Hubei province, China announced that a new coronavirus disease was found in Wuhan, the capital of Hubei province [1]. This coronavirus disease quickly spread around the world, and was tentatively named by World Health Organization (WHO) as the 2019 novel coronavirus (COVID-19) [2, 3]. As of July 10, 2020, over 12.27 million people have infected worldwide. Hence, the COVID-19 is now in a state of pandemic worldwide. However, in some local areas, this virus has been completely controlled, such as Hubei, Guangdong and Zhejiang provinces in China. This shows that through strict nucleic acid detection and treatment measures, this virus can be completely controlled in local areas.

The aim of this paper is to introduce a new multigroup COVID-19 model to simulate the outbreak in different groups (country, city or community), and design some appropriate control strategies such that the COVID-19 dies out. To this end, we will first recall the references about multigroup epidemic model. Multigroup epidemic model has been attracting much interest from researchers in various fields because the transmission of many infectious diseases can be modeled in this form, such as mumps, gonorrhea, measles, HIV/AIDS. In reality, the group can describe country, city or community. The study of multigroup epidemic model can be traced back to the reference [4], in which a class of SIS multigroup model has been presented. Various multigroup epidemic models have been investigated since then, such as [5-10]. For each multigroup epidemic model, the first task is to study the stability of disease-free equilibrium and endemic equilibrium. If the disease-free equilibrium is globally asymptotically stable, then it means that the disease will eventually die out. However, if the endemic equilibrium is globally asymptotically stable, then it means that the disease will eventually become a pandemic. When the basic reproduction number is less than or equal to one, the disease-free equilibrium is globally asymptotically stable [8, 9]. The stability of endemic equilibrium is more complex, and various methods are presented for dealing with it, such as Lyapunov function [4] and graph-theoretic approach [8]. The above-mentioned multigroup epidemic models generally do not take into account the effects of immunity. For an infectious disease, if there exists an effective vaccine, then this infectious disease can be effectively controlled by expanding the proportion of vaccinated people [11-15]. The treatment and vaccination are two effective ways to control the spread of virus [16, 17]. For multigroup epidemic model, there were some references which studied the impact of vaccination, such as multigroup SVEIR epidemic model [18], and multigroup SVIR model with age of vaccination [18, 19].

Recently, there have been some studies on multigroup COVID-19 model [20-22]. In [20], the authors introduced a multi-group SEIRA model for the spread of COVID-19, and gave a numerical simulation for this model. In order to sustain the point that social distancing is not enough to control the outbreak of COVID-19, a simple multi-group SEIR model was presented in [21]. A multi-group SIR model was presented in [22], including some parameters for presenting different hospitalization and fatality rates among different age groups. Nevertheless, the above different multigroup COVID-19 models do not consider the effective control strategies for the spread of COVID-19, such as strict nucleic acid detection and treatment measures.

In this paper, we propose a new multigroup COVID-19 model with immunity. For each group, the total population divides into five compartments, that is, a susceptible compartment, an exposed compartment, an infective compartment, an infective compartment in treatment and a recovered compartment. Assume that part of the population entering a group already has antibodies. Our model can be seen as a generalization of the model in [12, 13], and consider the cross infection among different groups. Although vaccination is a very useful method to control the infectious diseases, there are some epidemic viruses that have no vaccine, such as SARS and COVID-19 [1–3, 23]. Moreover, it is unrealistic for people to get group immunity naturally. Hence, for an infectious disease without a vaccine, in order to effectively control the spread of this disease, large-scale testing and treatment may be effective methods. For example, for COVID-19 the currently effective control methods include expanding the scale of nucleic acid detection and treatment. Hence, for our multigroup COVID-19 model we consider the control problem and use two strategies, that is, expanding the scale of nucleic acid detection and treatment.

For the proposed multigroup COVID-19 model we study the stability of disease-free equilibrium, and obtain that the disease-free equilibrium is globally asymptotically stable when the basic reproduction number is less than or equal to one, and the infection graph is strongly connected. When the disease-free equilibrium is unstable, two control strategies are provided to make the disease-free equilibrium globally asymptotically stable, namely linear state feedback and finite-time controllers. Finally, we use the multigroup COVID-19 model to simulate the COVID-19 outbreak in Hubei, Guangdong and Zhejiang provinces in China.

The paper is organized as follows: Sect. 2 presents some mathematical notations, some facts about digraph, the problem statement, and two lemmas. Section 3.1 considers the stability of disease-free equilibrium. Section 3.2 studies the stability of endemic equilibrium. Section 3.3 presents some control strategies for the proposed multigroup COVID-19 model. Section 4 provides the simulation of the COVID-19 outbreak in Hubei, Guangdong and Zhejiang provinces in China. Section 5 collects our conclusions and ideas for future work.

Preliminaries

Notations

Let denote the set of real n-dimensional column vectors. For , let . For a real square matrix A, stands for the spectral radius of matrix A, which is the smallest real positive number such that , where are the eigenvalues of A. A matrix A is irreducible if and only if it cannot be transformed into a block upper-triangular form by simultaneous row/column permutations.

Some facts about digraph

Let be a weighted digraph (directed graph), where is the set of nodes, is a set of edges, and is an adjacency matrix which assigns real numbers to the edges. Let denote a graph corresponding to adjacency matrix A. A directed path (length ) is a sequence of directed edges of the form with distinct nodes. A digraph is said to be strongly connected if there exists a directed path between any two distinct nodes. For an undirected graph, the adjacency matrix is symmetric. An undirected graph is said to be connected if there exists a path between any two distinct vertices of the graph. The adjacency matrix is irreducible if and only if the digraph (undirected graph) is strongly connected (connected).

Problem statements

In order to formulate a multigroup COVID-19 model with immunity, we partition the total population into N groups, and divide the total population of ith group into five compartments, that is, a susceptible compartment , an exposed compartment , an infective compartment , an infective compartment in treatment and a recovered compartment . Figure 1 depicts the flow diagram of multigroup COVID-19 model.

Fig. 1

Flow diagram of the multigroup COVID-19 model

The model is presented as the following set of ODEs:where . All variables and parameters are described in Table 1. Let denote the infection graph of multigroup COVID-19 model, where . In model (1), we have the following important assumptions:

Table 1

Variables and parameters of multigroup COVID-19 model

Variable	Explanation
S_i	Susceptible population of ith group
E_i	Exposed population of ith group
I_i	Infective population of ith group
T_i	Infective population in treatment of ith group
R_i	Recovered population of ith group
W_i	The total population in ith group
b_i	The rate at which new individual enter the population of ith group
p_i	A fraction of b_i has antibodies (0≤p_i<1)
α_ij	The cross infection coefficient from ith group to jth group
σ	The parameter describing the psychological effect of the general public toward the infectives
μiS	The death rate of populations occurring in S_i compartment
μiE	The death rate of populations occurring in E_i compartment
μiI	The death rate of populations occurring in I_i compartment
μiT	The death rate of populations occurring in T_i compartment
μiR	The death rate of populations occurring in R_i compartment
δ_1i	The disease-related death rate of infective population in ith group who are not in treatment
δ_2i	The disease-related death rate of infective population in ith group who are in treatment
β_i	The transmission coefficient from exposed population to infective population
γ_i	The transmission coefficient from infective population to treatment
σ_i	The recovery rate of infective population
ρ_i	The proportion of recovered population who become susceptible again to the disease

The population in susceptible compartment enter into the exposed compartment by adequate contact with infective individuals not in treatment.

Infective individuals not in treatment are infectious to the susceptible population.

For different groups cross infection occurs between susceptible compartment and infective compartment.

Part of the recovered compartment become susceptible again to the disease.

The assumptions (1)–(4) are reasonable for many epidemic diseases, such as measles, severe acute respiratory syndrome (SARS), mumps, gonorrhea [13]. Because the spread of COVID-19 and SARS is very similar, all assumptions (1)–(4) are reasonable for COVID-19. For a single group, the model (1) is investigated in [13], and the incidence rate is in bilinear form . However, this bilinear incidence rate may not be suitable for some epidemic diseases, such as SARS and COVID-19. As pointed out in [24], when SARS broke out early, the incidence rate is increasing when the infective population is small. When the number of infected people grows quickly, since the government has implemented aggressive measures and policies, such as border inspections, nucleic acid testing, mask wearing, isolation, etc., the incidence rate is decreasing as the infective population is large. In order to describe this phenomenon, the authors provided the following incidence ratewhere describes the infection force of the disease and describes the effects of aggressive measures and policies.

Similar to the 2003 SARS outbreak, when COVID-19 broke out, China enforced strict lockdown measures, such as losing all non-essential businesses and encouraging or mandating that people stay in their home. As the COVID-19 spreads in other countries, several countries have also implemented strict lockdown measures. These policies will lead to a decrease in the incidence rate due to the increase in the number of infections. Hence, the bilinear incidence rate is inappropriate for COVID-19. In order to accurately describe the characteristics of the incidence rate of COVID-19, we use the following incidence rate among different groups:Note that although the parameter can describe the effects of aggressive measures and policies, the basic reproduction number has nothing to do with parameter (see Sect. 3.1 for details). Hence, in order to completely control COVID-19, we need to design appropriate control strategies (see Sect. 4 for details). Although COVID-19 has become a pandemic worldwide, it has been died out in some areas, such as Hubei, Guangdong and Zhejiang provinces in China. Thus, COVID-19 can be completely controlled in some local areas.

Generally speaking, the model (1) exists a disease-free equilibrium and an endemic equilibrium. When basic reproduction number is less than or equal to one, we will show that the disease-free equilibrium is globally asymptotically stable and the disease dies out. When the basic reproduction number is larger than one, the disease may persist in the population. Hence, when the basic reproduction number is larger than one, this paper aims to design appropriate feedback control such that the disease-free equilibrium is globally asymptotically stable.

Some lemmas

The following lemmas will be used in the Proof of Theorem 4.

Lemma 1

[25] Consider the system , where . If there exist a positive definite function V(x) ( and ), numbers and such that . Then, V(x) approaches 0 in finite time, and the setting time is estimated by .

Lemma 2

[26] If are positive numbers and , then

Main results

This section will study the property of model (1) and design appropriate control when the disease-free equilibrium is unstable. First, we will investigate the stability of the disease-free equilibrium. Second, we will study the stability of the endemic equilibrium. Third, when the disease-free equilibrium is unstable, we will design feedback control such that the disease-free equilibrium of closed-loop system is globally asymptotically stable.

Stability analysis of disease-free equilibrium

For a new model, we first need to study its well-posed problem. The following Proposition shows that the model (1) is well-posed epidemiologically and mathematically.

Proposition 1

Define feasible region aswhere denotes the total number of high-risk human population in ith group, . The feasible region is positively invariant with respect to the model (1).

Proof

By adding the equations of the model (1), one can obtainOne can obtain from (2) thatWhen , .

In the sequel, we will show that the feasible region is positively invariant with respect to the model (1). One can obtain from system (1) thatSuppose that and . Then, one has and . It is easy to getHence,If , then .

One can obtain from the second equation of system (1) thatThen, it is easy to obtainIf , then . Using the similar method, one can obtain that are nonnegative when are nonnegative. Hence, the feasible region is positively invariant with respect to the model (1). This completes the Proof of Proposition 1.

It is easy to see that is a disease-free equilibrium of the model (1), where , . For the multigroup COVID-19 model (1), we take the infected compartments to be and . Let and denote the rate of secondary infections increase and the rate of disease progression, death, recovery decrease in the ith compartment respectively. ThenAs pointed out in [27], the next generation matrix of model (1) can be presented as , where , , . Then, the basic reproduction number of model (1) can be chosen as . By direct calculation, and (a block diagonal matrix), whereIt is easy to obtain that , whereHence, , whereLetThen, is the principal submatrix of , that is, there exists a permutation matrix P such thatwhere B and D are zero matrices with appropriate dimension. Hence, , and we define the basic reproduction number as .

The following theorem states that the disease-free equilibrium is globally stable under appropriate conditions.

Theorem 1

Suppose that the infection graph is strongly connected. If , then the disease-free equilibrium is globally asymptotically stable. If , then the disease-free equilibrium is unstable.

When the infection graph is strongly connected, the adjacency matrix is irreducible, and the nonnegative matrix is also irreducible. By Perron–Frobenius Theorem, the spectral radius is an algebraically simple eigenvalue of , and the corresponding left eigenvector is positive. DefineV(t) is nonnegative for all , and if and only if . Hence, V(t) is a good candidate for Lyapunov function. The time derivative of (5) along system (1) iswhere . If , then . If , then if and only if . If , then implies thatSincethe above equal sign holds if and only , . Hence, if , then implies that or , . It is easy to see that if and only if . Hence, by LaSalle’s Invariance Principle, the disease-free equilibrium is globally asymptotically stable.

One can see from the above proof that at . If and in a neighborhood of , then . Hence, the disease-free equilibrium is unstable when . This completes the Proof of Theorem 1.

Stability analysis of endemic equilibrium

An equilibrium in the interior of is called an endemic equilibrium, where satisfiesIf and , then system (1) has a unique endemic equilibrium when [13, see the Theorem 3]. However, when , the system (1) may have multiple endemic equilibriums. If the assumption (4) does not hold, that is, , then the endemic equilibrium is unique and globally asymptotically stable under appropriate conditions.

Theorem 2

Suppose that the infection graph is strongly connected, and is an arbitrary endemic equilibrium. If , , andthen there exists a unique endemic equilibrium , and is globally asymptotically stable in the interior of .

When , the variables do not exist in the first three equations of (1). Hence, we just need to consider the equations about . LetThen for all ,If the condition (7) holds, then through direct calculation, we haveHence, the conditions of Theorem 2.2 in [9] hold. LetSimilar to the Proof of Theorem 2.2 [9], we can obtain that there exists a unique endemic equilibrium , and is globally asymptotically stable in the interior of . This completes the Proof of Theorem 2.

If the assumption (4) holds, that is, , then it is hard to obtain that the endemic equilibrium is unique and globally asymptotically stable. There may exist multiple endemic equilibrium points, and the dynamic behavior of the system (1) will be very rich. For example, the system (1) may exist bifurcation phenomenon. In this case, it may be more difficult to control the spread of COVID-19. What is the impact of multiple endemic equilibrium points on the spread of COVID-19? We will try to solve it in future work.

Multigroup COVID-19 model with control

If , then the disease-free equilibrium is unstable. In this case, we want to design appropriate control such that the disease-free equilibrium of closed-loop system is globally asymptotically stable, that is, the disease dies out under appropriate control. In reality, the government can reduce the number of exposed people and increase the number of people treated such that the disease dies out. Suppose that partially exposed people can be tested to eliminate the risk of infection, such as temperature and pharynx swab test for COVID-19. Suppose that the excluded rate (the number of newly excluded exposed population per unit time) is given as a function of by , and the treatment rate (the number of newly added treatment population per unit time) is given as a function of by . Then, the epidemic model (1) with two controllers becomeswhere .

The most natural choice is linear state feedback. The Algorithm 1 provides a systematic method to identify the groups that need to be controlled.

The following theorem states that the disease-free equilibrium of closed-loop system (8) is globally asymptotically stable. Let , where

Theorem 3

Consider the epidemic model (8) with two controllers provided in Algorithm 1. Suppose that the infection graph is strongly connected, and . If andthen the disease-free equilibrium of closed-loop system (8) is globally asymptotically stable.

When the spectral radius and there does not exist controllers, by the Theorem 1, the disease-free equilibrium is unstable. For the epidemic model (8) with two controllers, the disease-free equilibrium is also . DefineThe time derivative of (10) along system (8) isSince the controllers are obtained by Algorithm 1, . By the Algorithm 1, . Then When the condition (9) holds, then and if and only if . Hence, by LaSalle’s Invariance Principle, the disease-free equilibrium of closed-loop system (8) is globally asymptotically stable. This completes the Proof of Theorem 3.

Although the linear state feedback controllers is a natural choice, there may take a long time to approach the disease-free equilibrium. If we design controllers as , , where and is a constant, then the time required for convergence is finite. We call these controllers as finite-time controllers. The following theorem states that under appropriate choice of the disease-free equilibrium of closed-loop system (8) is globally asymptotically stable, and the convergence can be achieved in finite time.

Theorem 4

Consider the epidemic model (8) with two controllers , , where and is a constant. Suppose that the infection graph is strongly connected, and . If andthen the disease-free equilibrium of closed-loop system (8) is globally asymptotically stable, and the exposed compartment and infective compartment can achieve zero in finite time.

If the spectral radius and there does not exist controllers, one can know from the Theorem 1 that the disease-free equilibrium is unstable. For the epidemic model (8) with two controllers , , the disease-free equilibrium is also . DefineThe time derivative of (12) along system (8) isBy the Algorithm 1, . When the condition (11) holds,There exists such that where the second inequality is obtain by Lemma 2. Hence,By Lemma 1, one can obtain that V(t) approaches 0 in finite time, and the setting time is estimated by . Hence, the exposed compartment and infective compartment can achieve zero in finite time, and the setting time is estimated by .

In addition, similar to the Proof of Theorem 3, one can obtain that the disease-free equilibrium of closed-loop system (8) is globally asymptotically stable, and the details are omitted here for brevity. This completes the Proof of Theorem 4.

Remark 1

Note that when the epidemic model (8) has two controllers , , where and is a constant, the closed-loop system does not satisfy Existence and Uniqueness Theorem. However, since the two controllers are continuous, one can know from Peano’s Existence Theorem and Extension Theorem that the system (8) exists at least one solution on for any initial state in .

Remark 2

Note that in the Proof of Theorem 4, the constant c can be chosen as , whereThe setting time is estimated by

Remark 3

In reality, one needs to consider the required energy cost of different controllers. Let , , where , and are constants. Then, it is easy to obtain that and . The functions f, g can be seen as cost of energy. Hence, when the gain are the same and the linear state feedback controllers need much energy compared to the finite-time controllers. If the gain are the same and the linear state feedback controllers need less energy compared to the finite-time controllers. If we define two new switching controllersthen the switching controllers (13) can minimize the amount of the required energy cost. One can see reference [28] for more details about the energy of control of complex networks.

Simulating the COVID-19 outbreak in Hubei, Guangdong and Zhejiang provinces

This section will provide a simulation of the COVID-19 outbreak in Hubei, Guangdong and Zhejiang provinces in China to illustrate the effectiveness of the obtained results. Consider a multigroup COVID-19 model (1) with three groups which represent Hubei, Guangdong and Zhejiang provinces, respectively. In order to apply the model (1) to stimulate the progression of COVID-19 in Hubei, Guangdong and Zhejiang provinces, we need to estimate the parameters in Table 1. The topology of cross infection among Hubei, Guangdong and Zhejiang provinces is provided in Fig. 2. Because there is no direct data about cross infection coefficient among three provinces, we set the cross infection coefficients as the migration ratio of the Spring Festival in 2020 which are presented in Fig. 2, and these data come from the Baidu migration map (https://qianxi.baidu.com/2020/).

Fig. 2

The topology of cross infection among Hubei, Guangdong and Zhejiang provinces in China

The matrix isOne can see that the matrix A is irreducible, and the directed graph is strongly connected. As pointed out in [3], the incubation period of the COVID-19 has been reported to be between 2 and 14 days. We take 5.2 days as time units which is the mean incubation period of the COVID-19 [29]. Hence, the rates of latent individuals becoming symptomatic (incubation rates) . One can obtain from the National Bureau of Statistics of China [30] that the annual death rate is equal to (number of deaths/annual average number of population), and the annual birth rate is equal to (number of births/average number of population). Hence, we can get that the death rate within 5.2 days is approximately equal to , and the birth rate within 5.2 days is approximately equal to . We assume that the death rates of population occurring in , , , compartments are equal to , that is, . According to the data in [31], for COVID-19 the death rate of population occurring in compartment is . We assume that the disease-related death rate of infective population who are not in treatment , and the transmission coefficient from infective population to treatment . Since the disease-related death rate of infective population who are in treatment should be less than the death rate of population occurring in compartment, we assume that . In addition, there are lack of real data on the recovery rate of infective population, the fraction of total population who already have antibodies, and the proportion of recovery population who become susceptible again. Hence, we assume that . We assume that the rates at which new individual enter the population of each province are equal to the birth rate, that is, .

With above values of parameters, the matrix isSince , one can know from Theorem 1 that the disease-free equilibrium is unstable. As pointed out in [32], the obtained basic reproduction number is reasonable for COVID-19 from China and overseas.

Since on January 23, Wuhan began to implement aggressive measures and policies, such as the closure of the city, nucleic acid testing, mask wearing, we set January 23, 2020 as time zero. We assume that . Since other parts of the China have implemented the closure of the city from January 23, the number of people who stay susceptible in each province is similar to that of its resident population. Hence, there are million in Hubei province, million in Guangdong province, and million in Zhejiang province. The number of people who stay exposed in each province is similar to that of its capital cities, that is, million in Hubei province, million in Guangdong province, and million in Zhejiang province. Table 2 depicts the daily reported cases in Hubei, Guangdong and Zhejiang provinces from January 23 to May 3 (see National Health Commission of the People’s Republic of China [33]). On January 23, one can obtain from Health Commission of Hubei province that the number of infected people in Hubei is , and the recovery population of Hubei is . On January 23, one can obtain from Health Commission of Guangdong province that the number of infected people in Guangdong is , and the recovery population of Guangdong is . On January 23, one can obtain from Health Commission of Zhejiang province that the number of infected people in Zhejiang is , and the recovery population of Zhejiang is .

Table 2

Daily reported cases in Hubei, Guangdong and Zhejiang provinces from January 23 to May 3

Date	1/23	1/27	1/31	2/4	2/8	2/12	2/16	2/20	2/24
Hubei	549	2567	6738	15,679	24,881	43,455	49,847	49,156	43,369
Guangdong	51	184	511	838	994	956	845	664	535
Zhejiang	43	172	584	832	902	818	701	555	410
Date	2/28	3/3	3/7	3/11	3/15	3/19	3/23	3/27	3/31
Hubei	34,715	25,905	19,710	14,427	9605	6287	4318	2526	1283
Guangdong	370	214	95	60	47	62	87	118	136
Zhejiang	217	105	45	5	17	15	18	25	30
Date	4/4	4/8	4/12	4/16	4/21	4/25	4/29	5/3
Hubei	648	401	244	129	97	12	0	0
Guangdong	129	108	93	90	63	36	23	11
Zhejiang	32	31	27	21	15	10	4	2

Based on the above data, we give the simulation of the COVID-19 outbreak in Hubei, Guangdong and Zhejiang provinces. Based on the model (1), Fig. 3 depicts the number of active infections in the three provinces under no control. One can see from Fig. 3 that the disease-free equilibrium is unstable and the number of COVID-19 infections in the three provinces will increase rapidly. Moreover, because Guangdong has a large population, its peak is higher than other provinces.

Fig. 3

The number of active infections in Hubei, Guangdong and Zhejiang provinces under no control

In order to make the disease-free equilibrium is globally asymptotically stable, we design two linear feedback controllers using Algorithm 1. By using Algorithm 1, we can obtain that the controlled set . If we choose , then the conditions of Theorem 3 hold. One can see from Fig. 4a that the disease-free equilibrium is globally asymptotically stable under proposed controllers. However, the number of active infections in Guangdong and Zhejiang provinces is much higher than the actual data (see Fig. 4b). The reason for this phenomenon is that we assume that there exists persistent mutual infections in three provinces. In fact, after the closure of Wuhan on January 23, the mutual infection between the three provinces should be very small. Hence, our simulation results can be regarded as the worst case that there are still many people moving among provinces.

Fig. 4

a The number of active infections in Hubei, Guangdong and Zhejiang provinces under linear control; b The actual number of active infections in Hubei, Guangdong and Zhejiang provinces

If we design controllers as , , where and , then the conditions of Theorem 4 hold. One can see from Fig. 5 that the disease-free equilibrium is globally asymptotically stable under proposed controllers. However, these controllers do not significantly deduce the convergence time. Moreover, the peak number of active infections is much higher than linear control and actual data (see Fig. 4b). The reason for this phenomenon may be that the finite time controllers consume higher energy than linear control. This is a very meaningful phenomenon and we will solve it in our future work.

Fig. 5

The number of active infections in Hubei, Guangdong and Zhejiang provinces under finite time control

Concluding remarks

We have introduced a new multigroup COVID-19 model with immunity, and studied the stability of disease-free equilibrium and endemic equilibrium. When the basic reproduction number is less than or equal to one, and the infection graph is strongly connected, the disease-free equilibrium is globally asymptotically stable. When the disease-free equilibrium is unstable, some effective control strategies were given for the proposed multigroup COVID-19 model. The proposed control strategy contains two methods, one is to increase the number of treatment, and the other is to increase the detection to exclude the virus-free people. Future research will consider the control problem for the case that the infection graph infection graph is not strongly connected.