Next Article in Journal
The Effect of Learning Rate on Fractal Image Coding Using Artificial Neural Networks
Previous Article in Journal
Financial Applications on Fractional Lévy Stochastic Processes
Previous Article in Special Issue
New Exact Solutions of Some Important Nonlinear Fractional Partial Differential Equations with Beta Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic Optimal Control Analysis of a Mathematical Model: Theory and Application to Non-Singular Kernels

1
Department of Mathematics, Sun Yat-Sen University, Guangzhou 510275, China
2
Department of Mathematics, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Submission received: 1 December 2021 / Revised: 20 December 2021 / Accepted: 27 December 2021 / Published: 23 May 2022

Abstract

:
Some researchers believe fractional differential operators should not have a non-singular kernel, while others strongly believe that due to the complexity of nature, fractional differential operators can have either singular or non-singular kernels. This contradiction in thoughts has led to the publication of a few papers that are against differential operators with non-singular kernels, causing some negative impacts. Thus, publishers and some Editors-in-Chief are concerned about the future of fractional calculus, which has generally brought confusion among the vibrant and innovative young researchers who desire to apply fractional calculus within their respective fields. Thus, the present work aims to develop a model based on a stochastic process that could be utilized to portray the effect of arbitrary-order derivatives. A nonlinear perturbation is used to study the proposed stochastic model with the help of white noises. The required condition(s) for the existence of an ergodic stationary distribution is obtained via Lyapunov functional theory. The finding of the study indicated that the proposed noises have a remarkable impact on the dynamics of the system. To reduce the spread of a disease, we imposed some control measures on the stochastic model, and the optimal system was achieved. The models both with and without control were coded in MATLAB, and at the conclusion of the research, numerical solutions are provided.

1. Introduction to the Problem and Model Formulation

Different types of epidemic models exist, including fractional, deterministic, stochastic, and age-structure models [1,2,3,4,5,6,7,8,9,10]. The stochastic differential equation models are more suited for mathematical modeling than the deterministic ones [1,3], because they may give a higher level of realism than their deterministic counterparts. When we run a stochastic model multiple times, we may build up a distribution of the anticipated outcomes, such as the size of infectious compartments at time t, which is more helpful than when we run a deterministic model. A deterministic model, on the other hand, will only provide a single predicted value.
Because of their extensive use, fractional integration and differentiation are rapidly growing subjects that have piqued the interest of a variety of researchers. Since then, this issue has piqued the interest of researchers from many areas [11,12,13,14,15,16,17,18,19,20,21]. The problem arose from a query posed by L’Hopital to Leibniz, which originally highlighted the issue of exponential function differentiation. At later stages, Liouville proposed a name for the derivative with a fractional index, while also cautioning that the characteristics of this new operator should not be confused with those of the integer-order derivative [11,12,13,14,15,16,22,23]. Following that, a collaboration between Liouville and Riemann produced the well-known integral of arbitrary order and, later, an arbitrary-order derivative. A derivative of a function that is continuous and a power-law function is their differential operator t α Γ ( 1 α ) . The Laplace transform of this derivative produces peculiar starting circumstances that would be hard to compute otherwise [11,12,13,14,15,16].
As a result, a novel set of differential operators was required, with a few variations proposed. The Caputo–Fabrizio (CF) derivatives and the Atangana–Baleanu (AB) derivative with the Mittag–Leffler and exponential kernels are two examples. One of these operator’s peculiarities is that its kernels are non-singular [9,11,12,13,14,15,16] and that it does not satisfy several classical derivative characteristics such as the index law. These fractional derivatives do not have the same characteristics as fractional or classical derivatives with the power-law kernel because of their continuous kernels. This shows that in nature, exponential, power-law, and Mittag–Leffler functions do not all perform the same role [24,25]. There are issues in nature that follow power-law procedures; a few follow decaying procedures; many others have degradation and afterwards power-law procedures. Because the fractional derivative with a singular kernel and the fractional derivative with non-singular kernels are two different subsets of fractional calculus, both ideas may have been explored individually in a normal and constructive field. As a result, they should be addressed separately.
Nevertheless, rather than examining the characteristics of non-singular kernel fractional derivatives separately, some writers have chosen to write some ill-structured papers that combine the two sub-sets; while other research articles are based on incorrect analyses, which are the result of a lack of understanding of the entire concept. As a result, this paper will go over all of the criticisms leveled about non-singular kernels’ differential and integral operators. The following fractional calculus chaotic survival characterized by non-singular kernels was suggested in [26]:
d F c = Λ d F c ( t ) ϕ 2 F c ( t ) β F c ( t ) I ( t ) N d t + η 1 F c ( t ) d B 1 ( t ) , d I ( t ) = β F c ( t ) I ( t ) N ( d + τ + δ + ψ ) I ( t ) d t + η 2 I ( t ) d B 2 ( t ) , d I P ( t ) = τ I ( t ) ( d + γ 1 + ϕ 1 ) I P ( t ) d t + η 3 I P ( t ) d B 3 ( t ) , d I N ( t ) = ψ I I N ( t ) ( d + σ + γ 3 ) d t + η 4 I N ( t ) d B 4 ( t ) , d R ( t ) = ϕ 1 I P ( t ) d R ( t ) + δ I ( t ) + γ 3 I N ( t ) γ 2 R ( t ) d t + η 5 R ( t ) d B 5 ( t ) , d D ( t ) = γ 2 R ( t ) + γ 1 I P ( t ) + ϕ 2 F c ( t ) + σ I N ( t ) d D ( t ) d t + D ( t ) η 6 d B 6 ( t ) ,
where F c ( t ) is the class of researchers that use the arbitrary-order operators of the differential type with non-singular kernels vulnerable to being misguided by critics (harmful) and I ( t ) stands for the class of those researchers that are the affectees of harmful research articles. The notion I P ( t ) denotes those researchers that are affected and yet they have a constructive approach towards the non-singular kernel. I N ( t ) represents those researchers that are affected and also have a bad opinion about derivatives that have a non-singular kernel. R ( t ) are those researchers that overcame the division of criticism, and finally, the researchers that are dead, retired, or have left the research of fractional derivatives and integrals are denoted by D ( t ) . The topic is extremely appealing because the notion is widely applicable; therefore, each time a new researcher joins the area, we use Λ as the recruiting rate in this model. Naturally, some researchers working in this area may die, and as a result, they will no longer be working in the field; as a result, we will choose d as the death rate. Retirement can be explained by d as well. β is interpreted as a contact coefficient, or the likelihood that a researcher will read a work concerning non-singular kernel derivatives with negative content. Similarly, the rate at which a fractional calculus researcher joins class D owing to division is represented by the parameter ϕ 2 . τ is the percentage of people who are impacted by criticism, yet still believe in non-singular kernel variants. The term δ is the rate of recovery. Due to division, ψ is the rate at which researchers join class I N , and γ 1 is the rate at which persons with favorable opinions join class D. The recovery rate of the I P class is ϕ 1 . The rate at which class I N joins class D owing to division is called σ . The recovery rate of class I is γ 3 . The recovery rate in class D is shown by γ 2 , whereas B i ( t ) such that i = 1 , 2 , 3 , 4 , 5 , 6 denotes the independent Brownian motion, and η i 2 for i = 1 , 2 , 3 , 4 , 5 , 6 stands for the white noise intensity.
In this work, the researchers were interested in formulating a stochastic model for the proposed social phenomenon and studying its possible control policy, besides the stability of the model. It is also of great concern to study the underlying deterministic model of the developed model, and thus, we analyzed its corresponding ODE model. The main goal of the control section of the study is to characterize the infection rate (being a control measure). This control reduces the amount of the infected population (over a given time interval), while keeping the expense of the control strategies in mind. We minimize the expected value when it comes to stochastic models. For the deterministic approach, the Hamilton–Jacobi–Bellman equation is utilized, but for the control theory of stochastic problems, the “Hamilton–Jacobi–Bellman equation” is employed. Hence, this paper reviews all issues raised against non-singular kernels’ differential and integral operators. A stochastic mathematical model depicting the survival of fractional calculus based on non-singular kernels is herein proposed and analyzed.
The rest of the article is constructed in the following fashion. In Section 2, we show the existence and uniqueness of global positive solutions of System (1) having a starting approximation. In Section 3, we give some necessary axioms for the extinction of the epidemic. In Section 4, we derive the existence of the solution of a stationary distribution by using the method of Khasminskii and making a proper Lyapunov operator. Likewise, we look into a stochastic control problem for optimality in Section 5. Some parallels were found in both the deterministic and stochastic control models. In Section 6, a few numerical simulations are presented to validate the required scheme. The main results of this paper end with the conclusion in Section 7.

