Next Article in Journal
Estimating the Expected Discounted Penalty Function in a Compound Poisson Insurance Risk Model with Mixed Premium Income
Next Article in Special Issue
Optimal Impulse Vaccination Approach for an SIR Control Model with Short-Term Immunity
Previous Article in Journal
Extended Fuzzy Metrics and Fixed Point Theorems
Previous Article in Special Issue
Use of Enumerative Combinatorics for Proving the Applicability of an Asymptotic Stability Result on Discrete-Time SIS Epidemics in Complex Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Role of Media and Effects of Infodemics and Escapes in the Spatial Spread of Epidemics: A Stochastic Multi-Region Model with Optimal Control Approach

Department of Mathematics and Computer Sciences, Faculty of Sciences Ben M’Sik, Hassan II University of Casablanca, Casablanca 20000, Morocco
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(3), 304; https://doi.org/10.3390/math7030304
Submission received: 23 January 2019 / Revised: 19 February 2019 / Accepted: 6 March 2019 / Published: 25 March 2019
(This article belongs to the Special Issue Mathematical Models in Epidemiology )

Abstract

:
Mass vaccination campaigns play major roles in the war against epidemics. Such prevention strategies cannot always reach their goals significantly without the help of media and awareness campaigns used to prevent contacts between susceptible and infected people. Feelings of fear, infodemics, and misconception could lead to some fluctuations of such policies. In addition to the vaccination strategy, the movement restriction approach is essential because of the factor of mobility or travel. However, anti-epidemic border measures may also be disturbed if some infected travelers manage to escape and infiltrate into a safer region. In this paper, we aim to study infection dynamics related to the spatial spread of an epidemic in interconnected regions in the presence of random perturbations caused by the three above-mentioned reasons. Therefore, we devise a stochastic multi-region epidemic model in which contacts between susceptible and infected populations, vaccination-based and movement restriction optimal control approaches are all assumed to be unpredictable, and then, we discuss the effectiveness of such policies. In order to reach our goal, we employ a stochastic maximum principle version for noised systems, state and prove the sufficient and necessary conditions of optimality, and finally provide the numerical results obtained using a stochastic progressive-regressive schemes method.

1. Introduction

Media plays a tremendous role in mounting awareness among susceptible populations in an attempt to reduce their contact with infection. In fact, it has the potential of generating a psychological impact on the social conduct, as explained in [1,2,3]. Then, many modelers of epidemics saw it was also important to introduce and discuss the effect of awareness through media in the outbreaks of diseases; see studies in [4,5,6,7,8].
On the other hand, in times of recently serious epidemic outbreaks, infodemics or epidemics of rumors may appear very quickly in the virtual world. Because of contradictory views and unreliable and misleading information diffused by some Internet users, people become confused and fearful. In such circumstances, the role of media is essential, and journalists with scientists are obliged to report and exhibit concrete and convincing proof = about the nature of the epidemic and should explain the reasons for control interventions led by health authorities [9]. Here, we try to focus on rumors that can prevent a large portion of the population from being informed about the necessity of following an urgent anti-epidemic measure, namely vaccination.
Recently, multi-region epidemic models as in [10,11] have been interested in the study of the spatial spread of some epidemics using metapopulation-like differential equations. Other papers as in [12,13,14,15,16,17,18,19,20,21] chose to model the same phenomenon, namely the regional spread of epidemics, using the framework of difference equations. All these last references treated dynamics of epidemics in the presence of a so-called travel-blocking control strategy. In fact, the mobility factor is very important in such considered systems, and then, it is important to discuss the effectiveness of movement restrictions between regions; otherwise, a vaccination control policy alone would not seem sufficient, especially when we talk about an infection spreading in or around a large geographical territory. In order to discuss other missing considerations on this subject, we try to propose here a stochastic version of those epidemic modeling approaches based on a susceptible-infected-removed (SIR) multi-region stochastic model that describes the spatial-temporal spread of an epidemic in the presence of optimal vaccination and movement restriction strategies under perturbations. In the analysis of the optimization of these control interventions, we focus in this paper on providing existence results of the sought controls in a proposition and theorem of sufficient and necessary conditions, and thereafter, we discuss the numerical results.
In the last decade, many applied mathematicians have been interested in the study of infection dynamics under perturbations based on stochastic compartmental models; see for example [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37]. The main idea here is that contacts between susceptible and infected populations are unpredictable. As an example, media coverage alone could lead to random spread of an epidemic because some susceptible individuals would avoid meeting infected people at any time once they receive alerts. We investigate in this case the stochastic dynamics of infection when it occurs in regions that are neighbors and interconnected by any kind of anthropological movement. The model is devised here for the study of the spread of an epidemic in a domain Ω and that has the form of a stochastic control differential state equation written at a time t as
x ˙ Ω ( t ) = f ( t , x Ω ( t ) , u Ω ( t ) ) + g ( t , x Ω ( t ) , u Ω ( t ) ) d W Ω ( t ) d t
with t [ 0 , T ] .
In the following sections, we define our stochastic multi-region SIR model and apply thereafter a stochastic maximum principle for characterizing the sought optimal control functions and that is associated with the mass vaccination strategy and movement restriction policies.

2. Model Description and Definitions

Presentation of the Stochastic Model without Control

First of all and based on the assumptions of the deterministic modeling approach proposed in [10,11], we assume that there are p geographical regions denoted Ω j (domains) of the domain studied. Ω = p j = 1 Ω j . Let N Ω j ( t ) be the population of domain Ω j at time t, presenting the number of individuals who are physically present in Ω j , both residents and travelers, and let the host population of Ω j be grouped into three epidemiological compartments. Let S Ω j ( t ) , I Ω j ( t ) , and R Ω j ( t ) be the number of individuals in the susceptible, infective, and removed compartments of Ω j at time t, respectively.
The stochastic disease transmission in a given domain Ω j at time t is modeled using a perturbed standard incidence, which we present by:
k = 1 p ρ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t )
where the stochastic disease transmission coefficient ρ j k ( t ) is the stochastic proportion of adequate contacts in domain Ω j between a susceptible from Ω j ( j = 1 , , p ) and an infective from another domain Ω k at a time t, and which we define by:
ρ j k ( t ) = β j k + σ j d W Ω j ( t ) d t
where β j k > 0 is the equivalent disease transmission coefficient to ρ j k , but in the deterministic case, which had been defined also as the proportion of adequate contacts in the study cases of [10,11], σ j ( j = 1 , , p ) are real constants and represent the intensities of fluctuations caused by media, and { W Ω j ( t ) } t [ 0 , T ] is an independent random variable composed of continuous white noises independent of F t F and that is a standard Brownian motion supposed to be caused in this part without control, by media coverage diffused in region Ω j .
When there is no control introduced yet in Ω j , the stochastic multi-regional continuous-time SIR model associated with Ω j is presented as follows:
S ˙ Ω j ( t ) = - k = 1 p ρ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + N Ω j ( t ) - S Ω j ( t ) d j
I ˙ Ω j ( t ) = k = 1 p ρ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) - γ j I Ω j ( t ) - d j I Ω j ( t )
R ˙ Ω j ( t ) = γ j I Ω j ( t ) - d j R Ω j ( t )
where d j is the birth and death rate and γ j is the recovery rate. The biological background requires that all parameters be non-negative.
N Ω j ( t ) = S Ω j ( t ) + I Ω j ( t ) + R Ω j ( t ) is the population size corresponding to domain Ω j at time t. The population size remains constant for all t [ 0 , T ] ; in fact:
N ˙ Ω j ( t ) = S ˙ Ω j ( t ) + I ˙ Ω j ( t ) + R ˙ Ω j ( t ) = 0
Therefore, as a function of the deterministic proportion of adequate contacts β j k and the continuous Wiener process W Ω j ( t ) , the stochastic system (6) becomes the Itô multi-region stochastic differential equations (SDEs) model:
S ˙ Ω j ( t ) = - k = 1 p β j k I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + N Ω j ( t ) - S Ω j ( t ) d j - σ j k = 1 p I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) d W Ω j ( t ) d t
I ˙ Ω j ( t ) = k = 1 p β j k I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) - γ j I Ω j ( t ) - d j I Ω j ( t ) + σ j k = 1 p I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) d W Ω j ( t ) d t
R ˙ Ω j ( t ) = γ j I Ω j ( t ) - d j R Ω j ( t )

3. The Model with Vaccination

3.1. Presentation of the Control Model

In this section, we introduce a function θ Ω j ( t ) as a perturbation of the control variable denoted by u Ω j ( t ) and that characterizes the effectiveness of vaccination in the above-mentioned model (4)–(6) as in [11]. This perturbation with disturbances can be caused by infodemics and rumors and also ideas from feelings of fear and misconception. Then, for a given domain Ω j targeted by vaccination, the model is given by the following equations:
S ˙ Ω j ( t ) = - k = 1 p ρ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + N Ω j ( t ) - S Ω j ( t ) d j - θ Ω j ( t ) S Ω j ( t )
I ˙ Ω j ( t ) = k = 1 p ρ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) - γ j I Ω j ( t ) - d j I Ω j ( t )
R ˙ Ω j ( t ) = γ j I Ω j ( t ) - d j R Ω j ( t ) + θ Ω j ( t ) S Ω j ( t )
with:
θ Ω j ( t ) = u Ω j ( t ) + δ j d W Ω j ( t ) d t .
Thus, in a more general form, it refers to the stochastic control equation written at a time t as
x ˙ Ω ( t ) = f ( t , x Ω ( t ) , u Ω ( t ) ) + g ( t , x Ω ( t ) , u Ω ( t ) ) d W Ω ( t ) d t
where at time t and for j = 1 , , p :
x Ω ( t ) = x Ω j ( t ) = S Ω j ( t ) I Ω j ( t ) R Ω j ( t )
u Ω ( t ) = u Ω j ( t ) ,
f ( t , x Ω ( t ) , u Ω ( t ) )
= - k = 1 p β j k I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + N Ω j ( t ) - S Ω j ( t ) d j - u Ω j ( t ) S Ω j ( t ) k = 1 p β j k I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) - γ j I Ω j ( t ) - d j I Ω j ( t ) γ j I Ω j ( t ) - d j R Ω j ( t ) + u Ω j ( t ) S Ω j ( t )
and:
g ( t , x Ω ( t ) , u Ω ( t ) ) = - σ j k = 1 p I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + δ j S Ω j ( t ) σ j k = 1 p I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) δ j S Ω j ( t )
Our goal is to try to minimize the population of the infected group and the cost of vaccination, while maximizing the number of removed people in all regions. Our control functions take values between u m i n Ω j and u m a x Ω j , where u m i n Ω k , u m a x Ω k 0 , 1 , k = 1 , , p .

3.2. A Stochastic Optimal Control Approach

3.2.1. Optimal Control Characterization and Necessary Conditions

We devise in this paper an optimal control approach that aims to minimize the number of the infected people and maximize the ones in the removed category for all regions, while minimizing the cost of vaccination.
The effort expended in preventing the epidemic in each region Ω j is proportional to u Ω j , i.e., achieving a higher value of control u Ω j means that more money, equipment, personnel, and resources must be brought to Ω j . This implies that the cost of vaccination in region Ω j , lets say a function C ( u Ω j ) , may not be linear; that is, doubling u Ω j may require more than double the expenditure [38].
Then, we are interested in minimizing the functional
J ( u Ω j ) = IE 0 T f 0 ( t , x Ω j ( t ) , u Ω j ( t ) ) d t
with:
f 0 ( t , x Ω j ( t ) , u Ω j ( t ) ) = α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2
where A j > 0 , α j I > 0 , α j R > 0 are the weight constants of control, the infected, and the removed in region Ω j , respectively.
In [11], the authors studied the special case of the minimization problem of the cost functional J when there were no perturbations.
Here, our goal is to minimize the number of infected people and minimize the systemic costs attempting to increase the number of removed people in each Ω j . In other words, we are seeking an optimal control u Ω j * such that:
J ( u Ω j * ) = min { J ( u Ω j ) / u Ω j U j }
where U j is the control set defined by:
U j ( [ 0 , T ] ) = { u Ω j ( t ) F t - p r o g r e s s i v e l y m e a s u r a b l e | u m i n Ω j u Ω j ( t ) u m a x Ω j , t [ 0 , T ] }
Let us define the Hamiltonian function H by:
H ( x Ω , u Ω , μ Ω , ν Ω ) = f 0 ( x Ω , u Ω ) + f ( x Ω , u Ω ) , μ Ω + t r ν Ω T g ( x Ω , u Ω )
At time t [ 0 , T ] and for j = 1 , . . . , p , it can be rewritten as:
H ( x Ω j ( t ) , u Ω j ( t ) , μ Ω j ( t ) , ν Ω j ( t ) ) = f 0 ( x Ω j ( t ) , u Ω j ( t ) ) + f ( x Ω j ( t ) , u Ω j ( t ) ) , μ Ω j ( t ) + l = 1 3 g l T ( x Ω j ( t ) , u Ω j ( t ) ) ν Ω j l ( t )
Here, . T means the transposition, while in a domain Ω , ( μ ( t ) , ν ( t ) ) ) is a pair of adjoint variables satisfying the following adjoint BSDE (backward stochastic differential equation):
d μ Ω ( t ) = - [ f x T ( t , x Ω ( t ) , u Ω ( t ) ) μ Ω ( t ) + l = 1 3 g x Ω l T ( t , x Ω ( t ) , u Ω ( t ) ) ν Ω l ( t ) + f 0 x Ω ( t , x Ω ( t ) , u Ω ( t ) ) ] d t + ν Ω ( t ) d W Ω ( t ) , μ Ω ( T ) = 0 .
Using a stochastic version of Pontryagin’s maximum principle [39], we characterize the optimal control u in the following theorem to find its analytical formulation.
Theorem 1. 
(Stochastic maximum principle and characterization of u Ω * )
If there exists an optimal pair ( x Ω * , u Ω * ) and a pair of processes ( μ ( t ) , ν ( t ) ) satisfying (27), then for j = 1 , . . . , p , we have:
H ( x Ω j * ( t ) , u Ω j * ( t ) , μ Ω j ( t ) , ν Ω j ( t ) ) = min u Ω j U H ( x Ω j ( t ) , u Ω j * ( t ) , μ Ω j ( t ) , ν Ω j ( t ) ) .
Moreover, we obtain the bounded stochastic control:
u Ω j * = m i n ( m a x ( u m i n Ω j , - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j * A j ) , u m a x Ω j ) ,
the solution of the FBSDEs (forward-backward stochastic differential equations):
d x Ω j ( t ) = f ( t , x Ω j ( t ) , u Ω j ( t ) ) d t + g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) d μ Ω j ( t ) = - [ f x Ω j T ( t , x Ω j ( t ) , u Ω j ( t ) ) μ Ω j ( t ) + l = 1 3 g x Ω j l T ( t , x Ω j ( t ) , u Ω j ( t ) ) ν Ω j l ( t ) + f 0 x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) ] d t + ν Ω j ( t ) d W Ω j ( t ) , x Ω j ( 0 ) = ( S 0 Ω j , I 0 Ω j , R 0 Ω j ) μ Ω j ( T ) = 0 .
Proof. 
Since our control u Ω j is bounded, we then prove the previous theorem by using the following Lagrangian:
L ( x Ω j ( t ) , u Ω j ( t ) , μ Ω j ( t ) , ν Ω j ( t ) , ω 1 ( t ) , ω 2 ( t ) ) = α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 + μ Ω j T ( t ) f ( t , x Ω j ( t ) , u Ω j ( t ) ) + l = 1 3 g l T ( x Ω j ( t ) , u Ω j ( t ) ) ν Ω j l ( t ) + ω 1 ( t ) ( u m a x Ω j - u Ω j ( t ) ) + ω 2 ( t ) ( u Ω j ( t ) - u m i n Ω j )
where ω 1 , ω 2 0 , verifying at u Ω j = u Ω j * the two conditions:
ω 1 ( u m a x Ω j - u Ω j * ) = 0 and ω 2 ( u Ω j * - u m i n Ω j ) = 0 .
Owing to the condition of minimization, we define by:
L ( x Ω j * ( t ) , u Ω j * ( t ) , μ Ω j ( t ) , ν Ω j ( t ) , ω 1 ( t ) , ω 2 ( t ) ) = min u Ω j U j L ( x Ω j ( t ) , u Ω j ( t ) , μ Ω j ( t ) , ν Ω j ( t ) , ω 1 ( t ) , ω 2 ( t ) ) .
We differentiate the Lagrangian with respect to u Ω j on the set:
{ t | u m i n Ω j u Ω j ( t ) u m a x Ω j }
to obtain the optimality equation
d L d u Ω j ( x Ω j ( t ) , u Ω j ( t ) , μ Ω j ( t ) , ν Ω j ( t ) , ω 1 ( t ) , ω 2 ( t ) ) | u Ω j
= u Ω j * = A j u Ω j ( t ) + ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j - ω 1 ( t ) + ω 2 ( t ) = 0 .
Furthermore, we find u Ω j * ( t ) = - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j - ω 1 + ω 2 A j .
-
If u m i n Ω j < u Ω j * ( t ) < u m a x Ω j , then w 1 ( t ) = w 2 ( t ) = 0 ; therefore:
u Ω j * ( t ) = - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j .
-
If u Ω j * ( t ) = u m i n Ω j , then w 1 ( t ) = 0 ;, therefore, u m i n Ω j = - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j + w 2 ( t ) A j , implying that w 2 ( t ) = - ( ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j + A j u m i n Ω j ) .
Due to w 2 ( t ) 0 and A j > 0 , we obtain u Ω j * ( t ) - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j
-
If u Ω j * ( t ) = u m a x Ω j , then w 2 ( t ) = 0 ; thus, u m a x Ω j = - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j - w 1 ( t ) A j implying that w 1 ( t ) = ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j + A j u m a x Ω j .
In view of w 1 ( t ) 0 and A j > 0 , we get u Ω j * ( t ) - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j .
Using these standard optimality arguments, we characterize the control u Ω j * ( t ) by
u Ω j * ( t ) = - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j i f u m i n Ω j < - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j < u m a x Ω j u m i n Ω j i f - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j u m i n Ω j u m a x Ω j i f - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j A j u m a x Ω j or by a more reduced form, we can rewrite u Ω j * ( t ) = m i n ( m a x ( u m i n Ω j , - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j * A j ) , u m a x Ω j ) ).  □

