Literature DB >> 34925704

A revisit to the past plague epidemic (India) versus the present COVID-19 pandemic: fractional-order chaotic models and fuzzy logic control.

Manashita Borah1, Binoy Krishna Roy2, Tomasz Kapitaniak3, Karthikeyan Rajagopal4, Christos Volos5.   

Abstract

India is one of the worst hit regions by the second wave of COVID-19 pandemic and 'Black fungus' epidemic. This paper revisits the Bombay Plague epidemic of India and presents six fractional-order models (FOMs) of the epidemic based on observational data. The models reveal chaotic dispersion and interactive coupling between multiple species of rodents. Suitable controllers based on fuzzy logic concept are designed to stabilise chaos to an infection-free equilibrium as well as to synchronise a chaotic trajectory with a regular non-chaotic one so that the unpredictability dies out. An FOM of COVID-19 is also proposed that displays chaotic propagation similar to the plague models. The index of memory and heredity that characterise FOMs are found to be crucial parameters in understanding the progression of the epidemics, capture the behaviour of transmission more accurately and reveal enriched complex dynamics of periodic to chaotic evolution, which otherwise remain unobserved in the integral models. The theoretical analyses successfully validated by numerical simulations signify that the results of the past Plague epidemic can be a pathway to identify infected regions with the closest scenarios for the present second wave of Covid-19, forecast the course of the outbreak, and adopt necessary control measures to eliminate chaotic transmission of the pandemic.
© The Author(s), under exclusive licence to EDP Sciences, Springer-Verlag GmbH Germany, part of Springer Nature 2021.

Entities:  

Year:  2021        PMID： 34925704      PMCID： PMC8665904          DOI： 10.1140/epjs/s11734-021-00335-2

Source DB:  PubMed          Journal:  Eur Phys J Spec Top        ISSN： 1951-6355            Impact factor:   2.891


Introduction

Under the clutches of the second wave of COVID-19 pandemic, India is struggling to deal with the immediate crisis and emergency situation. The number of COVID-19-infected cases which was 18000 per day on 5th March, 2021, accelerated exponentially to 0.4 million in 2 months by 6th May in India, recording the highest ever total for a single day [1]. The unprecedented steep rate of propagation of the pandemic led to the over-burden of health infrastructure, saturation of medical facilities, shortage of vaccines, oxygen supply, hospital beds, etc., so much so that the number of deaths per day rose to a staggering 4454 by the end of May 2021. To worsen complications is the rising insurgent threat of another growing epidemic, ‘Black fungus’ or Mucormycosis, reported to have 9000 cases in India by May, 2021, particularly infecting post COVID victims with comorbidities such as Diabetes with a mortality rate of 50% [2]. Increase in testing measures and easing lockdown restrictions until March, 2021, were most likely the causes of the immediately following spikes in number of cases in the second wave of COVID in India by the month of April, 2021. These alarming figures signify the urgency of an efficient mathematical analysis and modelling of COVID-19 propagation so that its course of outbreak can be foreseen and controlled through a solution driven approach. This paper hereby proposes to study its complex dynamical evolution by revisiting a past epidemic that was the cause of numerous historically prominent pandemics: Plague. The outbreak of Bombay (now Mumbai) Plague epidemic (1896–1911), India, was first reported in the British medical Journal, 1896 [3]. The surge in the casualties and infections dropped in 1897, but the epidemic turned into an endemic and lasted for more than a decade, returning every year with a seasonal dissemination. This has a close similarity with the COVID-19 pandemic, the first wave of which hit India in January 2020, and returned with a more severe second wave in early 2021, though it is too early to say whether it is seasonal. The bubonic plague, called so because of the development of buboes or sore lymph nodes in the infected, did not stay confined to Bombay alone, it hit all over India [4]. It is caused by bacillus Yersinia pestis and develops symptoms such as fever, vomiting, bleeding, organ failure, headache, weakness, etc. which if untreated may prove fatal in the first 36 hours of being infected [5]. Previous mathematical models of the plague epidemic are demonstrated by compartmental models based on the three compartments: susceptible, infected, recovered (SIR) [6] and fail to consider the transition of the disease between rodents and fleas. Plague is generally transmitted to humans by infected rodents who had been bitten by a flea. In fact, numerous species of rodents such as black rats, marmots, prairie dogs, gerbils, and squirrels are known carriers of plague [7]. Mathematical models of the epidemic based on nonlinear dynamics to study the interaction between a single species of rat and flea is proposed in 2012 by Bacaër et al. [8]. In 2015, a new model of plague was proposed in [5] that took into account multiple epizootics of rats such as Mus decumanus and Mus rattus. It led to the finding that the progression of the plague epidemic model displays chaotic behaviour. But all these models based on chaos theory are confined only to integer-order models (IOMs) which often fail to encompass in-depth understanding of the dynamics of the epidemic such as involvement of asymptomatic carriers, irregular transmission, critical cases with comorbidities, complexity of interaction amongst species, etc. Fractional order models (FOMs), defined by fractional calculus, have an additional parameter of the non-integral derivative with an index of memory that incorporates various unmodelled dynamics of real physical systems resulting in a closer and more accurate interpretation of the model dynamics [9]. FOMs extract the underlying randomness in various engineering applications as switching circuits [10], jerk oscillators [11], hyperchaotic systems [12], homopolar disc dynamo [13], secure communication [14], etc. FOSs are known to produce various enriched complexities as period doubling route to chaos [15], oscillator dynamics [16], multistability [17], memristive dynamics [18], etc. In epidemiology, the memory and probability of transmission of FO chaotic models are found to be crucial parameters to capture the real dynamics of epidemic and infectious diseases [19]. For example, the FO-SIR chaotic compartmental model of Dengue epidemic is found to generate a closer fit to the observational data as validated from the Dengue cases of Malaysia [20]. The FO chaotic epidemic model of HIV comprising AIDS related cancer cells suggests that FOMs are better predictors and significantly impact the dynamic progress [21]. Also, nonlinear dynamic approach has been used to analyse cardiovascular diseases [22], a study especially relevant in the present COVID situation where heart patients are at greater risks of severe infections. Our work therefore, applies the superior benefits of fractional calculus and nonlinear dynamics approach to investigate the presence of chaos in the dispersal of plague epidemic and help in better understanding of the physiological functioning of infectious diseases. Chaos, due to its high sensitivity to initial conditions, aperiodicity, erratic permanence and unpredictability, is an undesirable phenomenon in epidemiology [23]. This is why control of chaotic motion to stabilise it to a disease-free equilibrium or synchronise it with regular motion is a relevant area of research. Fuzzy logic control (FLC) has been a very practical and effective technique of controlling and synchronising chaos in nonlinear systems [24]. They are advantageous over conventional controllers since they are flexible over a wide range of operating conditions due to their soft logic principle and can be applicable to real-world systems with uncertainties and external disturbances, thus considerably reducing the tracking error and control effort [25]. The authors use the Takagi–Sugeno Fuzzy logic control (TS-FLC) [26] to suppress the undesirable effects of chaotic oscillations in the epidemic models. The remainder of the paper is as follows: Section 2 presents the proposed FOMs of plague epidemic and relevant FO fundamentals, Sect. 3 depicts the design of controllers for stabilisation and synchronisation of the proposed models, and Sect. 4 presents the results and discussion. The paper is finally concluded in Sect. 5.

