key: cord-0855868-76swbf7v
authors: Veisi, Amir; Delavari, Hadi
title: Fractional‐order backstepping strategy for fractional‐order model of COVID‐19 outbreak
date: 2021-11-28
journal: Math Methods Appl Sci
DOI: 10.1002/mma.7994
sha: 46a1399dcf46adf6ab0403c860fb48750650cdd1
doc_id: 855868
cord_uid: 76swbf7v

The coronavirus disease (COVID‐19) pandemic has impacted many nations around the world. Recently, new variant of this virus has been identified that have a much higher rate of transmission. Although vaccine production and distribution are currently underway, non‐pharmacological interventions are still being implemented as an important and fundamental strategy to control the spread of the virus in countries around the world. To realize and forecast the transmission dynamics of this disease, mathematical models can be very effective. Various mathematical modeling methods have been proposed to investigate the transmission patterns of this new infection. In this paper, we utilized the fractional‐order dynamics of COVID‐19. The goal is to control the prevalence of the disease using non‐pharmacological interventions. In this paper, a novel fractional‐order backstepping sliding mode control (FOBSMC) is proposed for non‐pharmacological decisions. Recently, new variant of this virus have been identified that have a much higher rate of transmission, so finally the effectiveness of the proposed controller in the presence of new variant of COVID‐19 is investigated.

In late 2019, in Wuhan, China, the latest human viral pathogen, extreme acute respiratory syndrome coronavirus-2, the source of the coronavirus disease 2019 (COVID-19) pandemic, emerged. The corona virus is spreading around the world and has brought many problems for different countries. The corona virus (COVID 19) has infected and killed many people globally, and counting. 1 The virus enters humans by crossing the species barrier, continuing the infection through human-to-human transmission. 2 There are now ubiquitous mathematical models of infectious disease transmission dynamics. These models play a significant role in helping to measure possible control and prevention methods for infectious diseases. [3] [4] [5] The typical classical differential equations, such as the well-known and widely utilized mathematical models for the spread of infections, are the SI, SIS, SIR, SEIR, SIRD, and SEIRD models. Each variable in these models denotes the number of people in different groups. Since the discovery of COVID 2019, numerous models have been suggested to study its dynamics. [6] [7] [8] [9] [10] [11] According to the first information reported in China, Zhong et al 12 have developed a simple SIR model for predicting the corona virus. Given the environmental impact, Yang and Wang 13 proposed a broad SEIR model for COVID-19 with variable velocity. Liang 14 described the spread of the three epidemics, COVID-19, SARS, and MERS, utilizing mathematical models and found that the progress rate of COVID 19 was much higher than that of SARS and MERS. Fractional-order differential equations have recently been utilized to describe the behavior of epidemics. [15] [16] [17] [18] [19] [20] [21] Fractional derivatives, besides to the current state, also depend on historical states and therefore have memory properties. 16, 17 The theory of control has recently emerged as one of the suggested approaches for non-pharmacological interventions. 22, 23 Thus, in recent years, different kinds control techniques have emerged in the literature, such as government feedback control, 23 feedback linearization, 22, 24 and observer-based control. 25 There are many techniques for designing robust controllers, such as H ∞ , fuzzy, or neural network-based techniques. [26] [27] [28] In Gondim and Machado, 29 the structure of the period and the quarantine rate (percentage of quarantine practices) are examined. The sliding mode controller (SMC) is also intended to design vaccination strategies 30 and regulate infected individuals. 31, 32 The study of optimal control over the outbreak of the virus in Gondim and Machado 29 and Morato et al 33 has been studied. An optimal control method with the aim of reducing the socioeconomic cost function in quarantine decisions can be very efficient. 34 One of the major problems with epidemics is the over usage of hospital facilities. In Berger, 35 the bang-bang controller is used to manage hospital beds by knowing information such as the number of infected people. In the present paper, the fractional-order model is used. In Oud et al, 36 a fractional epidemic model in the Caputo sense with the consideration of quarantine, isolation, and environmental impacts to examine the dynamics of the COVID-19 outbreak is proposed. In Chu et al, 37 a fractional-order transmission model is considered to study its dynamical behavior using the real cases reported in Saudi Arabia. In Khan and Atangana, 38 the mathematical modeling and dynamics of a novel corona virus (2019-nCoV) is studied. The brief details of interaction among the bats and unknown hosts, then among the peoples and the infections reservoir (seafood market), are investigated. In khan et al, 39 a new mathematical model is developed to understand COVID-19 dynamics and possible control. The model first is formulated in the integer order, and then, the Atangana-Baleanu derivative concept with a non-singular kernel is used for its generalization.

