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
(with
denoting the transmission rate) [
7] and
(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
, where g is a non-negative function with
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
, 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
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.
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 (
) who are at risk of contracting the disease. The second compartment comprises the infected individuals (
) 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 (
). 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:
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
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
. 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
.
, 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:
In this scenario, the initial conditions, which are non-negative, are provided as
The designated control set is clearly outlined as
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
(where
). Therefore, we define the objective function
subject to the model (
37) with initial conditions, where
T represents the control period,
denotes the weight for infected individuals, and
and
serve as the weighting constants, each corresponding to the financial burden associated with implementing the respective control measures.
, and
are the implementation costs for these three measures, respectively. Our main objective is to determine the optimal control
for the system, guaranteeing that the objective function reaches its minimal value, that is,
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 such thatsubject 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
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
exhibits convexity within the set
U;
(v) There exist constants
and
such that the integrand
satisfies
□
To derive the optimal control solution, we introduce Lagrangian function
and Hamiltonian function
where
for
are the adjoint variables. By applying Pontryagin’s maximum principle, we derive the adjoint variable equation and the transversal condition as indicated below:
satisfying transversality conditions
To proceed, let us designate
,
,
as the optimal values corresponding to the variables
. Additionally, let us acknowledge
as the solutions to Equation (
41). With these clarifications in mind, we are now poised to present the subsequent theorem.
Theorem 8. Consider as the optimal value for the control parameter U; hence, can be expressed as the equations: Proof of Theorem 8. To derive the characteristic equation of optimal control
, proceed by solving the equations:
The optimal control solution
can be expressed as
given that the control variable
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:
and the initial populations in each compartment are set as
.
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
, the disease-free equilibrium point
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
, which corresponds to
, the endemic equilibrium point
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
, we can see that a higher
q leads to a lower value of
, 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, , and . In this case, the basic reproduction number of the model without vaccination is . We observe that , thus . This correlation highlights the influence that vaccination has on diminishing the basic reproduction number, making the disease more controllable.
It should be noted that
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., ).
Strategy B: Using vaccination and treatment control (i.e., and ).
Strategy C: Using isolation and treatment control control (i.e., and ).
Strategy D: Using isolation and vaccination control (i.e., and ).
Figure 8 illustrates the variation in the objective function
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
varying with time under different control strategies. Moreover, as seen in
Figure 9a, for control strategy A, control measure
should be implemented on a large scale in the initial days and then stopped after 13 days. Control measure
needs to maintain the control intensity at around 0.213 after the 22nd day, while control measure
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
in strategy B lasts for approximately 13 days. As seen in
Figure 9c, in comparison to strategy A, the control measures
and
in strategy C have both increased in terms of time and control intensity.
Figure 9d indicates that the control measures
and
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
surpasses that of
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
of strategy C maintains the highest control intensity,
Figure 10a reveals that strategy C has the largest susceptible population. Similarly, the control measure
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 and the endemic equilibrium , and calculate the basic reproduction number . Furthermore, we investigate the necessary conditions for ensuring the local and global stability of both and . We demonstrate the critical role of vaccination in disease elimination by investigating its influence on . Additionally, we develop an optimal control strategy to reduce the spread of the disease. In the model, we introduce three control variables, , , and , 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.