Next Article in Journal
AFB-GPSR: Adaptive Beaconing Strategy Based on Fuzzy Logic Scheme for Geographical Routing in a Mobile Ad Hoc Network (MANET)
Next Article in Special Issue
A Robust Deep Learning Approach for Accurate Segmentation of Cytoplasm and Nucleus in Noisy Pap Smear Images
Previous Article in Journal
Solving the Problem of Elasticity for a Layer with N Cylindrical Embedded Supports
Previous Article in Special Issue
Genomic Phylogeny Using the MaxwellTM Classifier Based on Burrows–Wheeler Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MPC Controllers in SIIR Epidemic Models

by
Nikita Kosyanov
*,†,
Elena Gubar
*,† and
Vladislav Taynitskiy
*,†
Faculty of Applied Mathematics and Control Processes, Saint Petersburg State University, Universitetsky Pr. 35., Peterhof, Saint Petersburg 198504, Russia
*
Authors to whom correspondence should be addressed.
All authors contributed equally to this work.
Computation 2023, 11(9), 173; https://doi.org/10.3390/computation11090173
Submission received: 6 July 2023 / Revised: 9 August 2023 / Accepted: 22 August 2023 / Published: 4 September 2023
(This article belongs to the Special Issue 10th Anniversary of Computation—Computational Biology)

Abstract

:
Infectious diseases are one of the most important problems of the modern world, for example, the periodic outbreaks of coronavirus infections caused by COVID-19, influenza, and many other respiratory diseases have significantly affected the economics of many countries. Hence, it is therefore important to minimize the economic damage, which includes both loss of work and treatment costs, quarantine costs, etc. Recent studies have presented many different models describing the dynamics of virus spread, which help to analyze the epidemic outbreaks. In the current work we focus on finding solutions that are robust to noise and take into account the dynamics of future changes in the process. We extend previous results by using a nonlinear model-predictive-control (MPC) controller to find effective controls. MPC is a computational mathematical method used in dynamically controlled systems with observations to find effective controls.

1. Introduction

Infectious diseases are some of the modern world’s most important problems, e.g., influenza, the periodic outbreaks of coronavirus infections caused by COVID-19, and many other respiratory diseases have significantly affected the economies of many countries (lockdowns, masks, and reduced production efficiency) [1,2]. It is therefore important to minimize the economic damage, which includes both loss of work and treatment costs. Recent studies have presented many different models describing the dynamics of virus spread, which help to analyze and explain epidemic outbreaks [3,4]. Previously, it was shown that several infections can circulate in the population at the same time. In this case, the controlled Susceptible–Infected–Recovered model with multiple infected subgroups can be used [5]. The solution of this model is usually found by using classical optimality criteria, e.g., the Pontryagin’s maximum principle [6]. In this study, we focus on finding solutions that are robust to noise and take into account the dynamics of future changes in the process. We extend previous results by using a nonlinear model-predictive-control (MPC) controller to find effective controls [7]. MPC is a computational mathematical method used in dynamically controlled systems with observations to find effective controls.
It has been shown that several virus species, or several different virus strains of the same species, can circulate simultaneously in a population. In this case, the subpopulation of infected individuals can be divided into several subgroups. In our study, two types of virus were considered to simplify the calculations, but the results obtained can be generalized to a larger number of circulating strains.
Experience with influenza epidemics and the SARS-CoV-2 pandemic has shown that there can be several waves of illness in a population during a single epidemic season. To reduce the morbidity of an epidemic wave, the government and health authorities apply various measures to prevent and contain the epidemic. The use of different preventive and curative measures during a single wave reduces the number of cases and the rate at which the virus spreads in the population. Thus, the rise of the disease in the next wave of the epidemic begins with new initial states of the system. A review of the literature led to the idea of using predictive control models to develop effective control that takes into account changes in the environment. MPC controllers allow for the resulting control to be adapted by dividing the entire time interval into smaller sub-intervals. At each interval, the control system is checked to see if the constraints are satisfied, and if not, the control action is adjusted.
The main applications of this method are in the petrochemical [8], woodworking [9], and energy [10,11] industries, where, among other things, the stability of the obtained controllers and the possibility of their automatic adjustment are important. The MPC controller finds a solution in the form of piecemeal functions on the control interval, but the optimization problem is solved over a larger time interval—the forecast horizon. The MPC controller finds the solution in the form of piecewise constant functions on the control interval, but the optimization problem is solved on the forecast horizon.
As the main contribution of this work, we formulated a control problem and solved it using an MPC controller, a series of numerical experiments are carried out to confirm the effectiveness of the solutions obtained by the MPC controller for the SIIR model.
The paper is organized as follows. Section 2 presents the epidemic models and modified infection rate parameter. Section 3 shows the application of MPC-controllers to a SIIR epidemic model. In Section 4, we use numerical simulation to illustrate our results. The paper is concluded in Section 5.

2. Epidemic Process in Total Urban Population

