Next Article in Journal
A Multi-Branch Feature Extraction Residual Network for Lightweight Image Super-Resolution
Previous Article in Journal
Unified Modeling and Multi-Objective Optimization for Disassembly Line Balancing with Distinct Station Configurations
Previous Article in Special Issue
Advancing Fractional Riesz Derivatives through Dunkl Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Dynamics and Optimal Control of a Fractional-Order SIV Epidemic Model with Nonmonotonic Occurrence Rate

School of Mathematics and Computer Science, Yunnan Minzu University, Yuehua Street No. 2929, Chenggong District, Kunming 650500, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2024, 12(17), 2735; https://doi.org/10.3390/math12172735
Submission received: 3 August 2024 / Revised: 27 August 2024 / Accepted: 29 August 2024 / Published: 1 September 2024
(This article belongs to the Special Issue Fractional-Order Systems: Control, Modeling and Applications)

Abstract

:
This paper performs a detailed analysis and explores optimal control strategies for a fractional-order SIV epidemic model, incorporating a nonmonotonic incidence rate. In this paper, the population of vaccinated individuals is included in the disease dynamics model. After proving the non-negative boundedness of the fractional-order SIV model, we focus on analyzing the equilibrium point characteristics of the model, delving into its existence, uniqueness, and stability analysis. In addition, our research includes formulating optimal control strategies specifically aimed at minimizing the number of infections while keeping costs as low as possible. To validate the theoretical findings and uncover the practical efficacy and prospects of control measures in mitigating epidemic spread, numerical simulations are performed.

1. Introduction

Throughout the early 21st century, the globe has been beset by recurrent outbreaks of novel infectious diseases, which constitute a formidable hazard to human well-being and the very essence of life’s safety [1]. Instances include the severe acute respiratory syndrome (SARS) epidemic in 2003 [2,3] and the recent novel coronavirus pneumonia (COVID-19) outbreak towards the conclusion of 2019 [4,5,6], both of which underscore the ongoing global challenge posed by infectious diseases. The outbreak and spread of these epidemics have brought tremendous challenges to the global public health system and prompted researchers to pay more attention to the study of infectious disease dynamics. Among them, the epidemic research method based on mathematical models has received widespread attention because it is capable of quantitatively portraying the transmission dynamics of diseases, enabling forecasts regarding the evolution and progression of epidemic outbreaks [7].
Vaccination is crucial for eliminating infectious diseases, such as hepatitis B, measles, smallpox, poliomyelitis, and so on [8,9]. Since the pioneering work on smallpox of Edward Jenner (a doctor who worked in the UK and noticed that people infected with cowpox rarely contracted smallpox) [10], the practice of safeguarding individuals against infections via vaccination has evolved into a standard and routine procedure. However, achieving full vaccination coverage among all susceptible individuals within a given community remains a formidable challenge, particularly in nations where accessibility to such vaccines is hindered [11]. Certain clinical findings suggest that vaccines may confer only a temporary degree of protection against certain diseases. Hence, when a vaccine’s effectiveness diminishes within a vaccinated individual’s system, they once again become vulnerable to contracting the disease [12]. Consequently, identifying the optimal level of vaccination needed to eliminate a disease becomes paramount.
The incidence rate of a disease, which signifies the quantity of fresh cases arising within a specified time frame, holds significant weight in the realm of mathematical epidemiology research. Conventional epidemic models often utilize standard and bilinear incidence rates, like β S I (with β denoting the transmission rate) [7] and β S I N (where N represents the total population) [5], to depict the occurrence of new infections within a population. However, it is increasingly acknowledged that the utilization of linear incidence rates may fall short in precisely capturing the intricate progression patterns of infectious diseases within a population. In order to overcome this limitation, scientists and investigators have explored nonlinear incidence rates as a potential solution. For example, in their research of cholera outbreaks in Paris, Capasso and Serio [13] first proposed a nonlinear incidence rate with saturation effects. They further introduced the concept of saturated incidence rate into infectious disease models, specifically in the form of g ( I ) S , where g is a non-negative function with g ( 0 ) = 0 that ultimately tends to a saturation level. As research on infectious disease models has progressed, the form of saturated incidence rates has also become more complex. W. Liu et al. [14] proposed a general incidence rate in the form of β S I a 1 + λ I b , where a and b are positive constants and λ is a positive parameter, used to consider the influence of changes in the behavior of susceptible individuals. Furthermore, λ determines the magnitude of psychological or inhibitory effects, while 1 + λ I b represents how behavioral changes slow disease spread by reducing effective transmission as infected cases I increase, or what is known as the crowding effect of I.
Researchers have constructed and analyzed various mathematical models to investigate the influence of vaccines on the transmission of infectious diseases. Ref. [15] explored the SIS model in the context of partially effective vaccination, assuming that the vaccination strategy is primarily targeted at susceptible individuals while the total population size remains constant and that the spread of disease adheres to a standard incidence rate. In contrast, ref. [16] considered the same model in an environment where the overall population size is subject to change and assumed that the vaccine is fully effective for all susceptible individuals as well as newly arrived members of the population. Meanwhile, refs. [17,18], respectively, studied similar models with differing assumptions regarding incidence rates and population sizes. Ref. [19] presented a discrete-time SIS model that includes a vaccination scheme, with disease transmission occurring at a saturated incidence rate and the assumption that the population size fluctuates and that infected individuals acquire temporary immunity upon recovery. A deterministic and stochastic SIS system with continuous vaccination strategies is proposed in [20]. At the same time, it is assumed that all vaccines are effective and the disease spreads with bilinear incidence.
In models formulated through ordinary differential equations, the velocity of variation exclusively relies on the instantaneous state, implying a localized influence [21,22]. In reality, however, changes in individual behavior also depend on historical states, i.e., states prior to the current state. One method to simulate this phenomenon is through fractional-order derivatives. Fractional calculus generalizes the integer-order derivative and integral to arbitrary orders. Compared to integer-order calculus, the nonlocal characteristics of fractional calculus are more capable of accurately describing dynamic systems with features such as time dependency, path dependency, memory, and heredity [23,24]. This methodology offers a refined comprehension of the intricate interplay among diverse factors in disease transmission, proving instrumental for researchers across numerous physical contexts and facilitating successful implementations of Caputo derivatives [25,26].
Given the continuous advancements in the field of epidemiological models, this study focuses on analyzing the equilibrium point characteristics of the fractional-order SIV model after proving its non-negative boundedness and delves into its existence, uniqueness, and stability analysis. Additionally, we comprehensively consider core epidemiological factors such as fractional derivatives, neonatal immunization coverage, vaccination rates among susceptible populations, vaccine efficacy, and public behavior effects to explore how vaccination strategies impact the dynamics of disease transmission. By meticulously examining the interactions between these parameters, we aim to uncover the actual effectiveness and potential of control strategies in containing the spread of epidemics.
Furthermore, our research involves formulating optimal control strategies specifically aimed at minimizing the number of infections while keeping costs as low as possible through control measures such as modifying the behavior of susceptible populations, optimizing vaccination rates, and enhancing treatment for infected individuals. This strategy formulation process considers the practical constraints and desired goals of disease control, striving to achieve the most optimal allocation and efficient utilization of vaccination resources under limited resources.
The organization of the rest of this paper is as follows: Section 2 introduces pertinent definitions and lemmas essential for subsequent analyses. Section 4 delves into the dynamic behaviors of model (10). Section 5 examines the stability of model (10) at its equilibrium point. Section 6 formulates the fractional-order optimal control problem and outlines the optimal control strategy’s architecture. Section 7 employs numerical simulations to corroborate the theoretical findings. Lastly, Section 8 summarizes the conclusions and outlines directions for future research.