3.2.2. Existence of Solutions and Sufficient Conditions

Note that f : [ 0 , T ] × R 3 × U j R 3 , g : [ 0 , T ] × R 3 × U j R 3 R 3 and f 0 : [ 0 , T ] × R 3 × U j R are measurable such that f ( t , x , . ) : U j R 3 and f 0 ( t , x , . ) : U j R 3 are continuous, f, g, and f 0 are bounded, and there exists a constant K > 0 such that for all t [ 0 , T ] and for all x Ω j , x ^ Ω j R 3 , the following properties are checked [40,41]:
| f ( t , x Ω j ( t ) , u Ω j ( t ) ) - f ( t , x ^ Ω j ( t ) , u Ω j ( t ) ) | + | | g ( t , x Ω j ( t ) , u Ω j ) - g ( t , x ^ Ω j ( t ) , u Ω j ) | | K | x Ω j - x ^ Ω j |
In fact, if we suppose x Ω j = ( S Ω j , I Ω j , R Ω j ) , x ^ Ω j = ( S ^ Ω j , I ^ Ω j , R ^ Ω j ) Γ R 3 with Γ = ] a 1 , a 2 [ × ] b 1 , b 2 [ × ] c 1 , c 2 [ ( a i = 1 , 2 , b i = 1 , 2 and c i = 1 , 2 are positive constants). In order to prove (12), we can take constants
K 11 = δ j + σ j p b 2 a 1 + b 1 + c 1 ,
K 21 = σ j p a 2 a 1 + b 1 + c 1 ,
K 31 = δ j and K 1 = max ( K 11 , K 21 , K 31 ) to obtain:
| | g ( t , x Ω j ( t ) , u Ω j ( t ) ) - g ( t , x ^ Ω j ( t ) , u Ω j ( t ) ) | | K 1 | x Ω j - x ^ Ω j |
while taking also
K 2 = u m a x Ω j + ( b 2 + c 2 ) d j + k = 1 p β j k b 2 a 1 + b 1 + c 1 ,
K 3 = d j + γ j + k = 1 p β j k a 2 a 1 + b 1 + c 1
and K 4 = u m a x Ω j + d j + b 2 γ j ,
Thus, using the maximum norm, we have
| f 1 ( t , x Ω j ( t ) , u Ω j ( t ) ) - f 1 ( t , x ^ Ω j ( t ) , u Ω j ( t ) ) | K 2 | x Ω j - x ^ Ω j | ,
| f 2 ( t , x Ω j ( t ) , u Ω j ( t ) ) - f 2 ( t , x ^ Ω j ( t ) , u Ω j ( t ) ) | K 3 | x Ω j - x ^ Ω j | ,
and | f 3 ( t , x Ω j ( t ) , u Ω j ( t ) ) - f 3 ( t , x Ω j ( t ) , u Ω j ( t ) ) | K 4 | x Ω j - x ^ Ω j |
with f 1 , f 2 , and f 3 , the right-hand sides of differential equations in (24), (25), and (26). Thus, by setting f = ( f 1 , f 2 , f 3 ) , the property (12) is checked by considering
K = K 1 + max ( K 2 , K 3 , K 4 ) ,
and we should also have:
| f 0 ( t , x Ω j ( t ) , u Ω j ( t ) ) - f 0 ( t , x ^ Ω j ( t ) , u Ω j ( t ) ) | K | x Ω j - x ^ Ω j | .
which can easily be checked for the particular integrand supposed here.
Proposition 1.
The optimal control problem identified by the objective functional J and linked to the state system (24)–(26) that satisfies the properties (12) and (13) admits an optimal control pair.
Proof. 
At first, a backward stochastic differential equation (BSDE) with a terminal condition is introduced:
d Y ( t ) = - f 0 ( t , x Ω j ( t ) , u Ω j ( t ) ) d t + Z ( t ) d W Ω j ( t ) , t [ 0 , T ] Y ( T ) = 0 g i v e n ;
where:
f 0 ( t , x Ω j ( t ) , u Ω j ) = α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2
Note that under the previous observations, (12) and (13), for ( x 0 Ω j , u Ω j ) R 3 × U , the backward stochastic differential Equation (14) becomes linear. Therefore, the state system admits a unique strong solution [39] that is written in the following form x Ω j ( . ) x Ω j ( . ; x 0 Ω j , u Ω j ) .
Moreover, for a given ( x Ω j ( . ) , u Ω j ) , the backward stochastic differential equation (14) admits a unique adapted solution ( Y ( . ) , Z ( . ) ) ( Y ( . ; x 0 Ω j , u Ω j ) , Z ( . ; x 0 Ω j , u Ω j ) ) , which depends on ( x 0 Ω j , u Ω j ) through ( x Ω j ( . ) , u Ω j ) .
Notice that for t [ 0 , T ] , the process Y ( t ) , the solution of (14) has the general form:
Y ( t ) = Y ( T ) + t T f 0 ( t , x Ω j ( t ) , u Ω j ( t ) ) d t - t T Z ( t ) d W Ω j ( t ) ,
Denote that F t is the natural filtration of the Brownian motion W ( t ) . Here, Y ( t ) is ( F t ) t 0 -adapted [39], implying that:
Y ( t ) = E [ Y ( t ) | F t ] , t [ 0 , T ] .
Thus,
Y ( 0 ) = J ( u Ω j ) = E [ 0 T ( α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 ) d t ]
Now, the objective functional J can be rewritten as:
J ( u Ω j ) = Y ( 0 ; x 0 Ω j , u Ω j )
Considering the new formulation of the objective function (18), a forward backward stochastic differential equations (FBSDE) system is introduced:
d x Ω j ( t ) = ( f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ) d t + g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) , d Y ( t ) = - ( α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 ) d t + Z ( t ) d W Ω j ( t ) , x Ω j ( 0 ) = x 0 Ω j , Y ( T ) = 0 given .
In the following, we will use an appropriate approach in order to use a parabolic formulation and to establish the existence result.
In fact, if ( x Ω j , Y , Z ) is an adapted solution of (19), then there exists an appropriate function θ such that the following relationship is verified [42]:
Y ( t ) = θ ( t , x Ω j ( t ) , u Ω j ( t ) ) , t [ 0 , T ] a . s . P ,
where the undetermined function θ is assumed to be of C 1 , 2 , 0 ( [ 0 , T ] × R 3 × U j ) .
The control u Ω j is sought to minimize the cost functional Y ( 0 ) = J ( u Ω j ) . From (19), we have:
d Y ( t ) = - ( α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 ) d t + Z ( t ) d W Ω j ( t ) .
To apply Ito’s formula to θ ( t , x Ω j ( t ) , u Ω j ( t ) ) , start by using Taylor’s polynomial:
d θ ( t , x Ω j ( t ) , u Ω j ( t ) ) = θ t ( t , x Ω j ( t ) , u Ω j ( t ) ) d t + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) d x Ω j + 1 2 θ x Ω j x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) ( d x Ω j ) 2 ,
and recall that:
d x Ω j ( t ) = ( f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ) d t + g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) .
Replacing the formulation of d x Ω j ( t ) in d θ ( t ) :
d θ ( t , x Ω j ( t ) , u Ω j ( t ) ) = θ t ( t , x Ω j ( t ) , u Ω j ( t ) ) d t + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) [ f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ] d t + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) + 1 2 θ x Ω j x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) [ ( ( f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ) d t ) 2 + ( g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) ) 2 + 2 g ( t , x Ω j ( t ) , u Ω j ( t ) ) ( f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ) d t d W Ω j ( t ) ] ,
and using Itô’s multiplication, to obtain the final formulation of d θ ( t ) :
d θ ( t , x Ω j ( t ) , u Ω j ( t ) ) = [ θ t ( t , x Ω j ( t ) , u Ω j ( t ) ) + 1 2 θ x Ω j x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) g 2 ( t , x Ω j ( t ) , u Ω j ( t ) ) + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) ( f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ) ] d t + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) .
According to (19), by equating the corresponding drift terms d Y ( t ) and d θ ( t ) , the following parabolic PDE is obtained:
d θ ( t , x Ω j ( t ) , u Ω j ( t ) ) = [ θ t ( t , x Ω j ( t ) , u Ω j ( t ) ) + 1 2 θ x Ω j x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) g 2 ( t , x Ω j ( t ) , u Ω j ( t ) ) + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) ( f ^ ( t , x Ω j ( t ) ) + f 1 ( x Ω j ( t ) ) u Ω j ( t ) ] d t + θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) g ( t , x Ω j ( t ) , u Ω j ( t ) ) d W Ω j ( t ) . = - ( α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 ) d t + Z ( t ) d W Ω j ( t ) .
Substitute Z ( t ) = θ x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) ) g ( t , x Ω j ( t ) , u Ω j ( t ) ) into the above parabolic PDE to obtain:
0 = θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j ( f ^ + f 1 u Ω j ) + α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 = θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) + f 1 θ x Ω j u Ω j + A j 2 ( u Ω j ( t ) ) 2 = θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) + f 1 θ x Ω j u Ω j + A j 2 ( u Ω j ( t ) ) 2 + 1 2 A j | f 1 θ x Ω j | 2 - 1 2 A j | f 1 θ x Ω j | 2 = θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) - 1 2 A j | f 1 θ x Ω j | 2 + A j 2 [ u Ω j 2 + | f 1 θ x Ω j | 2 A j 2 + 2 f 1 θ x Ω j u Ω j A j ] = θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 [ u Ω j + f 1 θ x Ω j A j ] 2 - 1 2 A j | f 1 θ x Ω j | 2 .
By rearranging terms, the above equation can be written in the form of a backward parabolic PDE with a terminal condition on θ :
θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) - 1 2 A j | f 1 T θ x Ω j | 2 + A j 2 [ u Ω j + f 1 T θ x Ω j A j ] 2 = 0 θ ( T , x Ω j ( T ) ) = 0 , x Ω j R 3 .
Using a standard maximum principle for parabolic partial differential equations (PPDE) [43], the smallest θ solution should be the one of the following PPDE problem:
θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) - 1 2 A j | f 1 T θ x Ω j | 2 = 0 θ ( T , x Ω j ( T ) ) = 0 , x Ω j R 3 .
For the resolution of (22), results from Ladyzenskaja et al. [43] show that there exists a unique classical solution for comparable types of parabolic partial differential equations. The discussion of the solvability of this type of quasi-linear parabolic PDE in the presence of a quadratic gradient term can be found in [44]. Thus, there exists a unique classical solution to this equation. In this case, the control u Ω j is sought in order to minimize the objective function J such that:
min u Ω j U j J ( u Ω j ) = J ( u Ω j * ) ,
Thus:
min u Ω j U j θ ( 0 ; x Ω j ( 0 ) , u Ω j ) = θ ( 0 ; x Ω j ( 0 ) , u Ω j * ) ,
Then, use (15), (20), and (21) to obtain:
0 T f 0 ( t , x Ω j ( t ) , u Ω j ( t ) ) d t = θ ( 0 , x Ω j ( 0 ) ) - θ ( T , x Ω j ( T ) ) + 0 T Z ( t ) d W Ω j ( t ) + θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) - 1 2 A j | f 1 T θ x Ω j | 2 + A j 2 [ u Ω j + f 1 T θ x Ω j A j ] 2 .
Therefore,
J ( u Ω j ) = θ ( 0 , x Ω j ( 0 ) ) - E { θ ( T , x Ω j ( T ) ) } + E 0 T θ t + 1 2 θ x Ω j x Ω j g 2 + θ x Ω j f ^ + α j I I Ω j ( t ) - α j R R Ω j ( t ) - 1 2 A j | f 1 T θ x Ω j | 2 + A j 2 [ u Ω j + f 1 T θ x Ω j A j ] 2 d t = θ ( 0 , x Ω j ( 0 ) , u Ω j ) + A j 2 E 0 T [ u Ω j + f 1 T θ x Ω j A j ] 2 d t θ ( 0 , x Ω j ( 0 ) , u Ω j * ) = J ( u Ω j * ) .
Consequently, it is finally concluded that u Ω j * = - f 1 T θ x Ω j is an optimal control.

3.2.3. Numerical Results

In this part, we suppose we have three interconnected regions, and we investigate numerically using the stochastic progressive-regressive schemes method presented in [45] the effectiveness of the optimal control approach presented above. In our simulations, we note that:
-
We study three regions, denoted by Ω 1 , Ω 2 , and Ω 3 and that are all assumed to be infected.
-
β j k are replaced by β k just to avoid more complications in the program code. In other words, we assume that the probability to be infected does not depend on the source location of susceptibility, but on the source location of infectivity only, namely Ω k . More explicitly, β 11 = β 21 = β 31 ; they are represented by β 1 , β 12 = β 22 = β 32 ; they are represented by β 2 and, finally, β 13 = β 23 = β 33 = β 3 ; and they are represented by β 3 .
-
The unit of the parameters d j , γ j and β j k is d a y s - 1 ∀ fixed j and mobile k.
-
i t i + 1 = t i + h , Δ W i h N ( 0 , 1 ) with h the time step.
-
The coefficients in diffusions δ j , ς j k , and σ j are assumed to be all equal to 0.125. Larger values can be considered; however, they only increase the level of stochasticity, and this is not very interesting here.
Figure 1 depicts the dynamics of S Ω j , I Ω j , and R Ω j , j = 1 , 2 , 3 in the absence and presence of optimal controls u Ω 1 * , u Ω 2 * , and u Ω 3 * . As we can observe in the upper part of this figure and because of a strong immunity of a fraction of the population, there is a natural recovery from the infection and that causes an increase of the number of R Ω j in all three regions, but without exceeding the number of 40,000 people.
In the same part, we can also see that the number of infected people, namely I Ω j j , exceeds 70,000 people, and this proves the necessity of an urgent control measure to fight against the epidemic. In the lower part of this figure, we reach our goal that is summarized in maximizing the number of removed people while minimizing the number of infected people, and we can observe that the number of removed individuals has reached now a number that is near 90,000 people in all regions, while the number of infected people tended to zero values after 30, 40, and 48 days in regions Ω 1 , Ω 2 , and Ω 3 , respectively. We note that the optimal controls u Ω j * , j take the same values along the optimal control strategy period and that are all equal to one, and they tend to zero at the final time T. This is also the reason behind our avoidance of exhibiting simulations of the three optimal controls separately.
In Figure 2, we deduce that when we follow the optimal vaccination control strategy in the two regions Ω 2 and Ω 3 , the number R Ω j , j = 2 , 3 exceeded 100,000 individuals in the first five days, but decreased after towards values around 80,000 people, which represents an important value also, but they were smaller values than the ones obtained in the case treated in the first figure. In addition, the number I Ω j , j = 2 , 3 decreased to values around 2000 people, which shows also the effectiveness of the followed vaccination control strategy, but they were larger values than the ones obtained in the previous case. In parallel, despite the cancellation of vaccination in region Ω 1 , we can see that the number of removed people approached values around 50,000, which shows an increase by 10,000 people compared to the case when controls were not introduced in all regions, and this proves some influence of the optimal control policy followed in other regions in a region that is not yet controlled directly. We note also here that the optimal controls u Ω j * , j = 2 , 3 take the same values along the optimal control strategy period, which were all equal to one with some perturbations after 5 days, and they tended to zero at the final time T.
In Figure 3, we study the dynamics of S Ω j , I Ω j , and R Ω j , j = 1 , 2 , 3 when we introduce an optimal vaccination control in region Ω 1 only, and we can deduce some negative impact of the two regions Ω 2 and Ω 3 on the first region due to the absence of a control strategy in them. In fact, as they were at higher risk of infection than region Ω 1 , the optimal control u Ω 1 * seemed not sufficient for eradicating infection, as shown in first figure when controls were introduced in all regions, but this can at least increase the number of removed people in this region to an important value that exceeds 90,000 individuals. We note that the optimal control u Ω 1 * takes one as the maximal value.