2. Qualitative Analysis of the Global Non-Negative Solution

In order to investigate the dynamics of the proposed system, firstly, we are concerned about the solution being non-local and non-negative. In the following part, we make a proper “Lyapunov function” to verify the qualitative analysis of globally positive solutions of the problem (1).
Theorem 1.
For some initial approximation ( F c , I , I P , I N , R , D ) ( 0 ) R + 6 , a one positive root X = ( F c , I , I P , I N , R , D ) of model (1) exists on interval [ t , ) , and solution R + 4 has one chance of occurrence, namely the solution X R + 4 [ t , ) a.s.
Proof. 
For the proof of this theorem, see Theorem 5 in [26]. □

3. Extinction of the Proposed Model

We develop the conditions for model extinction in this section of the research. We include some notations and definitions below for convenience. Assume that:
X ( t ) = 1 t 0 t X ( r ) d r .
Lemma 1.
Let us assume that X represents a solution of Model (1) with given initial data ( F c , I , I P , I N , R , D ) ( 0 ) R + 6 , then:
lim t F c ( t ) + I P ( t ) + I ( t ) + I N ( t ) + D ( t ) + R ( t ) t = 0 , a . s .
Further, if d > ( η 1 2 η 2 2 η 3 2 η 4 2 η 5 2 η 6 2 ) 2 , then:
lim t 1 t 0 t F c ( r ) d B 1 ( r ) = 0 , lim t 1 t 0 t I ( r ) d B 2 ( r ) = 0 , lim t 1 t 0 t I P ( r ) d B 3 ( r ) = 0 , lim t 1 t 0 t I N ( r ) d B 4 ( r ) = 0 , lim t 1 t 0 t R ( r ) d B 5 ( r ) = 0 , lim t 1 t 0 t D ( r ) d B 6 ( r ) = 0 .
In the following, we calculate and define the threshold parameter R 0 for the model (1):
R 0 = β ( d + τ + δ + ψ + η 2 2 2 )
The following is the outcome of the illness extinction.
Theorem 2.
Assume that d > ( η 1 2 η 2 2 η 3 2 η 4 2 η 5 2 η 6 2 ) 2 and X represent a solution to the model (1) with subsidiary conditions ( F c , I , I P , I N , R , D ) ( 0 ) R 6 .
If R 0 < 1 , then:
lim t log I ( t ) t < 0 , lim t log I P ( t ) t < 0 a n d lim t log I N ( t ) t < 0 ,
a.s. This means that if we let I ( t ) 0 (exponentially a.s.), the affectees will be removed from the community with a probability of one. Besides:
lim t S ( t ) = Λ ( d + ϕ 2 ) , lim t I ( t ) = 0 , lim t I P ( t ) = 0 , lim t I N ( t ) = 0 , lim t R ( t ) = 0 , lim t D ( t ) = Λ ϕ 2 d ( Λ + ϕ 2 ) .
Proof. 
For the proof of this theorem, see Theorem 6 in [26]. □

4. The Stationary Distribution of the Disease

