Next Article in Journal
Prediction in Chaotic Environments Based on Weak Quadratic Classifiers
Next Article in Special Issue
Towards Dynamic Uncertain Causality Graphs for the Intelligent Diagnosis and Treatment of Hepatitis B
Previous Article in Journal
Risk Analysis and Control Factors Based on Excavation of a Large Underground Subway Station under Construction
Previous Article in Special Issue
A Geometric Accuracy Error Analysis Method for Turn-Milling Combined NC Machine Tool
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Stochastic Intelligent Computing with Neuro-Evolution Heuristics for Nonlinear SITR System of Novel COVID-19 Dynamics

by
Muhammad Umar
1,
Zulqurnain Sabir
1,
Muhammad Asif Zahoor Raja
2,3,
Muhammad Shoaib
4,
Manoj Gupta
5 and
Yolanda Guerrero Sánchez
6,*
1
Department of Mathematics and Statistics, Hazara University, Mansehra 21120, Pakistan
2
Future Technology Research Center, National Yunlin University of Science and Technology, 123 University Road, Section 3, Douliou, Yunlin 64002, Taiwan
3
Department Electrical and Computer Engineering, COMSATS University Islamabad, Attock Campus, Attock 43600, Pakistan
4
Department of Mathematics, COMSATS University Islamabad, Attock Campus, Attock 43600, Pakistan
5
Department of Electronics and Communication Engineering, JECRC University Jaipur, Jaipur 303905, Rajasthan, India
6
Department of Dermatology, Stomatology, Radiology and Physical Medicine, University of Murcia, 30100 Murcia, Spain
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(10), 1628; https://doi.org/10.3390/sym12101628
Submission received: 9 August 2020 / Revised: 13 September 2020 / Accepted: 17 September 2020 / Published: 2 October 2020

Abstract

:
The present study aims to design stochastic intelligent computational heuristics for the numerical treatment of a nonlinear SITR system representing the dynamics of novel coronavirus disease 2019 (COVID-19). The mathematical SITR system using fractal parameters for COVID-19 dynamics is divided into four classes; that is, susceptible (S), infected (I), treatment (T), and recovered (R). The comprehensive details of each class along with the explanation of every parameter are provided, and the dynamics of novel COVID-19 are represented by calculating the solution of the mathematical SITR system using feed-forward artificial neural networks (FF-ANNs) trained with global search genetic algorithms (GAs) and speedy fine tuning by sequential quadratic programming (SQP)—that is, an FF-ANN-GASQP scheme. In the proposed FF-ANN-GASQP method, the objective function is formulated in the mean squared error sense using the approximate differential mapping of FF-ANNs for the SITR model, and learning of the networks is proficiently conducted with the integrated capabilities of GA and SQP. The correctness, stability, and potential of the proposed FF-ANN-GASQP scheme for the four different cases are established through comparative assessment study from the results of numerical computing with Adams solver for single as well as multiple autonomous trials. The results of statistical evaluations further authenticate the convergence and prospective accuracy of the FF-ANN-GASQP method.

1. Introduction

Human life has faced various challenges and obstacles throughout time in the form of floods, earthquakes, and diseases, among others. We first mention some widespread diseases such as dengue fever, which is a mosquito-borne disease produced by the dengue virus. This disease spread extensively in the Africa, South and Central America, Eastern Mediterranean, Caribbean, Southeast and South Asia, and Oceania regions. It is a very serious disease, and almost half a million people infected with dengue virus are hospitalized every year, although most of the patients recover after a few days. Some major symptoms of this disease are high fever, vomiting, headache, and joint and muscle pains. Ebola, HIV, and Lassa are other infectious diseases of paramount interest to the research community concerning safeguarding humanity. Therefore, several medical utilities have been unified by humans to prevent the propagation of these diseases.
Recently, the world has suffered the dangerous and deathly disease of coronavirus disease 2019 (COVID-19), which is a respiratory infection spreading in an almost uncontrollable manner. COVID-19 is a disease that is highly transferable from one individual to another through droplets [1]. Many people have been infected with COVID-19, with the rate of death as a result of this virus increasing daily. The deathly COVID-19 was reported for the first time on 31 December 2019 in Wuhan city, which is in the Hubei province of China [2]. Currently, there is no particular treatment or vaccine for COVID-19, and it has become a big task for all countries to control it. Many people were killed in Wuhan, China by COVID-19, which led the Chinese government to lock down the infected cities as well as put the diseased people under quarantine. The number of infected cases is in the thousands in developing countries, and the reported deaths in advanced countries such as the United States of America, Italy, Spain, England, and Russia have overtaken those in China. The key symptoms of this virus are dry cough, fever, and tiredness. COVID-19 badly affects older people (around the age of 60 or older), as well as individuals facing serious medical diseases including cardiovascular disease, diabetes, chronic respiratory, and cancer. Those who have fever, breathing difficulty, and dry cough should seek medical care immediately. The spread ratio of COVID-19 is very problematic and has become a worry for the whole world. To present the mathematical model is another task, while predicting the outcomes of this disease is very frustrating as well. One of the main challenges is how to control the disease if it spreads to animals or birds. To prevent this virus from causing further damage, the proper vaccine should be discovered; otherwise, the world should be ready to face many new challenges such as human loss, shortage of food, unemployment, poverty, and many more.
The aim of the present investigation is to find the solution to a nonlinear SITR system using parameters for the dynamics of novel COVID-19, along with numerical analysis for a better understanding of the spread using stochastic approaches through the exploitation of artificial intelligence algorithms. The substantial potential of neuro-evolution-based computing is exploited/explored for linear and nonlinear systems by the functioning of neural networks trained with the competences of local and global search schemes [3,4,5,6,7]. A few recent studies using soft computing schemes include the singular three-point boundary value model [8], financial market forecasting [9], the Thomas–Fermi model [10], summer precipitation prediction for meteorological stations [11], nonlinear prey–predator models [12], nonlinear models of fuel ignition models on Bratu’s differential equation [13], the nonlinear Troesch’s problem [14], magneto hydrodynamic studies [15], the nonlinear functional differential model [16,17], the model of heartbeat dynamics [18], fluid dynamics [19], the corneal shape model [20], nonlinear models arising in electric circuits [21], multi-singularity based nonlinear models [22], the nonlinear reactive transport model [23], the singular periodic nonlinear model [24], the heat distribution model of human head based on a nonlinear singular system [25], the atomic physics model [26], the nonlinear Riccati equation [27], and the singular Lane–Emden model [28]. Recently, fractional order systems have also been solved with artificial intelligence (AI) techniques [29,30,31,32,33]. These provide the motivation to investigate AI methodologies for the numerical solution of the mathematical model of SITR representing the dynamics of COVID-19.
The SITR model for novel COVID-19 dynamics is designed on the mechanism of four classes: susceptible class (S), infected class (I), treatment class (T), and recovered class (R). The S is further segmented into two sub-classes S 1 ( χ ) and S 2 ( χ ) . The sub-class S 1 ( χ ) specifies the individuals not yet affected from COVID-19 disease and sub-class S 2 ( χ ) represents those not yet infected but who have some serious illnesses or older. I ( χ ) shows the class of those people who are infected by COVID-19 disease, T ( χ ) is the treatment class, and the last class R ( χ ) is the recovered class from COVID-19 disease. The SITR model governing the novel COVID-19 dynamics is given as follows [34]:
S 1 ( χ ) = B β I ( χ ) S 1 ( χ ) δ β T ( χ ) α S 1 ( χ ) , S 1 ( 0 ) = A 1 , S 2 ( χ ) = B β I ( χ ) S 2 ( χ ) δ β T ( χ ) α S 2 ( χ ) , S 2 ( 0 ) = A 2 , I ( χ ) = μ I ( χ ) + β I ( χ ) S 1 ( χ ) + S 2 ( χ ) α I ( χ ) + β δ T ( χ ) + σ I ( χ ) , I ( 0 ) = A 3 , T ( χ ) = μ I ( χ ) ρ T ( χ ) α T ( χ ) + ψ T ( χ ) + ε T ( χ ) , T ( 0 ) = A 4 , R ( χ ) = α R ( χ ) + ρ T ( χ ) , R ( 0 ) = A 5 , .
The objective of the proposed study is to find the solution of the system of equations given in set (1) by integrating intelligent neuro-evolution computing heuristics via feed-forward artificial neural networks (FF-ANNs), and its adjustable parameters are tuned with an optimization procedure of genetic algorithms (GAs) integrated with sequential quadratic programming (SQP), i.e., FF-ANN-GASQP solver. The parameter of interest with appropriate settings for particular descriptions of the dynamics are provided in Table 1 and Table 2, while the block diagram of the proposed study is given in Figure 1, which presents the system model, mathematical formulation, proposed methodology, evaluation criteria, and main results in pictorial form.
The major features of the suggested FF-ANN-GASQP scheme are briefly given as:
  • The solution of the mathematical expression for the nonlinear SITR model for novel COVID-19 dynamics is calculated viably by using the novel application of the intelligent neuro-evolution-based integrated computing paradigm, i.e., FF-ANN-GASQP.
  • Closely matching of the results of the proposed FF-ANN-GASQP solver with the solutions of the reference state of the art numerical procedure of Adams methods for variants of the nonlinear SITR-based mathematical model established the value and worth.
  • Authentication and verification of the performance through statistical assessments studies is proven on multiple implementations of FF-ANN-GASQP in terms of Theil’s inequality coefficient (TIC) as well as root mean square error (RMSE)-based indices.
  • In addition to the precise and accurate solutions for the SITR-based mathematical model of the COVID-19 pandemic, other valuable perks are that it is easy to understand the concepts, and it also has smooth operation, exhaustive applicability, consistency, and extendibility.
The remaining sections of the current study are organized as follows. The designed framework of FF-ANN-GASQP is given in Section 2, measures of the performance are listed in Section 3, a description of the results with necessary interpretations are provided in Section 4, and concluding remarks with related potential studies are given in the last section.

2. Proposed Methodology

ANN is an artificial intelligence-based computing approach that is used widely to effectively solve a variety of challenging stiff linear/nonlinear problems in terms of reliability, stability, and robustness. In this study, the strength of feed-forward ANN is exploited to solve the SITR model representing novel COVID-19. The weights of ANN are optimized with a genetic algorithms (GA)-based evolutionary metaheuristic along with sequential quadratic programming (SQP). The proposed FF-ANN-GASQP structure for getting the numerical solutions of the nonlinear SITR model based on COVID-19 is described in two phases:
  • To exploit the FF-ANN models, an error-based objective function is introduced.
  • Optimize the objective function for system (1) using the hybrid GA-SQ programming approach.

2.1. ANN Modeling

The mathematical measures using system (1) are given by the differential mapping of FF-ANN for approximate solutions S ^ 1 ( χ ) , S ^ 2 ( χ ) , I ^ ( χ ) , T ^ ( χ ) and R ^ ( χ ) along with their nth order derivatives, which are given as:
S ^ 1 ( χ ) , S ^ 2 ( χ ) , I ^ ( χ ) , T ^ ( χ ) , R ^ ( χ ) = i = 1 m ϕ S 1 , i P ( w S 1 , i χ + b S 1 , i ) , i = 1 m ϕ S 2 , i P ( w S 2 , i χ + b S 2 , i ) , i = 1 m ϕ I , i P ( w I , i χ + b I , i ) , i = 1 m ϕ T , i P ( w T , i χ + b T , i ) , i = 1 m ϕ R , i P ( w R , i χ + b R , i ) , S ^ 1 ( n ) ( χ ) , S ^ 2 ( n ) ( χ ) , I ^ ( n ) ( χ ) , T ^ ( n ) ( χ ) , R ^ ( n ) ( χ ) = i = 1 m ϕ S 1 , i P ( n ) ( w S 1 , i χ + b S 1 , i ) , i = 1 m ϕ S 2 , i P ( n ) ( w S 2 , i χ + b S 2 , i ) , i = 1 m ϕ I , i P ( n ) ( w I , i χ + b I , i ) , i = 1 m ϕ T , i P ( n ) ( w T , i χ + b T , i ) , i = 1 m ϕ R , i P ( n ) ( w R , i χ + b R , i ) ,
where the unknown weight vector is W and given as:
W = [ W S 1 , W S 2 , W I , W T , W R ] ,   for   W S 1 = [ ϕ S 1 , w S 1 , b S 1 ] ,   W S 2 = [ ϕ S 2 , w S 2 , b S 2 ] ,   W I = [ ϕ I , w I , b I ] ,   W T = [ ϕ T , w T , b T ]   and   W R = [ ϕ R , w R , b R ]
where
ϕ S 1 = [ ϕ S 1 , 1 , ϕ S 1 , 2 , ϕ S 1 , 3 , , ϕ S 1 , m ] , ϕ S 2 = [ ϕ S 2 , 1 , ϕ S 2 , 2 , ϕ S 2 , 3 , , ϕ S 2 , m ] , ϕ I = [ ϕ I , 1 , ϕ I , 2 , ϕ I , 3 , , ϕ I , m ] , ϕ T = [ ϕ T , 1 , ϕ T , 2 , ϕ T , 3 , , ϕ T , m ] , ϕ R = [ ϕ R , 1 , ϕ R , 2 , ϕ R , 3 , , ϕ R , m ] , ω S 1 = [ w S 1 , 1 , w S 1 , 2 , w S 1 , 3 , , w S 1 , m ] , ω S 2 = [ w S 2 , 1 , w S 2 , 2 , w S 2 , 3 , , w S 2 , m ] , ω I = [ w I , 1 , w I , 2 , w I , 3 , , w I , m ] , ω T = [ w T , 1 , w T , 2 , w T , 3 , , w T , m ] , ϕ R = [ w R , 1 , w R , 2 , w R , 3 , , w R , m ] , b S 1 = [ b S 1 , 1 , b S 1 , 2 , b S 1 , 3 , , b S 1 , m ] , b S 2 = [ b S 2 , 1 , b S 2 , 2 , b S 2 , 3 , , b S 2 , m ] , b I = [ b I , 1 , b I , 2 , b I , 3 , , b I , m ] , b T = [ b T , 1 , b T , 2 , b T , 3 , , b T , m ] , b R = [ b R , 1 , b R , 2 , b R , 3 , , b R , m ] .
Using the log-sigmoid, i.e., an function activation P ( χ ) = ( 1 + e χ ) 1 , the updated form of system (2) becomes:
S ^ 1 ( χ ) , S ^ 2 ( χ ) , I ^ ( χ ) , R ^ ( χ ) , T ^ ( χ ) = i = 1 m ϕ S 1 , i 1 + e w S 1 , i χ + b S 1 , i , i = 1 m ϕ S 2 , i 1 + e w S 2 , i χ + b S 2 , i , i = 1 m ϕ I , i 1 + e w I , i χ + b I , i , i = 1 m ϕ T , i 1 + e T I , i χ + b T , i , i = 1 m ϕ R , i 1 + e w R , i χ + b R , i , S ^ 1 ( χ ) , S ^ 2 ( χ ) , I ^ ( χ ) , R ^ ( χ ) , T ^ ( χ ) = i = 1 m w S 1 , i ϕ E , i e ρ S 1 , i χ + b S 1 , i 1 + e w S 1 , i χ + b S 1 , i 2 , i = 1 m w S 2 , i ϕ S 2 , i e w S 2 , i χ + b S 2 , i 1 + e w S 2 , i χ + b S 2 , i 2 , i = 1 m w I , i ϕ I , i e w I , i χ + b I , i 1 + e w I , i χ + b I , i 2 , i = 1 m w T , i ϕ T , i e w T , i χ + b T , i 1 + e w T , i χ + b T , i 2 , i = 1 m w R , i ϕ A r , i e w R , i χ + b R , i 1 + e w R , i χ + b R , i 2 .
Using the values of system (3), an error-based objective function e is given as:
e = e 1 + e 2 + e 3 + e 4 + e 5 + e 6 .
e 1 = 1 N m = 1 N ( S ^ 1 ) m B + β I m ( S ^ 1 ) m + δ β T m + α ( S ^ 1 ) m 2 ,
e 2 = 1 N m = 1 N ( S ^ 2 ) m B + β I m ( S ^ 2 ) m + δ β T m + α ( S ^ 2 ) m 2 ,
e 3 = 1 N m = 1 N I ^ m + μ I m β I m ( S ^ 1 ) m + ( S ^ 2 ) m + α I m β δ T m σ I m 2 ,
e 4 = 1 N m = 1 N T ^ m μ I m + ρ T ^ m + α T ^ m ψ T ^ m ε T ^ m 2 ,
e 5 = 1 N m = 1 N R ^ m + α R ^ m ρ T ^ m 2 ,
e 6 = 1 5 ( S ^ 1 ) 0 A 1 2 + ( S ^ 2 ) 0 A 2 2 + I ^ 0 A 3 2 + T ^ 0 A 4 2 + R ^ 0 I 5 2 ,
where N h = 1 , χ m = m h , ( S ^ 1 ) m = S ^ 1 χ m , ( S ^ 2 ) m = S ^ 2 χ m , I ^ m = I ^ χ m , T ^ m = T ^ χ m and R ^ m = R ^ χ m . The approximate results for both sub-classes of susceptible S 1 and S 2 , infected class I, treatment class T, and recovered class R are denoted as S ^ 1 , S ^ 2 , I ^ , T ^ , and R ^ . Accordingly, e 1 , e 2 , e 3 , e 4 , e 5 , and e 6 are the error-based functions for system dynamics and initial conditions of the SITR model representing COVID-19 as given in the set of Equation in (1). The proposed results are achieved from the adjustable weights for which the objective function as shown in system (4) tends to zero, i.e., mathematically e 0 . Hence, approximate solutions S ^ 1 ( χ ) , S ^ 2 ( χ ) , I ^ ( χ ) , T ^ ( χ ) , R ^ ( χ ) are matched with the exact results, i.e., [ S ^ 1 ( χ ) S ^ 1 ( χ ) ] , [ S ^ 2 ( χ ) S ^ 2 ( χ ) ] , [ I ^ ( χ ) I ( χ ) ] , [ T ^ ( χ ) T ( χ ) ] and [ R ^ ( χ ) R ( χ ) ] .