4. The Model with Vaccination Plus Movement Restriction

4.1. Presentation of the Control Model

Let I = 1 , . . . , p and I H I be the set of indices of regions at high-risk and then having the ability to spread the epidemic to other regions. Here, we study the case when a given region Ω j is under vaccination control u Ω j and at the same time under the threat of infection coming from other regions. For this, we add to the vaccination strategy another control denoted as v j Ω k to characterize the effectiveness of movement restriction operations, in order to prevent the infected of regions Ω k , k I H coming to the controlled region Ω j , where:
v j Ω k 0 k I H k j v j Ω k = 0 elsewhere
Then, the model (24)–(26) in the controlled region Ω j is rewritten as follows:
S ˙ Ω j ( t ) = - k = 1 p ϑ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + N Ω j ( t ) - S Ω j ( t ) d j - θ Ω j ( t ) S Ω j ( t )
I ˙ Ω j ( t ) = k = 1 p ϑ j k ( t ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) - γ j I Ω j ( t ) - d j I Ω j ( t )
R ˙ Ω j ( t ) = γ j I Ω j ( t ) - d j R Ω j ( t ) + θ Ω j ( t ) S Ω j ( t )
with the vaccination control defined as:
θ Ω j ( t ) = u Ω j ( t ) + δ j d W Ω j ( t ) d t .
and the function ϑ j k ( t ) defined as:
ϑ j k ( t ) = 1 - v j Ω k β j k + ( 1 - ς j k ) σ j d W Ω j ( t ) d t
where σ j and ς j k , k I H , ( j = 1 , , p ) are real constants and representing the intensities of fluctuations caused by media and escapes of infected people through borders between Ω k and Ω j , respectively.

4.2. A Stochastic Optimal Control Approach

Now, the problem (Section 3.2.1) is changed in this part to
J ( u Ω j , v j Ω ) = k I H E 0 T α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 + B j 2 ( v j Ω k ( t ) ) 2 d t where B j > 0 is the weight constant of the new control, while u Ω j U j and v j Ω = v j Ω k k I H belong to the control set V j I H defined as:
V j I H ( [ 0 , T ] ) = { v j Ω ( t ) F t - p r o g r e s s i v e l y m e a s u r a b l e | v m i n Ω j v j Ω k ( t ) v m a x Ω j , k I H , t [ 0 , T ] }
The Hamiltonian in this case is defined as:
H = k I H α j I I Ω j ( t ) - α j R R Ω j ( t ) + A j 2 ( u Ω j ( t ) ) 2 + B j 2 ( v j Ω k ( t ) ) 2 + μ Ω j T ( t ) f ( t , x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) + l = 1 3 g l T ( x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) ν Ω j l ( t )
where the state function f is defined as:
f ( t , x Ω ( t ) , u Ω ( t ) , v j Ω ( t ) )
= - k = 1 p 1 - v j Ω k β j k I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + N Ω j ( t ) - S Ω j ( t ) d j - u Ω j ( t ) S Ω j ( t ) k = 1 p 1 - v j Ω k β j k I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) - γ j I Ω j ( t ) - d j I Ω j ( t ) γ j I Ω j ( t ) - d j R Ω j ( t ) + u Ω j ( t ) S Ω j ( t )
and the diffusion matrix g is defined as:
g ( t , x Ω ( t ) , u Ω ( t ) , v j Ω ( t ) ) = - σ j k = 1 p ( 1 - ς j k ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) + δ j S Ω j ( t ) σ j k = 1 p ( 1 - ς j k ) I Ω k ( t ) N Ω j ( t ) S Ω j ( t ) δ j S Ω j ( t )
while in a domain Ω , ( μ ( t ) , ν ( t ) ) ) is a pair of adjoint variables satisfying the following adjoint BSDE (backward stochastic differential equation):
d μ Ω ( t ) = - [ f x T ( t , x Ω ( t ) , u Ω ( t ) , v j Ω ( t ) ) μ Ω ( t ) + l = 1 3 g x Ω l T ( t , x Ω ( t ) , u Ω ( t ) , v j Ω ( t ) ) ν Ω l ( t ) + f 0 x Ω ( t , x Ω ( t ) , u Ω ( t ) , v j Ω ( t ) ) ] d t + ν Ω ( t ) d W Ω ( t ) , μ Ω ( T ) = 0 .

4.2.1. Optimal Control Characterization and Necessary Conditions

Using the stochastic Pontryagin’s maximum principle as done for the first control strategy, we obtain the following optimal control characterization and necessary conditions:
Theorem 2.
If there exists an optimal pair ( x Ω * , u Ω * , v j Ω * ) and a pair of processes ( μ ( t ) , ν ( t ) ) satisfying (27), then for j = 1 , . . . , p , k I H , we have:
H ( t , x Ω j * ( t ) , u Ω j * ( t ) , v j Ω k * ( t ) , μ Ω j ( t ) , ν Ω j ( t ) )
= min ( u Ω j , v j Ω k ( t ) ) U j × V j I H H ( t , x Ω j ( t ) , u Ω j * ( t ) , v j Ω k * ( t ) , μ Ω j ( t ) , ν Ω j ( t ) ) .
Moreover, we obtain the bounded stochastic control:
u Ω j * = m i n ( m a x ( u m i n Ω j , - ( μ 3 Ω j ( t ) - μ 1 Ω j ( t ) ) S Ω j * A j ) , u m a x Ω j ) ,
v j Ω k * = m i n ( m a x ( v m i n Ω j , - ( μ 1 Ω j ( t ) - μ 2 Ω j ( t ) ) β j k I Ω k * S Ω j * B j ) , v m a x Ω j ) ,
and solutions of the FBSDEs (forward-backward stochastic differential equations):
d x Ω j ( t ) = f ( t , x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) d t + g ( t , x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) d W Ω j ( t ) d μ Ω j ( t ) = - [ f x Ω j T ( t , x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) μ Ω j ( t ) + l = 1 3 g x Ω j l T ( t , x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) ν Ω j l ( t ) + f 0 x Ω j ( t , x Ω j ( t ) , u Ω j ( t ) , v j Ω k ( t ) ) ] d t + ν Ω j ( t ) d W Ω j ( t ) , x Ω j ( 0 ) = ( S 0 Ω j , I 0 Ω j , R 0 Ω j ) μ Ω j ( T ) = 0 .

4.2.2. Numerical Results

In this part, we show the importance of following an optimal control strategy that is based on movement restrictions along with the presence of the vaccination policy discussed previously. For this, we take the example of preventing the epidemic from arriving at the region Ω 1 through infected travelers who come from other regions. This means we restrict movements of infected people coming from Ω 2 and Ω 3 either using v 1 Ω 2 * and/or v 1 Ω 3 * , respectively. Figure 4 presents the simulations of I Ω j and R Ω j , j = 1 , 2 , 3 in the presence of all optimal controls u Ω 1 * , u Ω 2 * , and u Ω 3 * , and we can see that when we restrict movement of individuals coming from region Ω 2 and aiming to enter Ω 1 , we can reach smaller values of infected people after only 15 days in all regions; at the same time, we deduce that the level of the number of removed people is approximately the same as in the absence of movement restrictions. We note that the optimal controls u Ω j * , j = 2 , 3 along with the optimal control v 1 Ω 2 * take the same values along the optimal control strategy period, which are all equal to one, and they tend to zero at the final time T.
In Figure 5, we can better see the importance of movement restrictions, especially when we prevent also individuals coming from region Ω 3 , which is at highest risk of infection. In fact, when we add the optimal control v 1 Ω 2 * to v 1 Ω 3 * , we can reach better results as in the previous case since the number of infected people has decreased more and the number of removed people has increased more in all regions compared to the case when optimal controls u Ω 1 * , u Ω 2 * , and u Ω 3 * are followed alone. We note also here that the optimal controls u Ω j * , j = 2 , 3 along with the optimal controls v 1 Ω k * , k = 2 , 3 take the same values along the optimal control strategy period, which are all equal to one, and they tend to zero at the final time T.

5. Discussion

In this section, we provide other simulations to compare between the results of the optimal control strategy proposed in Section 3 and the other one treated in Section 4, as well as the case when there was no control yet as defined in Section 1. We also discuss the effect of severity controls weight A j and B j and the initial conditions of infection on the optimal values of controls and states.
In Figure 6, we treat the case when there is no vaccination in region Ω 1 , and we compare it with the case when there is no control policy in all regions and with the case treated in Figure 1. From this figure, we can deduce that even in the absence of vaccination in the first region, there is an advantage of movement restrictions represented by optimal controls v 1 Ω k * , k = 2 , 3 , applied between region Ω 1 and other regions, in decreasing the number of infected people and increasing the number of removed people. This figure also shows the importance of vaccination when it is followed in all regions, and that can be strengthened more if we add v 1 Ω k * , k = 2 , 3 , as concluded in the previous figure. We note also here that the optimal controls u Ω j * , j = 1 , 2 , 3 along with the optimal controls v 1 Ω k * , k = 2 , 3 take the same values along the optimal control strategy period, which are all equal to one, with some perturbations after five days, and they tend to zero at the final time T.
As in Figure 7, we present the simulation of the function I Ω 1 when there is no control, and we compare it with the case with the control for different values of severity controls weights, namely A j and B j associated with u Ω j * and v j Ω k * , respectively, as noted previously. We observe that after introducing the five optimal controls u Ω j * , j = 1 , 2 , 3 , v 1 Ω 2 * , and v 1 Ω 3 * , the number of infected people has decreased towards the value of 400 individuals when A j = 1 , referring to the previous cases for obtaining u Ω j * , but with B j = 10 20 now. This number can be decreased earlier, as we can see below in the left-hand side part of this figure when we take A j = 10 4 . Defined as the denominator in the characterization of v j Ω k * , the value B j = 10 20 here implies a reduction of the optimal controls’ values, v 1 Ω 2 * and v 1 Ω 3 * , and the movement restriction is recommended for the first 22 and 25 days only. However, when we give a less important value to B j , 10 8 for example, we can observe that a long period of this strategy minimized the function I Ω 1 more significantly, and it can be minimized and even earlier than this when we take A j = B j = 10 4 , as seen below in the right-hand side of the figure.
As the last Figure 8, we are interested in doing the same as in the previous figure, but in the case of less important values for initial conditions I 0 Ω j , j = 1 , 2 , 3 and smaller values for optimal controls u Ω j * , j = 1 , 2 , 3 , v 1 Ω 2 * , and v 1 Ω 3 * . From this figure, we observe that when there is no control, the function I Ω 1 does not increase to larger values as observed in the first figures, and then, we understand that we will not need to vaccinate everyone and close all borders for all times of the control strategy period. As predicted, we observe that when we add the value 10 20 to both denominators A j and B j , implying an important reduction of all optimal controls, the number of infected people in Ω 1 decreased to very small values after 20 days only, and this can be done more effectively when vaccination is prolonged, as we can observe in the right-hand side of the figure, when we add 10 10 to A j .

6. Conclusions

In this paper, we have proposed a stochastic model for the study of infection dynamics when an epidemic emerges in regions that are connected by any kind of anthropological movement, considering perturbations that are due to media, infodemics, and escapes. We have also suggested two optimal control approaches related to vaccination and movement restriction policies based on a stochastic version of Pontryagin’s maximum principle.

Author Contributions

All authors contributed equally to this work. All authors read and approved the final manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank all Managing Editors and members of the Editorial Board who were responsible for dealing with this paper, and the anonymous referees for their valuable comments and suggestions, improving the content of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

All symbols used in this paper are presented in the following nomenclature.
ttime
Ttime horizon
Ω a large domain or region
x Ω R n random state vector associated with Ω
u Ω R m random control vector associated with Ω
( Θ , F , P ) probability space
W Ω vector-valued Wiener process associated with Ω all over ( Θ , F , P )
f R n vector-valued nonlinear function
gmatrix-valued nonlinear diffusion
Ω j a region, an area, or subdomain of Ω ( 1 j p )
S Ω j number of susceptibles in Ω j
I Ω j number of infectives in Ω j
R Ω j number of removed in Ω j
ρ j k stochastic proportion of adequate contacts in Ω j between a susceptible from Ω j and an infective from Ω k
β j k deterministic proportion of adequate contacts in Ω j between a susceptible from Ω j and an infective from Ω k
σ j intensities of fluctuations caused by media
d j birth and death rate
γ j recovery rate
N Ω j population size corresponding to Ω j
u Ω j vaccination control introduced in Ω j
θ Ω j perturbation of control function u Ω j associated with Ω j
u m i n Ω j minimal bound of u Ω j ( t )
u m a x Ω j maximal bound of u Ω j ( t )
Jobjective function
f 0 current gain function
α j I weight parameter associated with the number of infectives in Ω j
α j R weight parameter associated with the number of removed in Ω j
A j vaccination control severity weight in Ω j
B j movement restriction control severity weight in Ω j
U j vaccination control set
μ Ω adjoint state variable associated with x Ω
ν Ω adjoint matrix diffusion associated with g
I H set of indices of regions at a high-risk of infection
V j I H movement restriction control set
v j Ω k movement restriction control introduced in Ω j to prevent infection from Ω k
ϑ j k perturbation of control function v j Ω k associated to Ω j
ς j k intensities of fluctuations caused by escapes
v m i n Ω j minimal bound of v j Ω k ( t )
v m a x Ω j maximal bound of v j Ω k ( t )

References

  1. Liu, R.; Wu, J.; Zhu, H. Media/psychological impact on multiple outbreaks of emerging infectious diseases. Comput. Math. Methods Med. 2007, 8, 153–164. [Google Scholar] [CrossRef]
  2. Sahu, G.P.; Dhar, J. Dynamics of an SEQIHRS epidemic model with media coverage, quarantine and isolation in a community with pre-existing immunity. J. Math. Anal. Appl. 2015, 421, 1651–1672. [Google Scholar] [CrossRef]
  3. Wang, A.; Xiao, Y. A Filippov system describing media effects on the spread of infectious diseases. Nonlinear Anal. Hybrid Syst. 2014, 11, 84–97. [Google Scholar] [CrossRef]
  4. Liu, Y.; Cui, J.A. The impact of media coverage on the dynamics of infectious disease. Int. J. Biomath. 2008, 1, 65–74. [Google Scholar] [CrossRef]
  5. Samanta, S.; Rana, S.; Sharma, A.; Misra, A.K.; Chattopadhyay, J. Effect of awareness programs by media on the epidemic outbreaks: A mathematical model. Appl. Math. Comput. 2013, 219, 6965–6977. [Google Scholar] [CrossRef]
  6. Misra, A.K.; Sharma, A.; Shukla, J.B. Stability analysis and optimal control of an epidemic model with awareness programs by media. Biosystems 2015, 138, 53–62. [Google Scholar] [CrossRef]
  7. Yuan, X.; Xue, Y.; Liu, M. Analysis of an epidemic model with awareness programs by media on complex networks. Chaos Solitons Fractals 2013, 48, 1–11. [Google Scholar] [CrossRef]
  8. Wang, Y.; Cao, J.; Jin, Z.; Zhang, H.; Sun, G.Q. Impact of media coverage on epidemic spreading in complex networks. Phys. A Stat. Mech. Its Appl. 2013, 392, 5824–5835. [Google Scholar] [CrossRef]
  9. World Health Organization. Managing Epidemics: Key Facts About Major Deadly Diseases; World Health Organization: Geneva, Switzerland, 2018. [Google Scholar]
  10. Zakary, O.; Larrache, A.; Rachik, M.; Elmouki, I. Effect of awareness programs and travel-blocking operations in the control of HIV/AIDS outbreaks: A multi-domains SIR model. Adv. Differ. Equ. 2016, 2016, 169. [Google Scholar] [CrossRef]
  11. Zakary, O.; Rachik, M.; Elmouki, I. A multi-regional epidemic model for controlling the spread of Ebola: Awareness, treatment, and travel-blocking optimal control approaches. Math. Methods Appl. Sci. 2017, 40, 1265–1279. [Google Scholar] [CrossRef]
  12. Zakary, O.; Rachik, M.; Elmouki, I. On the analysis of a multi-regions discrete SIR epidemic model: An optimal control approach. Int. J. Dyn. Control 2017, 5, 917–930. [Google Scholar] [CrossRef]
  13. Zakary, O.; Rachik, M.; Elmouki, I. A new analysis of infection dynamics: Multi-regions discrete epidemic model with an extended optimal control approach. Int. J. Dyn. Control 2017, 5, 1010–1019. [Google Scholar] [CrossRef]
  14. Abouelkheir, I.; El Kihal, F.; Rachik, M.; Zakary, O.; Elmouki, I. A multi-regions SIRS discrete epidemic model with a travel-blocking vicinity optimal control approach on cells. Br. J. Math. Comput. Sci. 2017, 20, 1–16. [Google Scholar] [CrossRef] [PubMed]
  15. El Kihal, F.; Rachik, M.; Zakary, O.; Elmouki, I. A multi-regions SEIRS discrete epidemic model with a travel-blocking vicinity optimal control approach on cells. Int. J. Adv. Appl. Math. Mech. 2017, 4, 60–71. [Google Scholar]
  16. Abouelkheir, I.; Rachik, M.; Zakary, O.; Elmouki, I. A multi-regions SIS discrete influenza pandemic model with a travel-blocking vicinity optimal control approach on cells. Am. J. Comput. Appl. Math. 2017, 7, 37–45. [Google Scholar]
  17. Zakary, O.; Rachik, M.; Elmouki, I.; Lazaiz, S. A multi-regions discrete-time epidemic model with a travel-blocking vicinity optimal control approach on patches. Adv. Differ. Equ. 2017, 2017, 120. [Google Scholar] [CrossRef]
  18. Chouayakh, K.; Rachik, M.; Zakary, O.; Elmouki, I. A multi-regions SEIS discrete epidemic model with a travel-blocking vicinity optimal control approach on cells. J. Math. Comput. Sci. 2017, 7, 468–484. [Google Scholar]
  19. Zakary, O.; Rachik, M.; Elmouki, I. A new epidemic modeling approach: Multi-regions discrete-time model with travel-blocking vicinity optimal control strategy. Infect. Dis. Model. 2017, 2, 304–322. [Google Scholar] [CrossRef] [PubMed]
  20. Bidah, S.; Rachik, M.; Zakary, O.; Boutayeb, H.; Elmouki, I. Travel-blocking Optimal Control Policy on Borders of a Chain of Regions Subject to SIRS Discrete Epidemic Model. Asian J. Res. Infect. Dis. 2018, 1, 1–12. [Google Scholar]
  21. Zakary, O.; Bidah, S.; Rachik, M.; Elmouki, I. Cell and patch vicinity travel restrictions in a multi-regions SI discrete epidemic control model. Int. J. Adv. Appl. Math. Mech. 2019, 6, 30–41. [Google Scholar]
  22. Zhang, X.B.; Huo, H.F.; Xiang, H.; Meng, X.Y. Dynamics of the deterministic and stochastic SIQS epidemic model with non-linear incidence. Appl. Math. Comput. 2014, 243, 546–558. [Google Scholar] [CrossRef]
  23. Ji, C.; Jiang, D. Threshold behaviour of a stochastic SIR model. Appl. Math. Model. 2014, 38, 5067–5079. [Google Scholar] [CrossRef]
  24. Ji, C.; Jiang, D.; Shi, N. The behavior of an SIR epidemic model with stochastic perturbation. Stoch. Anal. Appl. 2012, 30, 755–773. [Google Scholar] [CrossRef]
  25. Jiang, D.; Ji, C.; Shi, N.; Yu, J. The long time behavior of DI SIR epidemic model with stochastic perturbation. J. Math. Anal. Appl. 2010, 372, 162–180. [Google Scholar] [CrossRef] [Green Version]
  26. Liu, W. A SIRS epidemic model incorporating media coverage with random perturbation. In Abstract and Applied Analysis; Hindawi: Cairo, Egypt, 2013; Volume 2013. [Google Scholar]
  27. Cai, Y.; Wang, X.; Wang, W.; Zhao, M. Stochastic dynamics of an SIRS epidemic model with ratio-dependent incidence rate. In Abstract and Applied Analysis; Hindawi: Cairo, Egypt, 2013; Volume 2013. [Google Scholar]
  28. Cai, Y.; Kang, Y.; Wang, W. A stochastic SIRS epidemic model with nonlinear incidence rate. Appl. Math. Comput. 2017, 305, 221–240. [Google Scholar] [CrossRef]
  29. Cai, Y.; Kang, Y.; Banerjee, M.; Wang, W. A stochastic SIRS epidemic model with infectious force under intervention strategies. J. Differ. Equ. 2015, 259, 7463–7502. [Google Scholar] [CrossRef]
  30. Zhao, Y.; Jiang, D. The threshold of a stochastic SIRS epidemic model with saturated incidence. Appl. Math. Lett. 2014, 34, 90–93. [Google Scholar] [CrossRef]
  31. Zhao, Y.; Jiang, D. Dynamics of stochastically perturbed SIS epidemic model with vaccination. In Abstract and Applied Analysis; Hindawi: Cairo, Egypt, 2013; Volume 2013. [Google Scholar]
  32. Zhao, Y.; Jiang, D.; O’Regan, D. The extinction and persistence of the stochastic SIS epidemic model with vaccination. Phys. A Stat. Mech. Its Appl. 2013, 392, 4916–4927. [Google Scholar] [CrossRef]
  33. Lin, Y.; Jiang, D.; Wang, S. Stationary distribution of a stochastic SIS epidemic model with vaccination. Phys. A Stat. Mech. Its Appl. 2014, 394, 187–197. [Google Scholar] [CrossRef]
  34. Witbooi, P.J. Stability of an SEIR epidemic model with independent stochastic perturbations. Phys. A Stat. Mech. Its Appl. 2013, 392, 4928–4936. [Google Scholar] [CrossRef] [Green Version]
  35. Zhou, Y.; Zhang, W.; Yuan, S. Survival and stationary distribution of a SIR epidemic model with stochastic perturbations. Appl. Math. Comput. 2014, 244, 118–131. [Google Scholar] [CrossRef]
  36. Gray, A.; Greenhalgh, D.; Hu, L.; Mao, X.; Pan, J. A stochastic differential equation SIS epidemic model. SIAM J. Appl. Math. 2011, 71, 876–902. [Google Scholar] [CrossRef]
  37. Meng, X.; Zhao, S.; Feng, T.; Zhang, T. Dynamics of a novel nonlinear stochastic SIS epidemic model with double epidemic hypothesis. J. Math. Anal. Appl. 2016, 433, 227–242. [Google Scholar] [CrossRef]
  38. Hethcote, H.W.; Waltman, P. Optimal vaccination schedules in a deterministic epidemic model. Math. Biosci. 1973, 18, 365–381. [Google Scholar] [CrossRef]
  39. Yong, J.; Zhou, X.Y. Stochastic Controls Hamiltonian Systems and HJB equations. In Application of Mathematics; Springer: New York, NY, USA, 1999. [Google Scholar]
  40. Bahlali, K.; Mezerdi, B.; Ouknine, Y. The maximum principle for optimal control of diffusions with non-smooth coeffcients. Stoch. Stoch. Rep. 1996, 57, 303–316. [Google Scholar] [CrossRef]
  41. Peng, S. A general stochastic maximum principle for optimal control problems. SIAM J. Control Optim. 1990, 28, 966–979. [Google Scholar] [CrossRef]
  42. Ma, J.; Protter, P.; Yong, J. Solving Forward-Backward Stochastic Differential Equations Explicitly—A Four Step Scheme. J. Probab. Theory. Relat. Fields 1994, 98, 339–3359. [Google Scholar] [CrossRef]
  43. Ladyz̆henskaia, O.A.; Solonnikov, V.A.; Ural’ceva, N.N. Linear and Quasi-Linear Equations of Parabolic Type; American Mathematical Soc.: Washington, DC, USA, 1968; Volume 23. [Google Scholar]
  44. Düring, B.; Jüngel, A. Existence and uniqueness of solutions to a quasilinear parabolic equation with quadratic gradients in financial markets. Nonlinear Anal. Theory Methods Appl. 2005, 62, 519–544. [Google Scholar] [CrossRef]
  45. Aboulaich, R.; Darouichi, A.; Elmouki, I.; Jraifi, A. A Stochastic Optimal Control Model for BCG Immunotherapy in Superficial Bladder Cancer. Math. Model. Nat. Phenom. 2017, 12, 99–119. [Google Scholar] [CrossRef]
Figure 1. (Top) S Ω j I Ω j R Ω j , j = 1 , 2 , 3 stochastic dynamics without controls. (Bottom) S Ω 1 I Ω 1 R Ω 1 stochastic dynamics in the presence of all optimal controls u Ω 1 * , u Ω 2 * , and u Ω 3 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Figure 1. (Top) S Ω j I Ω j R Ω j , j = 1 , 2 , 3 stochastic dynamics without controls. (Bottom) S Ω 1 I Ω 1 R Ω 1 stochastic dynamics in the presence of all optimal controls u Ω 1 * , u Ω 2 * , and u Ω 3 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Mathematics 07 00304 g001
Figure 2. S Ω j I Ω j R Ω j , j = 1 , 2 , 3 stochastic dynamics in the absence of optimal control u Ω 1 * and the presence of optimal controls u Ω 2 * and u Ω 3 * . Same parameters and initial conditions as in the first figure.
Figure 2. S Ω j I Ω j R Ω j , j = 1 , 2 , 3 stochastic dynamics in the absence of optimal control u Ω 1 * and the presence of optimal controls u Ω 2 * and u Ω 3 * . Same parameters and initial conditions as in the first figure.
Mathematics 07 00304 g002
Figure 3. S Ω j I Ω j R Ω j , j = 1 , 2 , 3 stochastic dynamics in the absence of optimal controls u Ω 2 * and u Ω 3 * and the presence of optimal control u Ω 1 * only. S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Figure 3. S Ω j I Ω j R Ω j , j = 1 , 2 , 3 stochastic dynamics in the absence of optimal controls u Ω 2 * and u Ω 3 * and the presence of optimal control u Ω 1 * only. S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Mathematics 07 00304 g003
Figure 4. (Red and green curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * . (Black curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * plus v 1 Ω 2 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Figure 4. (Red and green curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * . (Black curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * plus v 1 Ω 2 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Mathematics 07 00304 g004
Figure 5. (Red and green curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * . (Black curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * plus v 1 Ω 2 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Figure 5. (Red and green curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * . (Black curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * plus v 1 Ω 2 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Mathematics 07 00304 g005
Figure 6. (Blue curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics without controls. (Red and green curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * . (Black curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 2 * and u Ω 3 * plus v 1 Ω 2 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Figure 6. (Blue curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics without controls. (Red and green curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 1 * , u Ω 2 * and u Ω 3 * . (Black curves) I Ω j and R Ω j , j = 1 , 2 , 3 dynamics with optimal controls u Ω 2 * and u Ω 3 * plus v 1 Ω 2 * . S 0 Ω 1 = 90,000, I 0 Ω 1 = 1200 , S 0 Ω 2 = 89,000, I 2 Ω 2 = 1100 , S 0 Ω 3 = 88,000, I 0 Ω 3 = 1000 , R 0 Ω j = 0 j = 1 , 2 , 3 , d 1 = 0 . 06 , γ 1 = 0 . 04 , β 1 = 0 . 5 , d 2 = 0 . 05 , γ 2 = 0 . 03 , β 2 = 0 . 4 , d 3 = 0 . 04 , γ 3 = 0 . 02 , β 3 = 0 . 1 .
Mathematics 07 00304 g006
Figure 7. I Ω 1 and optimal controls u Ω j * , j = 1 , 2 , 3 , v 1 Ω 2 * and v 1 Ω 3 * for different values of A j and B j .
Figure 7. I Ω 1 and optimal controls u Ω j * , j = 1 , 2 , 3 , v 1 Ω 2 * and v 1 Ω 3 * for different values of A j and B j .
Mathematics 07 00304 g007
Figure 8. I Ω 1 and optimal controls u Ω j * , j = 1 , 2 , 3 , v 1 Ω 2 * , and v 1 Ω 3 * for different values of A j and B j and less important initial condition of infection.
Figure 8. I Ω 1 and optimal controls u Ω j * , j = 1 , 2 , 3 , v 1 Ω 2 * , and v 1 Ω 3 * for different values of A j and B j and less important initial condition of infection.
Mathematics 07 00304 g008

Share and Cite

MDPI and ACS Style

El Kihal, F.; Abouelkheir, I.; Rachik, M.; Elmouki, I. Role of Media and Effects of Infodemics and Escapes in the Spatial Spread of Epidemics: A Stochastic Multi-Region Model with Optimal Control Approach. Mathematics 2019, 7, 304. https://doi.org/10.3390/math7030304

AMA Style

El Kihal F, Abouelkheir I, Rachik M, Elmouki I. Role of Media and Effects of Infodemics and Escapes in the Spatial Spread of Epidemics: A Stochastic Multi-Region Model with Optimal Control Approach. Mathematics. 2019; 7(3):304. https://doi.org/10.3390/math7030304

Chicago/Turabian Style

El Kihal, Fadwa, Imane Abouelkheir, Mostafa Rachik, and Ilias Elmouki. 2019. "Role of Media and Effects of Infodemics and Escapes in the Spatial Spread of Epidemics: A Stochastic Multi-Region Model with Optimal Control Approach" Mathematics 7, no. 3: 304. https://doi.org/10.3390/math7030304

APA Style

El Kihal, F., Abouelkheir, I., Rachik, M., & Elmouki, I. (2019). Role of Media and Effects of Infodemics and Escapes in the Spatial Spread of Epidemics: A Stochastic Multi-Region Model with Optimal Control Approach. Mathematics, 7(3), 304. https://doi.org/10.3390/math7030304

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