2. Mathematical Preliminaries

Definition 1
([27]). Given a function f ( t ) defined on the interval ( a , b ) , with μ > 0 and n being the smallest integer greater than μ, the Caputo fractional derivative of order μ for f ( t ) is formulated as
D t μ a C f ( t ) = 1 Γ ( n μ ) a t f ( n ) ( ζ ) ( t ζ ) 1 + μ n d ζ , n 1 < μ < n d n d t n f ( t ) , μ = n .
In particular, when 0 < μ < 1 , we can obtain
D t μ a C f ( t ) = 1 Γ ( 1 α ) a t f ( ζ ) ( t ζ ) μ d ζ .
Definition 2
([28]). The Mittag–Leffler function with two parameters is represented as
E α , θ ( u ) = i = 0 u i Γ ( i α + θ ) ,
where α > 0 , θ > 0 , u C .
Definition 3
([27,29]). The Laplace transform, when applied to the Caputo fractional differentiation, can be expressed as follows:
L D t μ a C f ( t ) ; s = s μ F ( s ) i = 0 n 1 s μ i 1 f ( i ) ( 0 ) ,
where n 1 < μ n . In particular, when the initial value satisfies f i ( 0 ) = 0 , i = 0 , 1 , , n 1 , it can be obtained that
L D t μ a C f ( t ) ; s = s μ F ( s ) .
Lemma 1
([30]). Let h C 1 ( a , b ] and μ ( 0 , 1 ) :
(i) For ∀ t [ a , b ] , if D μ h ( t ) > 0 , then h ( t ) is monotonically increasing.
(ii) For ∀ t [ a , b ] , if D μ h ( t ) < 0 , then h ( t ) is monotonically decreasing.
Lemma 2
([31]). Suppose ψ : [ t 0 , ) R is a continuous function, and satisfy
D μ ψ ( t ) + λ ψ ( t ) ν , ψ ( t 0 ) = ψ 0 .
where μ ( 0 , 1 ] , λ , ν R , and the initial time is t 0 0 . Then there is the following inequality:
ψ ( t ) ψ ν λ E μ ( λ ( t t 0 ) μ ) + ν λ .
Lemma 3
([32]). Assuming h(t) is a continuous and differentiable function defined on the positive real numbers R + , then
D t μ a C h ( t ) h * h * l n h ( t ) h * ( 1 h ( t ) h * ) D t μ a C h ( t ) ,
where h * R + , μ ( 0 , 1 ] .

3. Model Derivation

In this section, we present a fractional-order SIV epidemic model with a non-monotonic incidence rate. This model explores the intricate interactions between susceptible, infected, and vaccinated individuals.
We initially categorize the population into three distinct compartments. The first compartment represents the susceptible individuals ( S ( t ) ) who are at risk of contracting the disease. The second compartment comprises the infected individuals ( I ( t ) ) who have already been infected and pose a risk of transmitting the disease. Lastly, the third compartment includes those who have been vaccinated and have gained immunity ( V ( t ) ). However, it is important to note that vaccinated individuals may experience a loss of immunity over time, which implies they may either become infected again or revert back to the susceptible compartment. In reality, vaccines rarely offer one hundred percent efficacy. As such, we incorporated the rate of immunity loss and the vaccine’s effectiveness into our model to better understand their impact on the disease dynamics. The flowchart of SIV prevalence is shown in Figure 1.
The SIV model incorporating vaccination strategies can be represented by the following interconnected nonlinear differential equations:
S ( t ) = Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S I ( t ) = β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I V ( t ) = Λ q + p S β σ V 1 + m I 2 I ( d + ε ) V
The meanings of the other parameters within the model are detailed in Table 1 below.
As we discussed earlier, compared to integer-order calculations, the nonlocal characteristics of fractional-order differential equations are more suitable for accurately describing dynamic systems with characteristics such as time dependence, path dependence, memory, and heredity. Therefore, the aforementioned ordinary differential model (9) can be further extended to the following θ -order fractional-order system. Using the definition of Caputo’s fractional derivative, system (9) can be rewritten as
D t θ 0 C S ( t ) = Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S D t θ 0 C I ( t ) = β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I D t θ 0 C V ( t ) = Λ q + p S β σ V 1 + m I 2 I ( d + ε ) V

4. Model Analysis

4.1. Non-Negativity and Boundedness of Solutions

Let A + = { ( S , I , V ) A : S , I , V R + } , where R + encompasses all positive real numbers.
Theorem 1.
Every vertex of system (10) lies within A + , maintains non-negativity, and is uniformly bounded over t 0 .
Proof of Theorem 1.
Non-negativity: Suppose the system has an equilibrium at X ( t i ) = ( S ( t i ) , I ( t i ) , V ( t i ) ) A + at the initial time. From the system (10), we have
D S θ | S t 0 = 0 = Λ ( 1 q ) + ε V + γ I 0 ,
D I θ | I t 0 = 0 = 0 ,
D V θ | V t 0 = 0 = Λ q + p S 0 .
By Lemma 1, for t t 0 , we obtain S ( t ) , I ( t ) , V ( t ) 0 . Therefore, under the initial conditions, all the vertices of the system are finally in A + .
Boundedness: View the overall population as the aggregated total of the three segments, i.e., N ( t ) = S ( t ) + I ( t ) + V ( t ) . Then
D θ N ( t ) = D θ S ( t ) + D θ I ( t ) + D θ V ( t ) = Λ d N ( t ) α I ( t ) Λ d N ( t )
D θ N ( t ) + d N ( t ) Λ .
By Lemma 2, it follows that N ( t ) N ( t 0 ) Λ d E θ d ( t t 0 ) θ + Λ d , when t , N ( t ) Λ d holds true. Thus, the proof of boundedness is complete. In conclusion, the proof of non-negative boundedness is complete. □

4.2. Existence and Uniqueness of Solutions