Since there are no endemic equilibria in stochastic systems, as a result, the stability analysis cannot be utilized to investigate the disease’s permanence. As a result, one must focus on the uniqueness/existence hypothesis of the stationary distribution, which, in some ways, will aid in the disease’s survival. We use a well-known finding of Khasminskii [27] to illustrate this point.
Let X ( t ) be a Markov process and regular time and homogeneous over R d :
d G ( t ) = b G d t + r = 1 k σ r ( G ) d B r ( t ) .
The diffusion matrix is defined as follows:
A ( x ) = ( b i k ( x ) ) , b i k ( t ) = r = 1 K σ r i ( t ) σ r k ( x ) .
Lemma 2
([1]). The Markov processes G ( t ) has an ergodic stationary distribution Π ( · ) in an open domain U R d that is unique, bounded, and has a regular boundary, while its closure U ¯ R d satisfies the following properties:
1. 
In the set U and its neighbor thereof, the eigenvalue of the diffusion matrix A ( t ) that has the smallest magnitude is bounded away from the origin;
2. 
If x is in R d \ U , the average time τ in which a path starts from x reaching U is < , and sup x K E x τ < K R n where K is compact. Besides, for an integrable function f ( · ) w.r.t measure Π, we have:
P { lim T 1 T 0 T f ( G X ( t ) ) d t = R d f ( x ) Π ( d x ) } = 1 .
Define a parameter:
R 0 s = μ β α ξ + μ + σ 2 2 2 α + μ + σ 3 2 2 δ + μ + σ 4 2 2 .
Theorem 3.
The solution X of System (1) is an ergodic one, and there must exist a stationary distribution π ( . ) when R 0 S > 1 .
Proof. 
To prove Condition (2) of Lemma 2, we need to formulate a positive C 2 -function V : R + 6 R + . To do so, we consider the following functions:
V 1 = F c + R + I P + I + D + I N c 1 ln F c c 2 ln I ,
where c 1 and c 2 are positive constants to be determined later. By the Itô’s formula and system (1), we can get
Ł ( S + V + E + I + R ) = Λ d ( F c + I + I P + I N + R + D ) , Ł ( ln F c ) = Λ F c + d + ϕ 2 + β I N + η 1 2 2 , Ł ( ln I ) = β F c N + ( d + τ + δ + ψ ) + η 2 2 2 , Ł ( ln I P ) = τ I ( t ) I P + ( d + γ 1 + ϕ 1 ) + η 3 2 2 , Ł ( ln I N ) = ψ I ( t ) I N + ( d + σ + γ 3 ) + η 4 2 2 . Ł ( ln R ) = ϕ 1 I P R + d δ I R γ 3 I N R + γ 2 + η 5 2 2 , Ł ( ln D ) = γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + d + η 6 2 2 ,
Therefore, we have
Ł V 1 = Λ d ( F c + I P + I + R + I N + D ) c 1 Λ F c + c 1 d + c 1 ϕ 2 + c 1 β I N + c 1 η 1 2 2 c 2 β F c N + c 2 ( d + δ + τ + ψ ) + c 2 η 2 2 2 , = Λ d ( F c + I + I N + I P + D + R ) c 1 Λ F c + c 1 β I N + c 1 ( d + ϕ 2 + η 1 2 2 ) c 2 β F c N + c 2 ( d + τ + δ + ψ + η 2 2 2 ) . 3 d ( F c + I P + I + R + I N + D ) c 1 Λ F c c 2 β F c N 1 4 + c 1 β I N + Λ + c 1 ( d + ϕ 2 + η 1 2 2 ) + c 2 ( τ + d + ψ + δ + η 2 2 2 ) ,
Let
c 1 ( d + ϕ 2 + η 1 2 2 ) = c 2 ( d + τ + δ + ψ + η 2 2 2 ) = Λ
Namely,
c 1 = Λ d + ϕ 2 + η 1 2 2 c 2 = Λ d + τ + δ + ψ + η 2 2 2
Consequently:
L V 1 3 d c 1 c 2 Λ β d + ϕ 2 + η 1 2 2 d + τ + δ + ψ + η 2 2 2 1 3 Λ + c 1 β I N , L V 1 4 Λ ( R 0 S ) 1 / 4 1 + c 1 β I N
In addition, we obtain:
V 2 = c 3 ( F c + I P + I + R + I N + D c 1 ln F c c 2 ln I ) ln F c ln I P ln I N ln D ln R + F c + I P + I + R + I N + D = ( c 3 + 1 ) ( F c + I P + I + R + I N + D ) ( c 1 c 3 + 1 ) ln F c c 2 c 3 ln I ln I P ln I N ln D ln R .
where c 3 > 0 is constant which will be determined later. It is helpful to show that:
lim inf ( F c , I , I P , I N , R , D ) R + 6 \ U k V 2 ( X ) = + , as k ,
where U k = ( 1 k , k ) × ( 1 k , k ) × ( 1 k , k ) . Next, we must demonstrate that V 2 ( F c , I , I P , I N , R , D ) has a unique minimum value V 2 ( F c ( 0 ) , I ( 0 ) , I P ( 0 ) , I N ( 0 ) , R ( 0 ) , D ( 0 ) ) .
The partial derivative of V 2 ( F c , I , I P , I N , R , D ) with respect to F c , I , I P , I N , R , and D is as follows:
V 2 ( X ) F c = 1 + c 3 1 + c 1 c 3 F c , V 2 ( X ) I = 1 + c 3 c 2 c 3 I , V 2 ( X ) I P = 1 + c 3 1 I P , V 2 ( X ) I N = 1 + c 3 1 I N , V 2 ( X ) R = 1 + c 3 1 R , V 2 ( X ) D = 1 + c 3 1 D .
One can notice that V 2 has one and only one point of stagnation:
( F c , I , I P , I N , R , D ) ( 0 ) = 1 + c 3 c 1 1 + c 3 , 1 c 3 + 1 , c 2 c 3 c 3 + 1 , c 3 c 2 c 3 + 1 , 1 1 + c 3 , 1 c 3 + 1 .
Furthermore, the Hessian matrix of V 2 ( X ) at ( F c , I , I P , I N , R , D ) ( 0 ) is:
B = 1 + c 1 c 3 ( F c ) 2 ( 0 ) 0 0 0 0 0 0 c 2 c 3 I 2 ( 0 ) 0 0 0 0 0 0 1 ( I P ) 2 ( 0 ) 0 0 0 0 0 0 1 ( I N ) 2 ( 0 ) 0 0 0 0 0 0 1 R 2 ( 0 ) 0 0 0 0 0 0 1 D 2 ( 0 ) .
It is clear that the Hessian matrix is positive definite. Hence, V 2 ( X ) has a minimum value V 2 ( X ( 0 ) ) . Further, by using Equation (11) and utilizing the fact that V 2 ( X ) is continuous, we can say that V 2 ( X ) has one and only one minimum value V 2 ( X ( 0 ) ) inside R + 6 . Furthermore, we have to chose another positive C 2 -function V : R + 6 R + as follows:
V ( X ) = V 2 ( X ) V 2 ( X ( 0 ) ) .
Applying the I t o s formula of model (1), we can get and keeping in view the model, we have:
V c 3 3 Λ ( R s 0 ) 1 / 3 1 + c 1 β I N γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + d + η 6 2 2 + Λ d N Λ F c + d + ϕ 2 + β I N + η 1 2 2 τ I ( t ) I P ( d + γ 1 + ϕ 1 ) + η 3 2 2 ψ I ( t ) I N ( d + σ + γ 3 ) + η 4 2 2 ϕ 1 I P R + d δ I R γ 3 I N R + γ 2 + η 5 2 2 ,
A direct consequence of the above equation is given as:
Ł V c 3 c 4 + c 1 β I N γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + Λ Λ F c + 5 d + ϕ 2 + β I N τ I ( t ) I P + γ 1 + ϕ 1 ψ I ( t ) I N + σ + γ 3 ϕ 1 I P R δ I R γ 3 I N R + γ 2 + η 1 2 + η 3 2 + η 4 2 + η 5 2 + η 6 2 2 d N
where:
c 4 = { 3 Λ ( R 0 s ) 1 / 3 1 > 0 .
To move further, we define a set in the form:
D = ϵ 1 < F c < 1 ϵ 2 , ϵ 1 < I < 1 ϵ 2 , ϵ 1 < I P < 1 ϵ 2 , ϵ 1 < I N < 1 ϵ 2 , ϵ 1 < R < 1 ϵ 2 , ϵ 1 < D < 1 ϵ 2 ,
for ϵ i > 0 where ( i = 1 , 2 , , 10 ) are very small positive real numbers to be calculated at later stages. To make the process simpler, we make a partition of the space R + 4 \ D in the form of:
D 1 = X R + 6 , 0 < F c ϵ 1 , D 2 = X R + 6 , 0 < I ϵ 1 , S > ϵ 2 , D 3 = X R + 6 , 0 < I P ϵ 1 , V > ϵ 2 , D 4 = X R + 6 , 0 < I N ϵ 1 , E > ϵ 2 , D 5 = X R + 6 , 0 < R ϵ 1 , I > ϵ 2 , D 6 = X R + 6 , 0 < D ϵ 1 , I > ϵ 2 , D 7 = X R + 6 , F c 1 ϵ 2 , D 8 = X R + 6 , I 1 ϵ 2 , D 9 = X R + 6 , I P 1 ϵ 2 , D 10 = X R + 6 , I N 1 ϵ 2 , D 11 = X R + 6 , R 1 ϵ 2 , D 12 = X R + 6 , D 1 ϵ 2 ,
Next, it is worth showing that V ( F c , I , I P , I N , R , D ) < 0 on R + 6 \ D is equivalent to proving the result in the above ten sub-regions.
Case 1.
For any X D 1 , then by Equation (15), we obtain:
Ł V c 3 c 4 + c 1 β I N γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + Λ Λ F c + 5 d + ϕ 2 + β I N τ I ( t ) I P + γ 1 + ϕ 1 ψ I ( t ) I N + σ + γ 3 ϕ 1 I P R δ I R γ 3 I N R + γ 2 + η 1 2 + η 3 2 + η 4 2 + η 5 2 + η 6 2 2 d N , c 3 c 4 + c 1 β I N γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + Λ Λ F c + 5 d + ϕ 2 + β I N τ I ( t ) I P + γ 1 + ϕ 1 ψ I ( t ) I N + σ + γ 3 + γ 2 + η 1 2 + η 3 2 + η 4 2 + η 5 2 + η 6 2 2 d ϵ 1 .
Choosing ϵ 1 > 0 , we obtain L V ˜ 1 for every ( F c , I , I P , I N , R , D ) D 1 .
Moreover, by using the inequalities in Equation (15), with the same method used for the proof of Case (1), we conclude that L V 1 for all X D 2 , X D 3 , X D 4 , X D 5 , and X D 6 .
Case 2.
For any X D 7 , then by Equation (15), we obtain:
Ł V c 3 c 4 + c 1 β I N γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + Λ Λ F c + 5 d + ϕ 2 + β I N τ I ( t ) I P + γ 1 + ϕ 1 ψ I ( t ) I N + σ + γ 3 ϕ 1 I P R δ I R γ 3 I N R + γ 2 + η 1 2 + η 3 2 + η 4 2 + η 5 2 + η 6 2 2 d N , Ł V c 3 c 4 + c 1 β I N γ 2 R D γ 1 I P D ϕ 2 F c D σ I N D + Λ Λ F c + 5 d + ϕ 2 + β I N τ I ( t ) I P + γ 1 + ϕ 1 + σ + γ 3 + γ 2 + η 1 2 + η 3 2 + η 4 2 + η 5 2 + η 6 2 2 Λ ϵ 2 ,
Assuming the smallest possible value of ϵ 2 > 0 , L V ˜ 1 for every X D 7 .
Moreover, with a similar approach used to prove Case (2) and using Equation (15), we obtain L V 1 , which can be derived for each X D 8 , X D 9 , X D 10 , X D 11 , X D 11 .
In conclusion, we arrive at the final result that there must exist W > 0 , such that
V ( X ) < W < 0 X R + 6 \ D .
Hence
d V ( X ) < W d t + [ ( 1 + c 3 ) F c ( 1 + c 3 c 1 ) η 1 ] d B 1 ( t ) + [ ( 1 + c 3 ) I c 2 c 3 η 2 ] d B 2 ( t ) + [ ( c 3 + 1 ) I p η 3 ] d B 3 ( t ) + [ ( c 4 + 1 ) I N η 4 ] d B 4 ( t ) + [ ( c 3 + 1 ) R η 5 ] d B 5 ( t ) + [ ( c 3 + 1 ) D η 6 ] d B 6 ( t ) .
Assume that ( ( F c , I , I P , I N , R , D ) ( 0 ) ) = ( x 1 , x 2 , x 3 , x 4 , x 5 , x 6 ) = x R + 6 \ D , and τ x is the time at which a path starting from x approaches set D,
τ n = i n f { t : n = | X ( t ) | } and τ ( n ) ( t ) = min { τ x , t , τ n } .
Integrating Relation (18) with limits zero and τ ( n ) ( t ) , assuming the concept expectation, and finally, applying Dynkin’s formula, the following result my be obtained:
E V ( X ( τ ( n ) ( t ) ) ) V ( x ) = E 0 τ ( n ) ( t ) L V ( F c ( u ) , I ( u ) , I P ( u ) , I P ( u ) , R ( u ) , D ( u ) ) d u , E 0 τ ( n ) ( t ) W d u = W E τ ( n ) ( t ) .
As the function V ( x ) is positive, therefore:
E τ ( n ) ( t ) 1 W V ( x ) .
By considering the results of Theorem 3, we have P { τ e = } = 1 . Besides, from the regularity of System (1), we have that, by letting n , t , we have τ ( n ) ( t ) τ x almost surely.
Utilizing the well-known lemma of Fatou yields:
E τ ( n ) ( t ) 1 W V ( x ) <
Clearly, sup x K E τ x < , where K is compact, as discussed already. Consequently, this proves Condition (2) of Lemma 2.
It also worth mentioning that the diffusion matrix of Model (1) is of the form:
B = η 1 2 ( F c ) 2 0 0 0 0 0 0 η 2 2 i 2 0 0 0 0 0 0 η 3 2 ( I P ) 2 0 0 0 0 0 0 η 4 2 ( I N ) 2 0 0 0 0 0 0 η 5 2 R 2 0 0 0 0 0 0 η 6 2 D 2
Choosing M = min X D ¯ R + 6 { η 1 2 ( F c ) 2 , η 2 2 I 2 , η 3 2 ( I P ) 2 , η 4 2 ( I N ) 2 , η 5 2 R 2 , η 6 2 D 2 } , we obtain:
i , j = 1 5 a i j X ξ i ξ j = η 1 2 ( F c ) 2 ξ 2 + η 2 2 I 2 ξ 2 2 + η 3 2 ( I P ) 2 ξ 2 + η 4 2 ( I N ) 2 ξ 4 2 I 2 η 5 2 R 2 ξ 5 + η 6 2 D 2 M | ξ | 2 , X D ¯ ,
where ξ = ( ξ 1 , ξ 2 , ξ 3 , ξ 4 , ξ 5 , ξ 6 ) R + 6 . This means that, the criteria (1) of Lemma 2 hold. According to the above-mentioned results, one can easily conclude that Lemma 2 shows the ergodicity of Model (1), as well as its stationary distribution.

5. Stochastic Optimal Control

The use of optimal control theory techniques is becoming increasingly important, especially in the field of mathematical modeling. These techniques may be utilized to discover an effective control approach by evaluating appropriate control functions [1,3]. Control theory includes a large body of literature, and these approaches are widely used in economics, mathematical physics, and dynamical system theory [3]. The readers are advised to consult the book for the derivation of different optimality conditions and other well-known characteristics of control theory [1,3].
This section discusses the proposed model’s comparable control system. Taking into account the model’s specified hypotheses at the start of the study, as well as the control measures u 1 ( t ) , u 2 ( t ) , u 3 ( t ) and u 4 ( t ) , we create the randomized control version process that operates with it.
The control measures to be included in System (1) have the properties given below:
  • The control measure u 1 ( t ) shows physically the discussion about non-singular and singular kernels on social research forums such as Google scholar and Researchgate for example "https://pubpeer.com/";
  • The variable u 2 ( t ) describes the size of papers published/accepted, the books, etc., about non-singular and singular kernels;
  • The variable u 3 stands for the qualitative aspects such as the fairness of publishers and the editorial board;
  • The measure u 4 ( t ) denotes a conference presenting good talks on the subject of fractional derivatives and integrals.
The study’s main goal is to figure out a proper control approach for limiting the number of infectious and vulnerable persons while increasing the number of those who survive. The equivalent stochastic control variant of System (1) assumes the following form when the control variables are taken into account:
d F c = Λ d F c ( t ) u 1 ( t ) F c ( t ) u 3 ( t ) F c ( t ) ϕ 2 F c ( t ) β F c ( t ) I ( t ) N d t + η 1 F c ( t ) d B 1 ( t ) , d I ( t ) = β F c ( t ) I ( t ) N u 4 ( t ) I ( t ) ( d + τ + δ + ψ ) I ( t ) d t + η 2 I ( t ) d B 2 ( t ) , d I P ( t ) = τ I ( t ) ( d + ϕ 1 + γ 1 + u 3 ( t ) ) I P ( t ) d t + η 3 I P ( t ) d B 3 ( t ) , d I N ( t ) = ψ I ( t ) u 2 ( t ) I N ( t ) ( d + σ + γ 3 ) I N ( t ) d t + η 4 I N ( t ) d B 4 ( t ) , d R ( t ) = [ ϕ 1 I P ( t ) + u 3 ( t ) I P ( t ) d R ( t ) + δ I ( t ) + u 4 ( t ) I ( t ) + u 1 ( t ) F c ( t ) + u 2 ( t ) I N ( t ) , + γ 3 I N ( t ) γ 2 R ( t ) ] d t + η 5 R ( t ) d B 5 ( t ) , d D ( t ) = γ 2 R ( t ) + γ 1 I P ( t ) + ϕ 2 F c ( t ) + σ I N ( t ) d D ( t ) d t + η 6 d B 6 ( t ) .
where the conditions at time t = 0 are:
F c ( 0 ) > 0 , I ( 0 ) 0 , 0 I P ( 0 ) , I ( 0 ) 0 , 0 < I N ( 0 ) , R ( 0 ) > 0 , D ( 0 ) > 0 .
To keep things simple, we use a vector of the type:
x ( t ) = [ x 1 , x 2 , x 3 , x 4 , x 5 , x 6 ] , u ( t ) = [ u 1 , u 2 , u 3 , u 4 ] .
and:
d x ( t ) = f ( x , u ) d t + g ( x ) d w ( t ) ;
here, the functions x i and u i are both time dependent. The conditions at t = 0 in terms of x i may be written as:
x ( 0 ) = [ ( x 1 , x 2 , x 3 , x 4 , x 5 , x 6 ) ( 0 ) ] = x 0 .
The functions in Relation (21) (i.e., f and g) are vectors having components as:
f 1 = Λ d F c ( t ) u 1 ( t ) F c ( t ) ϕ 2 F c ( t ) β F c ( t ) I ( t ) N d t + η 1 F c ( t ) d B 1 ( t ) , f 2 = β F c ( t ) I ( t ) N u 4 ( t ) I ( t ) ( d + τ + δ + ψ ) I ( t ) d t + η 2 I ( t ) d B 2 ( t ) , f 3 = τ I ( t ) ( d + γ 1 + ϕ 1 + u 3 ( t ) ) I P ( t ) d t + η 3 I P ( t ) d B 3 ( t ) , f 4 = ψ I u 2 ( t ) I N ( t ) ( d + σ + γ 3 ) I N ( t ) d t + η 4 I N ( t ) d B 4 ( t ) , f 5 = [ ϕ I P ( t ) d R ( t ) + δ I ( t ) + u 4 ( t ) I ( t ) + u 1 ( t ) F c ( t ) + u 2 ( t ) I N ( t ) , + u 3 ( t ) I P ( t ) + γ 3 I N ( t ) γ 2 R ( t ) ] d t + η 5 R ( t ) d B 5 ( t ) , f 6 = γ 2 R ( t ) + γ 1 I P ( t ) + ϕ 2 F c ( t ) + σ I N ( t ) d D ( t ) d t + η 6 d B 6 ( t ) .
g 1 = η 1 F c , g 2 = η 2 I , g 3 = η 3 i P , g 4 = η 4 I N , g 5 = η 5 R , g 6 = η 5 D . Keeping in view the usefulness of the quadratic terms in the objective functional, we assumed the following functional:
G ( u ) = 1 2 E { 0 t f S 1 F c + S 2 I + S 3 I N + S 4 I P S 5 R + X 2 u 1 2 ( t ) 2 + X 2 u 2 2 ( t ) 2 + X 3 u 3 2 ( t ) 2 + X 4 u 4 2 ( t ) 2 d t + k 1 2 ( F c ) 2 + k 2 2 I 2 + k 3 2 ( I P ) 2 k 4 2 ( I N ) 2 + k 5 2 R 2 + k 6 2 D 2 } ,
where S i , X j , and k i for i = 1 , , 5 and j = 1 , , 4 are constants and strictly greater than zero. The objective of the present work is to find a control vector u * ( t ) = ( u 2 * ( t ) , u 1 * ( t ) , u 3 * ( t ) , u 4 * ( t ) ) that has the following property:
J ( u * ) J ( u ) , for all u U ;
here, the set U denotes the control set (admissible) and is defined as:
U = { u i ( t ) | u i ( t ) [ 0 , u i max ] , t ( 0 , t f ] , u i L 2 [ 0 , t f ] where i = 1 , , 4 } ,
where the terms u i m a x are positive real numbers. To use the stochastic maximal criterion, we must first determine the Hamiltonian H m ( p , q , u , x ) for the problem in operation:
H ( x , u , p , q ) = g ( x ) , q l ( x , u ) + f ( x , u ) , p ;
here, · , · stand for the inner product in the Euclidean space, whereas p = [ p 1 , p 2 , p 3 , p 4 , p 5 , p 6 ] and q = [ q 1 , q 2 , q 3 , q 4 , q 5 , , q 6 ] describe the two independent sets of adjoint variables. Following a similar approach of maximum criteria, we have:
d x * ( t ) = H ( x * , u * , p , q ) p d t + g ( x * ( t ) ) d W ( t ) .
d p * ( t ) = q ( t ) d W ( t ) x H ( x * , u * , p , q ) d t .
H m ( x * , u * , p , q ) = min u U H m ( x * , u * , p . q ) .
Here, the state x * ( t ) represents the optimal path followed by the state variable x ( t ) . The condition at t = 0 on the state and the final condition of the adjoint variables of Equations (24) and (25) are:
x * ( 0 ) = x 0 ,
p ( t f ) = x h ( x * ( t f ) ) ,
respectively. As Equation (26) shows, the optimal value of the control variable u * ( t ) consists of adjoint variables p ( t ) , q ( t ) and state x * ( t ) of the form:
u * ( t ) = ϕ ( x * , q , p ) ;
here, ϕ is computed by Equation (26). Thus, Equations (24) and (25) can be written as:
d x * ( t ) = H ( x * , u * , p , q ) p d t + g ( x * ( t ) ) d W ( t ) .
d p ( t ) = H ( x * , u * , p , q ) x d t + q ( t ) d W ( t ) .
Therefore, the given “Hamiltonian” is:
H = S 1 F c + S 2 I + S 3 I N + S 4 I P S 5 R + X 2 u 1 2 ( t ) 2 + X 2 u 2 2 ( t ) 2 + X 3 u 3 2 ( t ) 2 + X 4 u 4 2 ( t ) 2 d t + k 1 2 ( F c ) 2 + k 2 2 I 2 + k 3 2 ( I P ) 2 k 4 2 ( I N ) 2 + k 5 2 R 2 + k 6 2 D 2 ) + p 1 Λ d F c ( t ) u 1 ( t ) F c ( t ) ϕ 2 F c ( t ) β F c ( t ) I ( t ) N + p 2 β F c ( t ) I ( t ) N u 4 ( t ) I ( t ) ( d + τ + δ + ψ ) I ( t ) + p 3 τ I ( t ) ( d + γ 1 + ϕ 1 + u 3 ( t ) ) I P ( t ) + p 4 ψ I u 2 ( t ) I N ( t ) ( d + σ + γ 3 ) I N ( t ) + p 5 ( ϕ 1 I P ( t ) + u 3 ( t ) I P ( t ) d R ( t ) + δ I ( t ) + u 4 ( t ) I ( t ) + u 1 ( t ) F c ( t ) + u 2 ( t ) I N ( t ) + γ 3 I N ( t ) γ 2 R ( t ) ) + p 6 γ 2 R ( t ) + γ 1 I P ( t ) + ϕ 2 F c ( t ) + d I N ( t ) d D ( t ) + η 1 F c q 1 + η 2 I ( t ) q 2 + η 3 I P ( t ) q 3 + η 4 I N q 4 + η 5 R q 5 + η 6 D q 5 .
The maximum principle of Pontryagin gives us the following adjoint system:
d p 1 d t = S 1 p 1 d + ϕ 2 + u 1 + β I N + p 2 β I N + p 5 u 1 + p 6 ϕ 2 + η 1 q 1 , d p 2 d t = S 2 + p 2 p 1 β F c N p 2 d + τ + δ + ψ + u 4 + p 3 τ + p 4 ψ + p 5 δ + u 4 + η 2 q 2 , d p 3 d t = S 4 p 3 d + γ 1 + ϕ 1 + u 3 + p 5 ϕ 1 + p 5 u 3 + p 6 γ 1 + η 3 q 4 , d p 4 d t = S 3 + p 4 d + σ + γ 3 + u 2 + p 5 γ 3 + u 2 + η 4 q 4 , d p 5 d t = S 5 p 5 d + p 5 γ 2 + p 6 γ 2 + η 5 q 5 , d p 6 d t = d p 6 + η 6 q 6 ,
where p 1 ( t f ) = k 1 F c , p 2 ( t f ) = k 2 I , p 3 ( t f ) = k 3 I N , p 4 ( t f ) = k 4 I P , p 5 ( t f ) = k 5 R , p 6 ( t f ) = k 5 D .
Along with the condition at final time p k t f = 0 for k = 1 , 2 , 3 , 4 , 5 , 6 . Similarly, the supplementary initial and final conditions are:
( F c ) * ( 0 ) = F c ^ , I * ( 0 ) = I ^ , ( I N ) * ( 0 ) = I N ^ , ( I P ) * ( 0 ) = I P ^ , R * ( 0 ) = R ^ D * ( 0 ) = D ^ , p ( t f ) = h ( x * ( t f ) ) x ,
and:
h ( F c , I , I N , I P , R , D ) = k 1 2 ( F c ) 2 + k 2 2 I 2 + k 3 2 ( I N ) 2 k 4 2 ( I P ) 2 + k 5 2 R 2 + k 6 2 D 6 ,
Considering the derivative of function H with respect to the control variables gives us the characterization of the control:
u 1 = F c ( t ) p 1 p 5 2 X 1 , u 2 = I N ( t ) p 4 p 5 2 X 2 , u 3 = I P ( t ) p 3 p 5 2 X 3 , u 4 = I ( t ) p 2 p 5 2 X 4 .
Assuming the lower and upper bounds of the control variables, we have:
u 1 * = min u ˜ 1 , max 0 , F c ( t ) p 1 p 5 2 X 1 , u 2 * = min u ˜ 2 , max 0 , I N ( t ) p 4 p 5 2 X 2 , u 3 * = min u ˜ 3 , max 0 , I P ( t ) p 3 p 5 2 X 3 , u 4 * = min u ˜ 4 , max 0 , I ( t ) p 2 p 5 2 X 4 .

6. Numerical Simulations

We give some numerical simulations in order to illustrate the theoretical results obtained in this paper. We devised a technique for System (1) using the stochastic Runge–Kutta method of fourth order. The discretization transformation takes the form:
F i + 1 c = F i + 1 c + Λ d F i c ϕ 2 F i c β F i c I i N t + η 1 F i c t ζ 1 , i + η 1 2 2 F i c ( ζ 1 , i 2 1 ) t , I i + 1 = I i + 1 + β F i c I i N ( d + τ + δ + ψ ) I i t + η 2 I i t ζ 2 , i + η 2 2 2 I i ( ζ 2 , i 2 1 ) t , I i + 1 P = I i + 1 P + t a u I i ( d + γ 1 + ϕ 1 ) I i P t + η 3 I i P t ζ 3 , i + η 3 2 2 I i P ( ζ 3 , i 2 1 ) t , I i + 1 N = I i + 1 N + ψ I i ( d + σ + γ 3 ) I i N t + η 4 I i N t ζ 4 , i + η 4 2 2 I i N ( ζ 4 , i 2 1 ) t , R i + 1 = R i + 1 + ϕ 1 I i P d R i + δ I i + γ 3 I i N γ 2 R i t + η 5 R i t ζ 5 , i + η 5 2 2 R i ( ζ 5 , i 2 1 ) t , D i + 1 = D i + 1 + γ 2 R i + γ 1 I i P + ϕ 2 F i c + σ I N ( t ) d D i t + η 6 D i t η 5 , i + α 5 2 2 D i ( ζ 5 , i 2 1 ) t .
In the numerical scheme, the terms ζ i , j for i , j = 1 , , 5 are reserved for the Gaussian disturbances, and it has the property N ( 0 , 1 ) where Δ t is a uniform time step size. Let us choose the initial value and the parameter values from Table 1. The desired time period for simulation is considered to be [ 0 , 140 ] .

6.1. Numerical Simulations for the Stationary Distribution and Extinction

To show the numerical illustration of Theorem 2, we chose the values in Table 1 ( V 1 ). Thus, the theorem suggests that if one keeps R 0 < 1 , then the disease will be eventually eliminated from the population with probability one. These facts are supported by Figure 1 as the parameter values in this case fulfill the hypothesis of the theorem. That is to say, large noises can lead the disease to extinction, which is a different phenomenon from its corresponding deterministic model (1).
To numerically show the results suggested by Theorem 3, we take into consideration another set of parameter values in Table 1 ( V 2 ). This set of parameters guarantees that R 0 s > 1 . When the model’s prediction is plotted (see Figure 2), it clearly shows the persistence of the infection with respect to the subject matter. By simulating the model nearly one thousand times, we calculated the average infection extinction time. It was observed that varying the noises’ intensity yielded different times of extinction. Further, it is noted that there exists an inverse relationship between the noise intensity and extinction time. We present the histograms depicting the probability of each class in Figure 3, which shows the stationary distribution as density functions, which are close to each other.
Example 1.
From the assumed values of Table 1 ( V 1 ), we calculated R 0 , which is less than one, and thus, by the statement of Theorem 2, the solution of the system has the following properties:
lim t sup log I ( t ) t 0 , a . s .
and:
lim t sup log I P ( t ) t 0 , a . s .
lim t sup log I N ( t ) t 0 , a . s .
These relations indicate that the dislike of fractional calculus will be eventually eliminated out of the population, and the same is validated via Figure 1.
Example 2.
By following the same procedure, we assumed the values of the parameters from Table 1 ( V 2 ), which satisfy the hypothesis of Theorem 3, that is R 0 s greater than unity. The theorem conclusions, as well as Figure 2 and Figure 3 show that under such circumstances, the dislike of the subject matter will persist in the community.

