Complex Dynamics and Fractional-Order Optimal Control of an Epidemic Model with Saturated Treatment and Incidence

In this study, we have developed a novel SIR epidemic model by incorporating fractional-order differential equations and utilizing saturated-type functions to describe both disease incidence and treatment. The intricate dynamical characteristics of the proposed model, encompassing the determination of the conditions for the existence of all possible feasible equilibria with their local and global stability criteria, are investigated thoroughly. The model undergoes backward bifurcation with respect to the parameter representing the side effects due to treatment. This phenomenon emphasizes the critical role of treatment control parameters in shaping epidemic outcomes. In addition, to understand the optimal role of the treatment in mitigating the disease prevalence and minimizing the associated cost, we investigated a fractional-order optimal control problem. To further visualize the analytical results, we have conducted simulation works considering feasible parameter values for the model. Finally, we have employed local and global sensitivity analysis techniques to identify the factors that have the greatest potential to reduce the impact of the disease.