Our model is based on the extended Susceptible–Infected–Recovered model with two types of viruses, which circulate simultaneously in an urban population [6]. A size of population is equal to N. According to these assumptions, we can consider four groups: the Susceptible, the Infected by virus V 1 , the Infected by virus V 2 , and the Recovered group. Susceptible is a group of people who are not infected by viruses, but can be infected by one or both types of viruses and who have no immunity to the viruses. Depending on the viral strength, we observe that the number of people infected by virus 1 or by virus 2 can be different, and people infected by virus V 1 or V 2 belong to the subgroup Infected. The Recovered subgroup consists of people who have recovered from the infection. The mixing of populations allows for viruses to spread quickly, and each person in the population is assumed to have an equal probability of coming into contact with others. Therefore, when an infected person interacts with a susceptible person, the virus can spread. A virus with higher virulence, we assume, will have a higher probability of success in spreading when an encounter occurs between the infected and the susceptible. The scheme of process is represented in the Figure 1.
We model the spread of the virus in the population based on the classical SIR epidemiological model, where a system of differential equations is used to describe the proportion of the population as a function of time. Then, at time t, n S ( t ) , n I 1 ( t ) , n I 2 ( t ) , and n R ( t ) correspond to the fractions of the population that are Susceptible, Infected by virus v 1 , Infected by virus v 2 , and Recovered, respectively, and for all t the condition N = n S ( t ) + n I 1 ( t ) + n I 2 ( t ) + n R ( t ) is justified. Define S ( t ) = n S ( t ) N , I 1 ( t ) = n I 1 ( t ) N , I 2 ( t ) = n I 2 ( t ) N , R ( t ) = n R ( t ) N as proportions of Susceptible, Infected, and Recovered. At the beginning of the epidemic t = 0 , most people in the population belong to the susceptible subpopulation, a small group in the total population is infected and other people are in the recovered subpopulation. Therefore, initial states are 0 < S ( 0 ) = S 0 < 1 , 0 < I 1 ( 0 ) = I 1 0 < 1 , 0 < I 2 ( 0 ) = I 2 0 < 1 , R ( 0 ) = 1 S 0 I 1 0 I 2 0 .
We have extended the simple SIR model introduced by [12] to describe the situation with two virus types:
S ˙ ( t ) = δ 1 S ( t ) I 1 ( t ) δ 2 S ( t ) I 2 ( t ) ; I 1 ˙ ( t ) = ( δ 1 S ( t ) σ 1 u 1 ( t ) ) I 1 ( t ) ; I 2 ˙ ( t ) = ( δ 2 S ( t ) σ 2 u 2 ( t ) ) I 2 ( t ) ; R ˙ ( t ) = ( σ 1 + u 1 ( t ) ) I 1 ( t ) + ( σ 2 + u 2 ( t ) ) I 2 ( t ) ;
where δ i are infection rates for virus V i , i = 1 , 2 , σ i are recovery rates (see Figure 1). Infection rate is defined as a product of the contact rate l and transmissibility of infection, i.e., probability of transmission infection during the contact, δ i 0
δ i = l δ i 0 ( n I i N ) .
Infection rate is integrated into the evolution of mutation process to epidemics in the urban population.
Medical treatment or quarantine reduces the number of the infected individuals in the urban population. These prevention measures can be interpreted as control parameters in the system defined as u = ( u 1 , u 2 ) ; here, u i is a fraction of the infected which are quarantined or under intensive medical treatment. Recovered rates are inversely proportional to disease duration T ¯ , hence σ i = 1 T ¯ .
The objective function: We will minimize the overall cost in time interval [ 0 , T ] [13,14,15,16]. At any given t following costs exist in the system:
  • f i ( I i ( t ) ) , i = 1 , 2 ¯ is a function measuring the direct costs of infection, with a non-decreasing, twice differentiable function, such that f i ( 0 ) = 0 , f i ( I i ( t ) ) > 0 for I i ( t ) > 0 . These costs can be interpreted as the price the government or medical organization pays for the health system to treat the infected and increase the probability of recovery for the treated individuals. Infected individuals who are not detected because they are asymptomatic are not treated, and therefore there are no direct costs for them.
  • The treatment or medical cure of those infected produces costs for society. The functions h i ( u i ( t ) ) represents the costs of treating the infected to keep them out of the labor market until recovery. The function h i ( u i ( t ) ) is increasing on the argument u i ( t ) , is twice differentiable, and is such that h i ( u i ( 0 ) ) = 0 .
  • The function g ( R ( t ) ) is a benefit rate. These are the bonuses that society receives from Recovered individuals who have returned to work and started contributing to society. The function g ( t ) is non-decreasing and differentiable and g ( 0 ) = 0 .
  • Constants e I 1 , e I 2 , e R represent the costs and benefit for Infected and Recovered in the end of the epidemic.
Aggregated system costs is given by
J = 0 T f 1 ( I 1 ( t ) ) + f 2 ( I 2 ( t ) ) g ( R ( t ) ) + h 1 ( u 1 ( t ) ) + h 2 ( u 2 ( t ) ) d t + e I 1 I 1 ( T ) + e I 1 I 2 ( T ) e R R ( T ) .
and the optimal control problem is to minimize the cost, i.e.,
min { u 1 , u 2 } J
To simplify the analysis, we consider the case where e I 1 = e I 1 = e R = 0 .
Here, we assume that as a control parameters which help to reduce numbers of infected in human population, medical treatment or isolation can be considered, then define variable u = ( u 1 , u 2 ) as control 0 u 1 ( t ) 1 , 0 u 2 ( t ) 1 , for all t.

2.1. SIIRS Models

The SIIR models discussed in the previous section are best suited to describe epidemic processes in which viruses mutate slowly. Because of the small number of mutations, individuals can develop lifelong immunity to the virus. However, there are viruses that mutate rapidly and for these, SIIRS (Susceptible–Infected–Recovered–Susceptible) models have been developed (see Figure 2) to account for the decline of immunity in the population over time.
The system of differential equations describing the dynamics of state change is equal to
S ˙ ( t ) = γ R ( t ) δ 1 S ( t ) I 1 ( t ) δ 2 S ( t ) I 2 ( t ) , I 1 ˙ ( t ) = ( δ 1 S ( t ) σ 1 u 1 ) I 1 ( t ) , I 2 ˙ ( t ) = ( δ 2 S ( t ) σ 2 u 2 ) I 2 ( t ) , R ˙ ( t ) = ( σ 1 + u 1 ) I 1 ( t ) + ( σ 2 + u 2 ) I 2 ( t ) γ R ( t ) .
Here, as in the SIIR model, the parameters δ i , σ i , i = 1 , 2 ¯ are responsible for the intensity of infection and cure, u i —the proportion of infected people on medication, and γ —the intensity of immune decline in the population. The objective function for such tasks is no different from (2).
Note that due to the gradual decline of immunity, the epidemic process in the SIIRS problem can continue indefinitely. For this reason, when solving this problems we should choose large values for T. However, in such a case there is a risk (due to the high mutation rate of viruses) that the model parameters will no longer be relevant. Therefore, methods that can adapt to changes in the process by predicting its behavior can be effective in solving such problems.

2.2. SWIIRS Model

In 2020, the SWIIRS (Susceptible–Warned–Infected–Recovered–Susceptible) model was proposed, which takes into account citizens’ awareness of the epidemic [17]. This model divides the population into five groups instead of four. An additional group of individuals is the group of warned W individuals who are aware of the virus (see Figure 3). This approach allows a more accurate description of epidemic processes in large countries, as the speed of spread of the virus often depends not only on the characteristics of the virus itself, but also on the behavior of members of the population.
The changes in population groups for the controlled SWIIRS model are described by a system of five differential equations:
S ˙ ( t ) = γ R ( t ) δ 1 S S ( t ) I 1 ( t ) δ 2 S S ( t ) I 2 ( t ) ( η W ( t ) + u 3 ) S ( t ) , W ˙ ( t ) = ( η W ( t ) + u 3 ) S ( t ) ( δ 1 W I 1 ( t ) + δ 2 W I 2 ( t ) + σ 3 ) W ( t ) , I 1 ˙ ( t ) = ( δ 1 S S ( t ) + δ 1 W W ( t ) σ 1 u 1 ) I 1 ( t ) , I 2 ˙ ( t ) = ( δ 2 S S ( t ) + δ 2 W W ( t ) σ 2 u 2 ) I 2 ( t ) , R ˙ ( t ) = ( σ 1 + u 1 ) I 1 ( t ) + ( σ 2 + u 2 ) I 2 ( t ) + σ 3 W ( t ) γ R ( t ) .
In the system (4), the constants δ i S and δ i W describe the infection rates for groups S and W, respectively, and σ i – the cure rate i = 1 , 2 ¯ . The value σ 3 is responsible for the probability that a knowledgeable individual will acquire immunity to viruses (e.g., go to a health facility to get a vaccine). The intensity of the decline in immunity is described by γ . Two parameters are responsible for informing (warning) members of the population:
  • Parameter u 3 —percentage of susceptible individuals informed by the media (become Warned);
  • Parameter η —probability of informing the susceptible individual when communicating with the warned individual.

