Next Article in Journal
Graphs Defined on Rings: A Review
Next Article in Special Issue
Employing a Fractional Basis Set to Solve Nonlinear Multidimensional Fractional Differential Equations
Previous Article in Journal
Multivariate Bayesian Semiparametric Regression Model for Forecasting and Mapping HIV and TB Risks in West Java, Indonesia
Previous Article in Special Issue
Modelling Fractional Advection–Diffusion Processes via the Adomian Decomposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetic Behavior and Optimal Control of a Fractional-Order Hepatitis B Model

School of Mathematics and Physics, Xinjiang Institute of Engineering, Urumqi 830023, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(17), 3642; https://doi.org/10.3390/math11173642
Submission received: 21 May 2023 / Revised: 3 August 2023 / Accepted: 22 August 2023 / Published: 23 August 2023
(This article belongs to the Special Issue Fractional Modeling, Control, Analysis and Applications)

Abstract

:
The fractional-order calculus model is suitable for describing real-world problems that contain non-local effects and have memory genetic effects. Based on the definition of the Caputo derivative, the article proposes a class of fractional hepatitis B epidemic model with a general incidence rate. Firstly, the existence, uniqueness, positivity and boundedness of model solutions, basic reproduction number, equilibrium points, and local stability of equilibrium points are studied employing fractional differential equation theory, stability theory, and infectious disease dynamics theory. Secondly, the fractional necessary optimality conditions for fractional optimal control problems are derived by applying the Pontryagin maximum principle. Finally, the optimization simulation results of fractional optimal control problem are discussed. To control the spread of the hepatitis B virus, three control variables (isolation, treatment, and vaccination) are applied, and the optimal control theory is used to formulate the optimal control strategy. Specifically, by isolating infected and non-infected people, treating patients, and vaccinating susceptible people at the same time, the number of hepatitis B patients can be minimized, the number of recovered people can be increased, and the purpose of ultimately eliminating the transmission of hepatitis B virus can be achieved.
MSC:
92D25; 92D30; 93D20; 34A08; 34D23; 92B05

1. Introduction

Hepatitis B causes liver lesions, which is an infectious disease caused by the hepatitis B virus. Each outbreak will have a profound impact on human life and health, global social stability, and economic development. Globally, approximately 887,000 people die each year from various diseases associated with hepatitis B infection, with cirrhosis and primary hepatocellular carcinoma accounting for 52 percent and 38 percent of deaths, respectively. Hepatocellular carcinoma can occur in 10 to 25 percent of hepatitis B virus carriers.
Mathematical models are powerful tools for studying the transmission patterns of infectious diseases and finding ways to control their prevalence. It can help us know the dynamics of infectious disease spread and provide some new disease control strategies. Here, we focus on the mathematical model of hepatitis B transmission. Considering that the transmission of hepatitis B is affected by many factors, scholars have established various dynamic models of hepatitis B [1,2,3]. For example, in [4], Khan et al. developed the following hepatitis B model with bilinear incidence:
d S ( t ) d t = Λ β S ( t ) I ( t ) μ 0 + ν S ( t ) , d I ( t ) d t = β S ( t ) I ( t ) μ 0 + μ 1 + γ I ( t ) , d R ( t ) d t = γ I ( t ) + ν S ( t ) μ 0 R ( t ) ,
where β is the transmission rate. In [5,6,7,8], the authors proposed several types of stochastic models of hepatitis B with standard incidence rates influenced by environmental fluctuations. The extinction, persistence, and ergodic stability distributions of the hepatitis B model were obtained using stochastic Lyapunov general function theory. Khan et al. [9] established a hepatitis B model with saturation incidence in the following form:
d S ( t ) d t = Λ β S ( t ) I ( t ) 1 + η I ( t ) μ 0 + ν S ( t ) , d I ( t ) d t = β S ( t ) I ( t ) 1 + η I ( t ) μ 0 + μ 1 + γ I ( t ) , d R ( t ) d t = γ I ( t ) + ν S ( t ) μ 0 R ( t ) ,
where β and η , respectively, represent the transmission rate of hepatitis B and saturation rates. Stability theory and central manifold theory were used to study the stability behavior and bifurcation phenomena of the model. Furthermore, the optimal control strategy of hepatitis B transmission was studied by applying the optimal control theory.
In the biomedical field, many biological phenomena, such as biomolecular or cellular interactions, population interactions, microbial cultures, cell growth processes, and human immune processes, exhibit fractal geometry, and global correlation, and have memory genetic effects, so that their kinetic behavior cannot be accurately represented by traditional mathematical models of integer-order microequations. Fractional-order derivatives are closely related to fractal and fractional dimensions and have some unique properties, including global correlation, memorability, and heritability, which make fractional-order differential equation models well suited to describe some complex behaviors and phenomena in nature. Another characteristic of fractional differential equations is that their stability region is larger than that of integer differential equations. This implies that the instability of integer-order differential systems may be stable in fractional differential systems. The fractional-order differential model requires only a few parameters to achieve very good results, solving the defect that the theory of the integer-order differential model does not match well with the experimental results. Therefore, fractional differential equation theory has been widely used to study the modeling of various infectious diseases [10,11,12]. For example, in [13], the authors developed a fractional-order typhoid disease dynamics model. Using the Lyapunov function, the local and global stability of the model was studied. In [14], the authors developed a fractional mathematical model of COVID-19. The uniqueness, boundedness, and non-negativity of the model solutions were discussed. The kinetics and sensitivity analysis were also studied. However, some research methods in integer-order infectious disease dynamics models cannot be directly applied to fractional-order infectious disease dynamics models, and related extensions are being improved. At present, although few differential equations are using fractional-order derivatives to construct models of hepatitis B, more and more scholars have carried out research on fractional-order hepatitis B infectious disease models and achieved certain results. From reading the literature, it is clear that to model the hepatitis B transmission process, some scholars use bilinear incidence rates, standardized incidence rates, and saturation rates. For example, several types of fractional-order hepatitis B transmission models with bilinear incidence rates have been developed in the literature [15,16]. Simelane et al. [17] extended the work of the literature [9] on integer-order hepatitis B models and developed a class of hepatitis B models with Caputo-type fractional-order derivatives as follows:
0 C D t α S t = Λ β S ( t ) I ( t ) 1 + η I ( t ) μ 0 + ν S ( t ) , 0 C D t α I t = β S ( t ) I ( t ) 1 + η I ( t ) μ 0 + μ 1 + γ I ( t ) , 0 C D t α R t = γ I ( t ) + ν S ( t ) μ 0 R ( t ) ,
where 0 < α 1 , 0 C D t α is the left Caputo fractional derivative. The basic properties of the fractional-order model were studied by applying the generalized fractional-order median theorem. Then the stability of the model was studied by applying fractional-order Routh–Hurwitz stability criterion.
In this paper, we investigate a more general incidence than the one mentioned before. Furthermore, we will further extend the work of [17] by developing a class of fractional-order hepatitis B models with a general incidence rate, in the following form:
0 C D t α S t = Λ β S ( t ) I ( t ) f S ( t ) , I ( t ) μ 0 + ν S ( t ) , 0 C D t α I t = β S ( t ) I ( t ) f S ( t ) , I ( t ) μ 0 + μ 1 + γ I ( t ) , 0 C D t α R t = γ I ( t ) + ν S ( t ) μ 0 R ( t ) ,
with
S 0 = S 0 0 , I 0 = I 0 0 , R 0 = R 0 0 ,
where 0 < α 1 , S ( t ) , I ( t ) , and R ( t ) represent susceptible, infected, and recovered persons, respectively. f S , I = 1 + a 1 S + a 2 I + a 3 S I , a 1 , a 2 , a 3 0 . For convenience below, we define:
g S , I = β S I f S , I = β S I 1 + a 1 S + a 2 I + a 3 S I .
where g S , I is the incidence rate. In particular, g ( S , I ) is the general form of mutual interference between S and I, as follows:
(i)
When a 1 = a 2 = a 3 = 0 , g ( S , I ) becomes the bilinear incidence;
(ii)
When a 1 = a 3 = 0 or a 2 = a 3 = 0 , g ( S , I ) becomes the saturation incidence;
(iii)
When a 3 = 0 , g ( S , I ) becomes the Beddington–DeAngelis function;
(iv)
When a 3 = a 1 a 2 , g ( S , I ) becomes the Crowley–Martin function.
The meanings of the parameters in this article are as follows:
  • Λ : the birth rate;
  • β : infection rate from susceptible population to hepatitis B;
  • μ 0 : the natural mortality rate;
  • μ 1 : the death rate of hepatitis B;
  • γ : the recovery rate of hepatitis B
  • ν : hepatitis B vaccine coverage rate.
Next, we make the following assumptions about model (1):
  • ( A 1 ) .   N ( t ) represents the total population and is divided into three components: S ( t ) stands for susceptible person; I ( t ) stands for people infected with hepatitis B in the population at time t; R ( t ) represents individuals who recover from infection and have lifelong immunity.
  • ( A 2 ) . All parameters and state variables of the model are non-negative.
  • ( A 3 ) . Incidence is set as the nonlinear incidence rate.
  • ( A 4 ) . Once successfully vaccinated or cured by treatment, immunity is considered permanent.
The research motivation of this paper is to study a new fractional hepatitis B transmission model with a more general incidence than that studied in the literature above. The dynamic behavior is studied, and the optimal control problem of the fractional-order hepatitis B model is discussed, and thus, the optimal control strategy is given to eliminate the hepatitis B virus. The innovation of this paper lies in the further extension of the work of literature [17]. Ref. [17] studied fractional hepatitis B models with saturation incidence, and this paper studies a fractional hepatitis B model with a more general incidence. Specifically, the work of reference [17] is a special case of this article when a 1 = a 3 = 0 .

2. Fractional-Order Basic Concepts

This section introduces the knowledge of fractional-order calculus required for the paper.
Definition 1
([18]). The Riemann–Liouville type left-hand and right-hand fractional-order integrals are defined in the following form, respectively:
t 0 I t α u ( t ) = 1 Γ ( α ) t 0 t ( t s ) α 1 u ( s ) d s , t I T α u ( t ) = 1 Γ ( α ) t T ( s t ) α 1 u ( s ) d s , α > 0 .
Definition 2
([18]). The Riemann–Liouville type left and right fractional-order derivatives are defined in the following form, respectively:
t 0 D t γ u ( t ) = d d t t 0 D t γ 1 u ( t ) = 1 Γ ( 1 γ ) d d t t 0 t ( t s ) γ u ( s ) d s , 0 γ < 1 ,
t D T γ u ( t ) = d d t t D T γ 1 u ( t ) = 1 Γ ( 1 γ ) d d t t T ( s t ) γ u ( s ) d s , 0 γ < 1 .
Definition 3
([18]). The left-hand and right-hand fractional-order derivatives of the Caputo type are defined in the following form, respectively:
t 0 C D t γ u ( t ) = t 0 D t γ 1 u ( t ) = 1 Γ ( 1 γ ) t 0 t ( t s ) γ u ( s ) d s , 0 < γ < 1 ,
t C D T γ u ( t ) = t D T γ 1 u ( t ) = 1 Γ ( 1 γ ) t T ( s t ) γ u ( s ) d s , 0 < γ < 1 .
Lemma 1
([18]). Let t > 0 , n 1 < γ < n , n N . Then:
t 0 C D t γ u ( t ) = t 0 D t γ u ( t ) j = 0 n 1 u j t 0 Γ j γ + 1 t t 0 j γ ,
t C D T γ u ( t ) = t D T γ u ( t ) j = 0 n 1 u j T Γ j γ + 1 T t j γ .
Thus, if u ( t 0 ) = u ( t 0 ) = = u ( n 1 ) ( t 0 ) = 0 , then t 0 D t γ u ( t ) = t 0 C D t γ u ( t ) ; if u ( T ) = u ( T ) = = u ( n 1 ) ( T ) = 0 , then t D T γ u ( t ) = t C D T γ u ( t ) .
Lemma 2
([19]). For 0 < α 1 , we have t C D T α u ( t ) = 0 C D t α u T t .

3. Kinetic Analysis of Fractional-Order Hepatitis B Model

3.1. Positivity and Boundedness of Solutions

The following necessary citation is given first.
Lemma 3
([20]). If u ( t ) C [ t 0 , T ] and t 0 C D t α u ( t ) C [ t 0 , T ] , then:
u ( t ) = u ( t 0 ) + 1 Γ ( α ) t 0 C D t α u ( s ) t α , 0 < α 1 ,
where t 0 s t , t ( t 0 , T ] .
Lemma 4
([19]). Let 0 < α 1 . Suppose u ( t ) C [ t 0 , T ] and t 0 C D t α u ( t ) C [ t 0 , T ] . If t 0 C D t α u ( t ) 0 , t ( t 0 , T ) , then u ( t ) is a non-decreasing function for t [ t 0 , T ] ; If t 0 C D t α u ( t ) 0 , t ( t 0 , T ) , then u ( t ) is a non-increasing function for t [ t 0 , T ] .
Theorem 1.
Let ( S , I , R ) Ω be one solution of models (1)–(2), where:
Ω = ( S , I , R ) R + 3 : max | S | , | I | , | R | A
for sufficiently large A. Then, the closed set:
Π = S , I , R R + 3 : 0 < N t = S ( t ) + I ( I ) + R ( t ) Λ μ 0
is positively invariant for models (1)–(2).
Proof. 
Step 1. The existence and uniqueness of the solutions of (1) are investigate in the region Ω × ( 0 , + ) . Define a mapping H ( X ) = ( H 1 ( X ) , H 2 ( X ) , H 3 ( X ) ) , where:
H 1 ( X ) = Λ β S I f S , I μ 0 + ν S , H 2 ( X ) = β S I f S , I μ 0 + μ 1 + γ I , H 3 ( X ) = γ I + ν S μ 0 R .
Then:
H ( X ) , H ( X ¯ ) = H 1 ( X ) H 1 ( X ¯ ) + H 2 ( X ) H 2 ( X ¯ ) + H 3 ( X ) H 3 ( X ¯ ) = β S I f S , I μ 0 + ν S + β S ¯ I ¯ f S ¯ , I ¯ + μ 0 + ν S ¯ + β S I f S , I μ 0 + μ 1 + γ I β S ¯ I ¯ f S ¯ , I ¯ + μ 0 + μ 1 + γ I ¯ + γ I + ν S μ 0 R γ I ¯ ν S ¯ + μ 0 R ¯ μ 0 + 2 ν + β a 2 S S ¯ + μ 0 + μ 1 + 2 γ + β a 1 I I ¯ + μ 0 R R ¯ N 1 X X ¯ ,
where N 1 = max μ 0 + 2 ν + β a 2 , μ 0 + μ 1 + 2 γ + β a 1 , μ 0 . Thus, H ( X ) satisfies the Lipschitz condition. Hence, we obtain that the solution of model (1) on ( 0 , ) exists and is unique.
Step 2. we show that Π is one positively invariant set. From model (1), one has:
0 C D t α S t = Λ β S I f S , I μ 0 + ν S β S I 1 + a 1 S + a 2 I + a 3 S I μ 0 + ν S β S I a 2 I μ 0 + ν S β a 2 + μ 0 + ν S = C 1 S ,
where C 1 = β a 2 + μ 0 + ν > 0 . Thus, S ( t ) S ( 0 ) E α , 1 C 1 t α . Since S ( 0 ) 0 , S ( t ) 0 . From model (1), we obtain:
0 C D t α I t = β S I f S , I μ 0 + μ 1 + γ I μ 0 + μ 1 + γ I = C 2 I ,
where C 2 = μ 0 + μ 1 + γ > 0 . Since I ( 0 ) 0 , then I ( t ) I ( 0 ) E α , 1 C 2 t α 0 . From model (1), one has:
0 C D t α R t = γ I + ν S μ 0 R μ 0 R = C 3 R ,
where C 3 = μ 0 > 0 . Since R ( 0 ) 0 , then R ( t ) R ( 0 ) E α , 1 C 3 t α 0 . Thus, it has been proved that the solutions of (1) are non-negative. According to Lemma 4, the solution of model (1) is still in Π . From model (1), we get 0 C D t α N t = Λ μ 0 N μ 1 I Λ μ 0 N . Then:
0 N t N ( 0 ) E α μ 0 t α + Λ t α E α , α + 1 μ 0 t α , t 0 .
Thus, limsup t N ( t ) Λ μ 0 . Therefore, all feasible solutions of model (1) are in Π . That is, N ( t ) is bounded by Λ μ 0 . Thus, all solutions of model (1) are in Π for all t 0 . □

3.2. Equilibrium Points and Their Asymptotic Stability

3.2.1. Disease-Free Equilibrium Point E 0 and Stability

Let 0 C D t α S t = 0 , 0 C D t α I t = 0 , 0 C D t α R t = 0 , then:
E 0 = ( S 0 , I 0 , R 0 ) = Λ μ 0 + ν , 0 , ν Λ μ 0 + ν μ 0 .
Let X = ( S , I , R ) , then model (1) can be written in the form: d X d t = F ( X ) V ( X ) , where
F ( X ) = 0 β S I f S , I 0 , V ( X ) = Λ + β S I f S , I + μ 0 + ν S μ 0 + μ 1 + γ I γ I ν S + μ 0 R .
The Jacobian of Equation (4) around E 0 . Then, the basic reproduction number R 0 is the spectral radius of F V 1 , in the following form:
R 0 = β Λ μ 0 + ν + a 1 Λ μ 0 + μ 1 + γ .
Theorem 2.
If R 0 < 1 , then E 0 is locally asymptotically stable, and conversely if R 0 > 1 , then E 0 is unstable.
Proof. 
The Jacobian matrix of hepatitis B epidemic model (1) at E 0 has the following form:
J | E 0 = μ 0 + ν β Λ μ 0 + ν + a 1 Λ 0 0 β Λ μ 0 + ν + a 1 Λ μ 0 + μ 1 + γ 0 ν γ μ 0 .
After calculation, J | E 0 has three eigenvalues:
λ 1 = μ 0 + ν < 0 , λ 2 = μ 0 < 0 , λ 3 = μ 0 + μ 1 + γ 1 R 0 .
If R 0 < 1 , then λ 3 < 0 . Therefore, E 0 is locally asymptotically stable. On the contrary, if R 0 > 1 , then λ 3 is positive, so E 0 is an unstable. □

3.2.2. Endemic Equilibrium Point E * and Stability

The calculation leads to E * = S * , I * , R * in the following form:
S * = Λ t 1 I * t 2 , I * = Δ 1 + Δ 1 2 4 t 1 2 a 3 Δ 2 2 t 1 2 a 3 , R * = γ I * + ν S * μ 0 ,
where t 1 = μ 0 + μ 1 + γ , t 2 = μ 0 + ν , Δ 1 = t 1 2 a 1 Λ t 1 a 3 t 1 t 2 a 2 β t 1 , Δ 2 = β Λ t 1 t 2 Λ t 1 a 1 .
Lemma 5.
If R 0 > 1 , then E * exists.
Proof. 
If R 0 > 1 , we have:
Δ 1 = t 1 2 a 1 Λ t 1 a 3 t 1 t 2 a 2 β t 1 < μ 0 + μ 1 + γ μ 0 + ν μ 0 + μ 1 + γ Λ + a 2 Λ a 3 < 0 , Δ 2 = β Λ t 1 t 2 Λ t 1 a 1 = μ 0 + ν + a 1 Λ μ 0 + μ 1 + γ R 0 d 1 > 0 .
Furthermore, by calculation, we obtain Δ 1 2 4 t 1 2 a 3 Δ 2 > 0 . Thus, I * > 0 . Because 2 Λ t 1 a 3 + Δ 1 > Δ 1 2 4 t 1 2 a 3 Δ 2 , so Λ Δ 1 + Δ 1 2 4 t 1 2 a 3 Δ 2 2 t 1 a 3 > 0 . Then, Λ t 1 I * > 0 . Hence, we have S * > 0 . Thus, R * > 0 . Consequently, E * exists, if R 0 > 1 . □
Theorem 3.
If R 0 > 1 , then E * is locally asymptotically stable, and conversely, if R 0 < 1 , then E * is unstable.
Proof. 
The Jacobian matrix of hepatitis B epidemic model (1) at E * has the following form:
J | E * = μ 0 + ν β I * 1 + a 2 I * f 2 S * , I * β S * 1 + a 1 S * f 2 S * , I * 0 β I * 1 + a 2 I * f 2 S * , I * β S * 1 + a 1 S * f 2 S * , I * μ 0 + μ 1 + γ 0 ν γ μ 0 .
Obviously, λ 1 = μ 0 < 0 . The other two eigenvalues are calculated below. Let:
A = μ 0 + ν β I * 1 + a 2 I * f 2 S * , I * β S * 1 + a 1 S * f 2 S * , I * β I * 1 + a 2 I * f 2 S * , I * β S * 1 + a 1 S * f 2 S * , I * μ 0 + μ 1 + γ .
By d I d t = 0 , one has β S * f S * , I * = μ 0 + μ 1 + γ . Thus, β S * = μ 0 + μ 1 + γ f S * , I * . By the definition of f, we get β S * 1 + a 1 S * < μ 0 + μ 1 + γ f 2 S * , I * . Thus:
β S * 1 + a 1 S * f 2 S * , I * < μ 0 + μ 1 + γ .
Then, by (5), we obtain:
trac A = μ 0 + ν β I * 1 + a 2 I * f 2 S * , I * + β S * 1 + a 1 S * f 2 S * , I * μ 0 + μ 1 + γ < 0
and
det A = μ 0 + ν + β I * 1 + a 2 I * f 2 S * , I * μ 0 + μ 1 + γ β S * 1 + a 1 S * f 2 S * , I * + β S * 1 + a 1 S * f 2 S * , I * · β I * 1 + a 2 I * f 2 S * , I * > 0 .
Thus, trac A < 0 , det A > 0 when R 0 > 1 . According to the Routh–Hurwitz criterion, the two eigenvalues of the matrix A are negative, so E * is locally asymptotically stable when R 0 > 1 . □

3.3. Numerical Simulation

Example 1.
In model (1), let α = 0.98 , Λ = 0.8 , β = 0.01 , a 1 = 0.1 , a 2 = 0.2 , a 3 = 1.2 , μ 0 = 0.02 , μ 1 = 0.1 , ν = 0.01 , γ = 0.01 , S ( 0 ) = 200 , I 0 = 100 , R ( 0 ) = 100 . After the calculation, we obtain E 0 = ( 27 , 0 , 13 ) and R 0 = 0.56 < 1 . By observing Figure 1a, we can see that E 0 is locally asymptotically stable, which verifies the rationality of Theorem 2.
Example 2.
See Example 1 for the values of all parameters except β = 0.9 . After calculation, we obtain E * = ( 5 , 5 , 10 ) and R 0 = 50 > 1 . By observing Figure 1b, we can see that E * is locally asymptotically stable, which verifies the rationality of Theorem 3.

4. Optimal Control Analysis

Our aim is to find an effective control strategy to reduce hepatitis B transmission. Optimal control theory [21,22,23] is adopted in this section. The number of hepatitis B patients is reduced by three control variables u 1 ( t ) , u 2 ( t ) , and u 3 ( t ) . The specific explanation of the control variables is as follows:
(i).
u 1 ( t ) represents the isolation rate. Through this control variable, infected persons are isolated to avoid contact between infected persons and susceptible persons;
(ii).
u 2 ( t ) represents the cure rate. Through this control variable, the number of infected individuals can be reduced by using effective drugs to treat the infected;
(iii).
u 3 ( t ) represents vaccination rate. The spread of Hepatitis B can be reduced through vaccination.
To design a control strategy to eliminate hepatitis B, we will consider the optimal strategy of model (1). The control strategy is to minimize the following objective function:
J u 1 , u 2 , u 3 = 0 I T α ω 1 I t + 1 2 ω 2 u 1 2 t + ω 3 u 2 2 t + ω 4 u 3 2 t ,
subject to the control system:
0 C D t α S t = Λ β S I f S , I 1 u 1 t μ 0 + ν + u 3 t S , 0 C D t α I t = β S I f S , I 1 u 1 t μ 0 + μ 1 + γ + u 2 t + u 3 t I , 0 C D t α R t = γ + u 2 t + u 3 t I + ν + u 3 t S μ 0 R ,
with
S 0 0 , I 0 0 , R 0 0 .
In formula (6), ω 1 represents the weight constant of hepatitis B infection I ( t ) , ω 2 , ω 3 , and ω 4 denote the weight constants for isolation of infected and susceptible individuals, treatment of infected individuals and vaccination, respectively. 1 2 ω 2 u 1 2 t , 1 2 ω 3 u 2 2 t , and 1 2 ω 4 u 3 2 t represent the costs of isolation, treatment, and vaccination, respectively. Our goal is to search the control functions u 1 * , u 2 * , u 3 * so that:
J u 1 * , u 2 * , u 3 * = min J u 1 , u 2 , u 3 , u 1 , u 2 , u 3 U
subject to problem (7)–(8), where
U : = u 1 , u 2 , u 3 0 u i 1 , u i t is Lebesgue measurable on [ 0 , T ] , i = 1 , 3 ¯ .
By applying Pontryaginâs maximum principle in [24] and the results in [25], the following theorem can be obtained.
Theorem 4.
Let S * , I * , and R * be optimal state solution with associated optimal control variables u 1 * , u 2 * , u 3 * for the optimal control problems (6)–(7). Then, there are three adjoint variables λ 1 , λ 2 , and λ 3 such that:
0 C D t α λ 1 T t = β I * T t 1 + a 2 I * T t f 2 S * T t , I * T t × λ 2 T t λ 1 T t 1 u 1 * T t λ 1 T t λ 3 T t u 3 * T t λ 1 T t μ 0 + ν + λ 3 T t ν , 0 C D t α λ 2 T t = ω 1 + β S * T t 1 + a 1 S * T t f 2 S * T t , I * T t × λ 2 T t λ 1 T t 1 u 1 * T t λ 2 T t λ 3 T t u 2 * T t + u 3 * T t λ 2 T t μ 0 + μ 1 + γ + λ 3 T t γ , 0 C D t α λ 3 T t = λ 3 T t μ 0
with terminal conditions:
λ i ( T ) = 0 , i = 1 , 2 , 3 .
In addition, the optimal controls variables u 1 * , u 2 * , and u 3 * have the following specific forms:
u 1 * = max min λ 2 λ 1 β S * I * f S * , I * ω 2 , 1 , 0 ,
u 2 * = max min λ 2 λ 3 I * ω 3 , 1 , 0 ,
u 3 * = max min λ 1 λ 3 S * + λ 2 λ 3 I * ω 4 , 1 , 0 .
Proof. 
Let x = S , I , R , u = ( u 1 , u 2 , u 3 ) . Then, the Lagrangian L is defined as:
L x , u = ω 1 I + 1 2 ω 2 u 1 2 + ω 3 u 2 2 + ω 4 u 3 2 ,
and the related Hamiltonian function H is defined as:
H x , u , λ = L x , u + λ g x , u ,
where λ = λ 1 , λ 2 , λ 3 , and g x , u = g 1 x , u , g 2 x , u , g 3 x , u with:
g 1 x , u = Λ β S I f S , I 1 u 1 μ 0 + ν + u 3 S , g 2 x , u = β S I f S , I 1 u 1 μ 0 + μ 1 + γ + u 2 + u 3 I , g 3 x , u = γ + u 2 + u 3 I + ν + u 3 S μ 0 R .
Then, we calculate the following adjoint equations:
t D T α λ 1 = H S * = λ 1 β I * 1 + a 2 I * f 2 S * , I * 1 u 1 * μ 0 + ν + u 3 * + λ 2 β I * 1 + a 2 I * f 2 S * , I * 1 u 1 * + λ 3 ν + u 3 * , t D T α λ 2 = H I * = ω 1 + λ 1 β S * 1 + a 1 S * f 2 S * , I * 1 u 1 * + λ 2 β S * 1 + a 1 S * f 2 S * , I * 1 u 1 * μ 0 + μ 1 + γ + u 2 * + u 3 * + λ 3 γ + u 2 * + u 3 * , t D T α λ 3 = H R * = λ 3 μ 0 ,
with terminal conditions:
λ i ( T ) = 0 , i = 1 , 3 ¯
Using Lemma 1, the equivalent equation can be obtained in the following form:
t C D T α λ 1 t = β I * t 1 + a 2 I * t f 2 S * , I * λ 2 t λ 1 t 1 u 1 * t λ 1 t μ 0 + ν + u 3 * t + λ 3 t ν + u 3 * t , t C D T α λ 2 t = ω 1 + β S * t 1 + a 1 S * t f 2 S * , I * λ 2 t λ 1 t 1 u 1 * t λ 2 t μ 0 + μ 1 + γ + u 2 * t + u 3 * t + λ 3 t γ + u 2 * t + u 3 * t , t C D T α λ 3 t = λ 3 t μ 0 ,
where t C D T α is the fractional-order derivative of the right-hand side Caputo type. By Lemma 2, one has:
0 C D t α λ 1 T t = β I * T t 1 + a 2 I * T t f 2 S * T t , I * T t × λ 2 T t λ 1 T t 1 u 1 * T t λ 1 T t λ 3 T t u 3 * T t λ 1 T t μ 0 + ν + λ 3 T t ν , 0 C D t α λ 2 T t = ω 1 + β S * T t 1 + a 1 S * T t f 2 S * T t , I * T t × λ 2 T t λ 1 T t 1 u 1 * T t λ 2 T t λ 3 T t u 2 * T t + u 3 * T t λ 2 T t μ 0 + μ 1 + γ + λ 3 T t γ , 0 C D t α λ 3 T t = λ 3 T t μ 0 ,
where 0 C D t α is the fractional-order derivative of the left-hand side Caputo type. Let H u k = 0 , k = 1 , 3 ¯ , then we have:
H u 1 * = ω 2 u 1 * t + λ 1 t λ 2 t β S * I * f S * , I * = 0 , H u 2 * = ω 3 u 2 * t + λ 3 t λ 2 t I * = 0 , H u 3 * = ω 4 u 3 * t + λ 3 t λ 1 t S * + λ 3 t λ 2 t I * = 0 .
This achieves the optimality condition H x * , u * , λ = min u 0 , 1 H x * , u * , λ . By taking into account:
u ^ 1 * t = β S * I * λ 2 t λ 1 t f S * , I * ω 2 , u ^ 2 * t = λ 2 t λ 3 t I * ω 3 ,
u ^ 3 * t = λ 1 t λ 3 t S * + λ 2 t λ 3 t I * ω 4 ,
the optimal value of the objective function can be obtained:
u 1 * = max min u ^ 1 * t , 1 , 0 , u 2 * = max min u ^ 2 * t , 1 , 0 , u 3 * = max min u ^ 3 * t , 1 , 0 .

5. Numerical Simulation

To test the plausibility of the results, we simulate the hepatitis B model (1). Firstly, a set of discrete knots t j = t 0 + j h , h = T t 0 N ; j = 0 , 1 , , N are used to replace the time interval [ t 0 , T ] . By the fractional Euler method for model (7), we obtain the following discretized equations:
S t k = S t 0 + h α Γ α + 1 j = 0 k 1 c j , k Λ β S t j I t j f S t j , I t j 1 u 1 μ 0 + ν + u 3 S t j , I t k = I t 0 + h α Γ α + 1 j = 0 k 1 c j , k β S t j I t j f S t j , I t j 1 u 1 μ 0 + μ 1 + γ + u 2 + u 3 I t j , R t k = R t 0 + h α Γ α + 1 j = 0 k 1 c j , k γ + u 2 + u 3 I t j + ν + u 3 S t j μ 0 R t j ,
where c j , k = k j α k 1 j α , for k = 1 , N ¯ , u 1 ( t ) represents the isolation rate, u 2 ( t ) represents the cure rate, and u 3 ( t ) represents vaccination rate. We study the effect of α in the absence of isolation measures, vaccination, and treatment (i.e., u 1 = u 2 = u 3 = 0 ) on the kinetic behavior of the hepatitis B model (1) with a step size of = 0.1 . We plot all compartments for α = 0.9 ; 0.95 ; 0.98 ; 1 using the data from Example 1 (see Figure 2).
Figure 2 demonstrates that the curve more closely resembles the integer-order hepatitis B model. As the fractional order α increases, the number of people recovering increases faster and the disease becomes extinct faster.
We also solve numerically the fractional-order optimal control problems (6), (7), and (11) and compare them with the uncontrolled hepatitis B model (1) (see Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9).
Figure 3 describes the optimal control results of the fractional-order hepatitis B model under the measures of isolating infected persons and avoiding contact between infected and susceptible persons by controlling variable u 1 ( t ) . Figure 3a–c represent the dynamic behavior of susceptible, infected, and recovered persons under isolation measures, respectively. By observing Figure 3, it can be seen that the greater the isolation intensity (i.e., the greater the value of u 1 ( t ) ), the faster the decline in the number of infected patients.
Figure 4 depicts the optimal control results of the fraction-order hepatitis B model under the measure of effective drug treatment of infected persons through the control variable u 2 ( t ) . Figure 4a–c represent the dynamic behavior of susceptible, infected, and recovered persons under effective drug treatment, respectively. By observing Figure 4, it can be seen that the greater the intensity of effective drug treatment for infected patients (i.e., the greater the value of u 2 ( t ) ), the faster the number of infected patients will decline and the faster the recovery rate of recovered patients will be.
Figure 5 depicts the optimal control results of the fractional-order hepatitis B model using measures to vaccinate susceptible individuals by controlling variable u 3 ( t ) . Figure 5a–c represent the dynamic behavior of susceptible, infected, and recovered persons under vaccination measures, respectively. From the observation of Figure 5, it can be seen that the greater the intensity of vaccination for susceptible persons (i.e., the greater the value of u 3 ( t ) ), the faster the number of susceptible persons will decline, and the faster the number of infected persons will decline.
Figure 6 describes the optimal control results for the fractional hepatitis B model by controlling variables u 1 ( t ) , u 2 ( t ) . Figure 6a–c represent, respectively, the dynamic behavior of susceptible, infected, and recovered persons under simultaneous isolation and effective drug treatment of infected persons. It can be seen from Figure 6 that the greater the intensity of isolation and effective drug treatment for infected patients (i.e., the greater the value of u 1 ( t ) and u 2 ( t ) ), the faster the number of infected patients will decline, and the faster the recovery rate of recovered patients will be.
Figure 7 describes the optimal control results for the fractional hepatitis B model by controlling variables u 1 ( t ) , u 3 ( t ) . Figure 7a–c represent the dynamic behavior of susceptible, infected, and recovered persons, respectively, through simultaneous isolation of infected persons and vaccination of vulnerable populations. It can be seen from Figure 7 that the greater the intensity of isolation of infected persons and vaccination of susceptible persons (i.e., the greater the value of u 1 ( t ) and u 3 ( t ) ), the faster the number of susceptible persons will decline, and the faster the number of infected persons will decline.
Figure 8 describes the optimal control results for the fractional hepatitis B model by controlling variables u 2 ( t ) , u 3 ( t ) . Figure 8a–c represent the dynamic behavior of susceptible, infected, and recovered persons while effective drug treatment of infected persons and vaccination of susceptible persons, respectively. By observing Figure 8, it can be seen that the greater the intensity of effective drug treatment for infected persons and vaccination for susceptible persons (i.e., the greater the value of u 2 ( t ) and u 3 ( t ) ), the faster the number of susceptible persons will decline, and the faster the number of infected persons will decline.
Figure 9 depicts controlling the variables u 1 ( t ) , u 2 ( t ) , u 3 ( t ) for fractional-order optimal control results of the hepatitis B model. Figure 9a–c represent, respectively, the dynamic behavior of susceptible, infected, and recovered persons through simultaneous isolation of infected and susceptible persons, effective drug treatment of infected persons, and vaccination of susceptible persons. It can be seen from Figure 9 that the greater the intensity of isolation and effective drug treatment of infected persons and vaccination of susceptible persons (i.e., the greater the value of u 1 ( t ) , u 2 ( t ) , u 3 ( t ) ), the faster the number of susceptible persons will decline, and the faster the number of infected persons will decline. As can be seen from the observations in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9, the number of infected persons declines most rapidly through the simultaneous effect of isolating infected and susceptible persons, effective drug treatment of infected persons, and vaccination of susceptible persons. That is, the higher the value of u 1 ( t ) , u 2 ( t ) , u 3 ( t ) , the faster the extinction rate of hepatitis B virus.

6. Conclusions

In this paper, the transmission dynamics of the hepatitis B virus are studied, and the optimal control strategy is developed to control the transmission of the hepatitis B virus. Specific results are as follows:
(i)
A fractional-order hepatitis B transmission dynamics model with general incidence is proposed.
(ii)
The positive and boundedness of the model solutions are studied, and the basic regeneration number, equilibrium points, and stability of the model are obtained.
(iii)
In the fractional hepatitis B model, three control variables u 1 , u 2 , and u 3 are added, representing isolation, treatment, and vaccination, respectively. Based on the Pontryagin’s extreme value principle, the necessary optimality conditions for fractional-order optimal control problems are derived.
(iv)
Numerical simulation of fractional-order optimal control system is given. These numerical simulations show that by isolating infected and susceptible persons, ensuring effective drug treatment of infected persons, and ensuring vaccination of susceptible persons, hepatitis B virus transmission can be controlled and prevented.
In the future, we will study the impact of media coverage on curbing the spread of hepatitis B. At the same time, we will further investigate the mechanism by which environmental fluctuations affect the transmission dynamics of hepatitis B when considering intervention strategies.

Author Contributions

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

Funding

This work is supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region (Grant No. 2022D01A246, 2022D01A247, 2021D01B35, 2021D01A65, 2019D01B10), Natural Science Foundation of colleges and universities in Xinjiang Uygur Autonomous Region (Grant No. XJEDU2021Y048), and Foundation of Xinjiang Institute of Engineering (Grant No. 2015xgy161712).

Data Availability Statement

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

Acknowledgments

The authors sincerely thank the editors and anonymous reviewers for the careful reading of the original manuscript and valuable comments, which have improved the quality of our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, J.W.; Zhou, Y.Z.; Wang, Z.G.; Wang, H.H. Analysis and achievement for fractional optimal control of Hepatitis B with Caputo operator. Alex. Eng. J. 2023, 70, 601–611. [Google Scholar] [CrossRef]
  2. Xue, T.T.; Zhang, L.; Fan, X.L. Dynamic modeling and analysis of Hepatitis B epidemic with general incidence. Math. Biosci. Eng. 2023, 20, 10883–10908. [Google Scholar] [CrossRef]
  3. Yavuz, M.; Ozkose, F.; Susam, M.; Kalidass, M. A New Modeling of Fractional-Order and Sensitivity Analysis for Hepatitis-B Disease with Real Data. Fractal Fract. 2023, 7, 165. [Google Scholar] [CrossRef]
  4. Khan, T.; Khan, A.; Zaman, G. The extinction and persistence of the stochastic hepatitis B epidemic model. Chaos Solitons Fract. 2018, 108, 123–128. [Google Scholar] [CrossRef]
  5. Din, A.; Li, Y.J.; Khan, T.; Anwar, K.; Zaman, G. Stochastic dynamics of hepatitis B epidemics. Results Phys. 2021, 20, 103730. [Google Scholar] [CrossRef]
  6. Liu, P.J.; Din, A.; Huang, L.F.; Yusuf, A. Stochastic optimal control analysis for the hepatitis B epidemic model. Results Phys. 2021, 26, 104372. [Google Scholar] [CrossRef]
  7. Din, A.; Li, Y.J. Stationary distribution extinction and optimal control for the stochastic hepatitis B epidemic model with partial immunity. Phys. Scr. 2021, 96, 074005. [Google Scholar] [CrossRef]
  8. Din, A.; Li, Y.J.; Yusuf, A. Delayed hepatitis B epidemic model with stochastic analysis. Chaos Solitons Fract. 2021, 146, 110839. [Google Scholar] [CrossRef]
  9. Khan, T.; Ullah, Z.; Ali, Z.; Zaman, G. Modeling and control of the hepatitis B virus spreading using an epidemic model. Chaos Solitons Fract. 2019, 124, 1–9. [Google Scholar] [CrossRef]
  10. Haq, F.; Shah, K.; Rahman, G.U.; Shahzad, M. Numerical solution of fractional order smoking model via Laplace Adomian decomposition method. Alex. Eng. J. 2018, 57, 1061–1069. [Google Scholar] [CrossRef]
  11. Shah, K.; Khalil, H.; Khan, R.A. Investigation of positive solution to a coupled system of impulsive boundary value problems for nonlinear fractional order differential equations. Chaos Solitons Fract. 2015, 77, 240–246. [Google Scholar] [CrossRef]
  12. Ahmad, Z.; Bonanomi, G.; di Serafino, D.; Giannino, F. Transmission dynamics and sensitivity analysis of pine wilt disease with asymptomatic carriers via fractal-fractional differential operator of Mittag–Leffler kernel. Appl. Numer. Math. 2023, 185, 446–465. [Google Scholar] [CrossRef]
  13. Sinan, M.; Shah, K.; Kumam, P.; Mahariq, I.; Ansari, K.J.; Ahmad, Z.; Shah, Z. Fractional order mathematical modeling of typhoid fever disease. Results Phys. 2022, 32, 105044. [Google Scholar] [CrossRef]
  14. Malik, A.; Alkholief, M.; Aldakheel, F.M.; Khan, A.A.; Ahmad, Z.; Kamal, W.; Gatasheh, M.K.; Alshamsan, A. Sensitivity analysis of COVID-19 with quarantine and vaccination: A fractal-fractional model. Alex. Eng. J. 2022, 61, 8859–8874. [Google Scholar] [CrossRef]
  15. Din, A.; Li, Y.J.; Yusuf, A.; Ali, A.I. Caputo type fractional operator applied to Hepatitis B system. Fractals 2022, 30, 2240023. [Google Scholar] [CrossRef]
  16. Ucar, S. Analysis of hepatitis B disease with fractal-fractional Caputo derivative using real data from Turkey. J. Comput. Appl. Math. 2023, 419, 114692. [Google Scholar] [CrossRef]
  17. Simelane, S.M.; Dlamini, P.G. A fractional order differential equation model for Hepatitis B virus with saturated incidence. Results Phys. 2021, 24, 104114. [Google Scholar] [CrossRef]
  18. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  19. Ameen, I.; Baleanu, D.; Ali, H.M. An efficient algorithm for solving the fractional optimal control of SIRV epidemic model with a combination of vaccination and treatment. Chaos Solitons Fract. 2020, 137, 109892. [Google Scholar] [CrossRef]
  20. Odibat, Z.M.; Shawagfeh, N.T. Generalized Taylor’s formula. Appl. Math. Comput. 2007, 186, 286–293. [Google Scholar] [CrossRef]
  21. Zaman, G.; Kang, Y.H.; Jung, I.H. Optimal treatment of an SIR epidemic model with time delay. Biosystems 2009, 98, 43–50. [Google Scholar] [CrossRef] [PubMed]
  22. Xia, M.T.; Bottcher, L.; Chou, T. Controlling epidemics through optimal allocation of test kits and vaccine doses across networks. IEEE Trans. Netw. Sci. Eng. 2022, 9, 1422–1436. [Google Scholar] [CrossRef]
  23. Zaman, G.; Kang, Y.H.; Jung, I.H. Stability analysis and optimal vaccination of an SIR epidemic model. BioSystems 2008, 93, 240–249. [Google Scholar] [CrossRef] [PubMed]
  24. Pontryagin, S.; Boltyanskii, V.; Gamkrelidze, R.; Mishchenko, E. The Mathematical Theory of Optimal Processes; Gordon and Breach Science Publishers: London, UK, 1986. [Google Scholar]
  25. Ali, H.M.; Pereira, F.L.; Gama, S.M.A. A new approach to the pontryagin maximum principle for nonlinear fractional optimal control problems. Math. Meth. Appl. Sci. 2016, 39, 3640–3649. [Google Scholar] [CrossRef]
Figure 1. The plot results show the equilibrium points of the model (1).
Figure 1. The plot results show the equilibrium points of the model (1).
Mathematics 11 03642 g001
Figure 2. Transient behavior of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) at different fractional-order derivatives α = 0.9 ; 0.95 ; 0.98 ; 1 .
Figure 2. Transient behavior of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) at different fractional-order derivatives α = 0.9 ; 0.95 ; 0.98 ; 1 .
Mathematics 11 03642 g002
Figure 3. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control u 1 .
Figure 3. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control u 1 .
Mathematics 11 03642 g003
Figure 4. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control u 2 .
Figure 4. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control u 2 .
Mathematics 11 03642 g004
Figure 5. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control u 3 .
Figure 5. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control u 3 .
Mathematics 11 03642 g005
Figure 6. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 1 , u 2 .
Figure 6. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 1 , u 2 .
Mathematics 11 03642 g006
Figure 7. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 1 , u 3 .
Figure 7. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 1 , u 3 .
Mathematics 11 03642 g007
Figure 8. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 2 , u 3 .
Figure 8. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 2 , u 3 .
Mathematics 11 03642 g008
Figure 9. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 1 , u 2 , u 3 .
Figure 9. Simulations of ( S ( t ) , I ( t ) , R ( t ) ) in model (7) with different control pairs u 1 , u 2 , u 3 .
Mathematics 11 03642 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xue, T.; Fan, X.; Xu, Y. Kinetic Behavior and Optimal Control of a Fractional-Order Hepatitis B Model. Mathematics 2023, 11, 3642. https://doi.org/10.3390/math11173642

AMA Style

Xue T, Fan X, Xu Y. Kinetic Behavior and Optimal Control of a Fractional-Order Hepatitis B Model. Mathematics. 2023; 11(17):3642. https://doi.org/10.3390/math11173642

Chicago/Turabian Style

Xue, Tingting, Xiaolin Fan, and Yan Xu. 2023. "Kinetic Behavior and Optimal Control of a Fractional-Order Hepatitis B Model" Mathematics 11, no. 17: 3642. https://doi.org/10.3390/math11173642

APA Style

Xue, T., Fan, X., & Xu, Y. (2023). Kinetic Behavior and Optimal Control of a Fractional-Order Hepatitis B Model. Mathematics, 11(17), 3642. https://doi.org/10.3390/math11173642

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