Fundamentals of FO calculus and proposed FO plague models (FOPMs)

The modelling of the epidemics and design of controllers to suppress chaos are determined by FO stability theorems whose preliminaries are presented below.

Preliminaries of FO

The Caputo fractional derivative, of order of a continuous function f(t) is defined in (1).Let us define an FO nonlinear system (FONLS) as in (2),where the FOs lie in and

Theorem 1

[27] Let  be an equilibrium point for the FONLS (2), V(x) be a Lyapunov function, and  be functions of class K such that , where , then, the equilibrium point is asymptotically stable. For a quadratic positive definite function , where Q is a positive definite matrix, the condition (a) of Theorem 1 is always satisfied [28].

Lemma 1

If is a continuous and derivable function, then, for any time instant ,where is the Caputo fractional derivative of of FO . The Adams–Bashforth–Moulton (ABM) method established on the predictor–corrector technique [29] is used to solve the FODEs.

FO plague models (FOPMs)

Recently, IOMs of chaotic plague epidemic models are proposed in [5] derived from the observational data on the number of deaths due to plague infection from September 1896 to December 1911 and coupling between the epidemic and the epizootic of the two main species of rodents using a global modelling technique. Based on this, we explore the FO dynamics of the models and propose six corresponding FOPMs where, and represent the number of fatalities due to plague infection, infected rats of the M. decumanus species captured and infected rats of the M. rattus species captured, respectively. a) FOPM is a 10-term chaotic model as in (3),where tuning parameter . b) FOPM is a 11-term chaotic model given as in (4), where .c) FOPM is also a 11-term chaotic model given as in (5), where d) FOPM is a 12-term chaotic model given as in (6), where .e) FOPM is a 10-term periodic model given as in (7).This is the untuned model with, . f) FOPM is a 11-term periodic model given as in (8).This is the untuned model with . Table 1 describes the initial conditions (ICs), equilibrium points, eigen values of the proposed FOPMs and their Lyapunov exponents (LEs) calculated using the algorithm proposed by Danca et al. [30].
Table 1

Dynamical analyses of the FO plague epidemic models