2.3. Controllable Infection Rate Parameter

Consider a model (1) with one virus I i and let p i be the probability of virus transmission V i from a diseased individual to a susceptible one, and l be the average number of contacts of an individual (with other members of the population) per unit of time. Consider a random contact of a susceptible individual. Infection with the virus will only occur if two events occur simultaneously: the contact was with the infected individual, and the virus was transmitted on contact.
Define two incompatible events:
  • Event A— a susceptible individual becomes infected through one of an l contacts;
  • Event B—with a similar number of contacts, infection will not occur.
Obviously that events A and B form a complete group of events (see Figure 4), then
p ( A ) + p ( B ) = 1 .
The event B will only occur if there is no infection at any of the l contacts. In other words
p ( B ) = 1 p i I i ( t ) l .
From Equations (5) and (6) we see that the probability of transmission to a susceptible individual by l contact with infected individuals can be defined as
p ( A ) = 1 1 p i I i ( t ) l .
The same formula can be obtained in another way. The probability of the event A is given by
p ( A ) = p i I i ( t ) + 1 p i I i ( t ) p i I i ( t ) + 1 p i I i ( t ) 2 p i I i ( t ) + 1 p i I i ( t ) l p i I i ( t ) = p i I i ( t ) 1 p i I i ( t ) l 1 p i I i ( t ) = 1 1 p i I i ( t ) l .
This probability p ( A ) is called the modified infection rate parameter and it is denoted by δ ˜ i ( I i ( t ) , l ) .
Let us redefine u i ( · ) (intensive virus treatment) as u i 1 ( · ) . Let us introduce two types of control into the resulting system: isolation of infected individuals and reduction in the number of contacts in the whole system. Let u i 2 ( · ) be the fraction of infected individuals isolated, and u l ( · ) be the fraction of contacts prevented by reducing the mobility of individuals in the whole population. In this case we obtain
δ ˜ i ( I i ( t ) , l , u i 2 ( t ) , u l ( t ) ) = 1 p i I i ( t ) ( 1 u i 2 ( t ) ) l ( 1 u l ( t ) ) ,
where 0 u i 2 ( t ) , u l ( t ) 1 . Then the system of differential equations of the corresponding single-variable SIR problem is
S ˙ ( t ) = 1 p i I i ( t ) ( 1 u i 2 ( t ) ) l ( 1 u l ( t ) ) S ( t ) , I i ˙ ( t ) = 1 p i I i ( t ) ( 1 u i 2 ( t ) ) l ( 1 u l ( t ) ) S ( t ) ( σ 1 + u i 1 ( t ) ) I i ( t ) , R ˙ ( t ) = ( σ 1 + u i 1 ( t ) ) I 1 ( t ) .
If we perform similar transformations to (8) for the SIIR dual virus model, the calculated formulae for the controllable infection parameters are as follows
δ ˜ i = δ ˜ i ( I 1 ( t ) , I 2 ( t ) , l , u 1 2 , u 2 2 , u l ) = p i I i ( t ) ( 1 u i 2 ) 1 p 1 I 1 ( t ) ( 1 u 1 2 ) p 2 I 2 ( t ) ( 1 u 2 2 ) l ( 1 u l ) 1 p 1 I 1 ( t ) ( 1 u 1 2 ) p 2 I 2 ( t ) ( 1 u 2 2 ) ,
where i = 1 , 2 ¯ . If we substitute these parameters into the system (1), the differential equations take the following form:
S ˙ ( t ) = δ ˜ 1 S ( t ) δ ˜ 2 S ( t ) , I 1 ˙ ( t ) = δ ˜ 1 S ( t ) ( σ 1 + u 1 ) I 1 ( t ) , I 2 ˙ ( t ) = δ ˜ 2 S ( t ) ( σ 2 + u 2 ) I 2 ( t ) , R ˙ ( t ) = ( σ 1 + u 1 ) I 1 ( t ) + ( σ 2 + u 2 ) I 2 ( t ) .
For the introduced controls we also need to redefine the cost functions. Let h 1 2 ( u 1 2 , I 1 ) , h 2 2 ( u 2 2 , I 2 ) , h l ( u l ) be non-reducing functions corresponding to the costs of the controls u 1 2 , u 2 2 , and u l . Then the function (2) takes the form
J ˜ = 0 T [ f 1 ( I 1 ( t ) ) + f 2 ( I 2 ( t ) ) g ( R ( t ) ) + h l ( u l ( t ) ) + h 1 1 ( u 1 1 ( t ) ) + h 1 2 ( u 1 2 ( t ) + h 2 1 ( u 2 1 ( t ) ) + h 2 2 ( u 2 2 ( t ) ] d t .

3. MPC Controllers

Model predictive control (MPC) is a set of methods for finding effective control based on model prediction. The basic idea of MPC is to estimate the future behavior of the controlled system at a finite time interval and to compute an effective control signal that minimizes a given objective function, given the constraints on the system (see Figure 5).
The next algorithm of the model-predictive-control method includes the following steps:
1.
Modeling of an object or process using mathematical and computer modeling techniques;
2.
Defining the main objectives to be achieved by management, selecting the forecasting horizon;
3.
Predicting the behavior of the process under the influence of the control signal;
4.
Find an effective control signal, taking into account the full range of constraints imposed on the control and manipulated variables;
5.
Select the control horizon in which the effective control found is applied;
6.
Measuring the actual process states at the end of the control horizon, which are taken as the new initial conditions;
7.
Shift the forecast horizon by the size of the management horizon and repeat steps 3–6.

3.1. Linear Prediction Model

We consider a linear controllable and observable finite-difference system:
x ( t + 1 ) = A ( t ) x ( t ) + B ( t ) u ( t ) , y ( t ) = C ( t ) x ( t ) , t [ 0 , T ] , u ( t ) U .
where A ( t ) , B ( t ) , C ( t ) are some matrix functions of dimensions n × n , n × m , and n × k , respectively, U is set of valid controls (is a subset of piecemeal functions), and [ 0 , T ] is a given time interval.
Denote by x r , u r the reference values of state and control variables. The MPC is a computational tool that finds the minimum either of its built-in target functionality (15) or of a custom function that the user must describe. Efficient control is determined using numerical Knows-What-It-Knows algorithms (KWIK) [18].
0 T w x , ( x x r ) 2 + w u , ( u u r ) 2 d t , t [ 0 , T ] ,
where w x and w u weight vectors for variables x and u. However, the system (1) should be transformed to have a possibility to apply the MPC controller. In our case we find a control u ( t ) that provides the minimum to the objective function (2). The system (1) must be linearized before the MPC can be used. The predicted values are determined by the formula:
y ¯ = G x ( t ) + H v + F u ( t ) , where
  • x ( t ) , u ( t ) —the known current states and controls of the system,
  • y ¯ = y ¯ ( t + 1 ) y ¯ ( t + 2 ) . . . y ¯ ( t + N ) —model prediction,
  • v = v ( t + 1 ) v ( t + 2 ) . . . v ( t + N ) —planned control signals,
  • H = 0 0 . . . 0 C B 0 . . . 0 . . . . . . . . . . . . C A N 2 B C A N 3 B . . . 0 ,
  • G = C A C A 2 . . . C A N ,
  • F = C B C A B . . . C A N 1 B .
The solution of the obtained problem is carried out using standard methods for solving quadratic programming problems. Let
v * = v * ( t + 1 ) v * ( t + 2 ) . . . v * ( t + N )
solution of the problem (15), then v * ( t + 1 ) will be called the effective control for the original system (14) on the interval [ t , t + 1 ] obtained using MPC method. After applying the control, we obtain the new initial data x ( t + 1 ) and u ( t + 1 ) to make the following forecast and search for an effective control on the interval [ t + 1 , t + 2 ] .

3.2. Linearization of S I 1 I 2 R Problem

In the current section, we discuss the application of the MPC controller to the extended S I 1 I 2 R system [19].
Suppose, that the initial states for all states S ( 0 ) , I 1 ( 0 ) , I 2 ( 0 ) , R ( 0 ) are known. Otherwise, we can always measure it. Define S I I R ( 0 ) as a vector of initial states.
S I 1 I 2 R ( 0 ) = S 0 I 1 0 I 2 0 R 0 .
Another assumption is that the initial control u ( 0 ) = u 0 is defined. We divide the initial period of time [ 0 , T ] into k equal parts, then the objective function (2) by the integral property can be written as:
J = 0 T [ f 1 ( I 1 ( t ) ) + f 2 ( I 2 ( t ) ) g ( R ( t ) ) + h 1 ( u 1 ( t ) ) + h 2 ( u 2 ( t ) ) ] d t = J 1 + + J k , where J i = T i 1 T i [ f 1 ( I 1 ( t ) ) + f 2 ( I 2 ( t ) ) g ( R ( t ) ) + h 1 ( u 1 ( t ) ) + h 2 ( u 2 ( t ) ) ] d t .
Here, the total length of time interval is divided into several segments [ T i 1 , T i ] , where T i 1 = ( i 1 ) T k and T i = i T k , i = 1 , k ¯ .
Obviously, that min u U J 1 + min u U J 2 + . . . + min u U J k min u U J . Then, instead of the original problem, we consider k-th sub-interval with objective function of the form (20), taking into account the system (1):
J i min u U .
The initial conditions for the first iteration are defined in (18), for each subsequent problem they can be found from the previous task as:
S I 1 I 2 R ( T i ) = S ( T i ) I 1 ( T i ) I 2 ( T i ) R ( T i ) = S i I 1 i I 2 i R i , u i = u ( T i ) , i = 1 , k 1 ¯ iteration number .
Let us redefine S I 1 I 2 R ( t ) = S ( t ) I 1 ( t ) I 2 ( t ) R ( t ) as a vector x ( t ) = x 1 ( t ) x 2 ( t ) x 3 ( t ) x 4 ( t ) . The system (1) is also valid for initial values, i.e.,
x ˙ 1 i = f 1 ( x i , u i , t ) = ( δ 1 x 2 i + δ 2 x 3 i ) x 1 i ; x ˙ 2 i = f 2 ( x i , u i , t ) = ( δ 1 x 1 i σ 1 u 1 i ) x 2 i ; x ˙ 3 i = f 3 ( x i , u i , t ) = ( δ 2 x 1 i σ 2 u 2 i ) x 3 i ; x ˙ 4 i = f 4 ( x i , u i , t ) = ( σ 1 + u 1 i ) x 2 i + ( σ 2 + u 2 i ) x 3 i .
Consider some disturbance of this system:
  • x ˜ ( t ) = x i ( t ) + Δ x ( t ) , where Δ x ( t ) —state variable variation vector;
  • u ˜ ( t ) = u i ( t ) + Δ u ( t ) , where Δ u ( t ) —control variable variation vector.
Then, x ˙ i ( t ) + Δ x ˙ ( t ) = f ( x ˜ , u ˜ , t ) . The right hand sides of system (22) is decomposed into a Taylor row:
x ˙ i ( t ) + Δ x ˙ ( t ) = f ( x i ( t ) , u i ( t ) , t ) + d f d x | x i , u i Δ x ( t ) + d f d u | x i , u i Δ u ( t ) + O ( ( Δ x ) 2 , ( Δ u ) 2 ) .
Here, O ( ( Δ x ) 2 , ( Δ u ) 2 ) is the second-order residual term of smallness. Subtract from (22) and obtain
Δ x ˙ ( t ) = d f d x | x i , u i A i ( t ) Δ x ( t ) + d f d u | x i , u i B i ( t ) Δ u ( t ) .
As a result, the linearized controlled system is obtained. We substitute the solution ( x ¯ , u ¯ ) of the linearized system into the original model (1), (20) and obtain the following value of the objective function J i ( x ¯ , u ¯ ) = J ¯ i . Let ( x * , u * ) be a solution that provides a minimum to the functional (20), then we can denote by an error in the solution the value of the ε i , which is determined by the formula:
ε i = | J i ( x ¯ , u ¯ ) J i ( x * , u * ) | , i = 1 , k ¯ .
Going to the limit k , we have:
lim k x ˙ i = A i x + B i u = > lim k ( x ¯ , u ¯ ) = ( x * , u * ) = > lim k ε i = 0 .
This means that the greater the number of split points, the more accurate the solution will be.
We have conducted a series of numerical simulations to check whether the MPC controller solution converges to the optimal solution as the parameter k increases (the other parameters remain the same). It can be seen that all the costs (standard, efficient, and optimal) increase as k increases. This is due to the fact that the accuracy of the integral calculation in the function increases. In test no. 2, the value of the effective costs is extremely high, due to the fact that it is recommended to use the value k T . Starting from test no. 5, the values of effective costs and optimal costs differ insignificantly, so in practice we use the values k = 2 T or k = 4 T . However, the program execution time in the second case (test no. 9) is 10 times longer than in the first case (test no. 5).

3.3. Nonlinear MPC-Controller

Figure 6 shows the work of the nonlinear MPC controller [20]. The MPC controller construct the efficient control ( u * ) i and predictive values of the objective function using a vector of outputs y i of the system. The main difference between nonlinear and linear MPC controller is the linearization of the optimization problem before the application of the solver [21,22].
The received control is inputted to the dynamic system, in the next iteration MPC-controller obtains another output. The process of constructing a solution is complete when the control horizon is reached. The final effective control strategy is chosen to minimize the objective function on the predictive horizon. Since numerical algorithms are used to find the solution, the resulting effective value can be suboptimal.

4. Numerical Simulations

This section presents the results of a numerical simulation that corroborates our theoretical results. The numerical experiments are aimed at obtaining effective treatment strategies using nonlinear MPC controllers. We selected MATLAB as our platform for programming and numerical calculations for a series of experiments. MATLAB offers an extensive library of the model-predictive-control method, including various examples. Nonlinear-MPC controllers were employed in our study. These controllers utilize the linearization algorithm described in the article to find solutions. Nonlinear MPC is adequate to address these issues. However, larger systems may demand more precise predictive models based on machine learning and neural networks. For first three experiments, two variants of the SARS-CoV-2 virus were selected for study: Alpha [23] and Omicron [24]. The epidemic situation was considered for a period of 180 days (6 months). During this period, six MPC controllers were used sequentially with a forecast horizon of 30 days (1 month) and a control horizon of 21 days (3 weeks). Time dependence of MPC controller accuracy and performance on k parameter are shown in Table 1. The values of the main system parameters are shown in Table 2.
The serial application of multiple nonlinear MPC controllers allows the constraints, system parameters and objective function coefficients to be adjusted step by step. To simplify the calculations, we will use the same parameters in each iteration of the process.
Experiment 1. In the first experiment, we considered the epidemic process with a small number of contacts, which is possible in cases where most of the infected comply with the quarantine. It can be seen that applying control is having a positive impact on stopping the virus spreading (Figure 7). On the effective control graph (Figure 8), you can notice “jumps” due to the fact that a new MPC controller is connected each new month (30th day). However, the resulting control greatly reduces the damage from the epidemic compared to a situation where we do not control the epidemic process. So the cost in the uncontrolled case is J 0 = 3159 monetary units (m.u.) and with the effective strategy J e f f = 1016 m.u. The benefit of using control is therefore 2143 m.u.
Experiment 2. A second experiment (Figure 9 and Figure 10) was run with the average number of contacts per day k = 3 . This experiment is the most different from the others. The main wave of the epidemic does not occur in the first few days, but only at the beginning of the third month. The highest costs are in the same period. After t = 60 , the effective control is similar to that of experiments 1 and 3. With this strategy, it is possible to reduce the cost of the epidemic from J 0 = 3613 m.u. to J e f f = 1045 m.u., which saves 2568 monetary units.
It can also be seen that as the average number of contacts increased, so did the values of the objective functions. Just as importantly, the benefits of using MPC controllers have also increased.
Experiment 3. In the third simulation the number of contacts is k = 10 , which approximates the real situation. Figure 11 shows that with a large number of contacts, almost all agents are infected by t = 5 . However, when efficient control is used, the infection peaks are significantly reduced (from 0.63 to 0.29 and from 0.26 to 0.14 ). As in the previous experiments, the use of efficient control is economically beneficial, reducing costs from J 0 = 4282 m.u. to J e f f = 860 m.u. The benefit of using an effective strategy is 3422 m.u. (see Figure 12).
It should also be noted that despite the largest value of k the costs of using effective management in Experiment 3 turned out to be the smallest. This may be due to the fact that when solving the problem, they are linearized and the accuracy of the solution can vary from the parameters of the problem. At the same time, all three experiments showed that effective control can reduce the economic damage from the epidemic, which shows the practical importance of using the model predictive control in SIIR problem.
Experiment 4. In the fourth experiment we examined the SIIR model with the data taken from the paper “Optimal control of heterogeneous mutating viruses” [5]. Table 3 shows the parameter values used in the simulation.
Figure 13 shows the population states in the uncontrolled and controlled cases. It can be seen that the use of the control found with MPC significantly reduces the number of infections. In the uncontrolled case, the maximum proportion of the number of infections is I 1 m a x = 0.594 , I 2 m a x = 0.386 , and in the controlled case, I 1 m a x = 0.32 , I 2 m a x = 0.18 , corresponding to the initial proportion of infected in the population.
After the computer simulation, the solution of the formulated problem was obtained and the values of the functional were found (see Figure 14).
The strategies that increased to u 1 = 0.26 and u 2 = 0.24 during the first three days of the epidemic and then decreased to zero were found to be effective. The resulting control signals are similar in appearance to the continuous controls presented in the experiments of this paper. The aggregated costs are equal to 41.6 m.u. in the controlled case and 227.51 m.u. in the uncontrolled case. Thus, the benefit of using the control is 185.91 monetary units, which is slightly less than using the optimal control 189.17 m.u. In Figure 15, we can see the structure of the optimal control for the model and the data from this experiment. It can be observed that the control structure in the case of effective control is similar to the structure in the case of optimal control.
Experiment 4 showed that the developed software package finds near-optimal solutions to classical SIIR problems and can be used as an alternative to classical optimal control methods.
Experiment 5. As the software implementation of the MPC method showed good accuracy for standard tasks, Experiment 5 was conducted using a modified SIIRS model (with controlled infection rate parameter). Data for the experiment were obtained from sources such as the World Health Organization, Our World in Data, and Worldometer. The parameters of the experiment are shown in Table 4.
Using effective strategies (see Figure 16), it was possible to reduce the maximum proportion infected with V 2 virus from 0.6 to 0.3 , but the proportion infected with V 1 virus increased from 0.34 to 0.38 , which is due to the fact that the proportion of people (30% of the population) who were not infected with V 2 virus (due to the measures introduced) managed to become infected with V 1 virus. It can also be seen that from t = 15 , the proportion infected with both viruses does not exceed 1%.
By running the computer simulation program, effective controls for the formulated problem were obtained (see Figure 17). It can be observed that the treatment of individuals for both viruses takes a maximum value starting from t = 1 and decreases for virus V 1 after t = 21 and for virus V 1 after t = 18 . For V 1 , a small fraction ( u 1 2 = 0.1 ) of those infected during [ 5 , 8 ] are isolated, and after t = 15 the control u 1 2 turns on again, reaching a maximum value, and turns off after t = 24 . For V 2 , maximum isolation of the infected occurs at [ 0 , 4 ] , [ 5 , 8 ] , [ 12 , 21 ] , and partial isolation at [ 9 , 11 ] , [ 29 , 30 ] . Reducing the total number of contacts was most effective in periods [ 15 , 25 ] .
For problems with a modified infection rate parameter, no optimality conditions have yet been found, so the only way to assess the quality of the solution is the cost savings provided by the constructed efficient control. For this reason, the values of the target functional and the total control for the two viruses were calculated for this experiment (see Figure 18). With the combined effect of all controls it was possible to significantly reduce the value of the target function from 33,286 to 11,180 m.u. The benefit of using effective controls was therefore 22,106 m.u.
Experiment 6. The SWIIRS model was used for Experiment 6, and the main objective was to construct long term strategies ( T 90 ) using multiple MPC controllers. The experimental parameter values are available in Table 5.
Figure 19 and Figure 20 show the results of Experiment 6. The effective control plots show that when several MPC controllers are connected, there are pulses (switching points) in the control. They are related to the fact that the prediction model is recalculated to a new prediction interval. You may notice that from the third controller ( t = 60 ) the constructed controllers have a similar structure. This is reflected both in the behavior of the process and in the behavior of the total cost function (they become periodic).
By applying efficient management, the value of the target function was significantly reduced from 2918.8 to 177.8 m.u. Thus, the benefit of applying management was 2748 monetary units. It can be seen that model-predictive-control solutions remain effective over long time intervals, but require the connection of several MPC controllers.

5. Conclusions

In this paper, we have investigated the application of MPC controllers to an extended SIR model with two different viruses. Model predictive control is one of the ways to achieve efficient control in dynamically controlled systems. We linearized the system of differential equations describing the dynamics of the viruses and showed that as the linearization step is reduced, the accuracy of the solution increases. We carried out a series of numerical experiments in which we obtained effective strategies whose form is close to the optimal strategies obtained by Pontryagin’s maximum principle. The experiments also showed that the effective control has a similar form for different numbers of contacts. It can be seen that efficient control is much simpler and faster than optimal control and uses fewer resources. This method can also be used to assess the situation during the spread of epidemics, when time is of the essence.

Author Contributions

Conceptualization, N.K., V.T. and E.G.; methodology, N.K.; software, N.K.; validation, V.T., E.G. and N.K.; formal analysis, V.T.; investigation, N.K.; resources, E.G.; data curation, E.G.; writing—original draft preparation, N.K.; writing—review and editing, V.T. and E.G.; visualization, N.K.; supervision, E.G.; project administration, E.G.; funding acquisition, E.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by RFBR and DFG, project number 21-51-12007.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviation is used in this manuscript:
MPCModel Predictive Control

References

  1. Powell, M. Learning lessons from the COVID-19 pandemic. Health Econ. Policy Law 2022, 18, 1–16. [Google Scholar] [CrossRef]
  2. Roser, M.; Ritchie, H.; Ortiz-Ospina, E.; Hasell, J. Coronavirus Pandemic (COVID-19). Published Online at OurWorldInData.org. 2020. Available online: https://ourworldindata.org/coronavirus (accessed on 23 January 2023).
  3. Behncke, H. Optimal control of deterministic epidemics. Optim. Control. Appl. Methods 2000, 21, 269–285. [Google Scholar] [CrossRef]
  4. Khouzani, M.H.R.; Sarkar, S.; Altman, E. Maximum Damage Malware Attack in Mobile Wireless Networks. In Proceedings of the IEEE INFOCOM 2010, San Diego, CA, USA, 14–19 March 2010. [Google Scholar]
  5. Taynitskiy, V.; Gubar, E.; Zhu, Q. Optimal Control of Heterogeneous Mutating Viruses. Games 2018, 9, 103. [Google Scholar] [CrossRef]
  6. Gubar, E.; Quanyan, Z. Optimal Control of Influenza Epidemic Model with Virus Mutations. In Proceedings of the European Control Conference (ECC), Zurich, Switzerland, 17–19 July 2013. [Google Scholar] [CrossRef]
  7. Biegler, L.T. A Perspective on Nonlinear Model Predictive Control; Chemical Engineering Department, Carnegie Mellon University Pittsburgh: Pittsburgh, PA, USA, 2021. [Google Scholar] [CrossRef]
  8. Arkun, Y. Design of Steady-State Optimizing Control Structures for Chemical Processes; ProQuest Dissertations Publishing: Ann Arbor, MI, USA, 1979; pp. 76–80. [Google Scholar]
  9. Ricker, N.L.; Subrahmanian, T.; Sim, T. Case studies of model-predictive control in pulp and paper production. IFAC Proc. Vol. 1988, 21, 13–22. [Google Scholar] [CrossRef]
  10. Culter, C.R.; Huwkins, R.B. Constrained Multivariable Control of a Hydrocracker Reactor. In Proceedings of the 1987 American Control Conference, Minneapolis, MN, USA, 10–12 June 1987. [Google Scholar] [CrossRef]
  11. Sidorov, D. Preface to “Model Predictive Control and Optimization for Cyber-Physical Systems”. Mathematics 2023, 11, 1004. [Google Scholar] [CrossRef]
  12. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. 1927, 115, 700–721. [Google Scholar]
  13. Altman, E.; Khouzani, M.; Sarkar, S. Optimal control of epidemic evolution. In Proceedings of the 2011 Proceedings IEEE INFOCOM, Shanghai, China, 10–15 April 2011; pp. 1683–1691. [Google Scholar] [CrossRef]
  14. Gubar, E.; Policardo, L.; Sánchez Carrera, E.J.; Taynitskiy, V. On optimal lockdown policies while facing socioeconomic costs. Ann. Oper. Res. 2023. [Google Scholar] [CrossRef]
  15. Kissler, S.M.; Tedijanto, C.; Goldstein, E.; Grad, Y.H.; Lipsitch, M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science 2020, 368, 860–868. [Google Scholar] [CrossRef] [PubMed]
  16. Perkins, T.A.; España, G. Optimal Control of the COVID-19 Pandemic with Non-pharmaceutical Interventions. Bull. Math. Biol. 2020, 82, 118. [Google Scholar] [CrossRef] [PubMed]
  17. Gubar, E.; Taynitskiy, V.; Fedyanin, D.; Petrov, I.; Quanyan, Z. Optimal Control of Joint Multi-Virus Infection and Information Spreading. IFAC-PapersOnLine 2020, 53, 6650–6655. [Google Scholar]
  18. Li, L.; Littman, M.; Walsh, T. Knows What It Knows: A Framework for Self-Aware Learning; Department of Computer Science, Rutgers University: Piscataway, NJ, USA, 2010. [Google Scholar] [CrossRef]
  19. Kolodin, A.A.; Elshin, V.V. Research and Development of the Controller Based on the Model Predictive Control; Irkutsk National Research Technical University: Irkutsk, Russia, 2023. [Google Scholar] [CrossRef]
  20. Allgower, F.; Findeisen, R.; Nagy, Z.K. Nonlinear Model Predictive Control: From Theory to Application; Institute for Systems Theory in Engineering, University of Stuttgart: Stuttgart, Germany, 2004. [Google Scholar]
  21. MPC Controller. Available online: https://www.mathworks.com/help/mpc/ref/mpccontroller.html (accessed on 24 January 2023).
  22. Nonlinear MPC. Available online: https://www.mathworks.com/help/mpc/ug/nonlinear-mpc.html (accessed on 3 February 2023).
  23. Coronavirus (COVID-19) Vaccinations. Available online: https://ourworldindata.org/covid-vaccinations (accessed on 7 April 2023).
  24. SARS-CoV-2 B.1.1.529 (Omicron) Variant Transmission within Households— Four U.S. Jurisdictions, November 2021–February 2022. Available online: https://www.cdc.gov/mmwr/volumes/71/wr/mm7109e1.htm (accessed on 8 April 2023).
Figure 1. Diagram of the two virus propagation model.
Figure 1. Diagram of the two virus propagation model.
Computation 11 00173 g001
Figure 2. SIIRS model diagram.
Figure 2. SIIRS model diagram.
Computation 11 00173 g002
Figure 3. Diagram of the SWIIRS model.
Figure 3. Diagram of the SWIIRS model.
Computation 11 00173 g003
Figure 4. The graph of outcomes.
Figure 4. The graph of outcomes.
Computation 11 00173 g004
Figure 5. Diagram for Model-Predictive-Control algorithm.
Figure 5. Diagram for Model-Predictive-Control algorithm.
Computation 11 00173 g005
Figure 6. Diagram of the nonlinear MPC controller.
Figure 6. Diagram of the nonlinear MPC controller.
Computation 11 00173 g006
Figure 7. Experiment 1. Behavior of SIIR model with the nonlinear-MPC controller. Average number of contacts per day is equal to 1.
Figure 7. Experiment 1. Behavior of SIIR model with the nonlinear-MPC controller. Average number of contacts per day is equal to 1.
Computation 11 00173 g007
Figure 8. Experiment 1. Structure of the effective control and aggregated system costs in uncontrolled and controlled cases.
Figure 8. Experiment 1. Structure of the effective control and aggregated system costs in uncontrolled and controlled cases.
Computation 11 00173 g008
Figure 9. Experiment 2. Behavior of SIIR model with the nonlinear MPC controller. Average number of contacts per day is equal to 3.
Figure 9. Experiment 2. Behavior of SIIR model with the nonlinear MPC controller. Average number of contacts per day is equal to 3.
Computation 11 00173 g009
Figure 10. Experiment 2. Structure of the effective control and aggregated system costs in uncontrolled and controlled cases.
Figure 10. Experiment 2. Structure of the effective control and aggregated system costs in uncontrolled and controlled cases.
Computation 11 00173 g010
Figure 11. Experiment 3. Behavior of SIIR model with the nonlinear MPC controller. Average number of contacts per day is equal to 10.
Figure 11. Experiment 3. Behavior of SIIR model with the nonlinear MPC controller. Average number of contacts per day is equal to 10.
Computation 11 00173 g011
Figure 12. Experiment 3. Structure of the effective control and aggregated system costs in uncontrolled and controlled cases.
Figure 12. Experiment 3. Structure of the effective control and aggregated system costs in uncontrolled and controlled cases.
Computation 11 00173 g012
Figure 13. Experiment 4. Behavior of SIIR model with the MPC controller.
Figure 13. Experiment 4. Behavior of SIIR model with the MPC controller.
Computation 11 00173 g013
Figure 14. Experiment 4. Structure of the effective control and aggregated system costs in uncontrolled and controlled models.
Figure 14. Experiment 4. Structure of the effective control and aggregated system costs in uncontrolled and controlled models.
Computation 11 00173 g014
Figure 15. Experiment 4. Structure of the optimal control and comparison of the aggregated costs in controlled and uncontrolled cases.
Figure 15. Experiment 4. Structure of the optimal control and comparison of the aggregated costs in controlled and uncontrolled cases.
Computation 11 00173 g015
Figure 16. Experiment 5. Behavior of SIIRS model with the nonlinear-MPC controller.
Figure 16. Experiment 5. Behavior of SIIRS model with the nonlinear-MPC controller.
Computation 11 00173 g016
Figure 17. Experiment 5. Structure of the effective control.
Figure 17. Experiment 5. Structure of the effective control.
Computation 11 00173 g017
Figure 18. Experiment 5. Structure of the effective control for isolated individuals and aggregated system costs.
Figure 18. Experiment 5. Structure of the effective control for isolated individuals and aggregated system costs.
Computation 11 00173 g018
Figure 19. Experiment 6. Behavior of SWIIRS model with the nonlinear MPC controller.
Figure 19. Experiment 6. Behavior of SWIIRS model with the nonlinear MPC controller.
Computation 11 00173 g019
Figure 20. Experiment 6. Structure of the effective control and aggregated system costs.
Figure 20. Experiment 6. Structure of the effective control and aggregated system costs.
Computation 11 00173 g020
Table 1. Dependence of MPC controller accuracy and performance on k parameter.
Table 1. Dependence of MPC controller accuracy and performance on k parameter.
TestT kStandard CostsEffective CostsOptimal CostOperation Time
11806042648258035.31 s
218090427310178179.64 s
318012042788518399.17 s
4180180428286085615.00 s
518036042868668651 min 0.72 s
618045042868688671 min 58.87 s
718054042878708704 min 6.98 s
818063042878708706 min 20.63 s
9180720428787087010 min 16.43 s
Table 2. Values of system parameters and initial states for Experiments 1, 2, and 3.
Table 2. Values of system parameters and initial states for Experiments 1, 2, and 3.
ParameterValueCharacteristics
Dynamic System Settings
n4Number of system states
r4Number of system outputs
m2Number of control variables
Time parameters
T l 6Number of time intervals
T s 0.25Step of the controller
P H 30Prediction Horizon
C H 21Control Horizon
Parameters of the S I I R model
1 (Exp. 1)Average number of contacts per day
k3 (Exp. 2)(Exp. stands for Experiment)
10 (Exp. 3)
Virus 1 (SARS-CoV-2 Alpha)
δ 1 0.256Infection rate of Virus 1
σ 1 0.025Recovery rate of Virus 1
Virus 2 (SARS-CoV-2 Omicron)
δ 2 0.517Infection rate of Virus 2
σ 2 0.07Recovery rate of Virus 2
Cost function parameters
f 1 134.5Infection costs for virus 1
f 2 134.5Infection costs for virus 2
h 1 192.07Treatment costs for virus 1
h 2 192.07Treatment costs for virus 2
Table 3. Values of system parameters and initial states for Experiment 4.
Table 3. Values of system parameters and initial states for Experiment 4.
ParameterValueDescription
Characteristics of the virus
δ 1 0.4 Infection rate of V 1
δ 2 0.5 Infection rate of V 2
σ 1 0.001 Recovery rate from V 1
σ 2 0.002 Recovery rate from V 2
Cost functions
f 1 ( I 1 ) 5 I 1 Infection costs of V 1
f 2 ( I 2 ) 6 I 2 Infection costs of V 2
g ( R ) 0.1 R Benefit rate
h 1 ( u 1 ) 15 u 1 2 Cost of treatment u 1
h 2 ( u 2 ) 10 u 2 2 Cost of treatment u 2
Initial states
S ( 0 ) 0.50 Fraction of Susceptible
I 1 ( 0 ) 0.32 Fraction of Infected with V 1
I 2 ( 0 ) 0.18 Fraction of Infected with V 2
R ( 0 ) 0.00 Fraction of Recovered
Values of MPC controller
T l 1Number of time intervals
T s 0.5 Step of the controller
P H 90Prediction horizon
C H 90Control horizon
Table 4. Values of system parameters and initial states for Experiment 5.
Table 4. Values of system parameters and initial states for Experiment 5.
ParameterValueDescription
Characteristics of the virus
p 1 0.256 Infection probability with virus V 1 (single contact)
p 2 0.517 Infection probability with virus V 2 (single contact)
σ 1 0.025 Recovery rate from V 1
σ 2 0.007 Recovery rate from V 2
γ 1 / 120 0.0083 Decreasing immunity rate
l30Average number of contacts in the population
f 1 ( I 1 ) 135 I 1 Infection costs of V 1
f 2 ( I 2 ) 135 I 2 Infection costs of V 2
g ( R ) 0Benefit rate
h 1 1 ( u 1 1 ) 192 ( u 1 1 ) 2 Treatment cost u 1 1
h 1 2 ( u 1 2 , I 1 ) 94 u 1 2 I 1 Isolation cost u 1 2
h 2 1 ( u 2 1 ) 192 ( u 2 1 ) 2 Treatment cost u 2 1
h 2 2 ( u 2 2 , I 2 ) 94 u 2 2 I 2 Isolation cost u 2 2
h l ( u l ) 94 ( u l ) 2 Cost of reducing the number of contacts u l
Initial states
S ( 0 ) 0.8 Fraction of Susceptible
I 1 ( 0 ) 0.13 Fraction of Infected with V 1
I 2 ( 0 ) 0.07 Fraction of Infected with V 2
R ( 0 ) 0.00 Fraction of Recovered
Values of MPC controller
k1Number of controllers
T s 0.5 Step the controller
P H 45Forecast horizon
C H 30Control horizon
Table 5. Values of system parameters and initial states for Experiment 6.
Table 5. Values of system parameters and initial states for Experiment 6.
ParameterValueDescription
Characteristics of the virus
δ 1 S 0.25 Infection rate for Susceptible ( V 1 )
δ 2 S 0.3 Infection rate for Susceptible ( V 2 )
δ 1 W 0.2 Infection rate for Warned ( V 1 )
δ 2 W 0.25 Infection rate for Warned ( V 2 )
σ 1 0.3 Recovery rate from V 1
σ 2 0.4 Recovery rate from V 2
σ 3 0.3 Vaccination rate
γ 0.3 Decreasing immunity rate
η 0.15 Information spreading rate
Cost functions
f 1 ( I 1 ) 30 I 1 Infection costs of V 1
f 2 ( I 2 ) 40 I 2 Infection costs of V 2
g ( R ) 0Benefit rate
h 1 ( u 1 ) 20 u 1 2 Treatment cost u 1
h 2 ( u 2 ) 25 u 2 2 Treatment cost u 2
h 3 ( u 3 ) 10 u 3 2 Cost of spreading information u 3
Initial states
S ( 0 ) 0.85 Fraction of Susceptible
W ( 0 ) 0.00 Fraction of Warned
I 1 ( 0 ) 0.11 Fraction of Infected with V 1
I 2 ( 0 ) 0.04 Fraction of Infected with V 2
R ( 0 ) 0.00 Fraction of Recovered
Values of MPC controller
k6Number of controllers
T s 0.5 Step the controller
P H 30Forecast horizon
C H 20Control horizon
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kosyanov, N.; Gubar, E.; Taynitskiy, V. MPC Controllers in SIIR Epidemic Models. Computation 2023, 11, 173. https://doi.org/10.3390/computation11090173

AMA Style

Kosyanov N, Gubar E, Taynitskiy V. MPC Controllers in SIIR Epidemic Models. Computation. 2023; 11(9):173. https://doi.org/10.3390/computation11090173

Chicago/Turabian Style

Kosyanov, Nikita, Elena Gubar, and Vladislav Taynitskiy. 2023. "MPC Controllers in SIIR Epidemic Models" Computation 11, no. 9: 173. https://doi.org/10.3390/computation11090173

APA Style

Kosyanov, N., Gubar, E., & Taynitskiy, V. (2023). MPC Controllers in SIIR Epidemic Models. Computation, 11(9), 173. https://doi.org/10.3390/computation11090173

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop