Graphical abstract Highlights • A novel mathematical model of Corona virus with new hybrid fractional operator derivative are presented. • Three control variables are presented to minimize the number of infected population. • Necessary control conditions are derived. • Two numerical methods are constructed to study the behavior of the obtained fractional optimality system. • The stability of the proposed methods are proved. • Numerical simulations and comparative studies are given. Introduction Coronavirus COVID-19 pandemic is the defining global health crisis of our time and the greatest challenge we have faced since world war two. To describe this disease mathematically, we noted that COVID-19, due to uncertainties associated to the pandemic, ordinal derivatives and their associated integral operators show deficient. The fractional order differential equations models seem more consistent with this disease than the integer order models. This is due to the fact that fractional derivatives and integrals enable the description of the memory and hereditary properties inherent in various materials and processes. Hence there is a growing need to study and use the fractional order differential equations. Also, optimal control theory is very important topic to control the variables in mathematical models of infectious disease. Moreover, a hybrid fractional operator which may be expressed as a linear combination of the Caputo fractional derivative and the Riemann–Liouville fractional integral is recently introduced. This new operator is more general than the operator of Caputo’s fractional derivative. Numerical techniques are very important tool in this area of research because most fractional order problems do not have exact analytic solutions. Objectives A novel fractional order Coronavirus (2019-nCov) mathematical model with modified parameters will be presented. Optimal control of the suggested model is the main objective of this work. Three control variables are presented in this model to minimize the number of infected populations. Necessary control conditions will be derived. Methods The numerical methods used to study the fractional optimality system are the weighted average nonstandard finite difference method and the Grünwald-Letnikov nonstandard finite difference method. Results The proposed model with a new fractional operator is presented. We have successfully applied a kind of Pontryagin’s maximum principle and were able to reduce the number of infected people using the proposed numerical methods. The weighted average nonstandard finite difference method with the new operator derivative has the best results than Grünwald-Letnikov nonstandard finite difference method with the same operator. Moreover, the proposed methods with the new operator have the best results than the proposed methods with Caputo operator. Conclusions The combination of fractional order derivative and optimal control in the Coronavirus (2019-nCov) mathematical model improves the dynamics of the model. The new operator is more general and suitable to study the optimal control of the proposed model than the Caputo operator and could be more useful for the researchers and scientists.
【저자키워드】 coronavirus diseases, Fractional order optimal control problems, A proportional derivative, Weighted average nonstandard finite difference method, Grünwald-Letnikov nonstandard finite difference method,