2.2. Optimization Technique: Hybrid of GA with SQP

In this section, a brief overview of GASQP is given to optimize networks for system (4) representing COVID-19 dynamics.
GAs is an efficient global search optimization approach, which was introduced by Holand [30,31]. GA is one of the mathematical models of human genetics, which is implemented as a global search tool for linear/nonlinear optimization processes in assorted domains of applied sciences. Population in GAs represents the set of possible solutions as per the obligation of the optimizing systems, while the populations are modified through the exploitation of selection, crossover, elitism, and mutation such that an appropriate outcome is achieved. Some well-known applications of GA in various fields include multilayer piezoelectric transducer [32,33,34,35,36,37], Hammerstein controlled autoregressive models [38], nonlinear circuits [39], economic load dispatch [40], active noise control system [41], and power signal modeling [42].
The GA optimization procedure is further enhanced normally using the hybridization process with the speedy local search approach, so in this study, SQP is applied for a fine-tuning mechanism. The GA optimization procedure is further enhanced normally using the hybridization process with the speedy local search approach, so in this study, SQP is applied for nonlinear models in combustion theory [43], Bagley–Torvik systems [44], and system identification [45]. In the presented investigation, the integrated metaheuristics through GASQP is exploited to optimize the unknown weights of FF-ANN models, which are demonstrating the SITR model based on COVID-19. The pseudocode detail of FF-ANN-GASQP for a nonlinear SITR system is given in Table 3.

3. Performance Measures

The mathematical form of the performance grades using the root mean square error (RMSE) and Theil’s inequality coefficient (TIC) for SITR model based on COVID-19 is presented as:
RMS E S 1 , RMS E S 2 , RMS E I , RMS E T , RMS E R = 1 m i = 1 m ( S 1 ) i ( S ^ 1 ) i 2 , 1 m i = 1 m ( S 1 ) i ( S ^ 1 ) i 2 , 1 m i = 1 m I i I ^ i 2 , 1 m i = 1 m T i T ^ i 2 , 1 m i = 1 m R i R ^ i 2 ,
TI C S 1 , TI C S 2 , TI C I , TI C T , T I C R = 1 m i = 1 m ( S 1 ) i ( S ^ 1 ) i 2 1 m i = 1 m ( S 1 ) i 2 + 1 m i = 1 m ( S ^ 1 ) i 2 , 1 m i = 1 m ( S 2 ) i ( S ^ 2 ) i 2 1 m i = 1 m ( S 2 ) i 2 + 1 m i = 1 m ( S ^ 2 ) i 2 1 m i = 1 m I i I ^ i 2 1 m i = 1 m I i 2 + 1 m i = 1 m I ^ i 2 , 1 m i = 1 m T i T ^ i 2 1 m i = 1 m T i 2 + 1 m i = 1 m T ^ i 2 1 m i = 1 m R i R ^ i 2 1 m i = 1 m R i 2 + 1 m i = 1 m R ^ i 2 ,

4. Results and Discussion