The parameters of the pandemic model are always associated with unwanted uncertainties because the parameters of disease spread and transmission rates, mortality, and contact rates between different people in the community are not constant and may change depending on social conditions and disease variations. Also, new variants of COVID-19 have different parameters and behavior. So a robust control method is necessary for control of COVID-19 outbreak.

The combination of backstepping method and sliding mode control (BSMC) can effectively deal with the problems caused by parameters uncertainties in the model of COVID-19 disease and can improve the robustness of the control system. Also using the long-term memory properties of fractional-order operator in control signal design can enhance the robustness against model uncertainty. So a novel fractional-order backstepping sliding mode control (FOBSMC) is proposed to deal with parameter uncertainty in model of COVID-19 disease. Also, the fractional-order model of COVID-19 is used to verify the effectiveness of the proposed method. Some contributions and advantages of this paper are as follows:

• The new nonlinear fractional-order SEIARM model for COVID19 is utilized to design the proposed controller.

• A novel FOBSMC is proposed.

• The performance of the proposed FOBSMC and the conventional BSMC are investigated against the new variants of COVID-19 (parameter uncertainty). • The simulation results of the proposed method guarantee that the beds of hospital intensive care unit will never reach saturation level, even if the disease transmission rate increases by 10%. • The proposed controller can be designed for each country according to the population and hospital facilities of that country. • Stability of the proposed method is theoretically proved.

Fractional calculus is the generalized form of the integer calculation. In this approach, the fractional-order integration and differentiation are defined by the fundamental operator of a D α t α ∈ R ð ), in which a and t are the bounds of the operation. The continuous fractional-order integral and differential operators are defined as follows 40 :

Basically, three definitions are presented for fractional calculations as follows: Grünwald-Letnikov (GL), Riemann-Liouville (RL), and Caputo. GL derivative definition of α order is defined as follows:

where h is time increment, t À a ð Þ=h ½ is the integer part and

where Γ is the gamma function. RL derivative definition of α order is defined as follows:

where n is the first integer number bigger than α. Caputo derivative definition of α order is defined as follows:

To define the Laplace operator S in the fractional calculus, the integer-order transfer function is first approximated. Based on the Oustaloup approximation, the following transfer function is defined in the frequency domain 40 :

Then, Equation (6) can be approximated by Equation (7).

where K is the obtained gain from two sides of Equation (7) with a unit (1 rad/s), and the number of poles and zeros is 2 N + 1, and ω z,n and ω p,n are defined as follows.

where ω b and ω h are the estimation bounds of the upper and lower frequency, respectively. Generally, the following relation is satisfied:

The definition of fractional integral is as follows 40 :

Theorem 1. If the function x t ð Þ∈R is a continuous and real derivative function and μ ¼ m n ≥ 1, which is m > 0 odd numbers and n∈N þ , then Equation (11) for t ≥ a is established 41 

Proof: Details of proof of Theorem 1 are provided in Dai and Chen. 41 Theorem 2. Assume a nonlinear fractional-order system C a D α

as a derivative and continuous function that is relative to the variable x Lipschitz and holds in conditions (12) and (13):

If α ∈ 0, 1 ð Þ, x ∈  , t ≥ 0 and the constants α 1 , α 2 , α 3 , a, b are selected positively and arbitrarily, then x ¼ 0 is a Mittag-Leffler stable and if the relations (12) and (13) generally held in the spaceR n , the system is a global Mittag-Leffler stable.

Proof: Details of proof for Theorem 2 are provided in reference Li et al. 42 

In this section, dynamic equations using the Caputo fractional-order derivative will be briefly defined as Equation (14) .

In this fractional-order model, the total human population N(t) can be divided into five categories: Susceptible S(t), Exposed E(t), Infected (symptomatic) I(t), Asymptomatic infected A(t), and the recovered people R(t) and M(t) represents COVID-19 in a reservoir or place or seafood market. The fractional-order equations of COVID-19 disease are defined as follows 43 :

where 0 < α < 1 and u is the control signal. The system model parameters are given in Table 1 . The basic conditions are relevant

The parameters are defined for simplicity as follows:

The values in Table 1 are for the population of China and N(0) = 1.4 billion.

The disease has a high rate of transmission. If in a specific and short time, the number of infected people in the community exceeds the maximum health facilities in a community; the community will suffer adverse consequences such as saturation of hospital beds, reduced efficiency of medical staff, and public dissatisfaction with the lack of health facilities. As a result, the number of deaths may increase dramatically. Therefore, in different countries, different strategies have been implemented to reduce the rate of disease transmission so that medical facilities are not saturated for a certain period of time. The following constraints can be set accordingly:

where Imax ¼ 6:87e6 indicates the maximum level of health facilities, which determines the number of beds in the hospital intensive care unit. This amount has been announced by the health authorities of each country. The goal is to control the prevalence of the disease in order to not saturate the number of hospital intensive care beds. Therefore, the goal is to control the number of patients. In fact, the number of patients should not exceed the I max . The basic and important parameter for deciding on non-pharmacological interventions is u(t) in (14) . u(t) is based on percentage, which can be called the percentage of quarantine and social constraints. However, despite the challenges of quarantine in most countries, it is important to consider how society will be affected by any quarantine rate. In fact, the fixed quarantine rate has its limitations, because if the quarantine rate is high, it will be difficult to implement due to the economic costs incurred during the quarantine, both by the community and the authorities. Also, high-rate quarantine works only for a limited time, because the country is practically shut down after several months of quarantine, and irreparable damage will be done to society. On the other hand, quarantine at a low rate makes the disease prevalence and has virtually no effect on reducing the prevalence of the disease. Also, the parameters of the pandemic model are always associated with unwanted uncertainties, because the parameters of disease spread and transmission rates, mortality, and contact rates between different people in the community are not constant and may change depending on social conditions and disease variations. Therefore, the quarantine rate can be considered variable over time to decrease or increase if needed. The parameters of the pandemic model are always associated with unwanted uncertainties because the parameters of disease spread and transmission rates, mortality, and contact rates between different people in the community are not constant and may change depending on social conditions and disease variations. Also, new variants of COVID-19 have different parameters and behavior. So a robust control method is necessary for control of COVID-19 outbreak. Hence, in the next section, a novel fractional-order robust method is proposed.

Designing a robust controller can be very effective in reducing the prevalence of the disease. The controller can overcome the condition in Equation (17) . In this section, a novel FOBSMC will be suggested to control the prevalence of the disease. Based on the above discussion, the error can be defined as Equation (18) 

Then, the fractional derivative of error with respect to time can be given by (19) c

where c 0 D α t I max will go to zero. Define the stabilizing function

where δ > 0 is a stabilizing coefficient. An auxiliary error variable will be defined as Equation (21) 

The first Lyapunov function is defined as follows:

By taking c 0 D α t : ð Þ form both side of Equation (22), we can obtain

A fractional sliding surface is proposed as Equation (24) σ

where, C 1 , C 2 , C 3 , C 4 > 0 are the controller parameters and have constant value. By taking fractional derivative of Equation (24), we have

A second augmented Lyapunov function is defined as Equation (26) 

According to Theorem 1 and by taking fractional derivative of Equation (26) and utilizing Equation (23), we have

By putting Equation (21) in Equation (28), one can obtain

By substituting Equation (14) in Equation (29), the control law u t ð Þ ð

Þis obtained by the proposed FOBSMC as

where C 4 > 0 is the controller parameter and has constant value, and K 1 , K 2 are defined in Equation (16), and the other parameters and its value are described in Table 1 . The proposed control signal (30) is applied to the fractionalorder model of COVID-19 disease in Equation (14) to obtain the simulation results for the proposed FOBSMC method.

Consider the second augmented Lyapunov function in Equation (26) and according to Equations (26) (27) (28) (29) , and by using the fractional-order model of COVID in Equation (14), one can obtain

By applying the control signal, we have Equation (32) c

simplifying Equation (32) and using Theorem 2 results:

The stability of the closed-loop system is guaranteed by Equation (33).

The results in this article are based on statistics from China, and the controller parameters are δ ¼ 1:1, C 1 ¼ 10, C 2 ¼ 8, 

The percentage of quarantine and social constraints (non-pharmacological interventions, i.e., u(t)) can directly effect on the disease prevalence. However, despite the challenges of quarantine in most countries, it is important to consider how society will be affected by any quarantine rate. In fact, the fixed quarantine rate has its limitations, because if the quarantine rate is high (close to 100%), it will be difficult to implement due to the economic costs incurred during the quarantine, both by the community and the authorities. Also, high-rate quarantine works only for a limited time, because the country is practically shut down after several months of quarantine and irreparable damage will be done to society. On the other hand, quarantine at a low rate (close to 0%) makes the disease prevalence and has virtually no effect on reducing the prevalence of the disease. In this section, we examined the behavior of the outbreak model with open-loop controller. In this case, the control signal is assumed to be u(t) = [0 0.2 0.4 0. 5 1] . The behavior of fractional model dynamic is examined and compared in terms of 0%, 20%, 40%, 50%, and 100% quarantine (non-pharmacological interventions, i.e., u(t)).

The performance of susceptible and exposed behavior is shown in Figures 1 and 2 . As shown by increasing the amount of input signal (increasing the quarantine limit), the number of susceptible individuals increases. It can be concluded that fewer people have been infected with this disease. In the beginning, when the disease is spreading, everyone in a society is susceptible. As quarantine increases in a community, communication between individuals decreases, which naturally reduces disease transmission, and it can be concluded that the number of susceptible increases with increasing quarantine because fewer people became infected. As restrictions increase, the number of exposed people also decreases ( Figure 2) . Figures 3 and 4 show symptomatic and asymptomatic patients. Reducing the rate of disease transmission through place or seafood markets is shown in Figure 5 . Naturally, increasing the restrictions will reduce the number of symptomatic and asymptomatic patients, but the main and important point is that increasing the restrictions can cause adverse consequences such as social and economic crises and so on.

In open-loop systems, no information from the output will affect the amount of input, and in any case, without measuring the output of the model, decisions are made to adopt quarantine, and the amount of input is constant, that does not lead to good performance. The percentage of 0% to 100% quarantines by trial and error has been implemented on the system model, which according to the results, 50% was the best results for open-loop controllers.

In the closed-loop system according to the amount of I max , the proposed controller produces the proper control law. The performance of the proposed method (FOBSMC) is shown in Figures 6-10 . To obtain Figures 6 and 7, our proposed control signal in Equation (30) is applied to the fractional-order model of COVID-19 disease in Equation (14) . In this paper, the proposed FOBSMC method is compared with the conventional BSMC method. The proposed method has more precise and better performance as shown in Figures 6-10 . This method has been able to reduce the level of exposed people and patients in the community. The susceptible people are shown in Figure 6 . At the beginning of the outbreak, all people in a society are considered susceptible. Reducing the susceptible population to zero means removing people from the category of susceptible and being divided into exposed-infected groups with symptomsasymptomatic infections, and so forth. The proposed method (FOBSMC) puts a higher percentage of people in the category of susceptible, which means less infect of people in the community ( Figure 6 ). Fractional-order value α ¼ 0:8 is utilized.

The purpose of this article is to manage hospital intensive care beds. The amount of hospital beds is announced by governments (I max ). The proposed controller defines the best control law according to this I max . The exposed people are shown in Figure 7 . The increase and decrease of exposed people are directly related to the number of infected people. Reducing the number of exposed people in society has been the most important non-pharmacy intervention strategy of countries. As shown in Figure 7 , the proposed method has the lower level of exposed individuals, which in turn reduces the number of infected individuals. Also it can be seen from Figure 7 that the slope and the trend of increasing the number of exposed people is small, which makes it possible to use the hospital facilities without its saturated or overuse. The fewer exposed people naturally means fewer people are affected. However, the BSMC method has reached the exposed population more rapidly and to a higher threshold, which can disrupt the management of the intensive care unit.