6.2. Numerical Simulations for Stochastic Optimality

In this section, we give some numerical simulations in order to illustrate the theoretical result of the optimal control theory of the stochastic model. The method used is again the well-known RK4 method. For this simulation, the authors assumed the values of the parameters and noise intensity from Table 1 ( S 3 ). The figures that are of great concern are Figure 4 and Figure 5. In Figure 4a,c,e, we plot the solution curves F c ( t ) , I ( t ) , and I P ( t ) both in the presence and absence of control measures in the case of the associated deterministic system. The objective of the work is well explained by the figures. Figure 5 represents the graph of ( I N ( t ) , R ( t ) , D ( t ) ) both with and without controls. Figure 5a,c,e, represents the graphs of r ( I N ( t ) , R ( t ) , D ( t ) ) with and without the (optimal) control of the corresponding deterministic version of Model (1). Furthermore, Figure 4c shows the corresponding deterministic part of System (1). The effect of the control variables is clearly shown in the figures.
Besides, we used the stochastic RK4 method and simulated the proposed system (1). The numerical solution of the optimal system was obtained using the scheme that discretizes the state SDEs, the system of adjoint equations, while imposing the terminal conditions. The state system (1) was approximated with the help of RK4 techniques. After that, the same procedure was used to find a solution to the adjoint system (33) utilizing the final conditions (34). By using convex combinations of the already-obtained control measures and the formulas for the characterization of the control measures, the variables were updated. The process was repeated recursively until the desired accuracy was achieved. The figures suggest that if we impose all of the control, a drastic decrease in the I class may be observed; for instance, one can see the effect in Figure 4, Figure 5 and Figure 6.
The computed results for optimal control under the strategies are reflected in Figure 4b,d,f, and one can observe that the control strategy resulted in a decrease in the number of ( I ( t ) , I P ( t ) ) and the control strategy result of F C ( t ) ) increasing for Model (1), while in Figure 5b,d,f, we can clearly observe the differences between the two cases with and without controls for the ( I N ( t ) , R ( t ) , D ( t ) ) of Model (1). Furthermore, Figure 6 reflects the dynamics of the model’s best control variables (1). Obviously, there is a difference between these two instances.

7. Conclusions and Prediction

A divide arose as a result of the development of the fractional differential operator with non-singular kernels, which marked a significant milestone in this discipline. With the introduction of these novel integral and differential operators, potential ways for both the application and theory of fractional calculus to real-world phenomenon have opened up. Some scholars, however, anticipated differential operators with non-singular kernels to follow the characteristics of classical derivatives in the same way as fractional derivatives with singular kernels do, which is impossible because the two operators are not employed for the very same reasons. This article looked at the asymptotic behavior of a probabilistic model with nonlinear white noise disturbances. Initially, we calculated the threshold R s 0 and then observed that model (1) has a stationary distribution, which is naturally ergodic. This indicates that the illness will continue to exist. This shows that stochastic disturbance can prevent negativity from spreading.
We came up with some interesting findings on the stochastic version of the optimal system control (1). The study also included suggestions for several types of control techniques. These approaches can support us in eradicating hostility in the community. To confirm the obtained analytical results, many numerical simulations were performed using the stochastic Runge–Kutta technique of fourth order. This idea might serve as a solid foundation for research into comparable illnesses, with significant implications in biologicalsciences. In addition, our suggested theory may be used to study other communicable infections including HIV, COIVD-19, and tuberculosis (TB), to name a few.

Author Contributions

Conceptualization, A.D.; Methodology, A.D.; Validation, Q.T.A.; Formal analysis, A.D., Q.T.A.; Investigation, A.D., Q.T.A.; Resources, A.D.; Data curation, Q.T.A.; Writing—original draft preparation, A.D.; writing—review and editing, A.D., Q.T.A.; Software, A.D.; Visualization, A.D., Q.T.A.; Supervision, A.D.; Project administration, A.D.; Funding acquisition, A.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was sponsored by the Fundamental Research Funds for the Central Universities, Sun Yat-sen University (Grant No. 22qntd2802).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated nor analyzed during the current study.

Acknowledgments

The authors would like to acknowledge the sponsorship from the Fundamental Research Funds for the Central Universities, Sun Yat-sen University (Grant No. 22qntd2802).

Conflicts of Interest

No conflict of interest exists regarding the publications of this paper.

References

  1. Din, A.; Li, Y. Stationary distribution extinction and optimal control for the stochastic hepatitis B epidemic model with partial immunity. Phys. Scr. 2021, 96, 074005. [Google Scholar] [CrossRef]
  2. Li, Y.X.; Rauf, A.; Naeem, M.; Binyamin, M.A.; Aslam, A. Valency-based topological properties of linear hexagonal chain and hammer-like benzenoid. Complexity 2021, 2021, 9939469. [Google Scholar] [CrossRef]
  3. Agarwal, R.P.; Badshah, Q.; Rahman, G.; Islam, S. Optimal control and dynamical aspects of a stochastic pine wilt disease model. J. Frankl. Inst. 2019, 356, 3991–4025. [Google Scholar] [CrossRef]
  4. Chen, S.-B.; Soradi-Zeid, S.; Jahanshahi, H.; Alcaraz, R.; Gómez-Aguilar, J.F.; Bekiros, S.; Chu, Y.-M. Optimal control of time-delay fractional equations via a joint application of radial basis functions and collocation method. Entropy 2020, 22, 1213. [Google Scholar] [CrossRef] [PubMed]
  5. Li, X.-P.; al Bayatti, H.; Din, A.; Zeb, A. A vigorous study of fractional order COVID-19 model via ABC derivatives. Results Phys. 2021, 29, 104737. [Google Scholar] [CrossRef]
  6. Chen, S.-B.; Rajaee, F.; Yousefpour, A.; Alcaraz, R.; Chu, Y.-M.; Gómez-Aguilar, J.F.; Bekiros, S.; Aly, A.A.; Jahanshahi, H. Antiretroviral therapy of HIV infection using a novel optimal type-2 fuzzy control strategy. Alex. Eng. J. 2021, 60, 1545–1555. [Google Scholar] [CrossRef]
  7. Din, A.; Li, Y. Controlling heroin addiction via age-structured modeling. Adv. Differ. Equ. 2020, 2020, 521. [Google Scholar] [CrossRef]
  8. Chen, S.-B.; Rashid, S.; Hammouch, Z.; Noor, M.A.; Ashraf, R.; Chu, Y. Integral inequalities via Raina’s fractional integrals operator with respect to a monotone function. Adv. Differ. Equ. 2020, 2020, 647. [Google Scholar] [CrossRef]
  9. Ye, W.; Khan, M.A.; Alshahrani, M.Y.; Muhammad, T. A dynamical study of SARS-COV-2: A study of third wave. Results Phys. 2021, 29, 104705. [Google Scholar]
  10. Din, A.; Khan, A.; Zeb, A.; Ammi, M.R.S.; Tilioua, M.; Torres, D.F.M. Hybrid Method for Simulation of a Fractional COVID-19 Model with Real Case Application. Axioms 2021, 10, 290. [Google Scholar] [CrossRef]
  11. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel, Theory and Application to Heat Transfer Model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  12. Caputo, M.; Fabrizio, M. On the notion of fractional derivative and applications to the hysteresis phenomena. Meccanica 2017, 52, 3043–3052. [Google Scholar] [CrossRef]
  13. Caputo, M. Linear model of dissipation whoseQ is almost frequency independent. II. Geophys. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  14. Leibniz, G.W. Letter from Hanover, Germany to G.F.A. L’Hospital, September 30, 1695. Math. Schr. 1849, 2, 301–302. [Google Scholar]
  15. Abel, N.H. Oplösning af et par opgaver ved hjelp af bestemte integraler (Solution de quelques problèmes à l’aide d’int égrales définies, Solution of a couple of problems by means of definite integrals). Mag. Nat. 1823, 2, 55–68. [Google Scholar]
  16. Depnath, L. A brief historical introduction to fractional calculus. Int. J. Math. Educ. Sci. Technol. 2004, 35, 487–501. [Google Scholar] [CrossRef]
  17. Chen, S.-B.; Beigi, A.; Yousefpour, A.; Rajaee, F.; Jahanshahi, H.; Bekiros, S.; Martínez, R.A.; Chu, Y. Recurrent neural network-based robust non-singular sliding mode control with input saturation for a non-holonomic spherical robot. IEEE Access 2020, 8, 188441–188453. [Google Scholar] [CrossRef]
  18. Shen, Z.-H.; Chu, Y.; Khan, M.A.; Muhammad, S.; Al-Hartomy, O.A.; Higazy, M. Mathematical modeling and optimal control of the COVID-19 dynamics. Results Phys. 2021, 31, 105028. [Google Scholar] [CrossRef]
  19. Saima, R.; Noor, M.A.; Ashraf, R.; Chu, Y.-M. A new approach on fractional calculus and probability density function. AIMS Math. 2020, 5, 7041–7054. [Google Scholar]
  20. Rashid, S.; Sultana, S.; Karaca, Y.; Khalid, A.; Chu, Y.-M. Some further extensions considering discrete proportional fractional operators. Fractals 2021, 2240026. [Google Scholar] [CrossRef]
  21. Balzotti, C.; D’Ovidio, M.; Loreti, P. Fractional SIS epidemic models. Fractal Fract. 2020, 4, 44. [Google Scholar] [CrossRef]
  22. Almeida, R.; Qureshi, S. A fractional measles model having monotonic real statistical data for constant transmission rate of the disease. Fractal Fract. 2019, 3, 53. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, Y.; Chen, Y. Dynamic analysis of the viscoelastic pipeline conveying fluid with an improved variable fractional order model based on shifted Legendre polynomials. Fractal Fract. 2019, 3, 52. [Google Scholar] [CrossRef] [Green Version]
  24. Tateishi, A.A.; Ribeiro, H.V.; Lenzi, E.K. The role of fractional time-derivative operators on anomalous diffusion. Front. Phys. 2017, 5, 52. [Google Scholar] [CrossRef] [Green Version]
  25. Hristov, J. Steady-state heat conduction in a medium with spatial non-singular fading memory: Derivation of Caputo-Fabrizio space-fractional derivative from Cattaneo concept with Jeffreys Kernel and analytical solutions. Therm. Sci. 2017, 21, 827–839. [Google Scholar] [CrossRef] [Green Version]
  26. Atangana, A. Mathematical model of survival of fractional calculus, critics and their impact: How singular is our world? Results Phys. 2021, 31, 105208. [Google Scholar] [CrossRef]
  27. Khasminskii, R. Stochastic Stability of Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011; Volume 66. [Google Scholar]