The detailed description of the results with necessary interpretations are provided here to solve the SITR system for COVID-19 dynamics, as given in the system of Equation (1). The numerical results are also determined with Adams method, which are used in this work as reference solutions for comparison with the designed FF-ANN-GASQP solver. Moreover, statistical observations are also listed to evaluate the performance. The stability and reliability are verified by the results obtained through 30 independent runs achieving a higher level of accuracy by means of statistical measures based on TIC and RMSE, which further validate the worth of the proposed FF-ANN-GASQP solver. The convergence analysis based on the performance measures of statistical operators TIC and RMSE is provided for reliability and robustness. The graphs of absolute error (AE) using the proposed results and numerical reference solutions are presented for the validity of the proposed scheme.

Nonlinear SITR Model Based on COVID-19

The different suitable values of the epidemic parameters have been given in Table 4, and the updated form of the SITR model based on COVID-19 given in system (1) using these Table 4 values is given as:
S 1 ( χ ) = 0.3 0.3 I ( χ ) S 1 ( χ ) 0.9 T ( χ ) 0.25 S 1 ( χ ) , S 1 ( 0 ) = 0.65 , S 2 ( χ ) = 0.3 0.3 I ( χ ) S 2 ( χ ) 0.9 T ( χ ) 0.25 S 2 ( χ ) , S 2 ( 0 ) = 0.15 , I ( χ ) = 0.1 I ( χ ) + 0.3 I ( χ ) S 1 ( χ ) + S 2 ( χ ) 0.25 I ( χ ) + 0.9 T ( χ ) + 0.005 I ( χ ) , I ( 0 ) = 0.75 , T ( χ ) = 0.1 I ( χ ) 0.3 T ( χ ) 0.25 T ( χ ) + 0.2 T ( χ ) + 0.1 T ( χ ) , T ( 0 ) = 0.35 , R ( χ ) = 0.25 R ( χ ) + 0.3 T ( χ ) , R ( 0 ) = 0.1 , .
The objective function of set (13) is formulated as follows:
e = 1 N m = 1 N ( S ^ 1 ) m 0.3 + 0.3 I m ( S ^ 1 ) m + 0.9 T m + 0.25 ( S ^ 1 ) m 2 + ( S ^ 2 ) m 0.3 + 0.3 I m ( S ^ 2 ) m + 0.9 T m + 0.25 ( S ^ 2 ) m 2 + I ^ m + 0.1 I m 0.3 I m ( S ^ 1 ) m + ( S ^ 2 ) m + 0.25 I m 0.9 T m 0.005 I m 2 + T ^ m 0.1 I m + 0.3 T ^ m + 0.25 T ^ m + 0.2 T ^ m 0.1 T ^ m 2 + R ^ m + 0.25 R ^ m 0.3 T ^ m 2 + 1 5 ( S ^ 1 ) 0 0.65 2 + ( S ^ 2 ) 0 0.15 2 + I ^ 0 0.75 2 + T ^ 0 0.35 2 + R ^ 0 0.1 2 .
Optimization of the SITR system for COVID-19 dynamics is performed by the hybrid-computing structure of GA and SQP for 30 numbers of independent trials to achieve the parameters of FF-ANN having five neurons. The set of the weight vectors of FF-ANN by GASQP are provided in Figure 2 for the best fitness, and these weights are applied to achieve the numerical outcomes of the FF-ANN-GASQP solver for the SITR system. Values of different parameters against four cases are given in Table 5. The mathematical formulations of the approximate numerical outcomes of FF-ANN-GASQP is given as follows:
S ^ 1 ( χ ) = 1.4704 1 + e ( 0.2502 χ 1.1644 ) 0.2502 1 + e ( 1.9550 χ + 0.4603 ) 0.3483 1 + e (   1.1161 χ 2.8067 ) + 1.5199 1 + e ( 0.2957 χ + 1.1530 ) + 0.1591 1 + e (   0.7159 χ 0.3006 ) ,
S ^ 2 ( χ ) = 1.5969 1 + e ( 0.2577 χ + 0.2176 ) 0.7521 1 + e ( 0.5993 χ + 0.8183 ) 0.5261 1 + e (   0.7527 χ 0.7915 ) 1.1453 1 + e ( 0.7548 χ 2.4308 ) + 0.1292 1 + e (   0.7124 χ 0.6757 ) ,
I ^ ( χ ) = 0.00001 1 + e ( 0.4977 χ 1.0265 ) + 0.0815 1 + e ( 0.9180 χ + 0.0571 ) + 1.3974 1 + e (   2.0394 χ + 1.0937 ) + 0.1292 1 + e ( 0.2777 χ 1.2619 ) 1.0806 1 + e (   2.1174 χ 0.6641 ) ,
T ^ ( χ ) = 2.0966 1 + e ( 0.9354 χ 0.6124 ) 0.5464 1 + e ( 1.3223 χ 0.8716 ) 0.2651 1 + e (   1.4587 χ 0.0180 ) + 0.9895 1 + e ( 1.3325 χ + 1.0869 ) + 0.8115 1 + e (   0.2659 χ + 1.3133 ) ,
R ^ ( χ ) = 0.7228 1 + e ( 1.0915 χ 0.7499 ) + 0.0451 1 + e ( 1.5696 χ 1.9934 ) 0.2232 1 + e (   0.2699 χ 1.0638 ) + 0.4927 1 + e ( 0.7613 χ + 0.7571 ) + 0.2232 1 + e (   0.2270 χ 1.2859 ) ,
Figure 2a depicts that the rate of contact is the source of a significant increase in the number of susceptible persons primarily, but with the passage of time, it starts declining. This happens because the higher contact rate causes most of the people to get infected and transfer to an infected class; therefore, the number of people in the susceptible class is reduced. Figure 2b clearly reflects that the number of infected persons increases with the increase of contact rate. Reduction in the contact rate results in a slight rise in the number of infected persons, while an increase in the contact rate suddenly enhances the number of infected persons. Figure 2c shows the growing behavior of susceptible persons with higher recovery rate values. Figure 2d depicts that as the recovery rate upsurges, so does the number of persons in the infectious class. It can be noticed that when the recovery rate is low, then less persons are recovered from the virus, because of the death of infected persons.
Figure 3a demonstrates the increasing behavior of susceptible persons with various death rates. Figure 3b shows the variation of recovered persons with different death rate values. Figure 3 depicts that with a high death rate, the number of recovered persons are reduced abruptly. This is because when the death rate is high, then a large number of persons from the infected and recovered class die, which results in the reduction of persons in all of these classes. As the death rate becomes very high, the infectious and recovered persons are approximately wiped out.
The set of Equations (15)–(19) are used to get the solution of the nonlinear SITR model based on COVID-19 as given in system (13) using the FF-ANN-GASQP, and numerical results are provided in Figure 2, Figure 3 and Figure 4 using five neurons. The trained weights set for the five neurons-based FF-ANN model is shown in Figure 4a–e for all the classes of a nonlinear SITR system for COVID-19 dynamics. The graphs of absolute error (AE) using the proposed results and numerical Adams scheme are plotted in Figure 5. It is seen that most of the AE values lie around 10−4 to 10−5 for all the parameters of the nonlinear SITR model using five numbers of neurons.
The convergence analysis based on the performance measures of the statistical operators TIC and RMSE values is drawn in Figure 6. The first part of Figure 6 shows the performance measures of the TIC values, while the second part of Figure 6 indicates the performance of the RMSE values for all the parameters of the nonlinear SITR model of COVID-19. One may conclude from these obtainable results that most of independent runs achieved a higher level of accuracy.