ModelICsEquilibrium pointsEigen values \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\lambda _{i},\, i=1,2,3)$$\end{document}(λi,i=1,2,3)FOLEsDynamics
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{0}$$\end{document}P0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(166.8903, \\ 1070.193, \\ 138.7208)\end{array}$$\end{document}(166.8903,1070.193,138.7208)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}E_{1}=(-\, 2723.906,0, \\ -\, 459.002)\end{array}$$\end{document}E1=(-2723.906,0,-459.002)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (64.556, − 9.44,32.0185)0.98\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\mathrm {L}\mathrm {E}_{\mathrm {1}}=0.059548\\ \mathrm {L}\mathrm {E}_{2}=0\\ \mathrm {L}\mathrm {E}_{3}= -\, 3.0896\end{array}$$\end{document}LE1=0.059548LE2=0LE3=-3.0896Chaotic
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{2}=(0,0,0)$$\end{document}E2=(0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (− 7.5489, 0,0)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{3}=$$\end{document}E3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(428.4239,948.939, \\ 97.436)\end{array}$$\end{document}(428.4239,948.939,97.436)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 7.2477, \\ 2.8296 + 11.2178i,\\ 2.8296 - 11.2178i)\end{array}$$\end{document}(-7.2477,2.8296+11.2178i,2.8296-11.2178i)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{1}$$\end{document}P1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(895.9957, \\ 2391.688, \\ 312.7087)\end{array}$$\end{document}(895.9957,2391.688,312.7087)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{1}=(0,0,0)$$\end{document}E1=(0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (− 12.623, 0,0)0.979\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\mathrm {L}\mathrm {E}_{\mathrm {1}}=0.22158\\ \mathrm {L}\mathrm {E}_{2}=0\\ \mathrm {L}\mathrm {E}_{3}= -\, 0.45977\end{array}$$\end{document}LE1=0.22158LE2=0LE3=-0.45977Chaotic
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{2}=$$\end{document}E2= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 1112.367,0, \\ -\, 379.308)\end{array}$$\end{document}(-1112.367,0,-379.308)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (26.363, − 19.1324, 18.7156)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{3}=$$\end{document}E3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(533.355,1181.3568, \\ 200.333)\end{array}$$\end{document}(533.355,1181.3568,200.333)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 7.8815, \\ 1.4159 + 12.3523i,\\ 1.4159 -\, 12.3523i)\end{array}$$\end{document}(-7.8815,1.4159+12.3523i,1.4159-12.3523i)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{2}$$\end{document}P2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(836.9658, \\ 1721.620, \\ 328.9299)\end{array}$$\end{document}(836.9658,1721.620,328.9299)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{1}=(0,0,0)$$\end{document}E1=(0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (− 11.453,0,0)0.996\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\mathrm {L}\mathrm {E}_{\mathrm {1}}=0.32517\\ \mathrm {L}\mathrm {E}_{2}=0\\ \mathrm {L}\mathrm {E}_{3}= -\, 6.099\end{array}$$\end{document}LE1=0.32517LE2=0LE3=-6.099Chaotic
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{2}=$$\end{document}E2= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 2659.1148,0, \\ - 460.8589)\end{array}$$\end{document}(-2659.1148,0,-460.8589)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (63.0210, − 14.0915, 31.7721)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{3}=$$\end{document}E3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(512.1487,1134.3856, \\ 178.6772)\end{array}$$\end{document}(512.1487,1134.3856,178.6772)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 7.9728, \\ 2.4939+12.7008i, \\ 2.4939 - 12.7008i)\end{array}$$\end{document}(-7.9728,2.4939+12.7008i,2.4939-12.7008i)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{3}$$\end{document}P3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(344.0896, \\ 627.7882, \\ 101.5058)\end{array}$$\end{document}(344.0896,627.7882,101.5058)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{1}=$$\end{document}E1= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(810.5939,1786.4544, \\ 145.9542)\end{array}$$\end{document}(810.5939,1786.4544,145.9542)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 17.624, \\ 2.5635+19.2493i, \\ 2.5635 - 19.2493i)\end{array}$$\end{document}(-17.624,2.5635+19.2493i,2.5635-19.2493i)0.996\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\mathrm {L}\mathrm {E}_{1}=0.76966\\ \mathrm {L}\mathrm {E}_{2}=0\\ \mathrm {L}\mathrm {E}_{3}= - 15.0367\end{array}$$\end{document}LE1=0.76966LE2=0LE3=-15.0367Chaotic
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{2}=$$\end{document}E2= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(1125.8296,2481.1972, \\ 1006.1933)\end{array}$$\end{document}(1125.8296,2481.1972,1006.1933)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(- 6.9833 + 14.039i, \\ - 6.9833 - 14.039i, \\ 52.146)\end{array}$$\end{document}(-6.9833+14.039i,-6.9833-14.039i,52.146)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{3}=(0,0,0)$$\end{document}E3=(0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (- 11.403, 0, 0)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{4}$$\end{document}P4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(166.8903, \\ 1070.193, \\ 138.7208)\end{array}$$\end{document}(166.8903,1070.193,138.7208)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{1}=(0,0,0)$$\end{document}E1=(0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 12.623, \\ 0,0)\end{array}$$\end{document}(-12.623,0,0)0.938\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\mathrm {L}\mathrm {E}_{\mathrm {1}}=0\\ \mathrm {L}\mathrm {E}_{2}= -\, 0.33791\\ \mathrm {L}\mathrm {E}_{3}= - 6.4779\end{array}$$\end{document}LE1=0LE2=-0.33791LE3=-6.4779Periodic
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{2}=$$\end{document}E2= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 1911.3082,0, \\ -\, 497.203)\end{array}$$\end{document}(-1911.3082,0,-497.203)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(45.298, -\, 16.6145, \\ 21.3474)\end{array}$$\end{document}(45.298,-16.6145,21.3474)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{3}=$$\end{document}E3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(532.1154,1178.6108, \\ 200.9904)\end{array}$$\end{document}(532.1154,1178.6108,200.9904)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 6.1764, \\ 1.3416 + 12.8042i, \\ 1.3416 - 12.8042i)\end{array}$$\end{document}(-6.1764,1.3416+12.8042i,1.3416-12.8042i)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{5}$$\end{document}P5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(836.9658, \\ 1721.620, \\ 328.9299)\end{array}$$\end{document}(836.9658,1721.620,328.9299)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{1}=(0,0,0)$$\end{document}E1=(0,0,0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= (− 12.726, 0,0)0.962\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\mathrm {L}\mathrm {E}_{\mathrm {1}}=0\\ \mathrm {L}\mathrm {E}_{2}= -\, 0.059592\\ \mathrm {L}\mathrm {E}_{3}= -\, 6.4876\end{array}$$\end{document}LE1=0LE2=-0.059592LE3=-6.4876Periodic
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{2}=$$\end{document}E2= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(534.5436,1183.9892, \\ 203.8006)\end{array}$$\end{document}(534.5436,1183.9892,203.8006)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(-\, 7.4083, \\ 1.9474 + 13.0348i, \\ 1.9474 - 13.0348i)\end{array}$$\end{document}(-7.4083,1.9474+13.0348i,1.9474-13.0348i)
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{3}=$$\end{document}E3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(654.4082,1449.4837, \\ 486.6842)\end{array}$$\end{document}(654.4082,1449.4837,486.6842)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _{1,2,3}=$$\end{document}λ1,2,3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}(4.9346, \\ -\, 0.6295 + 19.7595i, \\ -\, 0.6295 - 19.7595i)\end{array}$$\end{document}(4.9346,-0.6295+19.7595i,-0.6295-19.7595i)
Dynamical analyses of the FO plague epidemic models

Design of controllers for FOPMs

This section is dedicated to design controllers for two types of control: (a) stabilisation of chaos and (b) synchronisation of chaos.

Design of controllers for stabilisation of chaos in FOPMs

Stabilisation controllers to suppress chaos are applied only to the chaotic FOPMs, i.e., , , and as discussed in Table 1. As represented in Fig. 1a, TS Fuzzy Logic Control (TSFLC) [26] is designed with a rule base in which the consequents of each rule are functions of inputs, the fuzzifier converts the crisp inputs to linguistic variables using the membership functions stored in the database, the inference engine enables decision making by adequate reasoning and the defuzzifier converts the fuzzy output values back into crisp values.
Fig. 1

Fuzzy Logic control system

Fuzzy Logic control system Let the controlled FONLS be defined as in (9).The control signal, , is calculated by weighted sum defuzzification method. The fuzzy control rule base for the TSFLC is defined as in (10),where are fuzzy sets describing the linguistic terms (LTs) of input variables, is the control input of rule i, function AND is a fuzzy operator and r is the total number of fuzzy rules. Each fuzzy rule defined in (10) generates a weight defined as in (11).From (10), it is assumed that for any in the input universe of discourse X, there exists at least one weight among all rules that is non-zero. Applying weighted sum defuzzification method, we obtain u as in (12).Since a positive definite quadratic function always satisfies Theorem 1, the Lyapunov candidate is chosen as in (13). is a positive definite matrix, which implies that V is positive definite, having continuous partial FO derivatives.

Theorem 2

[26] Let the origin be an equilibrium of the controlled epidemic (9) and there exists a Lyapunov function, , where Q is a positive definite matrix on domain X containing the origin of such that , for . Let . If except the trivial solution , no solution of (9) exists identically in S, then, the infection-free equilibrium at the origin is asymptotically stable in the domain X. Using Theorem 2 based on FO Lyapunov stability, the stability of an infection free equilibrium may be obtained. On adding TSFLC to the FOPM (4), we have (14),where are the TSFLCs. The fuzzification module of TSFLC in Fig. 1b) shows the triangular membership function of LTs of the chosen linguistic variable as antecedent. Notations P, Z and N for the LTs represent ‘positive’, ’zero’ and ‘negative’, respectively, and the parameters are chosen as Table 2 describes the set of fuzzy control rules as per the membership function.
Table 2

Rule base of TSFLCs for the chaotic FOPMs

RuleAntecedent \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_{1}$$\end{document}x1Consequent \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$u=\left[ {\begin{array}{*{20}c} u_{1}\\ u_{2}\\ u_{3}\\ \end{array} } \right] $$\end{document}u=u1u2u3
1P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}$$\end{document}U1
2N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}$$\end{document}U2
3Z\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}$$\end{document}U3
Rule base of TSFLCs for the chaotic FOPMs The universe of discourse is From (13), the Lyapunov function is chosen as (15).Using Lemma 1 and Caputo fractional derivative in (15), we have (16).From Theorem 1, we getUsing Theorem 2, which implies . Next, the fuzzy rules are derived as follows. so that , i.e. the stability requirement is satisfied. so that , i.e. the stability requirement is satisfied. so that , i.e. the stability requirement is satisfied. Fuzzy rule 1, : for antecedent as P, we have the consequent as Fuzzy rule 2, : for antecedent as N, we have the consequent as Fuzzy rule 3, : for antecedent as Z, we have the consequent as Since Theorem 2 is satisfied in the FO sense, thus the designed FLC converges the controlled epidemic to an infection-free asymptotically stable equilibrium. Similarly, the FLCs designed for the chaotic FOPMs ,, and fuzzy rules, , where are enlisted in Table 3.
Table 3

Proposed TSFLCs to stabilise chaos in FOPMs in antecedent

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{k}$$\end{document}Pk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {R}_{\mathbf {i}}$$\end{document}RiConsequent\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{k}$$\end{document}Pk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {R}_{\mathbf {i}}$$\end{document}RiConsequent
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{0}$$\end{document}P0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{1}$$\end{document}R1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} 0.1123\, x_{3}^{2}-0.0491x_{2}x_{3}+5.912x_{3}-\frac{x_{3}^{2}}{x_{1}} \\ -0.0107x_{2}^{2}+0.0237x_{1}x_{2}-x_{2}-\frac{0.0108x_{3}^{3}}{x_{2}}-1.4512x_{3} \\ 0 \\ \end{array}} \right] $$\end{document}U1=u1u2u3=0.1123x32-0.0491x2x3+5.912x3-x32x1-0.0107x22+0.0237x1x2-x2-0.0108x33x2-1.4512x30\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{2}$$\end{document}P2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{1}$$\end{document}R1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] $$\end{document}U1=u1u2u3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=\left[ {\begin{array}{c} 0.1581x_{3}^{2}-0.0688x_{2}x_{3}+5.912x_{3} \\ -0.0107x_{2}^{2}+0.0262x_{1}x_{2}-x_{2}-1.4512x_{3}-\frac{0.0108x_{3}^{3}}{x_{2}} \\ -0.0108x_{3}^{2}-x_{3} \\ \end{array}} \right] $$\end{document}=0.1581x32-0.0688x2x3+5.912x3-0.0107x22+0.0262x1x2-x2-1.4512x3-0.0108x33x2-0.0108x32-x3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{2}$$\end{document}R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} 0.1123\, x_{3}^{2}-0.0491x_{2}x_{3}+5.912x_{3}-3.76\frac{x_{3}^{2}}{x_{1}} \\ -0.0107x_{2}^{2}+0.0237x_{1}x_{2}-2.93x_{2}-\frac{0.0108x_{3}^{3}}{x_{2}}-1.4512x_{3} \\ 0 \\ \end{array}} \right] $$\end{document}U2=u1u2u3=0.1123x32-0.0491x2x3+5.912x3-3.76x32x1-0.0107x22+0.0237x1x2-2.93x2-0.0108x33x2-1.4512x30\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{2}$$\end{document}R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} 0.1581x_{3}^{2}-0.0688x_{2}x_{3}+5.912x_{3} \\ -0.0107x_{2}^{2}+0.0262x_{1}x_{2}-3.92x_{2}-1.4512x_{3}-\frac{0.0108x_{3}^{3}}{x_{2}} \\ -0.0108x_{3}^{2}-1.4271x_{3} \\ \end{array}} \right] $$\end{document}U2=u1u2u3=0.1581x32-0.0688x2x3+5.912x3-0.0107x22+0.0262x1x2-3.92x2-1.4512x3-0.0108x33x2-0.0108x32-1.4271x3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{3}$$\end{document}R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} 0.1123\, x_{3}^{2}-0.0491x_{2}x_{3}+5.912x_{3}-\frac{x_{3}^{2}}{x_{1}} \\ -0.0107x_{2}^{2}+0.0237x_{1}x_{2}-2.93x_{2}-\frac{0.0108x_{3}^{3}}{x_{2}}-1.4512x_{3} \\ -x_{3} \\ \end{array}} \right] $$\end{document}U3=u1u2u3=0.1123x32-0.0491x2x3+5.912x3-x32x1-0.0107x22+0.0237x1x2-2.93x2-0.0108x33x2-1.4512x3-x3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{3}$$\end{document}R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} 0.1581x_{3}^{2}-0.0688x_{2}x_{3}+5.912x_{3} \\ -0.0107x_{2}^{2}+0.0262x_{1}x_{2}-{5.98x}_{2}-1.4512x_{3}-\frac{0.0108x_{3}^{3}}{x_{2}} \\ -0.0108x_{3}^{2}-1.09374x_{3} \\ \end{array}} \right] $$\end{document}U3=u1u2u3=0.1581x32-0.0688x2x3+5.912x3-0.0107x22+0.0262x1x2-5.98x2-1.4512x3-0.0108x33x2-0.0108x32-1.09374x3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{1}$$\end{document}P1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{1}$$\end{document}R1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} 0.1231\, x_{3}^{2}-0.0528x_{2}x_{3}+6.9174x_{3} \\ -0.0107x_{2}^{2}+0.0237x_{1}x_{2}-x_{2}-\frac{0.0213x_{3}^{3}}{x_{2}}-2.0814x_{3} \\ 0.0013x_{2}^{2}-x_{3} \\ \end{array}} \right] $$\end{document}U1=u1u2u3=0.1231x32-0.0528x2x3+6.9174x3-0.0107x22+0.0237x1x2-x2-0.0213x33x2-2.0814x30.0013x22-x3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{3}$$\end{document}P3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{1}$$\end{document}R1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} \, 0.1291x_{3}^{2}-0.04906x_{2}x_{3}+0.004x_{1}x_{3}+6.007x_{3} \\ -0.0103x_{2}^{2}-1.7263x_{3}+0.0227x_{1}x_{2}+\left( 0.00007 \right) x_{2}x_{3}-x_{2} \\ -0.0318x_{3}^{2}-x_{3} \\ \end{array}} \right] $$\end{document}U1=u1u2u3=0.1291x32-0.04906x2x3+0.004x1x3+6.007x3-0.0103x22-1.7263x3+0.0227x1x2+0.00007x2x3-x2-0.0318x32-x3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{2}$$\end{document}R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =$$\end{document}U2=u1u2u3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left[ {\begin{array}{c} 0.1231\, x_{3}^{2}-0.0528x_{2}x_{3}+6.9174x_{3} \\ -0.0107x_{2}^{2}+0.0237x_{1}x_{2}-{4.77x}_{2}-\frac{0.0213x_{3}^{3}}{x_{2}}-2.0814x_{3} \\ 0.0013x_{2}^{2}-{2.12x}_{3} \\ \end{array}} \right] $$\end{document}0.1231x32-0.0528x2x3+6.9174x3-0.0107x22+0.0237x1x2-4.77x2-0.0213x33x2-2.0814x30.0013x22-2.12x3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{2}$$\end{document}R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} \, 0.1291x_{3}^{2}-0.04906x_{2}x_{3}+0.004x_{1}x_{3}+6.007x_{3} \\ -0.0103x_{2}^{2}-1.7263x_{3}+0.0227x_{1}x_{2}+\left( 0.00007 \right) x_{2}x_{3}-4.57x_{2} \\ -0.0318x_{3}^{2}-1.75x_{3} \\ \end{array}} \right] $$\end{document}U2=u1u2u3=0.1291x32-0.04906x2x3+0.004x1x3+6.007x3-0.0103x22-1.7263x3+0.0227x1x2+0.00007x2x3-4.57x2-0.0318x32-1.75x3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{3}$$\end{document}R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =$$\end{document}U3=u1u2u3= \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left[ {\begin{array}{c} 0.1231\, x_{3}^{2}-0.0528x_{2}x_{3}+6.9174x_{3} \\ -0.0107x_{2}^{2}+0.0237x_{1}x_{2}-1.3164x_{2}-\frac{0.0213x_{3}^{3}}{x_{2}}-2.0814x_{3} \\ 0.0013x_{2}^{2}-6.247x_{3} \\ \end{array}} \right] $$\end{document}0.1231x32-0.0528x2x3+6.9174x3-0.0107x22+0.0237x1x2-1.3164x2-0.0213x33x2-2.0814x30.0013x22-6.247x3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{3}$$\end{document}R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}=\left[ {\begin{array}{c} u_{1} \\ u_{2} \\ u_{3} \\ \end{array}} \right] =\left[ {\begin{array}{c} \, 0.1291x_{3}^{2}-0.04906x_{2}x_{3}+0.004x_{1}x_{3}+6.007x_{3} \\ -0.0103x_{2}^{2}-1.7263x_{3}+0.0227x_{1}x_{2}+\left( 0.00007 \right) x_{2}x_{3}-6.335x_{2} \\ -0.0318x_{3}^{2}-9.59x_{3} \\ \end{array}} \right] $$\end{document}U3=u1u2u3=0.1291x32-0.04906x2x3+0.004x1x3+6.007x3-0.0103x22-1.7263x3+0.0227x1x2+0.00007x2x3-6.335x2-0.0318x32-9.59x3
Proposed TSFLCs to stabilise chaos in FOPMs in antecedent

Design of TS-fuzzy logic synchronisation controllers (TSFLSCs) for synchronisation of chaos in FOPMs

The aim is to design a synchronisation controller u such that the dynamics of the slave system synchronises with the master system. Let the periodic and chaotic system as master and slave be defined in (18) and (19), respectively,where andwhere is the synchronisation controller added to synchronise the slave with the master. The error states between the state variables of the master and slave systems are defined as in (20).Thus, the FO error dynamics are obtained as in (21).Therefore, the goal is to design such that e converge to zero as i.e., the trajectory of the slave system asymptotically approaches the trajectory of the master system. The fuzzy control rule for TSFLSC is defined as in (22),where are fuzzy sets describing the LTs of input variables, is the fuzzy logic synchronisation control input of rule following the same defuzzification process described in Sect. 3.1. Each fuzzy rule defined in (22) generates a weight as in (23).From (23), it is assumed that, for any in the input universe of discourse E, there exists at least one weighted output among all rules that is non-zero. From the weighted sum of the defuzzification process, the output is as in (24),where r is the total number of fuzzy rules. A quadratic Lyapunov function is defined as in (25) which satisfies condition (a) of Theorem 1.From Theorem 2, it can be concluded that, if the origin is an equilibrium of the controlled system (19) and is a positive definite function on domain X containing the origin of such that , for and , and , and suppose no solution of (19) can stay identically in S except the trivial solution , then, the equilibrium at the origin is asymptotically stable in the domain E. We now apply TSFLSC to synchronise the chaotic FOPM (3) as slave with the periodic as master, which is described in (26).From (21), the error dynamics are given by (27).The parameters of Fig. 1b chosen are so that . From (25), the Lyapunov function is chosen as in (28).Using Caputo derivative in (28), we have (29).The fuzzy rules are derived as follows. then from (29), Then from (29), . Then from (29),. Fuzzy rule 1, : for antecedent as P, we have the consequent as Fuzzy rule 2, : for antecedent as N, we have the consequent as Fuzzy rule 3, : for antecedent as Z, we have the consequent as Since which implies . Thus, the errors, converge to zero and the slave synchronises with the master . Similarly, the TSFLSC designed to synchronise FOPMs, the chaotic with the periodic is as in Table 4.
Table 4

Proposed TSFLSC for synchronisation in FOPMs, (slave) with (master); and (slave) with (master) at antecedent