In Figure 8 , the proposed FOBSMC has a very good performance and has reduced the number of infected people in the community much more than BSMC. This gentle increase process allows the health facilities to not reach their saturation level. On the other hand, gentle growth prevents the accumulation of many patients in a certain period of time; hence, medical staff will not lose their energy for service work and can be most efficient. It can be concluded that the control of the prevalence of the disease and the management of the intensive care unit has been done better with the proposed method. The number of asymptomatic infected individuals is shown in Figure 9 . A decrease in the number of asymptomatic infected people indicates a decrease in exposed people. Naturally, if the number of exposed people decreases, people in the community will have less contact with each other and the rate of disease transmission between people will decrease. As shown in Figure 9 , the proposed method has a lower level of infected (asymptomatic) people that can greatly help control the disease outbreak. Figure 10 shows the rate of disease transmission through contaminated seafood. Food contamination in the community will decrease when infected people decrease. 

Recently, new variant of this virus has been identified that have a much higher rate of transmission. This section contains simulations which illustrates the effectiveness of the proposed controller in the presence of new variants of COVID-19 as a parameters uncertainty in model. The actual parameters of the model are given by contact rate U ¼ 0:134, disease transmission coefficient U w ¼ 0:0000000080082, while the nominal parameters are given by contact rate U ¼ 0:1474, disease transmission coefficient U w ¼ 8:8090e À 09.

The nominal parameters do not match with the actual ones. Moreover, there is no information on their possible bounds. Though, regardless of these uncertainties, the designed controller would be able to achieve the control goal.

The performance of BSMC and FOBSMC controllers with parameters uncertainty is compared and shown in Figures 16-18 . Figure 16 shows the performance of the two controllers against parameters uncertainty for the exposed population. The proposed controller has shown good behavior and has made small changes in increasing the number of exposed people, which shows that the proposed FOBSMC has more robust performance against parameters uncertainty. On the other hand, the BSMC_Un (BSMC in the presence of uncertainty) has increased the number of exposed people, which can increase the disease in the community.

As shown in Figure 17 , the proposed controller has the lowest number of patients. On the other hand, the rate of disease increase was slower than the BSMC method. As we know, if the number of patients people increases significantly over a period of time, the health care system will face challenges such as saturating of hospital beds, reducing the efficiency of medical staff, and so on. According to Figure 17 , the proposed controller has shown a satisfactory behavior. A 10% increase in the rate of disease transmission and contact of people in a community can pose a major challenge in the management of the health care system, but nevertheless, the proposed controller has shown very good behavior and has been robust to parameter uncertainty. The proposed controller behavior in the presence of uncertainty ensures that the controller is robust and prevents the saturation of hospital beds.

One of the biggest challenges of the corona virus is its incubation period. The incubation period means that the sufferer is unaware that he or she has the disease. From the time the disease enters the body until the person has symptoms, it takes a period of time. This group of infected people who do not know about their disease can be included in the category of asymptomatic infected people. The implementation of social restrictions such as quarantine and school closures should be decided in such a way as to reduce the contact of asymptomatic infected people in the community so that the disease is not transmitted from person to person and this cycle of human transmission is stopped. As shown in Figure 18 , the performance of the proposed controller was very good in the presence of a 10% increase in disease transmission and human contact, and its changes compared to BSMC-Un were very small, which according to the results we can ensure the controller is robust. To show the superiority of the proposed method, the amount of area below the plot diagram is shown in the Table 2 .

Results show the proposed method has a smaller numerical value than BSMC, meaning that fewer people get infected with FOBSMC method.

The parameters of the pandemic model are always associated with unwanted uncertainties. Also, new variants of COVID-19 have different parameters and behavior. So a novel robust control method is proposed by combination of fractional calculus, backstepping and sliding mode control (FOBSMC). The proposed controller is applied to a fractional-order model of COVID- 19 . At first open-loop controller for the percentage of 0% to 100% quarantines by trial and error has been implemented on the system model, which according to the results, 50% was the best results for openloop controller. But the fixed quarantine rate has its limitations and difficulty for implementation; also, the model parameters are always associated with unwanted uncertainties. So the robust proposed method (FOBSMC) is applied. The area under the curve (infected I(t)) Value Simulation results show that the rate of disease increase was slower by proposed FOSMC than the BSMC method; the proposed controller has the lowest number of infected people, and also, it has the lower level of exposed individuals. The fewer exposed people naturally means fewer people are affected. Also, the slope and the trend of increasing the number of exposed people with FOSMC is small, which makes it possible to use the hospital facilities without its saturated or overuse. Finally, the performance of the proposed controller in the presence of new variants of COVID-19 is compared with BSMC and open-loop controller. Simulation results ensure the proposed FOBSMC has more robust performance against parameters uncertainty and prevents the saturation of hospital intensive care beds.

There are no funders to report for this submission.

This work does not have any conflicts of interest.

AUTHOR CONTRIBUTIONS Amir Veisi: Conceptualization; data curation; formal analysis; funding acquisition; investigation; methodology; resources; software; validation; visualization. Hadi Delavari: Conceptualization; data curation; formal analysis; funding acquisition; investigation; methodology; project administration; resources; software; supervision; validation; visualization.

Hadi Delavari https://orcid.org/0000-0002-2538-8801

World Health Organization. WHO coronavirus (COVID-19) dashboard

A hybrid multi-scale model of COVID-19 transmission dynamics to assess the potential of non-pharmaceutical interventions

A stochastic SICA epidemic model for HIV transmission

Mathematical modeling of Zika disease in pregnant women and newborns with microcephaly in Brazil

Dynamics and optimal control of Ebola transmission

Early dynamics of transmission and control of COVID-19: a mathematical modelling study

Real-time forecasts of the COVID-19 epidemic in China from

The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application

Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions

Forecasting covid-19

Analysis and forecast of COVID-19 spreading in China, Italy and France

Early prediction of the 2019 novel coronavirus outbreak in the mainland China based on simple mathematical model

A mathematical model for the novel coronavirus epidemic in Wuhan, China

Mathematical model of infection kinetics and its analysis for COVID-19, SARS and MERS

A fractional order SEIR model with density dependent death rate

On a fractional order Ebola epidemic model

A fractional order epidemic model for the simulation of outbreaks of influenza A (H1N1)

Analysis of a fractional SEIR model with treatment

An epidemiological MSEIR model described by the Caputo fractional derivative

A fractional order SEIR model with vertical transmission

Fractional order model for the spread of leptospirosis

Analysis and control of an SEIR epidemic system with nonlinear transmission rate

On vaccination controls for the SEIR epidemic model

Feedback linearization-based vaccination control strategies for true-mass action type SEIR epidemic models

An observer-based vaccination control law for an SEIR epidemic model based on feedback linearization techniques for nonlinear systems

HDo robust controller design of media advance systems with time domain specifications

Decentralized robust control for uncertain TS fuzzy large-scale systems with time-delay

Adaptive neural control of pure-feedback nonlinear time-delay systems via dynamic surface technique

Optimal quarantine strategies for the COVID-19 pandemic in a population with a discrete age structure

Robust sliding control of SEIR epidemic models

Sliding mode control of outbreaks of emerging infectious diseases

Analysis and control of seir epedemic model via sliding mode control

An optimal predictive control strategy for COVID-19 (SARS-CoV-2) social distancing policies in Brazil

Modeling, state estimation, and optimal control for the US COVID-19 outbreak

Feedback control of the COVID-19 pandemic with guaranteed non-exceeding ICU capacity

A fractional order mathematical model for COVID-19 dynamics with quarantine, isolation, and environmental viral load

Dynamics of fractional order COVID-19 model with a case study of Saudi Arabia

Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative

A robust study on 2019-nCOV outbreaks through non-singular derivative

Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and Some of their Applications

New power law inequalities for fractional derivative and stability analysis of fractional order systems

Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability

Modeling and simulation of the novel coronavirus in Caputo derivative

How to cite this article: Veisi A, Delavari H. Fractional-order backstepping strategy for fractional-order model of COVID-19 outbreak