5. Conclusions

A new computational intelligent solver FF-ANN-GASQP is designed for the solution of a nonlinear SITR system for COVID-19 dynamics using the approximation capability of FF-ANN optimized globally with GAs and locally with SQP. The dynamics of the nonlinear SITR model of COVID-19 are proficiently assessed by the designed FF-ANN-GASQP approach with the single hidden layer formulation of FF-ANNs using five neurons. The exactness of the designed FF-ANN-GASQP technique is endorsed by matching of the outcomes with reasonable precision from the reference Adams numerical scheme to solve a nonlinear SITR system based on COVID-19. Moreover, statistical clarifications through 30 autonomous nuns by means of TIC and RMSE indices further authenticate the worth and value of the proposed FF-ANN-GASQP solver.
In future research, the presented FF-ANN-GASQP scheme algorithm can be implemented as an efficient/accurate stochastic numerical solver for the nonlinear dynamical measures of a computational fluid mechanics model [46,47,48,49,50], renewed fractional order processing systems [51,52,53,54], higher-order singular models [55], parameter estimation problems and electromagnetic waves [56,57], a prediction differential model [58], and bioinformatics [59,60,61].

Author Contributions

Formal analysis, M.U. and Z.S.; Investigation, M.A.Z.R. and M.S.; Validation and model design, M.G. and Y.G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This paper has been partially supported by Ministerio de Ciencia, Innovación y Universidades, grant number PGC2018-097198-B-I00, and by Fundación Séneca of Región de Murcia, grant number 20783/PI/18. These two grants were received by Yolanda Guerrero. No other author obtained any grant.

Conflicts of Interest

All the authors of the manuscript have no conflict of interest and have worked in an equal way to obtain the results presented in this paper.