Master systemSlave systemRule \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\varvec{i}}}$$\end{document}iConsequent synchronisation controller
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{4}$$\end{document}P4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{0}$$\end{document}P0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{1}$$\end{document}R1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}=\left[ {\begin{array}{l} 0.0976{(x}_{3}^{2}-y_{3}^{2})-0.045(x_{2}x_{3}-y_{2}y_{3})+7.5489726x_{1}-12.6237y_{1}\, -e_{1}+5.912e_{3} \\ -0.0107{(x}_{2}^{2}-y_{2}^{2})+0.0237(x_{1}x_{2}-y_{1}y_{2})-e_{2}-1.4512e_{3} \\ -0.0108\left( x_{3}^{2}-y_{3}^{2} \right) +0.0147{(x}_{1}x_{3}-y_{1}y_{3})-0.0041(x_{1}x_{2}-y_{1}y_{2})-e_{3} \\ \end{array}} \right] $$\end{document}U1=0.0976(x32-y32)-0.045(x2x3-y2y3)+7.5489726x1-12.6237y1-e1+5.912e3-0.0107(x22-y22)+0.0237(x1x2-y1y2)-e2-1.4512e3-0.0108x32-y32+0.0147(x1x3-y1y3)-0.0041(x1x2-y1y2)-e3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{2}$$\end{document}R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}=\left[ {\begin{array}{l} 0.0976{(x}_{3}^{2}-y_{3}^{2})-0.045(x_{2}x_{3}-y_{2}y_{3})+7.5489726x_{1}-12.6237y_{1}-6.5e_{1}+5.912e_{3} \\ -0.0107{(x}_{2}^{2}-y_{2}^{2})+0.0237(x_{1}x_{2}-y_{1}y_{2})\, -9.5e_{2}-1.4512e_{3} \\ -0.0108\left( x_{3}^{2}-y_{3}^{2} \right) +0.0147{(x}_{1}x_{3}-y_{1}y_{3})-0.0041(x_{1}x_{2}-y_{1}y_{2})-e_{3} \\ \end{array}} \right] $$\end{document}U2=0.0976(x32-y32)-0.045(x2x3-y2y3)+7.5489726x1-12.6237y1-6.5e1+5.912e3-0.0107(x22-y22)+0.0237(x1x2-y1y2)-9.5e2-1.4512e3-0.0108x32-y32+0.0147(x1x3-y1y3)-0.0041(x1x2-y1y2)-e3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{3}$$\end{document}R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}=\left[ {\begin{array}{l} 0.0976{(x}_{3}^{2}-y_{3}^{2})-0.045(x_{2}x_{3}-y_{2}y_{3})+7.5489726x_{1}-12.6237y_{1}-6.5e_{1}+5.912e_{3} \\ -0.0107{(x}_{2}^{2}-y_{2}^{2})+0.0237(x_{1}x_{2}-y_{1}y_{2})\, -9.5e_{2}-1.4512e_{3} \\ -0.0108\left( x_{3}^{2}-y_{3}^{2} \right) +0.0147{(x}_{1}x_{3}-y_{1}y_{3})-0.0041(x_{1}x_{2}-y_{1}y_{2})-13e_{3} \\ \end{array}} \right] $$\end{document}U3=0.0976(x32-y32)-0.045(x2x3-y2y3)+7.5489726x1-12.6237y1-6.5e1+5.912e3-0.0107(x22-y22)+0.0237(x1x2-y1y2)-9.5e2-1.4512e3-0.0108x32-y32+0.0147(x1x3-y1y3)-0.0041(x1x2-y1y2)-13e3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{5}$$\end{document}P5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{2}$$\end{document}P2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{1}$$\end{document}R1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{1}=\left[ {\begin{array}{l} 0.1434{(x}_{3}^{2}-y_{3}^{2})-0.0674(x_{2}x_{3}-y_{2}y_{3})+0.0025{(x}_{2}^{2}-y_{2}^{2})+11.45376x_{1}-12.7264y_{1}-e_{1}+5.912e_{3} \\ -0.0107{(x}_{2}^{2}-y_{2}^{2})+0.0237(x_{1}x_{2}-y_{1}y_{2})-e_{2}-1.4512e_{3} \\ -0.0108\left( x_{3}^{2}-y_{3}^{2} \right) +0.0147{(x}_{1}x_{3}-y_{1}y_{3})-0.0041(x_{1}x_{2}-y_{1}y_{2})-e_{3} \\ \end{array}} \right] $$\end{document}U1=0.1434(x32-y32)-0.0674(x2x3-y2y3)+0.0025(x22-y22)+11.45376x1-12.7264y1-e1+5.912e3-0.0107(x22-y22)+0.0237(x1x2-y1y2)-e2-1.4512e3-0.0108x32-y32+0.0147(x1x3-y1y3)-0.0041(x1x2-y1y2)-e3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{2}$$\end{document}R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{2}=\left[ {\begin{array}{l} {0.1434(x}_{3}^{2}-y_{3}^{2})-0.0674(x_{2}x_{3}-y_{2}y_{3})+0.0025{(x}_{2}^{2}-y_{2}^{2})+11.45376x_{1}-12.7264y_{1}-7.5e_{1}+5.912e_{3} \\ -0.0107{(x}_{2}^{2}-y_{2}^{2})+0.0237(x_{1}x_{2}-y_{1}y_{2})-10.5e_{2}-1.4512e_{3} \\ -0.0108\left( x_{3}^{2}-y_{3}^{2} \right) +0.0147{(x}_{1}x_{3}-y_{1}y_{3})-0.0041(x_{1}x_{2}-y_{1}y_{2})-e_{3} \\ \end{array}} \right] $$\end{document}U2=0.1434(x32-y32)-0.0674(x2x3-y2y3)+0.0025(x22-y22)+11.45376x1-12.7264y1-7.5e1+5.912e3-0.0107(x22-y22)+0.0237(x1x2-y1y2)-10.5e2-1.4512e3-0.0108x32-y32+0.0147(x1x3-y1y3)-0.0041(x1x2-y1y2)-e3
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{3}$$\end{document}R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_{3}=\left[ {\begin{array}{l} {0.1434(x}_{3}^{2}-y_{3}^{2})-0.0674(x_{2}x_{3}-y_{2}y_{3})+0.0025{(x}_{2}^{2}-y_{2}^{2})+11.45376x_{1}-12.7264y_{1}-7.5e_{1}+5.912e_{3} \\ -0.0107{(x}_{2}^{2}-y_{2}^{2})+0.0237(x_{1}x_{2}-y_{1}y_{2})-10.5e_{2}-1.4512e_{3} \\ -0.0108\left( x_{3}^{2}-y_{3}^{2} \right) +0.0147{(x}_{1}x_{3}-y_{1}y_{3})-0.0041(x_{1}x_{2}-y_{1}y_{2})-15e_{3} \\ \end{array}} \right] $$\end{document}U3=0.1434(x32-y32)-0.0674(x2x3-y2y3)+0.0025(x22-y22)+11.45376x1-12.7264y1-7.5e1+5.912e3-0.0107(x22-y22)+0.0237(x1x2-y1y2)-10.5e2-1.4512e3-0.0108x32-y32+0.0147(x1x3-y1y3)-0.0041(x1x2-y1y2)-15e3
Proposed TSFLSC for synchronisation in FOPMs, (slave) with (master); and (slave) with (master) at antecedent Attractors and bifurcation of period doubling route to chaos in FO plague model

Results and discussion

Results obtained from numerical simulation are discussed in the following subsections.

Dynamical analyses of FOPM

The bifurcation analysis of FOPM (4) is performed in Fig. 2a for a varying FO as the bifurcation parameter against that represents the number of captured infected rats of M. rattus species. It is revealed that plague follows a path of multiple period doublings to finally give rise to chaos at . The Period 1 attractor (P1A) that persists in the region is shown in Fig. 2b. P1A doubles its period to a Period 2 attractor (P2A) that persists in the region and is displayed in Fig. 2c. Similarly, a Period 4 attractor (P4A) existing in the range is illustrated in Fig. 2d and the final chaotic attractor from onwards in Fig. 3b.
Fig. 2

Attractors and bifurcation of period doubling route to chaos in FO plague model

Fig. 3

Attractors of proposed FO plague models

Figure 3 depicts the dynamics of attractors of the FOPMs as investigated in Table 1 of Sect. 2.2. It is observed that display chaotic dynamics and display periodic dynamics. Attractors of proposed FO plague models

Stabilisation of chaos in FOPMs using TSFLC

Since, only the FOPMs, display chaotic dynamics, hence the controllers designed using TSFLC concept in Sect. 3.1 are added to the corresponding plague models at 10 s. It is observed from the plots of Fig. 4 that chaos dies out within 4 s once the controller is applied.
Fig. 4

Stabilisation of chaos in FO plague epidemic models when controller is added at 10 s

Stabilisation of chaos in FO plague epidemic models when controller is added at 10 s Synchronisation control between (master) and (slave) plague epidemic models using TSFLSC Synchronisation control between (master) and (slave) plague epidemic models using TSFLSC

Synchronisation of chaos in FOPM using TSFLSC

Another way of suppressing chaotic dynamics in plague epidemic model is using synchronisation controller. The significance lies in the fact that the uncontrolled chaotic dynamics of the slave epidemic model after controlled synchronisation becomes asymptotically periodic which means that the motion of its trajectories repeats after fixed intervals of time with a constant period and is imperturbable to minute changes in initial states of physical parameters. In real situations, such definite repetitive dynamics makes prediction of the epidemic model easier. The numerical results of synchronisation controllers designed in Sect. 3.2 are plotted in Figs. 5 and 6, where it is clear that on adding the TSFLSCs at 10 s, the initially chaotic slaves and display periodic dynamics once they are synchronised with the masters and , respectively.
Fig. 5

Synchronisation control between (master) and (slave) plague epidemic models using TSFLSC

Fig. 6

Synchronisation control between (master) and (slave) plague epidemic models using TSFLSC

Comparative analyses between the plague epidemic and COVID-19 pandemic models

The IOM of chaotic COVID-19 pandemic is modelled using global modelling technique [31] obtained from two real data sets: the official data from the National Health Commission of the People’s Republic of China [32] and the data from the Johns Hopkins University [33]. The FO COVID-19 model (FOCM) is proposed in (30),where are the numbers of confirmed cases, critical cases under intensive care and cumulative number of fatalities per day, respectively. In terms of dynamical analysis, it is observed that the FOCM (30) undergoes period doubling bifurcation to generate chaotic transmission with the gradual increase in FO; similar to that of the plague model. The bifurcation diagram plotted in Fig. 7 for the trajectory corresponding to cumulative number of COVID-19 infected deaths per day, reveals that it follows a rich complex dynamical path through successive cascades of period doubling which otherwise remain unobserved in the IO chaotic model. The FOCM displays a P1A in , which gets converted to a P2A in [0.96, 0.97]. It splits further to transform into a P4A that gives birth to chaotic attractor after infinite doublings from onwards.
Fig. 7

Period doubling bifurcation of FO COVID-19 pandemic model

The corresponding attractors of FOCM are plotted in Fig. 8 that shows asymptotically period doubling route to chaos.
Fig. 8

Phase portraits of FO COVID-19 pandemic model

The above results and analyses lead to the following qualitative comparisons between the plague epidemic and the COVID-19 pandemic situation in India: Both plague and COVID-19 propagation dynamics are governed by chaotic behaviour where an apparently trivial event may trigger a series of divergent events, indicating that a long-term prediction of the spread of the epidemic is not possible. One of the major shortcomings of recent mathematical models to exactly predict the course of COVID-19 is that they do not take into account the chaotic behaviour of the spread [34]. COVID-19 is a multi-strain pandemic, a characteristic similar to plague. More than 100 strains of COVID-19 have been discovered as reported in Brazil [35], the United Kingdom [36], South Africa [37], etc. The Covid-19 virus has undergone multiple mutations of which the variants named as ‘Kappa’ and ‘Delta’ by WHO [38], first discovered in India, in October, 2020 are far more contagious than their predecessors having a transmission rate of 2.5 times faster. The plague epidemic sustained repeated seasonal waves lasting for over 15 years (1896–1911), while the ongoing COVID-19 pandemic too has repeated waves of resurgence pointing towards the possibility of a third wave in India [39]. Both plague and COVID-19 dispersal in India have displayed close similarities demographically as well as geographically. The Bombay Plague first hit Mumbai in 1896 and then spread all over India [5]. Similarly, Mumbai was one of the first and worst hit cities in India in both the first and second waves of Covid-19 which gradually disseminated with a rapid transmission rate throughout the country. Period doubling bifurcation of FO COVID-19 pandemic model The control framework in real scenario may reflect cost effective treatments, social mobility restrictions, injection of vaccines and prophylactic drugs, delays, effective planning on medical infrastructure as hospital beds, ventilators, oxygen cylinders, ICU berths, ambulance, etc. and providing other health facilities. In our future scope of work, we shall explore the design of the controllers based on these actual control measures. Sensing and control strategies can be applied for close monitoring of the epidemic activities using real time technology such as remote sensors, web-based tools for surveillance, tracking individual movements, e- government initiatives such as the Aarogya Setu [40], high throughput sequencing technology, etc. [41]. In fact, lower-middle income countries like India need judicious planning on executing control policies to evaluate the risk and cost of prophylactic strategies and public health policies depending on the short-term estimation of the chaotic course of the pandemic based on past experience. Phase portraits of FO COVID-19 pandemic model

Conclusions

This paper investigates the complex dynamics of fractional-order models of the past plague epidemic which hit India over a hundred years ago and lasted for more than a decade. We also analyse the FOM of COVID-19 pandemic and found that both epidemics are multi-strain models, have probable seasonality, drastic transmission rate in high population density areas of India, manifest chaotic wave of propagation, all striking close similarities demographically and geographically. Both plague and COVID-19 are found to exhibit chaotic behaviour, which means a small change in a trivial event may activate a series of divergent unpredictable events due to their hyped sensitivity to initial states of driving parameters. Two types of controllers are proposed to suppress chaos in the epidemic: first, stabilisation control, where a TS fuzzy logic controller is designed that successfully stabilises the chaotic dynamics of the epidemic by minimizing the numbers of infections, critical cases and fatalities to a disease-free equilibrium, and second, synchronisation control, where the proposed controller efficiently synchronises the chaotic dynamics of the slave with the regular periodic dynamics of the master. It is significant since after synchronisation, the originally chaotic slave epidemic model transforms into a periodic one which means that its flow repeats after fixed intervals of time with a constant period and is no longer sensitive to small perturbations of parameters. In real situations, such repeated regular dynamics facilitates prediction of the course of outbreak, since the next surge or spike in infections occurring at a periodic regular interval can be estimated beforehand and an early alarm can be set. The proposed FOMs, possessing the features of flexibility, heredity and memory, can provide more accurate models of the chaotic pandemics by incorporating unmodelled dynamics that often fail to be encompassed by integer-order models. It was discovered that in both plague and COVID-19 models, with stronger index of memory (lower values of fractional-order), chaos subsides into regular periodic dynamics, while as memory weakens and fades away towards an integral order, the route to chaos in transmission dynamics is evident. In our future scope of work, we shall explore the design of the controllers based on real control measures Thus, the authors are hopeful that the results on the chaotic dynamics of propagation of the epidemic will be helpful to recognise regions with high risk of contamination and initiate timely action on crisis assessment, hazard management, policy making, emergency planning and economic control measures, i.e. learning from the knowledge of the past.
  11 in total

1.  XXII. Epidemiological observations in Bombay City.

Authors: 
Journal:  J Hyg (Lond)       Date:  1907-12

2.  The Reported Appearance of Plague in Bombay.

Authors: 
Journal:  Br Med J       Date:  1896-10-03

3.  The model of Kermack and McKendrick for the plague epidemic in Bombay and the type reproduction number with seasonality.

Authors:  Nicolas Bacaër
Journal:  J Math Biol       Date:  2011-03-15       Impact factor: 2.259

4.  Optimal vaccination and treatment policies for regional approximate controllability of the time-fractional reaction-diffusion SIR epidemic systems.

Authors:  Fudong Ge; YangQuan Chen
Journal:  ISA Trans       Date:  2021-01-15       Impact factor: 5.468

5.  Entire functions sharing a small function with their two difference operators.

Authors:  Feng Lü; Yanfeng Wang; Junfeng Xu
Journal:  Adv Differ Equ       Date:  2017-08-01

6.  Global dynamics of a fractional order model for the transmission of HIV epidemic with optimal control.

Authors:  Parvaiz Ahmad Naik; Jian Zu; Kolade M Owolabi
Journal:  Chaos Solitons Fractals       Date:  2020-06-18       Impact factor: 5.944

7.  Chaos theory applied to the outbreak of COVID-19: an ancillary approach to decision making in pandemic context.

Authors:  S Mangiarotti; M Peyre; Y Zhang; M Huc; F Roger; Y Kerr
Journal:  Epidemiol Infect       Date:  2020-05-08       Impact factor: 2.451

8.  A hybrid stochastic fractional order Coronavirus (2019-nCov) mathematical model.

Authors:  N H Sweilam; S M Al-Mekhlafi; D Baleanu
Journal:  Chaos Solitons Fractals       Date:  2021-02-10       Impact factor: 5.944

9.  Development of a multiplex RT-qPCR using the drop out strategy to screen the SARS-CoV-2 South African 501Y.V2 variant.

Authors:  François Maclot; Sébastien Bontems; Cécile Meex; Maria Artesi; Pablo Beckers; Vincent Bours; Keith Durkin; Marie-Pierre Hayette
Journal:  J Infect       Date:  2021-05-03       Impact factor: 6.072

10.  A revisit to the past plague epidemic (India) versus the present COVID-19 pandemic: fractional-order chaotic models and fuzzy logic control.

Authors:  Manashita Borah; Binoy Krishna Roy; Tomasz Kapitaniak; Karthikeyan Rajagopal; Christos Volos
Journal:  Eur Phys J Spec Top       Date:  2021-12-12       Impact factor: 2.891
View more
  2 in total

1.  A revisit to the past plague epidemic (India) versus the present COVID-19 pandemic: fractional-order chaotic models and fuzzy logic control.

Authors:  Manashita Borah; Binoy Krishna Roy; Tomasz Kapitaniak; Karthikeyan Rajagopal; Christos Volos
Journal:  Eur Phys J Spec Top       Date:  2021-12-12       Impact factor: 2.891

2.  Complex bio rhythms.

Authors:  Akif Akgul; Marcelo A Savi; Mustafa Zahid Yildiz; Miguel A F Sanjuan; Jun Ma
Journal:  Eur Phys J Spec Top       Date:  2022-04-20       Impact factor: 2.891
  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.