Figure 1. Simulations of ( F c ( t ) , I ( t ) , I P ( t ) , I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1) with its corresponding deterministic version. (a) F C ( t ) : researchers using fractional differential operators with non-singular kernels. (b) I ( t ) : researchers that have been affected by the harmful papers. (c) I P ( t ) : researchers that have been affected, but still have positive opinions about non-singular kernel derivatives. (d) I N ( t ) : researchers that have been affected and have negative opinions about non-singular kernel derivatives. (e) R ( t ) : researchers that have overcome divisive criticism. (f) D ( t ) : researchers that die, or retire.
Figure 1. Simulations of ( F c ( t ) , I ( t ) , I P ( t ) , I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1) with its corresponding deterministic version. (a) F C ( t ) : researchers using fractional differential operators with non-singular kernels. (b) I ( t ) : researchers that have been affected by the harmful papers. (c) I P ( t ) : researchers that have been affected, but still have positive opinions about non-singular kernel derivatives. (d) I N ( t ) : researchers that have been affected and have negative opinions about non-singular kernel derivatives. (e) R ( t ) : researchers that have overcome divisive criticism. (f) D ( t ) : researchers that die, or retire.
Fractalfract 06 00279 g001
Figure 2. Simulations of ( F c ( t ) , I ( t ) , I P ( t ) , I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1) with its corresponding deterministic version. (a) F C ( t ) : researchers using fractional differential operators with non-singular kernels. (b) I ( t ) : researchers that have been affected by the harmful papers. (c) I P ( t ) : researchers that have been affected, but still have positive opinions about non-singular kernel derivatives. (d) I N ( t ) : researchers that have been affected and have negative opinions about non-singular kernel derivatives. (e) R ( t ) : researchers that have overcome divisive criticism. (f) D ( t ) : researchers that die, or retire.
Figure 2. Simulations of ( F c ( t ) , I ( t ) , I P ( t ) , I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1) with its corresponding deterministic version. (a) F C ( t ) : researchers using fractional differential operators with non-singular kernels. (b) I ( t ) : researchers that have been affected by the harmful papers. (c) I P ( t ) : researchers that have been affected, but still have positive opinions about non-singular kernel derivatives. (d) I N ( t ) : researchers that have been affected and have negative opinions about non-singular kernel derivatives. (e) R ( t ) : researchers that have overcome divisive criticism. (f) D ( t ) : researchers that die, or retire.
Fractalfract 06 00279 g002
Figure 3. The probability distribution histogram of ( F c ( t ) , I ( t ) , I P ( t ) , I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1). (a) F C ( t ) : researchers using fractional differential operators with non-singular kernels. (b) I ( t ) : researchers that have been affected by the harmful papers. (c) I P ( t ) : researchers that have been affected, but still have positive opinions about non-singular kernel derivatives. (d) I N ( t ) : researchers that have been affected and have negative opinions about non-singular kernel derivatives. (e) R ( t ) : researchers that have overcome divisive criticism. (f) D ( t ) : researchers that die, or retire.
Figure 3. The probability distribution histogram of ( F c ( t ) , I ( t ) , I P ( t ) , I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1). (a) F C ( t ) : researchers using fractional differential operators with non-singular kernels. (b) I ( t ) : researchers that have been affected by the harmful papers. (c) I P ( t ) : researchers that have been affected, but still have positive opinions about non-singular kernel derivatives. (d) I N ( t ) : researchers that have been affected and have negative opinions about non-singular kernel derivatives. (e) R ( t ) : researchers that have overcome divisive criticism. (f) D ( t ) : researchers that die, or retire.
Fractalfract 06 00279 g003
Figure 4. Simulations of ( F c ( t ) , I ( t ) , I P ( t ) ) for the stochastic model (1) with its corresponding deterministic version, both with and without control. (a) F c ( t ) in the case of the deterministic model both with and without control. (b) F c ( t ) in the case of the stochastic model both with and without control. (c) I ( t ) in the case of the deterministic model both with and without control. (d) I ( t ) in the case of the stochastic model both with and without control. (e) I P ( t ) in the case of the deterministic model both with and without control. (f) I P ( t ) in the case of the stochastic model both with and without control.
Figure 4. Simulations of ( F c ( t ) , I ( t ) , I P ( t ) ) for the stochastic model (1) with its corresponding deterministic version, both with and without control. (a) F c ( t ) in the case of the deterministic model both with and without control. (b) F c ( t ) in the case of the stochastic model both with and without control. (c) I ( t ) in the case of the deterministic model both with and without control. (d) I ( t ) in the case of the stochastic model both with and without control. (e) I P ( t ) in the case of the deterministic model both with and without control. (f) I P ( t ) in the case of the stochastic model both with and without control.
Fractalfract 06 00279 g004
Figure 5. Simulations of ( I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1) with its corresponding deterministic version, both with and without controls. (a) I N ( t ) in the case of the deterministic model both with and without control. (b) I N ( t ) in the case of the stochastic model both with and without control. (c) R ( t ) in the case of the deterministic model both with and without control. (d) R ( t ) in the case of the stochastic model both with and without control. (e) D ( t ) in the case of the stochastic model both with and without control. (f) D ( t ) in the case of the stochastic model both with and without control.
Figure 5. Simulations of ( I N ( t ) , R ( t ) , D ( t ) ) for the stochastic model (1) with its corresponding deterministic version, both with and without controls. (a) I N ( t ) in the case of the deterministic model both with and without control. (b) I N ( t ) in the case of the stochastic model both with and without control. (c) R ( t ) in the case of the deterministic model both with and without control. (d) R ( t ) in the case of the stochastic model both with and without control. (e) D ( t ) in the case of the stochastic model both with and without control. (f) D ( t ) in the case of the stochastic model both with and without control.
Fractalfract 06 00279 g005
Figure 6. Simulations of ( u 1 ( t ) , u 2 ( t ) , u 3 ( t ) , u 4 ( t ) ) for the stochastic model (1) with its corresponding deterministic version. (a) Trajectories of the optimal control of the deterministic system with and without control. (b) Trajectories of the optimal control of the stochastic system with and without control.
Figure 6. Simulations of ( u 1 ( t ) , u 2 ( t ) , u 3 ( t ) , u 4 ( t ) ) for the stochastic model (1) with its corresponding deterministic version. (a) Trajectories of the optimal control of the deterministic system with and without control. (b) Trajectories of the optimal control of the stochastic system with and without control.
Fractalfract 06 00279 g006
Table 1. Values of the parameters for simulating Model (1).
Table 1. Values of the parameters for simulating Model (1).
ParametersUnite V 1 V 2 V 3
Λ Per week0.052.52.8
dPer week1/170.3651/70.3651/70.365
β Per week0.080.90.8
γ 1 Per week0.010.080.8
γ 2 Per week0.040.070.5
γ 3 Per week0.0040.070.004
ϕ 2 Per week0.020.040.04
τ Per week0.0040.0040.005
δ Per week0.010.030.03
σ Per week0.030.10.1
ψ vPer week0.010.010.01
ϕ 1 Per week0.020.0010.01
η 1 Noise intensity0.40.42
η 2 Noise intensity0.30.20.3
η 3 Noise intensity0.40.50.4
η 4 Noise intensity0.50.50.5
η 5 Noise intensity0.60.60.1
η 6 Noise intensity0.20.10.2
F c ( 0 ) Initial value10060430
I ( 0 ) Initial value105010
I P ( 0 ) Initial value064030
I N ( 0 ) Initial value033520
R ( 0 ) Initial value001510
D ( 0 ) Initial value021510
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Din, A.; Ain, Q.T. Stochastic Optimal Control Analysis of a Mathematical Model: Theory and Application to Non-Singular Kernels. Fractal Fract. 2022, 6, 279. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6050279

AMA Style

Din A, Ain QT. Stochastic Optimal Control Analysis of a Mathematical Model: Theory and Application to Non-Singular Kernels. Fractal and Fractional. 2022; 6(5):279. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6050279

Chicago/Turabian Style

Din, Anwarud, and Qura Tul Ain. 2022. "Stochastic Optimal Control Analysis of a Mathematical Model: Theory and Application to Non-Singular Kernels" Fractal and Fractional 6, no. 5: 279. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6050279

Article Metrics

Back to TopTop