The subsequent lemma serves as the foundation for verifying the existence and uniqueness of the solution.
Lemma 4
([32]). Consider the fractional-order system presented below:
D μ h ( t ) = g ( t , y ) , y ( t 0 ) = y 0 , and t 0 > 0 .
where μ ( 0 , 1 ] , g : t 0 , + × A R n . If g ( t , x ) satisfies the Lipschitz condition with respect to y, that is, there exists a constant L 0 such that g ( t , y 1 ) g ( t , y 2 ) L y 1 y 2 for all ( t , y 1 ) , ( t , y 2 ) N , then for the system there exists a unique solution in [ t 0 , ) × A .
Theorem 2.
Given any initial condition X t 0 = ( S t 0 , I t 0 , V t 0 ) within the positive space A + , the system (10) possesses a unique solution X ( t ) = ( S ¯ ( t ) , I ¯ ( t ) , V ¯ ( t ) ) that persists within A + for all later times t > t 0 .
Proof of Theorem 2.
Consider the region A + × [ t 0 , t 1 ] , where A + = { ( S , I , V ) R 3 : max S , I , V k } . Let t 1 be a finite number. Denote X ( t ) = ( S ( t ) , I ( t ) , V ( t ) ) and X ¯ ( t ) = S ¯ ( t ) , I ¯ ( t ) , V ¯ ( t ) , abbreviated as X ( t ) = X and X ¯ ( t ) = X ¯ .
Suppose P ( X ) = P 1 ( X ) , P 2 ( X ) , P 3 ( X ) , where
P 1 ( X ) = Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S ,
P 2 ( X ) = β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I ,
P 3 ( X ) = Λ q + p S β σ I 1 + m I 2 V ( d + ε ) V .
For X , X ¯ A + ,
P ( X ) P ( X ¯ ) = P 1 ( X ) P 1 ( X ¯ ) + P 2 ( X ) P 2 ( X ¯ ) + P 3 ( X ) P 3 ( X ¯ ) = Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S Λ ( 1 q ) + ε V ¯ + γ I ¯ ( β I ¯ 1 + m I ¯ 2 + p + d ) S ¯ + β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I β ( S ¯ + σ V ¯ ) 1 + m I ¯ 2 I ¯ ( γ + d + α ) I ¯ + Λ q + p S β σ I 1 + m I 2 V ( d + ε ) V Λ q + p S ¯ β σ I ¯ 1 + m I ¯ 2 V ¯ ( d + ε ) V ¯ V V ¯ + γ I I ¯ + ( γ + d + α ) I I ¯ + p S S ¯ + ( d + ε ) V V ¯ + 2 β S I 1 + m I 2 S ¯ I ¯ 1 + m I ¯ 2 + ( p + d ) S S ¯ + 2 β σ I V 1 + m I 2 σ I ¯ V ¯ 1 + m I ¯ 2 ( 2 p + d + 2 β k ) S S ¯ + ( 2 γ + d + α + 2 β σ k ) I I ¯ + ( 2 ε + d ) V V ¯ M 1 S S ¯ + M 2 I I ¯ + M 3 V V ¯ M X X ¯
where
M = max { M 1 , M 2 , M 3 } , M 1 = 2 p + d + 2 β k , M 2 = 2 γ + d + α + 2 β σ k , M 3 = 2 ε + d .
Hence, based on Lemma 4, the uniqueness of the system’s solution is made more unique by the fact that P ( X ) satisfies Lipschitz’s condition. □

4.3. Existence of Equilibria

To ascertain the equilibrium state of system (10), we set the expressions on the right side of the equations to null, effectively converting the system (10) into a form:
Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S = 0 β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I = 0 Λ q + p S β σ I 1 + m I 2 V ( d + ε ) V = 0
After simplification and calculation, two equilibrium points are obtained: the disease-free equilibrium point E 0 = ( Λ ε + d ( 1 q ) d ( d + p + ε ) , 0 , Λ ( d q + p ) d ( d + p + ε ) ) and the endemic equilibrium point E e = ( S e , I e , V e ) , where
S e = d ( 1 + m I e 2 ) ( γ + d + α ) β σ Λ I e ( d + α ) β d ( 1 σ ) ,
V e = β ( Λ I e ( d + α ) ) d ( γ 1 + d + α ) ( 1 + m I e 2 ) β d ( 1 σ ) ,
and I e can represent it with the equation
A I e 4 + B I e 3 + C I e 2 + D I e + E = 0 ,
where
A = d m 2 ( γ + d + α ) ( p + d + ε ) , B = β m ( d + α ) d + ( p + d ) σ + ε + d γ σ , C = ( d + α ) β 2 σ + 2 d m ( p + d + ε ) + 2 d γ m ( p + d + ε ) β Λ m ( p + d ) σ + d ( 1 q ) ( 1 σ ) + ε , D = β ( d + α ) d + ( p + d ) σ + ε β σ Λ + d γ σ , E = d ( γ + d + α ) ( p + d + ε ) β Λ ( p + d ) σ + d ( 1 q ) ( 1 σ ) + ε .

4.4. Basic Reproduction Number

The basic reproduction number is a parameter that describes the transmission capacity of a certain infectious disease within a susceptible population [33]. It signifies the anticipated quantity of additional infections that can arise from a single infected person throughout their contagious duration at the very beginning of a disease outbreak when all members of the population are susceptible. This paper calculates the basic reproduction number through the next-generation matrix method, drawing on the pioneering approach established by Driessche and Watmough [34,35].
In this system, F signifies the growth rate of secondary infections within the infected compartment, and V represents the reduction rate due to disease progression, death, and recovery within the infected compartment.
F = β ( S + σ V ) 1 + m I 2 I , V = ( γ + d + α ) I .
Hence, the Jacobian matrices of F and V evaluated at the disease-free equilibrium point E 0 are presented as follows:
F = β Λ ε + d ( 1 q ) d ( d + ε + p ) + σ Λ ( d q + p ) d ( d + ε + p ) , V = γ + d + α .
Therefore,
R 0 = ρ 0 ( F V 1 ) = β Λ ε + d ( 1 σ ) d q + σ p d ( d + ε + p ) ( γ + d + α ) .

5. Stability Analysis of the Equilibria

This section delves into the stability analysis of both the disease-free equilibrium and the endemic equilibrium, examining them from two distinct perspectives: locally asymptotically stable and globally asymptotically stable. First, the Jacobian matrix of the system (10) at an arbitrary point X = ( S , I , V ) is constructed as follows:
J ( S , I , V ) = ( β I 1 + m I 2 + p + d ) γ ε β I 1 + m I 2 β ( S + σ V ) ( 1 m I 2 ) ( 1 + m I 2 ) 2 ( γ + d + α ) β σ I 1 + m I 2 p β σ V ( 1 m I 2 ) ( 1 + m I 2 ) 2 β σ I 1 + m I 2 d ε
Theorem 3.
If R 0 < 1 , then the disease-free equilibrium point E 0 of the system (10) is locally asymptotically stable.
Proof of Theorem 3.
At the disease-free equilibrium E 0 of the system (10), the Jacobian matrix is
J ( S 0 , I 0 , V 0 ) = ( p + d ) γ ε 0 β ( S 0 + σ V 0 ) ( ( γ + d + α ) 0 p β σ V 0 d ε = ( p + d ) γ ε 0 β Λ [ ε + d + σ p ( 1 σ ) d q ] d ( p + d + ε ) ( ( γ + d + α ) 0 p β σ Λ ( d q + p ) d ( p + d + ε ) d ε
The corresponding characteristic equation at E 0 is
λ E J ( E 0 ) = λ β Λ ε + d + σ p ( 1 σ ) d q d ( p + d + ε ) + γ + d + α ( λ + p + d ) ( λ + d + ε ) ε p = 0 .
Through observation of the characteristic equation, one can find that at E 0 , one of the eigenvalues is
λ 1 = β Λ ε + d + σ p ( 1 σ ) d q d ( p + d + ε ) ( γ + d + α ) = ( γ + d + α ) ( R 0 1 ) .
When R 0 < 1 , it is evident that λ 1 < 0 . It can be found that λ 2 and λ 3 satisfy the equation:
( λ + p + d ) ( λ + d + ε ) ε p = 0 .
After simplification, it becomes
λ 2 + ( p + 2 d + ε ) λ + p d + d 2 + d ε = 0 .
It is easy to know p + 2 d + ε > 0 and p d + d 2 + d ε > 0 ; then the Routh–Hurwitz conditions are satisfied, confirming that under the condition R 0 < 1 , the point E 0 is locally asymptotically stable. □
Theorem 4.
If R 0 > 1 , then the endemic equilibrium point E e of the system (10) is locally asymptotically stable if a 2 > 0 , a 1 > 0 , a 0 > 0 and a 2 a 1 a 0 > 0 , where
a 2 = γ + 3 d + ε + p + α ( 1 m I e 2 ) ( γ + d + α ) β I e ( σ 1 ) 1 + m I e 2 , a 1 = ( d + p + β I e 1 + m I e 2 ) ( d + ε + β I e 1 + m I e 2 ) ε p + d + α β S e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 2 d + ε + p + β I e ( σ + 1 ) 1 + m I e 2 β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ ( p + d + ε ) + β σ I e 1 + m I e 2 β V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ + β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ , a 0 = d + α β S e ( 1 m I e 2 ) ( 1 m I e 2 ) 2 ( d + p + β I e 1 + m I e 2 ) ( d + ε + β σ I e 1 + m I e 2 ) ε p β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ ( p + d + ε ) d β σ I e 1 + m I e 2 β V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ d .
Proof of Theorem 4.
At the endemic equilibrium point E e of the system (10), the Jacobian matrix is
J ( S e , I e , V e ) = ( β I e 1 + m I e 2 + p + d ) γ ε β I e 1 + m I e 2 β ( S e + σ V e ) ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 ( γ + d + α ) β σ I e 1 + m I e 2 p β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 β σ I 1 + m I e 2 d ε .
Subsequently, the characteristic equation of the aforementioned matrix is expressed as
λ 3 + a 2 λ 2 + a 1 λ + a 0 = 0 ,
where
a 2 = γ + 3 d + ε + p + α ( 1 m I e 2 ) ( γ + d + α ) β I e ( σ 1 ) 1 + m I e 2 ; a 1 = ( d + p + β I e 1 + m I e 2 ) ( d + ε + β I e 1 + m I e 2 ) ε p + d + α β S e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 2 d + ε + p + β I e ( σ + 1 ) 1 + m I e 2 β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ ( p + d + ε ) + β σ I e 1 + m I e 2 β V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ + β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ ; a 0 = d + α β S e ( 1 m I e 2 ) ( 1 m I e 2 ) 2 ( d + p + β I e 1 + m I e 2 ) ( d + ε + β σ I e 1 + m I e 2 ) ε p β σ V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ ( p + d + ε ) d β σ I e 1 + m I e 2 β V e ( 1 m I e 2 ) ( 1 + m I e 2 ) 2 γ d .
According to the Routh–Hurwitz stability criteria, if a 2 > 0 , a 1 > 0 , a 0 > 0 and a 2 a 1 a 0 > 0 , the endemic equilibrium point E e is locally asymptotically stable. □
Theorem 5.
If R 0 < 1 , d ε > d q + p ε d ( 1 q ) , and d p > ε d ( 1 q ) d q + p , then the disease-free equilibrium point E 0 of the system (10) is globally asymptotically stable.
Proof of Theorem 5.
Let us consider the Lyapunov function as
L 1 ( S , I , V ) = S S 0 S 0 l n S S 0 + I + V V 0 V 0 l n V V 0 .
Then, by Lemma 3, we have
D θ L 1 ( 1 S 0 S ) D θ S + D θ I + ( 1 V 0 V ) D θ V = ( 1 S S 0 ) Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S + β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I + ( 1 V V 0 ) Λ q + p S β σ I 1 + m I 2 V ( d + ε ) V = Λ d ( S + I + V ) α I S 0 S Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S V 0 V Λ q + p S β σ I 1 + m I 2 V ( d + ε ) V
For R 0 < 1 , the steady-state equation at the equilibrium point is
Λ ( 1 q ) = ( p + d ) S 0 ε V 0 ,
Λ q = ( d + ε ) V 0 p S 0 ,
Then
D θ L 1 d S 0 + d V 0 d ( S + I + V ) α I S 0 S γ I + β I 1 + m I 2 S 0 + β σ I 1 + m I 2 V 0 S 0 S ε V + ( p + d ) S 0 V 0 V p S + ( d + ε ) V 0 ( p + d ) S 0 2 S + ε S 0 V 0 S ( d + ε ) V 0 2 V + p S 0 V 0 V ( d S 0 ε V 0 ) ( 2 S S 0 S 0 S ) + ( d V 0 p S 0 ) ( 2 V V 0 V 0 V ) + p S 0 ( 3 S 0 S V V 0 V 0 S V S 0 ) + ε V 0 ( 3 S S 0 V 0 V V S 0 V 0 S )
Since the arithmetic mean consistently exceeds the geometric mean, then D θ L 1 0 holds if and only if d S 0 ε V 0 > 0 and d V 0 p S 0 > 0 , i.e., d ε > d q + p ε d ( 1 q ) and d p > ε d ( 1 q ) d q + p , where D θ L 1 = 0 if and only if S = S 0 , I = I 0 , V = V 0 . In line with the Lyapunov–LaSalle asymptotic stability theorem [36,37], E 0 is globally asymptotically stable. □
Theorem 6.
If R 0 > 1 , d S e ε V e > 0 and d V e p S e > 0 , then the endemic equilibrium point E e of the system (10) is globally asymptotically stable.
Proof of Theorem 6.
Let us consider the Lyapunov function as
L 2 ( S , I , V ) = S S e S e l n S S e + I I e I e l n I I e + V V e V e l n V V e .
Then, by Lemma 3, we have
D θ L 2 ( 1 S e S ) D θ S + ( 1 I e I ) D θ I + ( 1 V e V ) D θ V = Λ d ( S + I + V ) α I S e S Λ ( 1 q ) + ε V + γ I ( β I 1 + m I 2 + p + d ) S I e I β ( S + σ V ) 1 + m I 2 I ( γ + d + α ) I V e V Λ q + p S β σ I 1 + m I 2 V ( d + ε ) V
For R 0 > 1 , the steady-state equation at the equilibrium point is
Λ ( 1 q ) = ( β I e 1 + m I e 2 + p + d ) S e ε V e γ I e ,
γ = β ( S e + σ V e ) 1 + m I e 2 ( d + α ) ,
Λ q = β σ I e 1 + m I e 2 V e + ( d + ε ) V e p S e .
Then
D θ L 2 d ( S e + I e + V e ) + α I e d ( S + I + V ) α I ( p + d ) S e 2 S + ε S e V e S + β σ V e I e 1 m I e 2 S e S ( d + α ) S e I e S ε V S e S β ( S e + σ V e ) 1 + m I e 2 ( d + α ) S e I S + ( β I 1 + m I 2 ) S e β ( S + σ V ) 1 + m I 2 I e + β ( S e + σ V e ) 1 + m I e 2 I e ( β σ I e 1 + m I 2 + d + ε ) V e 2 V + p V e S e V p V e S V + ( β σ I 1 + m I 2 + p + d ) V e ( d S e ε V e ) ( 2 S S e S e S ) + ( d V e p S e ) ( 2 V V e V e V ) + p S e ( 3 S e S V V e V e S V S e ) + ε V e ( 3 S S e V e V V S e V e S )
Since the arithmetic mean consistently exceeds the geometric mean, then D θ L 2 0 holds if and only if d S e ε V e > 0 and d V e p S e > 0 , where D θ L 2 = 0 if and only if S = S e , I = I e , V = V e . In line with the Lyapunov–LaSalle asymptotic stability theorem [36,37], E e is globally asymptotically stable. □

6. Fractional Optimal Control

In this section, a comprehensive control strategy is proposed for the aforementioned infectious disease model, striving to reduce the disease’s effects while also maintaining control expenditures at a minimum. The goal of these control measures is to lower the incidence of infections within the population, which necessitates the creation of an optimal control problem to realize the desired outcome. During the epidemic transmission process, by changing the behaviors of susceptible individuals, such as maintaining social distance, improving awareness, and prevention of the disease, the rate at which susceptible individuals become infected can be decreased, which is represented by the control variable u 1 ( t ) . By increasing the vaccination rate among susceptible individuals, for instance, through promoting vaccination campaigns, offering vaccine subsidies, or expanding vaccine supply, we can significantly reduce their risk of contracting the disease. This is represented by the control variable u 2 ( t ) . u 3 ( t ) , on the other hand, signifies the extent of therapeutic measures undertaken for infected individuals (those who have already contracted the disease) at time t. These measures may encompass pharmacological interventions (such as medication and vaccine therapy) and physical segregation practices (like isolation wards and home quarantine). Upon enhancing the aforementioned three control strategies, System (10) undergoes the following modification:
D t θ 0 C S ( t ) = Λ ( 1 q ) + ε V + ( γ + u 3 ( t ) ) I ( 1 u 1 ( t ) ) β S I 1 + m I 2 ( p + u 2 ( t ) + d ) S D t θ 0 C I ( t ) = ( 1 u 1 ( t ) ) β S I 1 + m I 2 + β σ V I 1 + m I 2 ( γ + u 3 ( t ) + d + α ) I D t θ 0 C V ( t ) = Λ q + ( p + u 2 ( t ) ) S β σ V 1 + m I 2 I ( d + ε ) V
In this scenario, the initial conditions, which are non-negative, are provided as
S ( 0 ) = S 0 , I ( 0 ) = I 0 , V ( 0 ) = V 0 .
The designated control set is clearly outlined as
U = { ( u 1 , u 2 , u 3 ) | u i i s L e b e s g u e m e a s u r a b l e o n [ 0 , 1 ] , i = 1 , 2 , 3 } .
Our overarching control objective involves minimizing both the proportion of infected individuals within the population and the associated expenses incurred through the utilization of control variables u i ( t ) (where i = 1 , 2 , 3 ). Therefore, we define the objective function
J ( u 1 , u 2 , u 3 ) = 0 T ( B 1 I + 1 2 C 1 u 1 2 + 1 2 C 2 u 2 2 + 1 2 C 3 u 3 2 ) d t
subject to the model (37) with initial conditions, where T represents the control period, B 1 denotes the weight for infected individuals, and C 1 , C 2 , and C 3 serve as the weighting constants, each corresponding to the financial burden associated with implementing the respective control measures. 1 2 C 1 u 1 2 , 1 2 C 2 u 2 2 , and 1 2 C 3 u 3 2 are the implementation costs for these three measures, respectively. Our main objective is to determine the optimal control U * = ( u 1 * , u 2 * , u 3 * ) U for the system, guaranteeing that the objective function reaches its minimal value, that is,
J ( u 1 * , u 2 * , u 3 * ) = min u U J ( u 1 , u 2 , u 3 )
Next, we apply Pontryagin’s maximum principle [38] to solve the optimal control problem. First of all, we need to prove that the optimal solution exists for the system (37).
Theorem 7.
There exists an optimal control solution ( u 1 * , u 2 * , u 3 * ) such that
J ( u 1 * , u 2 * , u 3 * ) = min u U J ( u 1 , u 2 , u 3 ) ,
subject to the fractional system (37).
Proof of Theorem 7.
From the aforementioned analysis of the model (10) and the results in [39], the existence of the optimal control U * is contingent upon the fulfillment of the following conditions:
(i) The set of control U and the corresponding set of state variables are non-empty;
(ii) U is closed and convex;
(iii) The right side of the system (37) is bounded by a linear function with the control and state variables;
(iv) The integrand of the objective function
L ( I , u 1 , u 2 , u 3 ) = B 1 I + 1 2 C 1 u 1 2 + 1 2 C 2 u 2 2 + 1 2 C 3 u 3 2
exhibits convexity within the set U;
(v) There exist constants K 1 , K 2 > 0 and Q > 1 such that the integrand L ( I , u 1 , u 2 , u 3 ) satisfies
L ( I , u 1 , u 2 , u 3 ) K 1 ( | u 1 | 2 + | u 2 | 2 + | u 3 | 2 ) Q 2 K 2 .
 □
To derive the optimal control solution, we introduce Lagrangian function
L ( I , u 1 , u 2 , u 3 ) = B 1 I + 1 2 C 1 u 1 2 + 1 2 C 2 u 2 2 + 1 2 C 3 u 3 2
and Hamiltonian function
H ( S , I , V , u , λ ) = L ( I , u 1 , u 2 , u 3 ) + λ 1 D t θ 0 C S ( t ) + λ 2 D t θ 0 C I ( t ) + λ 3 D t θ 0 C V ( t ) ,
where λ i ( t ) for i = 1 , 2 , 3 are the adjoint variables. By applying Pontryagin’s maximum principle, we derive the adjoint variable equation and the transversal condition as indicated below:
D θ λ 1 ( t ) = H S = ( 1 u 1 ( t ) ) β I 1 + m I 2 + ( p + u 2 ( t ) + d ) λ 2 ( 1 u 1 ( t ) ) β I 1 + m I 2 λ 2 ( p + u 2 ( t ) ) λ 3 D θ λ 2 ( t ) = H I = B 1 γ + u 3 ( t ) ( 1 u 1 ( t ) ) β S ( 1 m I 2 ) ( 1 + m I 2 ) 2 λ 1 ( 1 u 1 ( t ) ) β S + β σ V ( 1 m I 2 ) ( 1 + m I 2 ) 2 ( γ + u 3 ( t ) + d + α ) λ 2 + β σ V ( 1 m I 2 ) ( 1 + m I 2 ) 2 λ 3 D θ λ 3 ( t ) = H V = ε λ 1 β σ I 1 + m I 2 λ 2 + ( β σ I 1 + m I 2 d ε ) λ 3
satisfying transversality conditions
λ i ( T ) = 0 , i = 1 , 2 , 3 .
To proceed, let us designate S ¯ , I ¯ , V ¯ as the optimal values corresponding to the variables S , I , V . Additionally, let us acknowledge λ ¯ 1 , λ ¯ 2 , λ ¯ 3 as the solutions to Equation (41). With these clarifications in mind, we are now poised to present the subsequent theorem.
Theorem 8.
Consider U * as the optimal value for the control parameter U; hence, U * can be expressed as the equations:
u 1 * = max { 0 , min { β S ¯ I ¯ ( λ ¯ 2 λ ¯ 1 ) c 1 ( 1 + m I ¯ 2 ) , 1 } } u 2 * = max { 0 , min { S ¯ ( λ ¯ 1 λ ¯ 3 ) c 2 , 1 } } u 3 * = max { 0 , min { I ¯ ( λ ¯ 2 λ ¯ 1 ) c 3 , 1 } }
Proof of Theorem 8.
To derive the characteristic equation of optimal control U * , proceed by solving the equations:
H u 1 = c 1 u 1 * + β S ¯ I ¯ 1 + m I ¯ 2 λ ¯ 1 β S ¯ I ¯ 1 + m I ¯ 2 λ ¯ 2 = 0 H u 2 = c 2 u 2 * S ¯ λ ¯ 1 + S ¯ λ ¯ 3 = 0 H u 3 = c 3 u 3 * + I ¯ λ ¯ 1 I ¯ λ ¯ 2 = 0
Therefore,
u 1 * = β S ¯ I ¯ ( λ ¯ 2 λ ¯ 1 ) c 1 ( 1 + m I ¯ 2 ) ; u 2 * = S ¯ ( λ ¯ 1 λ ¯ 3 ) c 2 ; u 3 * = I ¯ ( λ ¯ 2 λ ¯ 1 ) c 3 .
The optimal control solution U * can be expressed as
u 1 * = 0 , β S ¯ I ¯ ( λ ¯ 2 λ ¯ 1 ) c 1 ( 1 + m I ¯ 2 ) 0 , β S ¯ I ¯ ( λ ¯ 2 λ ¯ 1 ) c 1 ( 1 + m I ¯ 2 ) , 0 < β S ¯ I ¯ ( λ ¯ 2 λ ¯ 1 ) c 1 ( 1 + m I ¯ 2 ) < 1 , 1 , β S ¯ I ¯ ( λ ¯ 2 λ ¯ 1 ) c 1 ( 1 + m I ¯ 2 ) 1 . u 2 * = 0 , S ¯ ( λ ¯ 1 λ ¯ 3 ) c 2 0 , S ¯ ( λ ¯ 1 λ ¯ 3 ) c 2 , 0 < S ¯ ( λ ¯ 1 λ ¯ 3 ) c 2 < 1 , 1 , S ¯ ( λ ¯ 1 λ ¯ 3 ) c 2 1 . u 3 * = 0 , I ¯ ( λ ¯ 2 λ ¯ 1 ) c 3 0 , I ¯ ( λ ¯ 2 λ ¯ 1 ) c 3 , 0 < I ¯ ( λ ¯ 2 λ ¯ 1 ) c 3 < 1 , 1 , I ¯ ( λ ¯ 2 λ ¯ 1 ) c 3 1 .
given that the control variable ( U ) is bounded. This concludes the demonstration of Theorem 8. □

7. Numerical Simulation

Given the intricacy of obtaining analytical solutions, the use of numerical methods becomes essential. This section utilizes the fractional Adams–Bashforth–Moulton method [40] and particle swarm optimization [41] for numerical experimentation via MATLAB R2023(b), aiming to validate some of the theoretical findings discussed earlier in the paper. Additionally, numerous numerical simulations have been carried out, encompassing a wide range of biologically feasible parameter spaces, yielding a comprehensive array of dynamic outcomes across various scenarios. With the scale for population set at 1 million individuals and the time scale at one day, the model’s parameters are defined as follows: α = 0.15 , β = 0.45 , d = 0.15 , γ = 0.25 , p = 0.2 , q = 0.4 , σ = 0.06 , Λ = 0.35 , ε = 0.05 , m = 0.3 , and the initial populations in each compartment are set as S ( 0 ) = 4.5 , I ( 0 ) = 0.2 , V ( 0 ) = 0.3 .

7.1. Numerical Simulation without Control Measures

In this subsection, we present numerical results for a fractional system. Figure 2 and Figure 3 showcase the dynamic behavior of the susceptible, infected, and vaccinated classes across various fractional orders θ . Figure 2 demonstrates that when the parameter takes the aforementioned value, which corresponds to R 0 = 0.7426 , the disease-free equilibrium point E 0 is locally asymptotically stable, validating Theorem 3. Additionally, it is evident that the fractional order θ holds a pivotal position in governing the dynamics of these three distinct compartments and influences the duration required for each to attain stability. Within the community, an augmentation in the value of θ is observed to be directly proportional to an elevation in the population of both susceptible and vaccinated individuals.
Figure 3 confirms that when β = 1.05 , d = 0.1 , p = 0.15 , which corresponds to R 0 = 2.9743 , the endemic equilibrium point E e is locally asymptotically stable. For the infected compartment, especially in the endemic state, as depicted in Figure 3b, a reduction in the value of θ leads to an enhancement in the system’s memory capacity, subsequently causing a deceleration in the rate of convergence. This means that it takes longer to eliminate the disease. Figure 4 illustrates the impact of the vaccination ratio q for newborns on the model. Here comes an interesting finding: the change in q significantly affects the variations in S and V. When the vaccination ratio escalates, the quantity of S diminishes, whereas the count of V augments. However, the change in q has a minor effect on I. This is because by analyzing the expression of R 0 , we can see that a higher q leads to a lower value of R 0 , making the disease more controllable.
Figure 5 demonstrates the influence of the vaccination ratio p among the susceptible population on the model. Clearly, when the p value decreases, infections accumulate at a gradual pace, eventually contributing to the widespread dissemination of the disease and its evolution into an endemic state. Figure 6 shows the effect of vaccine effectiveness σ on the model. A higher value of σ indicates lower vaccine effectiveness. It can be observed that the level of vaccine effectiveness has a significant impact on both the infected population and the vaccinated population. Higher effectiveness leads to a lower peak of infection, and the disease eventually subsides, along with a reduction in the number of vaccinations, thus saving on vaccination costs.
When considering a model that does not include vaccination, we set the parameters related to vaccination to zero, p = q = ε = 0 , and σ = 1 . In this case, the basic reproduction number of the model without vaccination is R 0 w v = β Λ d ( d + α + γ ) . We observe that R 0 = ε + d + σ p ( 1 σ ) d q d + ε + p R 0 w v , thus R 0 < R 0 w v . This correlation highlights the influence that vaccination has on diminishing the basic reproduction number, making the disease more controllable.
It should be noted that R 0 is independent of m, so it does not alter the stability at the equilibrium point, which is evident in Figure 7. Although m changes the peak value of I, it does not alter its equilibrium state.

7.2. Numerical Simulation with Control Measures

In this subsection, we introduce the numerical results of the model after the introduction of control measures. To more effectively manage the dissemination of diseases and mitigate the expenses associated with implementing preventive measures, four distinct control strategies were devised to combat the dissemination of the disease:
Strategy A: Using isolation, vaccination, and treatment control (i.e., u 1 , u 2 , u 3 0 ).
Strategy B: Using vaccination and treatment control (i.e., u 2 , u 3 0 and u 1 = 0 ).
Strategy C: Using isolation and treatment control control (i.e., u 1 , u 3 0 and u 2 = 0 ).
Strategy D: Using isolation and vaccination control (i.e., u 1 , u 2 0 and u 3 = 0 ).
Figure 8 illustrates the variation in the objective function J ( u ) in the four strategies. It is evident that strategy A yields the most effective control outcome. Figure 9 depicts the curves of the control variables u i , i = 1 , 2 , 3 . varying with time under different control strategies. Moreover, as seen in Figure 9a, for control strategy A, control measure u 1 should be implemented on a large scale in the initial days and then stopped after 13 days. Control measure u 2 needs to maintain the control intensity at around 0.213 after the 22nd day, while control measure u 3 should maintain the control intensity at approximately 0.096 after the 37th day. Figure 9b illustrates that when contrasted with control strategy A, although the final control intensities are not significantly different, the peak intensity of control measure u 3 in strategy B lasts for approximately 13 days. As seen in Figure 9c, in comparison to strategy A, the control measures u 1 and u 3 in strategy C have both increased in terms of time and control intensity. Figure 9d indicates that the control measures u 1 and u 2 in strategy D have also increased in terms of duration and the strength of control applied; however, the rise is not as significant as in control strategy C. As a result, the influence of the control measure u 2 surpasses that of u 3 in terms of effectiveness.
Figure 10 illustrates the variations in three variables, S, I, and V, in the model over time when there is no control and when different control strategies are used. Specifically, Figure 10b shows that after the introduction of control measures, the disease will eventually disappear, and when control strategy A is applied, the disease disappears first. Moreover, as the control measure u 3 of strategy C maintains the highest control intensity, Figure 10a reveals that strategy C has the largest susceptible population. Similarly, the control measure u 2 maintains the highest intensity in strategy D, and Figure 10c indicates that when strategy D is used, the vaccinated population reaches the highest level compared to other strategies.

8. Conclusions

In this paper, we propose a fractional-order SIV epidemic model featuring a nonmonotonic incidence rate and a continuous vaccination strategy. Initially, to maintain this study’s relevance in biology, we examine the model’s dynamic characteristics from a fractional-order perspective, which involves verifying the existence and uniqueness of solutions as well as confirming their non-negativity and boundedness.Then, we identify two equilibrium points in the model, namely, the disease-free equilibrium E 0 and the endemic equilibrium E e , and calculate the basic reproduction number R 0 . Furthermore, we investigate the necessary conditions for ensuring the local and global stability of both E 0 and E e . We demonstrate the critical role of vaccination in disease elimination by investigating its influence on R 0 . Additionally, we develop an optimal control strategy to reduce the spread of the disease. In the model, we introduce three control variables, u 1 , u 2 , and u 3 , representing isolation, vaccination, and treatment, respectively. Drawing upon Pontryagin’s maximum principle, we formulate the necessary optimality conditions for the fractional-order optimal control problem and introduce particle swarm optimization to enhance the effectiveness of disease control. Finally, the results of numerical simulations indicate that strategy A is superior in terms of curbing the spread of disease.
In future work, our model can take the isolation compartment into consideration, which may be an interesting extension. Bifurcation analysis can also be included in the subsequent model analysis. It is also a good choice to build a model in a stochastic system.

Author Contributions

J.Y.: methodology, software, formal analysis, writing—original draft, and writing—review and editing. W.W.: validation and writing—review and editing. Q.M.: project administration. X.T.: conceptualization, visualization, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Nos. 11361104, 12261104), the Youth Talent Program of Xingdian Talent Support Plan (XDYC-QNRC-2022-0514), the Yunnan Provincial Basic Research Program Project (No. 202301AT070016, No. 202401AT070036), and the Science Research Fund of Education Department of Yunnan Province (No. 2024Y468).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article as no data sets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that there is no conflicts of interest regarding the publication of this paper.

References

  1. Khan, T.; Ullah, Z.; Ali, N.; Zaman, G. Modeling and control of the hepatitis B virus spreading using an epidemic model. Chaos Solitons Fractals 2019, 124, 1–9. [Google Scholar] [CrossRef]
  2. Nill, F. Endemic oscillations for SARS-COV-2 Omicron—A SIRS model analysis. Chaos Solitons Fractals 2023, 173, 113678. [Google Scholar] [CrossRef] [PubMed]
  3. Girardi, P.; Gaetan, C. An SEIR Model with Time-Varying Coefficients for Analyzing the SARS-CoV-2 Epidemic. Risk Anal. 2023, 43, 144–155. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, G.; Li, Z.; Din, A.; Chen, T. Dynamic analysis and optimal control of a stochastic COVID-19 model. Math. Comput. Simul. 2024, 215, 498–517. [Google Scholar] [CrossRef]
  5. Vargas-De-León, C. On the global stability of SIS, SIR and SIRS epidemic models with standard incidence. Chaos Solitons Fractals 2011, 44, 1106–1110. [Google Scholar] [CrossRef]
  6. Memon, Z.; Qureshi, S.; Memon, B.R. Assessing the role of quarantine and isolation as control strategies for COVID-19 outbreak: A case study. Chaos Solitons Fractals 2021, 144, 110655. [Google Scholar] [CrossRef] [PubMed]
  7. Anderson, R.M.; May, R.M. Infectious Diseases of Humans: Dynamics and Control; Oxford University Press: Oxford, UK, 1991. [Google Scholar]
  8. Alexander, M.E.; Bowman, C.; Moghadas, S.M.; Summers, R.; Gumel, A.B.; Sahai, B.M. A vaccination model for transmission dynamics of influenza. SIAM J. Appl. Dyn. Syst. 2004, 3, 503–524. [Google Scholar] [CrossRef]
  9. Elbasha, E.H.; Gumel, A.B. Theoretical assessment of public health impact of imperfect prophylactic HIV-1 vaccines with therapeutic benefits. Bull. Math. Biol. 2006, 68, 577–614. [Google Scholar] [CrossRef]
  10. Kermack, W.O.; McKendrick, A.G. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. Ser. A Contain. Pap. Math. Phys. Character 1927, 115, 700–721. [Google Scholar]
  11. Brauer, F.; Van den Driessche, P.; Wu, J.; Allen, L.J. Mathematical Epidemiology; Springer: Berlin/Heidelberg, Germany, 2008; Volume 1945. [Google Scholar]
  12. Brauer, F.; Castillo-Chavez, C.; Castillo-Chavez, C. Mathematical Models in Population Biology and Epidemiology; Springer: Berlin/Heidelberg, Germany, 2012; Volume 2. [Google Scholar]
  13. Capasso, V.; Serio, G. A generalization of the Kermack-McKendrick deterministic epidemic model. Math. Biosci. 1978, 42, 43–61. [Google Scholar] [CrossRef]
  14. Liu, W.-m.; Levin, S.A.; Iwasa, Y. Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models. J. Math. Biol. 1986, 23, 187–204. [Google Scholar] [CrossRef] [PubMed]
  15. Jang, S.R.J. Backward bifurcation in a discrete SIS model with vaccination. J. Biol. Syst. 2008, 16, 479–494. [Google Scholar] [CrossRef]
  16. Parsamanesh, M.; Mehrshad, S. Stability of the equilibria in a discrete-time sivs epidemic model with standard incidence. Filomat 2019, 33, 2393–2408. [Google Scholar] [CrossRef]
  17. Farnoosh, R.; Parsamanesh, M. Disease extinction and persistence in a discrete-time SIS epidemic model with vaccination and varying population size. Filomat 2017, 31, 4735–4747. [Google Scholar] [CrossRef]
  18. Parsamanesh, M.; Erfanian, M.; Mehrshad, S. Stability and bifurcations in a discrete-time epidemic model with vaccination and vital dynamics. BMC Bioinform. 2020, 21, 1–15. [Google Scholar] [CrossRef] [PubMed]
  19. Parsamanesh, M.; Erfanian, M. Stability and bifurcations in a discrete-time SIVS model with saturated incidence rate. Chaos Solitons Fractals 2021, 150, 111178. [Google Scholar] [CrossRef]
  20. Elbaz, I.M.; Sohaly, M.A.; El-Metwally, H. Random dynamics of an SIV epidemic model. Commun. Nonlinear Sci. Numer. Simul. 2024, 131, 107779. [Google Scholar] [CrossRef]
  21. Baleanu, D.; Sajjadi, S.S.; Asad, J.H.; Jajarmi, A.; Estiri, E. Hyperchaotic behaviors, optimal control, and synchronization of a nonautonomous cardiac conduction system. Adv. Differ. Equations 2021, 2021, 157. [Google Scholar] [CrossRef]
  22. Jajarmi, A.; Baleanu, D.; Zarghami Vahid, K.; Mobayen, S. A general fractional formulation and tracking control for immunogenic tumor dynamics. Math. Methods Appl. Sci. 2022, 45, 667–680. [Google Scholar] [CrossRef]
  23. Huo, J.; Zhao, H.; Zhu, L. The effect of vaccines on backward bifurcation in a fractional order HIV model. Nonlinear Anal. Real World Appl. 2015, 26, 289–305. [Google Scholar] [CrossRef]
  24. Pinto, C.M.; Carvalho, A.R. A latency fractional order model for HIV dynamics. J. Comput. Appl. Math. 2017, 312, 240–256. [Google Scholar] [CrossRef]
  25. Arqub, O.A.; El-Ajou, A. Solution of the fractional epidemic model by homotopy analysis method. J. King Saud-Univ.-Sci. 2013, 25, 73–81. [Google Scholar] [CrossRef]
  26. Jana, S.; Mandal, M.; Nandi, S.K.; Kar, T.K. Analysis of a fractional-order SIS epidemic model with saturated treatment. Int. J. Model. Simul. Sci. Comput. 2021, 12, 2150004. [Google Scholar] [CrossRef]
  27. Podlubny, I. Fractional Differential Equations; Academic: San Diego, CA, USA, 1999. [Google Scholar]
  28. Dai, L.; Liu, X.; Chen, Y. Global dynamics of a fractional-order SIS epidemic model with media coverage. Nonlinear Dyn. 2023, 111, 19513–19526. [Google Scholar] [CrossRef]
  29. Phukan, J.; Dutta, H. Dynamic analysis of a fractional order SIR model with specific functional response and Holling type II treatment rate. Chaos Solitons Fractals 2023, 175, 114005. [Google Scholar] [CrossRef]
  30. 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 Fractals 2020, 137, 109892. [Google Scholar] [CrossRef]
  31. Majee, S.; Jana, S.; Das, D.K.; Kar, T. Global dynamics of a fractional-order HFMD model incorporating optimal treatment and stochastic stability. Chaos Solitons Fractals 2022, 161, 112291. [Google Scholar] [CrossRef]
  32. Li, Y.; Chen, Y.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag–Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef]
  33. Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A.J. On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef]
  34. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef]
  35. Van den Driessche, P.; Watmough, J. Further notes on the basic reproduction number. In Mathematical Epidemiology; Springer: Berlin/Heidelberg, Germany, 2008; pp. 159–178. [Google Scholar]
  36. La Salle, J.P. The Stability of Dynamical Systems; SIAM: Bangkok, Thailand, 1976. [Google Scholar]
  37. Lyapunov, A.M. The general problem of the stability of motion. Int. J. Control 1992, 55, 531–534. [Google Scholar] [CrossRef]
  38. Kopp, R.E. Pontryagin maximum principle. In Mathematics in Science and Engineering; Elsevier: Amsterdam, The Netherlands, 1962; Volume 5, pp. 255–279. [Google Scholar]
  39. Lukes, D.L. Differential equations: Classical to controlled. Am. Math. Mon. 1982, 92, 223–225. [Google Scholar]
  40. Diethelm, K. Efficient solution of multi-term fractional differential equations using P (EC) m E methods. Computing 2003, 71, 305–319. [Google Scholar] [CrossRef]
  41. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; IEEE: Piscataway, NJ, USA, 1995; Volume 4, pp. 1942–1948. [Google Scholar]
Figure 1. Disease transmission flow chart.
Figure 1. Disease transmission flow chart.
Mathematics 12 02735 g001
Figure 2. (ac) Effects of θ on S, I, and V when R 0 < 1 .
Figure 2. (ac) Effects of θ on S, I, and V when R 0 < 1 .
Mathematics 12 02735 g002
Figure 3. (ac) Effects of θ on S, I, and V when R 0 > 1 .
Figure 3. (ac) Effects of θ on S, I, and V when R 0 > 1 .
Mathematics 12 02735 g003
Figure 4. (ac) Effect of q on S, I, and V.
Figure 4. (ac) Effect of q on S, I, and V.
Mathematics 12 02735 g004
Figure 5. (ac) Effect of p on S, I, and V.
Figure 5. (ac) Effect of p on S, I, and V.
Mathematics 12 02735 g005
Figure 6. (ac) Effect of σ on S, I, and V.
Figure 6. (ac) Effect of σ on S, I, and V.
Mathematics 12 02735 g006
Figure 7. Effect of m on I.
Figure 7. Effect of m on I.
Mathematics 12 02735 g007
Figure 8. Changes in the objective function under different strategies ( θ = 0.95 ).
Figure 8. Changes in the objective function under different strategies ( θ = 0.95 ).
Mathematics 12 02735 g008
Figure 9. Time series plot of control variables under different strategies.
Figure 9. Time series plot of control variables under different strategies.
Mathematics 12 02735 g009
Figure 10. (ac) The impact of different strategies on S, I, and V.
Figure 10. (ac) The impact of different strategies on S, I, and V.
Mathematics 12 02735 g010
Table 1. Model (10) parameter description.
Table 1. Model (10) parameter description.
ParameterDescription
Λ Population supplement rate
pVaccination rate of susceptible populations
ε Immune loss rate
γ Recovery rate of infected populations
β Contact rate
qVaccination rate of newborns
dNatural mortality rate
α Mortality rate due to disease
σ Effectiveness of the vaccine
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

Yan, J.; Wu, W.; Miao, Q.; Tan, X. Global Dynamics and Optimal Control of a Fractional-Order SIV Epidemic Model with Nonmonotonic Occurrence Rate. Mathematics 2024, 12, 2735. https://doi.org/10.3390/math12172735

AMA Style

Yan J, Wu W, Miao Q, Tan X. Global Dynamics and Optimal Control of a Fractional-Order SIV Epidemic Model with Nonmonotonic Occurrence Rate. Mathematics. 2024; 12(17):2735. https://doi.org/10.3390/math12172735

Chicago/Turabian Style

Yan, Juhui, Wanqin Wu, Qing Miao, and Xuewen Tan. 2024. "Global Dynamics and Optimal Control of a Fractional-Order SIV Epidemic Model with Nonmonotonic Occurrence Rate" Mathematics 12, no. 17: 2735. https://doi.org/10.3390/math12172735

APA Style

Yan, J., Wu, W., Miao, Q., & Tan, X. (2024). Global Dynamics and Optimal Control of a Fractional-Order SIV Epidemic Model with Nonmonotonic Occurrence Rate. Mathematics, 12(17), 2735. https://doi.org/10.3390/math12172735

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