References

  1. China Virus Death Toll Rises to 41, more than 1300 Infected Worldwide. CNBC, 24 January 2020. Available online: https://www.cnbc.com/2020/01/24/chinas-hubei-province-confirms-15-more-deaths-due-to-coronavirus.html (accessed on 10 June 2020).
  2. Khan, M.A.; Atangana, A. Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative. Alex. Eng. J. 2020. [Google Scholar] [CrossRef]
  3. Sabir, Z.; Umar, M.; Guirao, J.L.; Shoaib, M.; Raja, M.A.Z. Integrated intelligent computing paradigm for nonlinear multi-singular third-order Emden–Fowler equation. Neural Comput. Appl. 2020, 1–20. [Google Scholar] [CrossRef]
  4. Zúñiga-Aguilar, C.; Ugalde, H.M.R.; Gómez-Aguilar, J.F.; Escobar-Jiménez, R.F.; Valtierra-Rodríguez, M. Solving fractional differential equations of variable-order involving operators with Mittag-Leffler kernel using artificial neural networks. Chaos Solitons Fractals 2017, 103, 382–403. [Google Scholar] [CrossRef]
  5. Sabir, Z.; Saoud, S.; Raja, M.A.Z.; Wahab, H.A.; Arbi, A. Heuristic computing technique for numerical solutions of nonlinear fourth order Emden–Fowler equation. Math. Comput. Simul. 2020, 178, 534–548. [Google Scholar] [CrossRef]
  6. Sabir, Z.; Wahab, H.A.; Umar, M.; Sakar, M.G.; Raja, M.A.Z. Novel design of Morlet wavelet neural network for solving second order Lane–Emden equation. Math. Comput. Simul. 2020, 172, 1–14. [Google Scholar] [CrossRef]
  7. Zuñiga-Aguilar, C.J.; Coronel-Escamilla, A.; Gómez-Aguilar, J.F.; Alvarado-Martínez, V.M.; Ugalde, H.M.R. New numerical approximation for solving fractional delay differential equations of variable order using artificial neural networks. Eur. Phys. J. Plus 2018, 133, 75. [Google Scholar] [CrossRef]
  8. Sabir, Z.; Baleanu, D.; Shoaib, M.; Raja, M.A.Z. Design of stochastic numerical solver for the solution of singular three-point second-order boundary value problems. Neural Comput. Appl. 2020, 1–17. [Google Scholar] [CrossRef]
  9. Bukhari, A.H.; Raja, M.A.Z.; Sulaiman, M.; Islam, S.; Shoaib, M.; Kumam, P. Fractional Neuro-Sequential ARFIMA-LSTM for Financial Market Forecasting. IEEE Access 2020, 8, 71326–71338. [Google Scholar] [CrossRef]
  10. Sabir, Z.; Manzar, M.A.; Raja, M.A.Z.; Sheraz, M.; Wazwaz, A.-M. Neuro-heuristics for nonlinear singular Thomas-Fermi systems. Appl. Soft Comput. 2018, 65, 152–169. [Google Scholar] [CrossRef]
  11. Bukhari, A.H.; Sulaiman, M.; Islam, S.; Shoaib, M.; Kumam, P.; Raja, M.A.Z. Neuro-fuzzy modeling and prediction of summer precipitation with application to different meteorological stations. Alex. Eng. J. 2020, 59, 101–116. [Google Scholar] [CrossRef]
  12. Umar, M.; Sabir, Z.; Raja, M.A.Z. Intelligent computing for numerical treatment of nonlinear prey–predator models. Appl. Soft Comput. 2019, 80, 506–524. [Google Scholar] [CrossRef]
  13. Hassan, A.; Ahmad, S.-U.-I.; Kamran, M.; Illahi, A.; Zahoor, R.M.A. Design of cascade artificial neural networks optimized with the memetic computing paradigm for solving the nonlinear Bratu system. Eur. Phys. J. Plus 2019, 134, 122. [Google Scholar] [CrossRef]
  14. Majeed, K.; Masood, Z.; Samar, R.; Raja, M.A.Z. A genetic algorithm optimized Morlet wavelet artificial neural network to study the dynamics of nonlinear Troesch’s system. Appl. Soft Comput. 2017, 56, 420–435. [Google Scholar] [CrossRef]
  15. Mehmood, A.; Zameer, A.; Ling, S.H.; Raja, M.A.Z. Design of neuro-computing paradigms for nonlinear nanofluidic systems of MHD Jeffery–Hamel flow. J. Taiwan Inst. Chem. Eng. 2018, 91, 57–85. [Google Scholar] [CrossRef]
  16. Sabir, Z.; Wahab, H.A.; Umar, M.; Erdogan, F. Stochastic numerical approach for solving second order nonlinear singular functional differential equation. Appl. Math. Comput. 2019, 363, 124605. [Google Scholar] [CrossRef]
  17. Sabir, Z.; Raja, M.A.Z.; Umar, M. Neuro-swarm intelligent computing to solve the second-order singular functional differential model. Eur. Phys. J. Plus. 2020, 135, 474. [Google Scholar] [CrossRef]
  18. Raja, M.A.Z.; Shah, F.H.; Alaidarous, E.S.; Syam, M.I. Design of bio-inspired heuristic technique integrated with interior-point algorithm to analyze the dynamics of heartbeat model. Appl. Soft Comput. 2017, 52, 605–629. [Google Scholar] [CrossRef]
  19. Mehmood, A.; Zameer, A.; Raja, M.A.Z. Intelligent computing to analyze the dynamics of Magnetohydrodynamic flow over stretchable rotating disk model. Appl. Soft Comput. 2018, 67, 8–28. [Google Scholar] [CrossRef]
  20. Umar, M.; Amin, F.; Wahab, H.A.; Baleanu, D. Unsupervised constrained neural network modeling of boundary value corneal model for eye surgery. Appl. Soft Comput. 2019, 85, 105826. [Google Scholar] [CrossRef]
  21. Mehmood, A.; Zameer, A.; Ling, S.H.; Rehman, A.U.; Raja, M.A.Z. Integrated computational intelligent paradigm for nonlinear electric circuit models using neural networks, genetic algorithms and sequential quadratic programming. Neural Comput. Appl. 2019, 1–21. [Google Scholar] [CrossRef]
  22. Raja, M.A.Z.; Mehmood, J.; Sabir, Z.; Nasab, A.K.; Manzar, M.A. Numerical solution of doubly singular nonlinear systems using neural networks-based integrated intelligent computing. Neural Comput. Appl. 2017, 31, 793–812. [Google Scholar] [CrossRef]
  23. Ahmad, I.; Ilyas, H.; Urooj, A.; Aslam, M.S.; Shoaib, M.; Raja, M.A.Z. Novel applications of intelligent computing paradigms for the analysis of nonlinear reactive transport model of the fluid in soft tissues and microvessels. Neural Comput. Appl. 2019, 31, 9041–9059. [Google Scholar] [CrossRef]
  24. Sabir, Z.; Raja, M.A.Z.; Guirao, J.L.G.; Shoaib, M. A Neuro-Swarming Intelligence-Based Computing for Second Order Singular Periodic Non-linear Boundary Value Problems. Front. Phys. 2020, 8. [Google Scholar] [CrossRef]
  25. Raja, M.A.Z.; Umar, M.; Sabir, Z.; Khan, J.A.; Baleanu, D. A new stochastic computing paradigm for the dynamics of nonlinear singular heat conduction model of the human head. Eur. Phys. J. Plus 2018, 133, 364. [Google Scholar] [CrossRef]
  26. Ahmad, S.U.I.; Faisal, F.; Shoaib, M.; Raja, M.A.Z. A new heuristic computational solver for nonlinear singular Thomas–Fermi system using evolutionary optimized cubic splines. Eur. Phys. J. Plus 2020, 135, 55. [Google Scholar] [CrossRef]
  27. Raja, M.A.Z.; Manzar, M.A.; Samar, R. An efficient computational intelligence approach for solving fractional order Riccati equations using ANN and SQP. Appl. Math. Model. 2015, 39, 3075–3093. [Google Scholar] [CrossRef]
  28. Sabir, Z.; Raja, M.A.Z.; Umar, M.; Shoaib, M. Design of neuro-swarming-based heuristics to solve the third-order nonlinear multi-singular Emden–Fowler equation. Eur. Phys. J. Plus 2020, 135, 410. [Google Scholar] [CrossRef]
  29. Lodhi, S.; Manzar, M.A.; Raja, M.A.Z. Fractional neural network models for nonlinear Riccati systems. Neural Comput. Appl. 2017, 31, 359–378. [Google Scholar] [CrossRef]
  30. Raja, M.A.Z.; Khan, J.A.; Bilal, M. A new stochastic approach for solution of Riccati differential equation of fractional order. Ann. Math. Artif. Intell. 2010, 60, 229–250. [Google Scholar] [CrossRef]
  31. Raja, M.A.Z.; Samar, R.; Manzar, M.A.; Shah, S.M. Design of unsupervised fractional neural network model optimized with interior point algorithm for solving Bagley–Torvik equation. Math. Comput. Simul. 2017, 132, 139–158. [Google Scholar] [CrossRef]
  32. Khan, N.; Hameed, T.; Khan, N.; Raja, M.A.Z. A heuristic optimization method of fractional convection reaction: An application to diffusion process. Therm. Sci. 2018, 22, 243–252. [Google Scholar] [CrossRef] [Green Version]
  33. Ara, A.; Alam Khan, N.; Razzaq, O.A.; Hameed, T.; Raja, M.A.Z. Wavelets optimization method for evaluation of fractional partial differential equations: An application to financial modelling. Adv. Differ. Equ. 2018, 2018, 8. [Google Scholar] [CrossRef]
  34. Sanchez, Y.G.; Sabir, Z.; Guirao, J.L.G. Design of a Nonlinear Sitr Fractal Model Based on the Dynamics of a Novel Coronavirus (COVID-19). Fractals 2020. [Google Scholar] [CrossRef]
  35. Srinivas, N.; Deb, K. Muilti-objective optimization using no dominated sorting in genetic algorithms. Evolut. Comput. 1994, 2, 221–248. [Google Scholar] [CrossRef]
  36. Holland, J.H. Genetic Algorithms. Sci. Am. 1992, 267, 66–72. [Google Scholar] [CrossRef]
  37. Zameer, A.; Majeed, M.; Mirza, S.M.; Raja, M.A.Z.; Khan, A.; Mirza, N.M. Bio-inspired heuristics for layer thickness optimization in multilayer piezoelectric transducer for broadband structures. Soft Comput. 2018, 23, 3449–3463. [Google Scholar] [CrossRef]
  38. Raja, M.A.Z.; Shah, A.A.; Mehmood, A.; Chaudhary, N.I.; Aslam, M.S. Bio-inspired computational heuristics for parameter estimation of nonlinear Hammerstein controlled autoregressive system. Neural Comput. Appl. 2016, 29, 1455–1474. [Google Scholar] [CrossRef]
  39. Tan, J.; Kerr, W.L. Determination of glass transitions in boiled candies by capacitance based thermal analysis (CTA) and genetic algorithm (GA). J. Food Eng. 2017, 193, 68–75. [Google Scholar] [CrossRef]
  40. Jamal, R.; Men, B.; Khan, N.H.; Raja, M.A.Z.; Men, K.R. Hybrid Bio-Inspired Computational Heuristic Paradigm for Integrated Load Dispatch Problems involving Stochastic Wind. Energies 2019, 12, 2568. [Google Scholar] [CrossRef] [Green Version]
  41. Khan, W.U.; Ye, Z.; Altaf, F.; Chaudhary, N.I.; Raja, M.A.Z. A novel application of fireworks heuristic paradigms for reliable treatment of nonlinear active noise control. Appl. Acoust. 2019, 146, 246–260. [Google Scholar] [CrossRef]
  42. Mehmood, A.; Chaudhary, N.I.; Zameer, A.; Raja, M.A.Z. Novel computing paradigms for parameter estimation in power signal models. Neural Comput. Appl. 2019, 32, 6253–6282. [Google Scholar] [CrossRef]
  43. Raja, M.A.Z. Solution of the one-dimensional Bratu equation arising in the fuel ignition model using ANN optimised with PSO and SQP. Connect. Sci. 2014, 26, 195–214. [Google Scholar] [CrossRef]
  44. Raja, M.A.Z.; Manzar, M.A.; Shah, S.M.; Chen, Y. Integrated Intelligence of Fractional Neural Networks and Sequential Quadratic Programming for Bagley–Torvik Systems Arising in Fluid Mechanics. J. Comput. Nonlinear Dyn. 2020, 15. [Google Scholar] [CrossRef]
  45. Khan, W.U.; Ye, Z.; Chaudhary, N.I.; Raja, M.A.Z. Backtracking search integrated with sequential quadratic programming for nonlinear active noise control systems. Appl. Soft Comput. 2018, 73, 666–683. [Google Scholar] [CrossRef]
  46. Umar, M.; Sabir, Z.; Imran, A.; Wahab, A.H.; Shoaib, M.; Raja, M.A.Z. The 3-D flow of Casson nanofluid over a stretched sheet with chemical reactions, velocity slip, thermal radiation and Brownian motion. Therm. Sci. 2020, 24, 2929–2939. [Google Scholar] [CrossRef] [Green Version]
  47. Sabir, Z.; Imran, A.; Umar, M.; Zeb, M.; Shoaib, M.; Raja, M.A.Z. A numerical approach for two-dimensional Sutterby fluid flow bounded at a stagnation point with an inclined magnetic field and thermal radiation impacts. Therm. Sci. 2020, 186. [Google Scholar] [CrossRef]
  48. Sabir, Z.; Akhtar, R.; Zhiyu, Z.; Umar, M.; Imran, A.; Wahab, H.A.; Shoaib, M.; Raja, M.A.Z. A Computational Analysis of Two-Phase Casson Nanofluid Passing a Stretching Sheet Using Chemical Reactions and Gyrotactic Microorganisms. Math. Probl. Eng. 2019, 2019, 1–12. [Google Scholar] [CrossRef] [Green Version]
  49. Sabir, Z.; Ayub, A.; Guirao, J.L.; Bhatti, S.; Shah, S.Z.H. The Effects of Activation Energy and Thermophoretic Diffusion of Nanoparticles on Steady Micropolar Fluid along with Brownian Motion. Adv. Mater. Sci. Eng. 2020, 2020, 1–12. [Google Scholar] [CrossRef]
  50. Umar, M.; Akhtar, R.; Sabir, Z.; Wahab, H.A.; Zhiyu, Z.; Imran, A.; Shoaib, M.; Raja, M.A.Z. Numerical Treatment for the Three-Dimensional Eyring-Powell Fluid Flow over a Stretching Sheet with Velocity Slip and Activation Energy. Adv. Math. Phys. 2019, 2019, 1–12. [Google Scholar] [CrossRef] [Green Version]
  51. Chaudhary, N.I.; Zubair, S.; Aslam, M.S.; Raja, M.A.Z.; Machado, J.A.T. Design of momentum fractional LMS for Hammerstein nonlinear system identification with application to electrically stimulated muscle model. Eur. Phys. J. Plus 2019, 134, 407. [Google Scholar] [CrossRef]
  52. Raja, M.A.Z.; Akhtar, R.; Chaudhary, N.I.; Zhiyu, Z.; Khan, Q.; Rehman, A.U.; Zaman, F. A new computing paradigm for the optimization of parameters in adaptive beamforming using fractional processing. Eur. Phys. J. Plus 2019, 134, 1–13. [Google Scholar] [CrossRef]
  53. Masood, Z.; Samar, R.; Raja, M.A.Z. Design of fractional order epidemic model for future generation tiny hardware implants. Futur. Gener. Comput. Syst. 2020, 106, 43–54. [Google Scholar] [CrossRef]
  54. Muhammad, Y.; Khan, R.; Raja, M.A.Z.; Ullah, F.; Chaudhary, N.I.; He, Y. Design of Fractional Swarm Intelligent Computing With Entropy Evolution for Optimal Power Flow Problems. IEEE Access 2020, 8, 111401–111419. [Google Scholar] [CrossRef]
  55. Sabir, Z.; Günerhan, H.; Guirao, J.L. On a New Model Based on Third-Order Nonlinear Multisingular Functional Differential Equations. Math. Probl. Eng. 2020, 2020, 1–9. [Google Scholar] [CrossRef] [Green Version]
  56. Akbar, S.; Zaman, F.; Asif, M.; Rehman, A.U.; Raja, M.A.Z. Novel application of FO-DPSO for 2-D parameter estimation of electromagnetic plane waves. Neural Comput. Appl. 2018, 31, 3681–3690. [Google Scholar] [CrossRef]
  57. Zaman, F.; Hassan, A.; Akbar, S.; Rehman, A.U.; Raja, M.A.Z. Novel Computational Heuristics for Wireless Parameters Estimation in Bistatic Radar systems. Wirel. Pers. Commun. 2019, 111, 909–927. [Google Scholar] [CrossRef]
  58. Sabir, Z.; Guirao, J.L.; Saeed, T.; Erdoğan, F. Design of a Novel Second-Order Prediction Differential Model Solved by Using Adams and Explicit Runge–Kutta Numerical Methods. Math. Probl. Eng. 2020, 2020, 9704968. [Google Scholar] [CrossRef]
  59. Umar, M.; Sabir, Z.; Amin, F.; Guirao, J.L.; Raja, M.A.Z. Stochastic numerical technique for solving HIV infection model of CD4+ T cells. Eur. Phys. J. Plus 2020, 135, 1–19. [Google Scholar] [CrossRef]
  60. Umar, M.; Raja, M.A.Z.; Sabir, Z.; Alwabli, A.S.; Shoaib, M. A stochastic computational intelligent solver for numerical treatment of mosquito dispersal model in a heterogeneous environment. Eur. Phys. J. Plus 2020, 135, 1–23. [Google Scholar] [CrossRef]
  61. Guerrero-Sánchez, Y.; Sabir, Z.; Günerhan, H.; Baskonus, H.M. Analytical and Approximate Solutions of a Novel Nervous Stomach Mathematical Model. Discret. Dyn. Nat. Soc. 2020, 2020, 5063271. [Google Scholar] [CrossRef]
Figure 1. Process workflow of proposed FF-ANN-GASQP structure for SITR model based on COVID-19. FF-ANN-GASQP: feed-forward artificial neural networks trained with global search genetic algorithms and speedy fine tuning by sequential quadratic programming.
Figure 1. Process workflow of proposed FF-ANN-GASQP structure for SITR model based on COVID-19. FF-ANN-GASQP: feed-forward artificial neural networks trained with global search genetic algorithms and speedy fine tuning by sequential quadratic programming.
Symmetry 12 01628 g001
Figure 2. Impact of contact rate and recovery rate on susceptible and infected classes.
Figure 2. Impact of contact rate and recovery rate on susceptible and infected classes.
Symmetry 12 01628 g002aSymmetry 12 01628 g002b
Figure 3. Impact of death rate on susceptible and recovered persons classes.
Figure 3. Impact of death rate on susceptible and recovered persons classes.
Symmetry 12 01628 g003
Figure 4. Set of best fitted trained weights/decision variables of FF-ANN-GASQP.
Figure 4. Set of best fitted trained weights/decision variables of FF-ANN-GASQP.
Symmetry 12 01628 g004aSymmetry 12 01628 g004b
Figure 5. Absolute error (AE) results for all the variables of the nonlinear SITR model based on COVID-19 using five neurons.
Figure 5. Absolute error (AE) results for all the variables of the nonlinear SITR model based on COVID-19 using five neurons.
Symmetry 12 01628 g005
Figure 6. Convergence analysis based on the performance measures of Theil’s inequality coefficient (TIC) and root mean square error (RMSE) values for all the variables of the nonlinear SITR model.
Figure 6. Convergence analysis based on the performance measures of Theil’s inequality coefficient (TIC) and root mean square error (RMSE) values for all the variables of the nonlinear SITR model.
Symmetry 12 01628 g006
Table 1. State variables for the dynamics of the SITR model.
Table 1. State variables for the dynamics of the SITR model.
VariableDescription
S 1 ( χ ) Non-infected individuals
S 2 ( χ ) Non-infected older or major diseased people
I ( χ ) Rate of infected from COVID-19
R ( χ ) Recovery rate from COVID-19
T ( χ ) Treatment
Table 2. Descriptions of the state variables based on the nonlinear SITR model.
Table 2. Descriptions of the state variables based on the nonlinear SITR model.
ParameterDescription
β Contact rate
B Rate of natural birth
δ Reduce infection from treatment
σ Fever, tiredness and dry cough rate
μ Recovery rate
α Death rate
ρ Rate of infection from treatment
ψ Healthy food rate
ε Sleep rate
A j , j = 1 , 2 , 3 , 4 , 5 Initial conditions
Table 3. Pseudocode of FF-ANN-GASQP for nonlinear SITR system for COVID-19 dynamics.
Table 3. Pseudocode of FF-ANN-GASQP for nonlinear SITR system for COVID-19 dynamics.
Start the GA process
Inputs: The individuals with genes equally representing the decision values of FF-ANN as:
   W = [ W S 1 , W S 2 , W I , W T , W R ] , where W S 1 = [ ϕ S 1 , w S 1 , b S 1 ] , W S 2 = [ ϕ S 2 , w S 2 , b S 2 ] ,
   W I = [ ϕ I , w I , b I ] , W T = [ ϕ T , w T , b T ] and W R = [ ϕ R , w R , b R ] as per the details provided
  in the system (3).
  Population: Number of chromosomes in a set define a population as:
   P = [ W 1 , W 2 , W 3 , , W n ] t , for ith component
   W i = [ W S 1 , i , W S 2 , i , W I , i , W T , i , W R , i ] with
   W S 1 , i = [ ϕ S 1 , i , w S 1 , i , b S 1 , i ] , W S 2 , i = [ ϕ S 2 , i , w S 2 , i , b S 2 , i ] , W I , i = [ ϕ I , i , ω I , i , b I , i ] ,
   W T , i = [ ϕ T , i , w T , i , b T , i ] , and W R , i = [ ϕ R , i , w R , i , b R , i ]
Output: The best global decision variables/trained weights of the ANN-GASQ programming scheme denoted as WGA-Best
Initialization: Generate chromosome W and P with pseudo random numbers.
  Initialization is performed for {GA} and {gaoptimset} routines with
  suitable declarations and settings.
Fitness evaluated: Calculate the fitness and its parts shown in Equations (5)–(10) for each {W} in P.
Termination: Terminate the procedure, when any requirements meets
   • {Achieved Fitness = 10−20},{Generations = 60}
   • {Tolerances: {TolFun= 10−20 and TolCon =10−21},
   • {StallGenLimit=100},{Population size = 300}
   • Others values: default.
  When the above conditions meet, go to storage
Ranking: Rank is proficient for every ‘W’ of ‘P’ indicates the attained
  fitness.
Reproduction:
   • {Selection: selectionuniform}
   • {Mutations: mutationadaptfeasible}.
   • {Crossover: crossoverheuristic}.
   • {Elitism: Transmit 5% individuals in P}
   Go to fitness assessment step.
Storage: Store the WGA-Best, i.e., the weight vector, fitness assessment,
  generations, time and count of functions for the present run of
  GAs.
End of GA
SQP Process Start
Inputs: Start point is WGA-Best
Output: GASQP best weights are denoted as WGASQ
Initialize: Set the limited constraints, iterations and other values of
  optimset.
Terminate: The SQ programming process terminates when one the criteria
  meets
  {Generations = 900}, {Fitness = 10−18},
  (TolFun = TolCon = TolX = 10−22) and {MaxFunEvals = 285000}.
While (Terminate)
  Fitness Evaluate: Compute fitness value of every W of P by using
  system (4) to (10).
  Adjustments: Fine-tune {fmincon} with SQ programming scheme to tune
  W and adjust again the fitness value by using systems (4) to (10).
  Accumulate Store fitness, time, WGASQ, function counts and generations for multiple trials of SQ programming.
End of the SQ programming scheme
Data Generations
  Repeat the procedure 30 times based on the GASQ programming to get a
  massive dataset of the optimization variables of ANNs for numerical
  solutions of the SITR model based on COVID-19
Table 4. Description and values for various involved parameters in SITR model.
Table 4. Description and values for various involved parameters in SITR model.
SymbolParameter DescriptionAssigned Value
β Contact rate0.3
B Natural birth rate0.3
δ Reduce infection from the treatment0.3
σ Fever, tiredness, and dry cough rate0.005
μ Recovery rate0.1
α Death rate0.25
ρ Rate of infection from the treatment0.3
ψ Healthy food rate0.2
ε Sleep rate0.1
Table 5. Variation of different involved parameters.
Table 5. Variation of different involved parameters.
ScenariosVariable ParameterCase ICase IICase IIICase IV
1Contact Rate β = 0.25 β = 0.30 β = 0.35 β = 0.40
2Recovery Rate μ = 0.08 μ = 0.10 μ = 0.12 μ = 0.14
3Death Rate α = 0.20 α = 0.25 α = 0.30 α = 0.35

Share and Cite

MDPI and ACS Style

Umar, M.; Sabir, Z.; Raja, M.A.Z.; Shoaib, M.; Gupta, M.; Sánchez, Y.G. A Stochastic Intelligent Computing with Neuro-Evolution Heuristics for Nonlinear SITR System of Novel COVID-19 Dynamics. Symmetry 2020, 12, 1628. https://doi.org/10.3390/sym12101628

AMA Style

Umar M, Sabir Z, Raja MAZ, Shoaib M, Gupta M, Sánchez YG. A Stochastic Intelligent Computing with Neuro-Evolution Heuristics for Nonlinear SITR System of Novel COVID-19 Dynamics. Symmetry. 2020; 12(10):1628. https://doi.org/10.3390/sym12101628

Chicago/Turabian Style

Umar, Muhammad, Zulqurnain Sabir, Muhammad Asif Zahoor Raja, Muhammad Shoaib, Manoj Gupta, and Yolanda Guerrero Sánchez. 2020. "A Stochastic Intelligent Computing with Neuro-Evolution Heuristics for Nonlinear SITR System of Novel COVID-19 Dynamics" Symmetry 12, no. 10: 1628. https://doi.org/10.3390/sym12101628

APA Style

Umar, M., Sabir, Z., Raja, M. A. Z., Shoaib, M., Gupta, M., & Sánchez, Y. G. (2020). A Stochastic Intelligent Computing with Neuro-Evolution Heuristics for Nonlinear SITR System of Novel COVID-19 Dynamics. Symmetry, 12(10), 1628. https://doi.org/10.3390/sym12101628

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop