Next Article in Journal
Strong Interacting Internal Waves in Rotating Ocean: Novel Fractional Approach
Next Article in Special Issue
Mathematical Analysis of a Fractional COVID-19 Model Applied to Wuhan, Spain and Portugal
Previous Article in Journal
Classical Partition Function for Non-Relativistic Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method

by
Maryamsadat Hedayati
1,
Reza Ezzati
1,* and
Samad Noeiaghdam
2,3
1
Department of Mathematics, Karaj Branch, Islamic Azad University, Karaj 3149968111, Iran
2
Industrial Mathematics Laboratory, Baikal School of BRICS, Irkutsk National Research Technical University, 664074 Irkutsk, Russia
3
Department of Applied Mathematics and Programming, South Ural State University, Lenin Prospect 76, 454080 Chelyabinsk, Russia
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(2), 122; https://doi.org/10.3390/axioms10020122
Submission received: 11 April 2021 / Revised: 21 May 2021 / Accepted: 7 June 2021 / Published: 16 June 2021
(This article belongs to the Special Issue Mathematics of the COVID-19)

Abstract

:
Coronaviruses are a group of RNA (ribonucleic acid) viruses with the capacity for rapid mutation and recombination. Coronaviruses are known to cause respiratory or intestinal infections in humans and animals. In this paper, a biologically compatible set of nonlinear fractional differential equations governing the outbreak of the novel coronavirus is suggested based on a model previously proposed in the literature. Then, this set is numerically solved utilizing two new methods employing sine–cosine and Bernoulli wavelets and their operational matrices. Moreover, the convergence of the solution is experimentally studied. Furthermore, the accuracy of the solution is proved via comparing the results with those obtained in previous research for the primary model. Furthermore, the computational costs are compared by measuring the CPU running time. Finally, the effects of the fractional orders on the outbreak of the COVID-19 are investigated.

1. Introduction

The viruses of the family Coronaviridae (Coronaviridae is a family of enveloped, positive-strand RNA viruses which infect amphibians, birds, and mammals) are ubiquitous in nature due to their existence in a wide spectrum of mammals and avian species, they cause respiratory or intestinal infections that can range from mild to lethal [1]. The elders and those with underlying medical problems including cardiovascular disease, diabetes, chronic respiratory disease, and cancer are more prone to develop a serious illness. The best way to prevent and diminish transmission is to become well-informed about COVID-19, how it is caused and how it spreads. It has been frequently recommended that an individual protect themselves and others from infection by sanitizing their hands and not touching their faces [2]. In December 2019, the World Health Organization (WHO) based in China identified the city of Wuhan, Hubei Province, as the origin of the novel coronavirus (2019-nCoV) outbreak (known as Acute Respiratory Syndrome) [3]. The virus spread rapidly around the world in 2020 and severely affected human life. It is reported that the virus might be bat origin, and the transmission of the virus might be related to a seafood market (Huanan Seafood Wholesale Market) exposure. Hence, it became of great importance to model the mechanism of the virus transmission from its source to the people. This drew some researchers’ attention to conducting studies focused on simulating the outbreak of the COVID-19.
According to the importance of mathematical modelling, Chen et al. [4] and Khan and Atangana [5] proposed the coronavirus models independently. Chen et al. [4] considered a transmission network consisting of four groups including bats as the source, hosts (probably some wild animals), seafood market considered as the reservoir of the virus and finally the people exposed to the market and consequently facing the risk of getting infected. In each group, the number of susceptible members, those exposed to the virus, infected ones as well as the number of members removed from the system either by death or by recovering from infection were considered. In addition, the birth and death rates along with the incubation and infection periods were considered. This model was finally presented as a set of nonlinear differential equations as follows:
d S p d t = Λ p m p S p β p S p I p + k A p β w S p W , d E p d t = β p S p I p + k A p + β w S p W 1 δ p ω p E p ω p E p m p E p , d I p d t = 1 δ p ω p E p γ p + m p I p , d A p d t = δ p ω p E p γ p + m p A p , d R p d t = γ p I p + γ p A p m p R p , d W p d t = μ p I p + μ p A p ε W .
where S p ( t ) , E p ( t ) , I p ( t ) , A p ( t ) , R p ( t ) and W ( t ) are, respectively, susceptible people, exposed people, symptomatic infected people, asymptomatic infected people and COVID-19 in reservoir in time t . Furthermore, Λ p = n p × N p where N p refers to the total number of people. Other parameters are listed and described in Table 1.
On the other hand, widespread attempts to simulate the outbreak of COVID-19 have inspired numerous applied mathematicians to present novel methods and solutions of these models. In [6], integer-order temporal derivatives of a set of differential equations, governing the transmission of the coronavirus, were replaced with fractional order derivatives in order that the memory effects be included in simulation resulting in effective modelling of the epidemic diseases. Furthermore, a new approach to solve uncertain SIR model was proposed in [7]. Then, parameter estimation and numerical solution of these equations were provided to be applied on development trend of COVID-19. Moreover, alternative Legendre polynomials (ALPs) were employed in [8] to solve the model introduced in [4] with integer-order temporal derivatives.
Since the biological systems have memory, their response at a certain moment depends on the history of their responses. Thus, the effect of the memory should be simulated via using temporal fractional derivatives instead of integer-order ones [9]. Hence, temporal derivatives have been replaced with fractional ones in several real-world modelling specially in medical applications. For instance, a competition model of breast cancer based on fractional differential equations was presented and numerically solved in [10]. Furthermore, fractional derivatives have been employed in various articles studying COVID-19 outbreak [11,12].
In recent years, remarkable research have been conducted in nonlinear fractional differential equations. For instance, the modified Riemann–Liouville fractional derivative rule and two kinds of fractional dual-function methods were combined to solve nonlinear fractional models [13]. Furthermore, the fractional mapping and fractional bi-function methods were employed to investigate nonlinear fractional partial differential Schrödinger equation [14]. Furthermore, the Jacobian elliptic function expansion method was modified by considering conformable fractional derivative to obtain the solution of conformable fractional discrete complex cubic Ginzburg–Landau equation [15]. In [16], fractional F-expansion method with the Hermit transformation were utilized to solve the wick-type stochastic fractional nonlinear Schrödinger equation. In addition, Hermit transformation, modified Riemann–Liouville fractional derivative rule and fractional mapping method were combined to analyze stochastic fractional solutions of a wick-type SFNLSE [17]. More mathematical modelling in medical sciences can be found in [18,19,20].
In this paper, integer-order derivatives in Equation (1) are replaced with Caputo fractional derivatives so that the obtained set of differential equations include memory effects and become biologically compatible with real-world evidences [21]:
D * υ 1 S p = Λ p m p S p β p S p ( I p + k A p ) β W S p W , D * υ 2 E p = β p S p ( I p + k A p ) + β W S P W ( 1 δ p ) ω p E p δ p ω p E p m p E p , D * υ 3 I p = ( 1 δ p ) ω p E p ( γ p + m p ) I p , D * υ 4 A p = δ p ω p E p ( γ p + m p ) A p , D * υ 5 R p = γ p I p + γ p A p m p R p , D * υ 6 W = μ p I p + μ p A p ε W .
where 0 < υ 1 , υ 2 , υ 3 , υ 4 , υ 5 , υ 6 1 and D * signifies the Caputo fractional derivative. Therefore, Equation (1) is a special case of Equation (2). In other words, in case the fractional orders equal to 1, Equation (2) changes to Equation (1).
Moreover, the aim of this paper is to solve nonlinear fractional differential Equation (2) using a collocation method based on the sine–cosine and Bernoulli wavelets.
This model of COVID-19 transmission can be employed to simulate any epidemic in a local region including any arbitrary country or state. In order to attain this goal, all parameters of the model including the fractional orders must be calculated in a way that the results of the simulation become best fitted to the real-world data. This necessitates employing a powerful numerical method as well as considering the effects of all parameters specifically fractional orders. Therefore, in this paper, the sine–cosine and Bernoulli wavelets are examined for rate of convergence, accuracy and computational cost to determine more powerful scheme for possible future modelling of any outbreak using this model. Thus, an experimental convergence analysis is conducted. In order to demonstrate the accuracy of the results and the capability of the methods, the results for the special case, Equation (1), are obtained and compared with those reported in [8] and with the solution attained using RK4 method. In addition, the CPU running time of these two methods are measured and compared with each other to evaluate the computational costs. Finally, the effects of the fractional orders on the COVID-19 spread are investigated.
The remainder of the paper has been organized as follows. In Section 2, some preliminaries and definitions related to the fractional calculus are presented. In Section 3, the properties of wavelets are briefly reviewed. This section contains two subsections presenting properties of sine–cosine and Bernoulli wavelets and their operational matrices for fractional integration as well as function approximations. In Section 4, numerical solutions of Equation (2) using the above-mentioned wavelets are presented. In Section 5, the results are presented and an extensive discussion on the convergence, accuracy of the solutions, computational costs as well as the effects of the fractional orders on the field variables are provided. Finally, in Section 6, concluding remarks are highlighted.

2. Preliminaries

Definition 1.
[22]. The Riemann-Liouville fractional integral operator of order α is defined as
I α u ( t ) = 1 Γ ( α ) 0 t ( t τ ) α 1 u ( τ ) d τ , α > 0 , u ( t ) = 0 , α = 0 ,
Definition 2.
[23]. The Caputo fractional derivative operator of order α is defined as
D α u ( t ) = 1 Γ ( n α ) 0 t ( t τ ) n α 1 u ( n ) ( τ ) d τ , α > 0 , n 1 < α < n , d ( n ) u ( t ) d t n , α = n ,
where t 0 and n .
The Riemann-Liouville integral operator and Caputo derivative operators have the following properties:
D α I α u ( t ) = u ( t )
and
I α D α u ( t ) = u ( t ) k = 0 n 1 u ( k ) ( 0 + ) k ! t k ,       t 0 , n 1 < α < n .

3. A Brief Review of Wavelets

Wavelet analysis is a new development in applied mathematics. Wavelets are special functions which exhibits oscillatory behavior in a short period and then vanishes.
Wavelets constitute a family of functions constructed from dilation and translation of a single function ψ ( x ) called the mother wavelet. When the dilation parameter α and the translation parameter b vary continuously we have the following family of continuous wavelets [24]:
ψ a , b ( t ) = a 1 2 ψ t b a ,     a , b , a 0 .
If we restrict the parameters a and b to discrete values by a = a 0 k , b = n b 0 a 0 k , a 0 > 1 , b 0 > 0 , we have the following family of discrete wavelets:
ψ k , n ( t ) = a 0 k 2 ψ ( a 0 k t n b 0 ) ,       k , n + ,
These wavelets for all integers k and n produce an orthogonal basis of L 2 ( ) . Specially, when a 0 = 2 and b 0 = 1 then ψ k , n ( t ) form an orthonormal basis.

3.1. Sine–Cosine Wavelets and Their Properties

3.1.1. Sine–Cosine Wavelets

Sine–cosine wavelets ψ n , m ( t ) = ψ ( k , n , m , t ) involve four arguments; k = 0 , 1 , 2 , , n = 0 , , 2 k 1 1 , the values of m are given in Equation (10) and t is the normalized time. They are defined on the interval 0 , 1 as follows [25,26]:
ψ n , m ( t ) = 2 k + 1 2 f m ( 2 k t n ) , n 2 k t < n + 1 2 k , 0 , e l s e w h e r e ,
where
f m ( t ) = 1 2 , m = 0 , cos ( 2 m π t ) , m = 1 , 2 , , L , sin ( 2 ( m L ) π t ) , m = L + 1 , L + 2 , , 2 L ,
In addition, L is any positive integer.
Sine–cosine wavelet have compact support, S u p p ( ψ n , m ( x ) ) = x : ψ n , m ( x ) 0 ¯ = n 1 2 k , n 2 k . The set of sine–cosine wavelet is an orthonormal set.

3.1.2. Function Approximation

Any function u ( t ) L 2 [ 0 , 1 ] can be expanded as
u ( t ) = m = 0 n = 0 c n , m ψ n , m ( t ) ,
where
c n , m = ( u ( t ) , ψ n , m ( t ) ) = 0 1 ψ n , m ( t ) u ( t ) d t ,
in which ( , ) denotes the inner product in L 2 0 , 1 . If the infinite series in Equation (11) is truncated, then it can be written as
u ( t ) m = 0 2 L n = 0 2 k 1 c n , m ψ n , m ( t ) = C T Ψ ( t ) ,
where C and Ψ ( t ) are m ^ × 1 , m ^ = 2 k ( 2 L + 1 ) , matrices given by
C = c 0 , 0 , c 0 , 1 , , c 0 , 2 L , c 1 , 0 , c 1 , 1 , , c 1 , 2 L , , c 2 k 1 , 0 , c 2 k 1 , 1 , , c 2 k 1 , 2 L T
and
Ψ ( t ) = ψ 0 , 0 ( t ) , ψ 0 , 1 ( t ) , , ψ 0 , 2 L ( t ) , ψ 1 , 0 ( t ) , ψ 1 , 1 ( t ) , , ψ 1 , 2 L ( t ) , , ψ 2 k 1 , 0 ( t ) , ψ 2 k 1 , 1 ( t ) , , ψ 2 k 1 , 2 L ( t ) T .
The collocation points for the sine–cosine wavelet are taken as t i = 2 i 1 2 m ^ , where i = 1 , 2 , , m ^ . The sine–cosine wavelet matrix Φ m ^ × m ^ is given as follows.
Φ m ^ × m ^ = Ψ 1 2 m ^ , Ψ 3 2 m ^ , , Ψ 2 m ^ 1 2 m ^ ,

3.1.3. Sine–Cosine Wavelet Operational Matrix of the Fractional Integration

In this subsection, the operational matrix of the sine–cosine wavelet for fractional integration is presented in the framework of collocation method.

Block Pulse Functions (BPFs)

The m ^ -set of block pulse functions (BPFs) on [0, T) is defined in the following form:
b i ( t ) = 1 , ( i 1 ) m ^ t < i m ^ , 0 , otherwise ,
where i = 1 , , m ^ . The functions b i ( t ) are disjoint and orthogonal. That is
b i ( t ) b j ( t ) = 0 ,       i j , b i ( t ) ,     i = j ,
0 1 b i ( τ ) b j ( τ ) d τ = 0 , i j , 1 m ^ , i = j .
Definition 3.
The tensor product of two vectors F m ^ = [ f i ] and G m ^ = [ g i ] is defined as
F G = ( f i × g i ) m ^ .
Similarly, for A = [ a i , j ] and B = [ b i , j ] , we define
A B = ( a i . j × b i , j ) m ^ × m ^ .
Lemma. 
Assume that f ( t ) and g ( t ) both are functions belonging to L 2 [ 0 , 1 ] , which can be expressed as
f ( t ) F T B m ^ ( t ) ,   g ( t ) G T B m ^ ( t ) ,
where F T = [ f 1 , f 2 , , f m ^ ] ,   G T = [ g 1 , g 2 , , g m ^ ] and B m ^ ( t ) = [ b 1 ( t ) , b 2 ( t ) , , b m ^ ( t ) ] T . Thus, from the properties of BPFs, we have
f ( t ) g ( t ) F T B m ^ ( t ) G T B m ^ ( t ) = ( F T G T ) B m ^ ( t ) ,
f ( t ) 2 ( F T B m ^ ( t ) ) 2 = ( F T ) 2 B m ^ ( t ) .
Proof. 
The proof is straightforward. □
The integration of the vector Ψ ( t ) defined in (15) can be obtained:
0 t Ψ m ^ ( τ ) d τ P Ψ m ^ ( t )
where P is the m ^ × m ^ operational matrix for integration. The sine–cosine operational matrix of integration has been derived in [27].
From the orthogonality property of BPF, the sine–cosine wavelet can be expanded into a m ^ —term BPF as
Ψ m ^ ( t ) = Φ m ^ × m ^ B m ^ ( t )
where
B m ^ ( t ) = [ b 1 ( t ) , b 2 ( t ) , , b m ^ ( t ) ] T ,
Kilicman [28] has presented the block pulse operational matrix of the fractional integration as
( I α B m ^ ) ( t ) F α B m ^ ( t ) ,
where
F α = 1 m ^ α 1 Γ ( α + 2 ) 1 ξ 1 ξ 2 ξ 3 ξ m ^ 1 0 1 ξ 1 ξ 2 ξ m ^ 2 0 0 1 ξ 1 ξ m ^ 3 0 0 0 1 ξ 1 0 0 0 0 1 ,
with ξ k = ( k + 1 ) α + 1 2 k α + 1 + ( k 1 ) α + 1 .
Next, the sine–cosine wavelet operational matrix of the fractional integration is derived as
I α Ψ m ^ ( t ) P m ^ × m ^ α Ψ m ^ ( t ) ,
where the matrix P m ^ × m ^ α is called the sine–cosine wavelet operational matrix of the fractional integration.
Using Equations (22) and (24), we have
( I α Ψ m ^ ) ( t ) ( I α Φ m ^ × m ^ B m ^ ) ( t ) = Φ m ^ × m ^ ( I α B m ^ ) ( t ) Φ m ^ × m ^ F α B m ^ ( t ) .
From Equations (26) and (27), we have
P m ^ × m ^ α Ψ m ^ ( t ) = P m ^ × m ^ α Φ m ^ × m ^ B m ^ ( t ) = Φ m ^ × m ^ F α B m ^ ( t ) .
So, with the help of (22), we get
B m ^ ( t ) = Φ m ^ × m ^ 1 Ψ m ^ ( t ) .
Consequently, from (27)–(29) the sine–cosine wavelet operational matrix of the fractional integration P m ^ × m ^ α is as follows:
P m ^ × m ^ α = Φ m ^ × m ^ F α Φ m ^ × m ^ 1 .
For k = 1 , L = 1 , α = 0.7 , the sine–cosine wavelet operational matrix of the fractional integration P m ^ × m ^ α is given by
P 6 × 6 0.7 = 0.3985 0.0119 0.1677 0.4977 0.0097 0.0529 0.0119 0.0609 0.0812 0.0097 0.0028 0.0073 0.1677 0.0812 0.0946 0.0529 0.0073 0.0246 0 0 0 0.3985 0.0119 0.1677 0 0 0 0.01193 0.0609 0.0812 0 0 0 0.1677 0.0812 0.0946

3.2. Bernoulli Wavelets and Their Properties

3.2.1. Bernoulli Wavelets

Bernoulli wavelets ψ n , m ( t ) = ψ ( k , n , m , t ) have four arguments, where n = 1 , , 2 k 1 , k can be any positive integer, m is the order for Bernoulli polynomial and t is the normalized time. We define them on the interval [ 0 , 1 ) as follows [29,30,31,32]:
ψ n , m ( t ) = 2 k 1 2 β ˜ m ( 2 k 1 t n + 1 ) , n 1 2 k 1 t < n 2 k 1 0 , o t h e r w i s e
With
β ˜ m ( t ) = 1 , m = 0 , 1 ( 1 ) m 1 ( m ! ) 2 ( 2 m ) ! α 2 m β m ( t ) , m > 0 ,
where m = 0 , 1 , , M 1 and n = 1 , 2 , , 2 k 1 . . The coefficient 1 ( 1 ) m 1 ( m ! ) 2 ( 2 m ) ! α 2 m guarantees orthonormality, the dilation parameter is a = 2 ( k 1 ) and translation parameter is b = ( n 1 ) 2 ( k 1 ) .
Here, β m ( t ) are the well-known m th order Bernoulli polynomials defined on the interval 0 , 1 by
β m ( t ) = i = 0 m m i α m i t i ,
where α i , i = 0 , 1 , , m are Bernoulli numbers. These numbers are a sequence of signed rational numbers arising in the series expansion of trigonometric functions and can be defined by the identity,
t e t 1 = i = 0 α i t i i !
The first few Bernoulli numbers are as follows.
α 0 = 1 , α 1 = 1 2 , α 2 = 1 6 , α 4 = 1 30 , α 6 = 1 42 , α 8 = 1 30 ,
with α 2 i + 1 = 0 , i = 1 , 2 , 3 , .
The first few Bernoulli polynomials are given as:
β 0 ( t ) = 1 , β 1 ( t ) = t 1 2 , β 2 ( t ) = t 2 t + 1 6 , β 3 ( t ) = t 3 3 2 t 2 + 1 2 t ,

3.2.2. Function Approximation

Suppose that ψ 10 ( t ) , ψ 11 ( t ) , , ψ 2 k 1 M 1 ( t ) L 2 [ 0 , 1 ] is the set of Bernoulli wavelets, Y = s p a n { ψ 10 ( t ) , ψ 11 ( t ) , , ψ 1 M 1 ( t ) , ψ 20 ( t ) , , ψ 2 M 1 ( t ) , , ψ 2 k 1 0 ( t ) , , ψ 2 k 1 M 1 ( t ) } , and f ( t ) be arbitrary elements of L 2 [ 0 , 1 ] . Since Y is a finite dimensional vector space, f ( t ) has the best approximation out of Y such as f 0 ( t ) Y , that is
y ( t ) Y , f ( t ) f 0 ( t ) f ( t ) y ( t ) .
Since f 0 ( t ) Y , there exists unique coefficients c 10 , c 11 , , c 2 k 1 M 1 such that
f ( t ) f 0 ( t ) = n = 1 2 k 1 m = 0 M 1 c n m ψ n m ( t ) = C T Ψ ( t )
where C and Ψ ( t ) are 2 k 1 M × 1 matrices given by
C = [ c 10 , c 11 , , c 1 M 1 , c 20 , , c 2 M 1 , , c 2 k 1 0 , , c 2 k 1 M 1 ] T ,
Ψ ( t ) = [ ψ 10 ( t ) , ψ 11 ( t ) , , ψ 1 M 1 ( t ) , ψ 20 ( t ) , , ψ 2 M 1 ( t ) , , ψ 2 k 1 0 ( t ) , , ψ 2 k 1 M 1 ( t ) ] T
Taking the collocation points as follows:
t i = 2 i 1 2 k M , i = 1 , 2 , , 2 k 1 M ,
the Bernoulli wavelet matrix Φ m ^ × m ^ is given as follows:
Φ m ^ × m ^ = Ψ 1 2 m ^ , Ψ 3 2 m ^ , , Ψ 2 m ^ 1 2 m ^ ,
where m ^ = 2 k 1 M .

3.2.3. Bernoulli Wavelet Operational Matrix of the Fractional Integration

The Bernoulli wavelet operational matrix of fractional integration P m ^ × m ^ α is obtained following the same procedure discussed in (Section 3.1.3) which leads to:
P m ^ × m ^ α = Φ m ^ × m ^ F α Φ m ^ × m ^ 1 .
In particular, for k = 2 , M = 3 , α = 0.5 , the Bernoulli wavelet operational matrix of fractional integration P m ^ × m ^ α is expressed as
P 6 × 6 0.7 = 0.528223 0.181881 0.0297821 0.443844 0.087099 0.0256378 0.14516 0.224295 0.132924 0.0798823 0.0449052 0.0198105 0.0598166 0.096441 0.168799 0.0417244 0.0001858 0.00286805 0 0 0 0.528223 0.181881 0.0297821 0 0 0 0.14516 0.224295 0.132924 0 0 0 0.0598166 0.096441 0.168799

4. Description of Numerical Method

In this section, sine–cosine and Bernoulli wavelets and their operational matrices are employed to obtain approximate solutions of Equation (2). For this aim, we rewrite the normalized Equation (2) using the following non-dimensional parameters:
s p = S p N p , e p = E p N p , i p = I p N p , a p = A p N p , r p = R p N p , w = ε W μ p N p , μ p = c μ p , b p = β p N p , b W = μ p β W N p ε .
under the initial conditions,
S p ( 0 ) = s 0 , e p ( 0 ) = e 0 , i p ( 0 ) = i 0 , a p ( 0 ) = a 0 , r p ( 0 ) = r 0 , w ( 0 ) = w 0 .
D * υ 1 s p = n p m p s p b p s p ( i p + k a p ) b W s p w , D * υ 2 e p = b p s p ( i p + k a p ) + b W s P w ( 1 δ p ) ω p e p δ p ω p e p m p e p , D * υ 3 i p = ( 1 δ p ) ω p e p ( γ p + m p ) i p , D * υ 4 a p = δ p ω p e p ( γ p + m p ) a p , D * υ 5 r p = γ p i p + γ p a p m p r p , D * υ 6 w = ε ( i p + c a p w ) .
Now, let
D * υ 1 s p = C 1 T Ψ m ^ ( t ) , D * υ 2 e p = C 2 T Ψ m ^ ( t ) , D * υ 3 i p = C 3 T Ψ m ^ ( t ) , D * υ 4 a p = C 4 T Ψ m ^ ( t ) , D * υ 5 r p = C 5 T Ψ m ^ ( t ) , D * υ 6 w = C 6 T Ψ m ^ ( t ) .
where
C 1 = [ c 1 , 1 , c 1 , 2 , , c 1 , m ^ ] T , C 2 = [ c 2 , 1 , c 2 , 2 , , c 2 , m ^ ] T , C 3 = [ c 3 , 1 , c 3 , 2 , , c 3 , m ^ ] T , C 4 = [ c 4 , 1 , c 4 , 2 , , c 4 , m ^ ] T , C 5 = [ c 5 , 1 , c 5 , 2 , , c 5 , m ^ ] T , C 6 = [ c 6 , 1 , c 6 , 2 , , c 6 , m ^ ] T .
Applying fractional integral operator on both sides of Equation (37) and then approximating the functions s p ( t ) , e p ( t ) , i p ( t ) , a p ( t ) , r p ( t ) and w ( t ) by sine–cosine wavelet operational matrix result in:
s p ( t ) = I υ 1 D * υ 1 s p ( t ) + s 0 C 1 T P m ^ × m ^ υ 1 Ψ m ^ ( t ) + s 0 e p ( t ) = I υ 2 D * υ 2 e p ( t ) + e 0 C 2 T P m ^ × m ^ υ 2 Ψ m ^ ( t ) + e 0 i p ( t ) = I υ 3 D * υ 3 i p ( t ) + i 0 C 3 T P m ^ × m ^ υ 3 Ψ m ^ ( t ) + i 0 a p ( t ) = I υ 4 D * υ 4 a p ( t ) + a 0 C 4 T P m ^ × m ^ υ 4 Ψ m ^ ( t ) + a 0 r p ( t ) = I υ 5 D * υ 5 r p ( t ) + r 0 C 5 T P m ^ × m ^ υ 5 Ψ m ^ ( t ) + r 0 w ( t ) = I υ 6 D * υ 6 w ( t ) + w 0 C 6 T P m ^ × m ^ υ 6 Ψ m ^ ( t ) + w 0
Collocating Equation (38) at the points t 2 i 1 2 m ^ i = 1 , , m ^ and then using Equation (22) lead to
s p ( t ) C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ B m ^ ( t ) + s 0 , , s 0 1 × m ^ B m ^ ( t ) , e p ( t ) C 2 T P m ^ × m ^ υ 2 Φ m ^ × m ^ B m ^ ( t ) + e 0 , , e 0 1 × m ^ B m ^ ( t ) , i p ( t ) C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ B m ^ ( t ) + i 0 , , i 0 1 × m ^ B m ^ ( t ) , a p ( t ) C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ B m ^ ( t ) + a 0 , , a 0 1 × m ^ B m ^ ( t ) , r p ( t ) C 5 T P m ^ × m ^ υ 5 Φ m ^ × m ^ B m ^ ( t ) + r 0 , , r 0 1 × m ^ B m ^ ( t ) , w ( t ) C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ B m ^ ( t ) + w 0 , , w 0 1 × m ^ B m ^ ( t ) ,
Then, the nonlinear terms in (36) can be expressed as
s p ( t ) i p ( t ) ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ ) B m ^ ( t ) + ( [ s 0 , , s 0 ] 1 × m ^ C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ ) B m ^ ( t ) + ( [ i 0 , , i 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ ) B m ^ ( t ) + [ s 0 i 0 , , s 0 i 0 ] 1 × m ^ B m ^ ( t ) ,
s p ( t ) a p ( t ) ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ ) B m ^ ( t ) + ( [ s 0 , , s 0 ] 1 × m ^ C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ ) B m ^ ( t ) + ( [ a 0 , , a 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ ) B m ^ ( t ) + [ s 0 a 0 , , s 0 a 0 ] 1 × m ^ B m ^ ( t ) ,
s p ( t ) w ( t ) ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ ) B m ^ ( t ) + ( [ s 0 , , s 0 ] 1 × m ^ C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ ) B m ^ ( t ) + ( [ w 0 , , w 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ ) B m ^ ( t ) + [ s 0 w 0 , , s 0 w 0 ] 1 × m ^ B m ^ ( t ) ,
Substituting Equations (37) and (39)–(42) into Equation (36), we have
C 1 T Φ m ^ × m ^ = [ n p , , n p ] 1 × m ^ m p C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 , , s 0 ] 1 × m ^ b p ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ + [ s 0 , , s 0 ] 1 × m ^ C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ + [ i 0 , , i 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 i 0 , s 0 i 0 ] 1 × m ^ ) k b p ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ s 0 , s 0 ] 1 × m ^ C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ a 0 , a 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 a 0 , s 0 a 0 ] 1 × m ^ ) b w ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ + [ s 0 , s 0 ] 1 × m ^ C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ + [ w 0 , w 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 w 0 , , s 0 w 0 ] 1 × m ^ ) , C 2 T Φ m ^ × m ^ = b p ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ ) + [ s 0 , , s 0 ] 1 × m ^ C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ + [ i 0 , , i 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 i 0 , , s 0 i 0 ] 1 × m ^ ) + k b p ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ s 0 , , s 0 ] 1 × m ^ C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ a 0 , , a 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 a 0 , , s 0 a 0 ] 1 × m ^ ) + b w ( C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ + [ s 0 , , s 0 ] 1 × m ^ C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ + [ w 0 , , w 0 ] 1 × m ^ C 1 T P m ^ × m ^ υ 1 Φ m ^ × m ^ + [ s 0 w 0 , , s 0 w 0 ] 1 × m ^ ) 1 δ p ω p C 2 T P m ^ × m ^ υ 2 Φ m ^ × m ^ + [ e 0 , , e 0 ] 1 × m ^ δ p ω p C 2 T P m ^ × m ^ υ 2 Φ m ^ × m ^ + [ e 0 , , e 0 ] 1 × m ^ m p C 2 T P m ^ × m ^ υ 2 Φ m ^ × m ^ + [ e 0 , , e 0 ] 1 × m ^ , C 3 T Φ m ^ × m ^ = ( 1 δ p ) ( ω p ) ( C 2 T P m ^ × m ^ υ 2 Φ m ^ × m ^ + [ e 0 , , e 0 ] 1 × m ^ ) ( γ p + m p ) ( C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ + [ i 0 , , i 0 ] 1 × m ^ ) , C 4 T Φ m ^ × m ^ = δ p ω p ( C 2 T P m ^ × m ^ υ 2 Φ m ^ × m ^ + [ e 0 , , e 0 ] 1 × m ^ ) ( γ p + m p ) ( C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ a 0 , , a 0 ] 1 × m ^ ) , C 5 T Φ m ^ × m ^ = γ p ( C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ + [ i 0 , , i 0 ] 1 × m ^ ) + γ p ( C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ a 0 , , a 0 ] 1 × m ^ ) m p ( C 5 T P m ^ × m ^ υ 5 Φ m ^ × m ^ + [ r 0 , , r 0 ] 1 × m ^ ) , C 6 T Φ m ^ × m ^ = ε ( ( C 3 T P m ^ × m ^ υ 3 Φ m ^ × m ^ + [ i 0 , , i 0 ] 1 × m ^ ) + c ( C 4 T P m ^ × m ^ υ 4 Φ m ^ × m ^ + [ a 0 , , a 0 ] 1 × m ^ ) ( C 6 T P m ^ × m ^ υ 6 Φ m ^ × m ^ + [ w 0 , , w 0 ] 1 × m ^ ) ) .
Finally, Equation (36) is converted into a system of nonlinear algebraic Equation (43). In order to solve this set, initial guesses are needed. To do this, nonlinear terms in Equation (43) are neglected to achieve a linear set of equations. The solutions of this linear set are employed as the initial guesses. Then, Equation (43) can be solved by a suitable numerical method, e.g., Newton’s method, to obtain the unknown coefficients c 1 , i ,   c 2 , i , c 3 , i ,   c 4 , i , c 5 , i and c 6 , i , for i = 1 , , m ^ . Substituting these coefficients into Equation (38), s p ( t ) , e p ( t ) , i p ( t ) , a p ( t ) , r p ( t ) and w ( t ) are attained.

5. Results and Discussion

In this section, Mathematica software was used to solve Equation (43) and to find an approximate solution of Equation (36). The initial values and parameters describing the model are given as [4]:
N p = 1000000000 , δ p = μ p = 0.5 , c = k = 0.5 , n p = m p = 0.0018 , w p = w p = 0.1923 , γ p = 0.1724 , β p = β w = γ p = 0.5 , ε = 0.1 .
under the following initial values:
s p 0 = 2 ,   e p 0 = 4 ,   i p 0 = 3 ,   a p 0 = 4 ,   r p 0 = 2 ,   w 0 = 3.5
.
Furthermore, it is assumed that: υ 1 = υ 2 = υ 3 = υ 4 = υ 5 = υ 6 = υ .
First, an experimental convergence analysis is conducted to guarantee a convergent solution for both sine–cosine and Bernoulli wavelets. Since the exact solution is not available, the solutions are verified by comparing them with the results obtained by the command NDSolve in software Mathematica using Runge-Kutta (RK4) method as well as the results obtained in [8] by ALP for υ = 1 . Finally, the effects of fractional orders on the COVID-19 outbreak are studied.

5.1. Convergence of the Solution

In order to obtain a convergent solution, Equation (43) is solved for υ = 0.5 and various values of k and M for Bernoulli wavelets as well as different k and L for sine–cosine wavelets. In addition, rate of convergence (ROC) is calculated to demonstrate how fast the solution converges as follows.
R O C m ^ = m ^ 4 = log f m ^ 4 f m ^ 3 f m ^ 3 f m ^ 2 log f m ^ 3 f m ^ 2 f m ^ 2 f m ^ 1
where, f stands for each of the calculated functions s p ,   e p ,   i p ,   a p ,   r p and w p at a certain time t = t i for four successive values of m ^ : m ^ 1 < m ^ 2 < m ^ 3 < m ^ 4 .
Table 2, Table 3, Table 4, Table 5, Table 6 and Table 7 demonstrate s p ,   e p ,   i p ,   a p ,   r p and w p evaluated at t = 0.2 ,   0.4 ,   0.6 ,   0.8 for υ = 0.5 , k = 2 and various values of M using Bernoulli wavelets. It is inferred that by increasing the order of approximations, the solutions converge with higher precision. The ROC reveals that the solutions converge more rapidly around starting and ending time ( t = 0.2 and t = 0.8 ). Furthermore, the ROC might be roughly estimated as R O C 1 . Moreover, it is clarified that e p ,   i p ,   a p ,   r p and w p express more rapid convergence in comparison with s p . It is evident that by using Bernoulli wavelets with k = 2 and M 8 , s p converges with at least two decimal places at t = 0.6 and t = 0.8 whereas e p ,   i p ,   a p ,   r p and w p converge with at least four decimal places at the mentioned time stages.
Furthermore, Table 8, Table 9, Table 10, Table 11, Table 12 and Table 13 show evaluated s p ,   e p ,   i p ,   a p ,   r p and w p at t = 0.2 ,   0.4 ,   0.6 ,   0.8 for υ = 0.5 , k = 0 and different values of L using sine–cosine wavelets. These tables clarify that higher orders of approximations result in higher precision. Moreover, the ROC proves that the solutions converge more rapidly around the mid time ( t = 0.4 and t = 0.6 ). Moreover, it can be inferred that e p ,   i p ,   a p ,   r p and w p are more convergent in comparison with s p . It is apparent that by employing sine–cosine wavelets with k = 0 and L 10 , s p converges with at least one decimal place at t = 0.4 and t = 0.6 whereas e p ,   i p ,   a p ,   r p and w p converge with at least two decimal places at those time stages.
Comparing two mentioned wavelets methods suggests that the Bernoulli wavelets method provides higher precision with lower order of approximation compared to the sine–cosine wavelets method. Furthermore, the sine–cosine wavelets method slowly converges around the starting and ending time stages.

5.2. Verification of the Solution

In this section, the accuracy of the solutions achieved by the two mentioned methods is considered. As the exact solution is not available, there is no absolute and definite criterion assuring us of the level of accuracy. However, comparing the results of four various methods including RK4, ALPs, Bernoulli and sine–cosine wavelets methods might relatively verify the obtained results.
Thus, the results attained via using Bernoulli wavelets with k = 2 and M = 8 ( m ^ = 16 ), sine–cosine wavelets with k = 0 and L = 8 ( m ^ = 17 ), ALPs with m ^ = 16 and RK4 are compared in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. It is evident that the results obtained by ALPs and sine–cosine wavelets remarkably deviate from RK4 solution around starting and ending time stages whereas the Bernoulli wavelets solutions agree well with RK4 results. Considering the evidence observed in Section 5.1 which states that the sine–cosine wavelets solutions slowly converge at starting and ending time stages, it can be concluded that Bernoulli wavelets solutions are more trustworthy compared to the sine–cosine wavelets solutions through whole domain.

5.3. Computational Cost

There are several ways of determining computational cost analytically in which the number of algebraic operations is calculated. However, the computational cost can be practically measured by calculating the CPU running time. Since the running time of the ALPs method was not reported in literature for this problem, only Bernoulli and sine–cosine wavelets methods are compared together. The CPU running time for these two procedures are tabulated in Table 14. It is obvious that the running time for the Bernoulli wavelets method is considerably low for lower orders of approximations, but it abruptly increases at M = 10 . In addition, it is slightly lower than that of the sine–cosine wavelets even with higher orders. It can be concluded that the Bernoulli wavelets method imposes less computational costs, yet it converges faster and provides higher level of accuracy in comparison with the sine–cosine wavelets. Hence, this method is recommended to be employed in solving the current model of COVID-19 outbreak in possible future studies. Specifically, it is utilized in Section 5.4 to study the effects of the fractional orders on the solution.

5.4. The Effects of Fractional Orders

As it is mentioned in Section 1, fractional derivatives include memory of the system leading to constructing a more biologically compatible model compared to the integer-order derivatives. However, it is essential to determine the proper values for the fractional orders.
The current model of COVID-19 outbreak might be used to simulate any epidemic in a local region including any arbitrary country or state. To achieve such a model, all parameters including the fractional orders must be evaluated in a way the results of the simulation become best adjusted to the real-world data. This might be accomplished if the effects of the fractional orders on the solutions are considered. In this section, these effects are studied.
Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 illustrate, respectively, s p ,   e p ,   i p ,   a p ,   r p and w p versus time for various fractional orders, υ = υ 1 = υ 2 = υ 3 = υ 4 = υ 5 = υ 6 = 0.2 ,   0.4 ,   0.6 ,   0.8 . It is observed that the functions s p ,   e p ,   i p and a p rise when the fractional orders increase. In contrast, higher values of the fractional orders cause the functions r p and w p to decrease. This proves that using integer-order derivatives instead of fractional ones may remarkably affect the prediction of the future outbreak whereas employing fractional orders might be more advantageous.

6. Conclusions

In this paper, the model of novel coronavirus is presented as a system of nonlinear fractional differential equations based on a model previously proposed in literature with integer-order derivatives. Then, two new procedures to approximate the solution of the model are proposed. The sine–cosine and Bernoulli wavelets and their operational matrices are employed to transform the studied system of nonlinear fractional differential equations into nonlinear set of algebraic equations. The two mentioned methods are examined for rate of convergence, accuracy and computational costs via conducting an experimental convergence analysis, solution verification and CPU running time comparison. The results are compared in special cases with those reported in the literature obtained by alternative Legendre polynomials (ALPs) and with the results attained using Runge-Kutta (RK4). It is proved that the Bernoulli wavelets method exhibits faster convergence, lower computational cost and higher level of accuracy compared to sine–cosine wavelets method. Furthermore, the superiority of this method over previously employed APL method is demonstrated. Thus, this method is recommended for possible future investigation of COVID-19 spread via using this model. Then, the effects of the fractional orders on the COVID-19 transmission are studied. Finally, possible application of this parametric study is clarified.

Author Contributions

Conceptualization, R.E. and M.H.; methodology, R.E., M.H. and S.N.; validation, R.E., M.H. and S.N.; software, M.H.; formal analysis, R.E., M.H. and S.N.; investigation, R.E., M.H. and S.N.; writing—review and editing, R.E., M.H. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

The research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alluwaimi, A.M.; Alshubaith, I.H.; Al-Ali, A.M.; Abohelaika, S. The Coronaviruses of Animals and Birds: Their Zoonosis, Vaccines, and Models for SARS-CoV and SARS-CoV2. Front. Vet. Sci. 2020, 7, 655. [Google Scholar] [CrossRef] [PubMed]
  2. Madjid, M.; Safavi-Naeini, P.; Solomon, S.D.; Vardeny, O. Potential effects of coronaviruses on the cardiovascular system: A review. JAMA Cardiol. 2020, 5, 840–931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. WHO. Novel Coronavirus (2019-nCoV): Situation Report 3; WHO: Geneva, Switzerland, 2020. [Google Scholar]
  4. Chen, T.M.; Rui, J.; Wang, Q.P.; Zhao, Z.Y.; Cui, J.A.; Yin, L. A mathematical model for simulating the phase-based transmissibility of a novel coronavirus. Infect. Dis. Poverty 2020, 9, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Khan, M.A.; Atangana, A. Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative. Alex. Eng. J. 2020, 59, 2379–2389. [Google Scholar] [CrossRef]
  6. Rajagopal, K.; Hasanzadeh, N.; Parastesh, F.; Hamarash, I.I.; Jafari, S.; Hussain, I. A fractional-order model for the novel coronavirus (COVID-19) outbreak. Nonlinear Dyn. 2020, 101, 711–718. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, X.; Li, J.; Xiao, C.; Yang, P. Numerical solution and parameter estimation for uncertain SIR model with application to COVID-19. Fuzzy Optim. Decis. Mak. 2021, 20, 189–208. [Google Scholar] [CrossRef]
  8. Hashemizadeh, E.; Ebadi, M.A. A numerical solution by alternative Legendre polynomials on a model for novel coronavirus (COVID-19). Adv. Differ. Equ. 2020, 2020, 527. [Google Scholar] [CrossRef]
  9. Kisela, T. Fractional Differential Equations and Their Applications; Faculty of Mechanical Engineering Institute of Mathematics: Brno-střed, Czech Republic, 2008. [Google Scholar]
  10. Solís-Pérez, J.E.; Gómez-Aguilar, J.F.; Atangana, A. A fractional mathematical model of breast cancer competition model. Chaos Solitons Fractals 2019, 127, 38–54. [Google Scholar] [CrossRef]
  11. Hussain, A.; Baleanu, D.; Adeel, M. Existence of solution and stability for the fractional order novel coronavirus (nCoV-2019) model. Adv. Differ. Equ. 2020, 2020, 1–9. [Google Scholar] [CrossRef]
  12. Atangana, A. Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Chaos Solitons Fractals 2020, 136, 109860. [Google Scholar] [CrossRef]
  13. Wang, B.H.; Wang, Y.Y.; Dai, C.Q.; Chen, Y.X. Dynamical characteristic of analytical fractional solitons for the space-time fractional Fokas-Lenells equation. Alex. Eng. J. 2020, 59, 4699–4707. [Google Scholar] [CrossRef]
  14. Yu, L.J.; Wu, G.Z.; Wang, Y.Y.; Chen, Y.X. Traveling wave solutions constructed by Mittag–Leffler function of a (2 + 1)-dimensional space-time fractional NLS equation. Results Phys. 2020, 17, 103156. [Google Scholar] [CrossRef]
  15. Fang, J.J.; Mou, D.S.; Wang, Y.Y.; Zhang, H.C.; Dai, C.Q.; Chen, Y.X. Soliton dynamics based on exact solutions of conformable fractional discrete complex cubic Ginzburg–Landau equation. Results Phys. 2021, 20, 103710. [Google Scholar] [CrossRef]
  16. Wang, B.-H.; Wang, Y.-Y. Fractional white noise functional soliton solutions of a wick-type stochastic fractional NLSE. Appl. Math. Lett. 2020, 110, 106583. [Google Scholar] [CrossRef]
  17. Wang, B.-H.; Wang, Y.-Y.; Dai, C.-Q. Fractional optical solitons with stochastic properties of a wick-type stochastic fractional NLSE driven by the Brownian motion. Waves Random Complex Media 2021, 1–14. [Google Scholar]
  18. Noeiaghdam, S. A novel technique to solve the modified epidemiological model of computer viruses. SEMA J. 2019, 76, 97–108. [Google Scholar] [CrossRef]
  19. Noeiaghdam, S.; Micula, S. Dynamical Strategy to Control the Accuracy of the Nonlinear Bio-Mathematical Model of Malaria Infection. Mathematics 2021, 9, 1031. [Google Scholar] [CrossRef]
  20. Noeiaghdam, S.; Suleman, M.; Budak, H. Solving a modified nonlinear epidemiological model of computer viruses by homotopy analysis method. Math. Sci. 2018, 12, 211–222. [Google Scholar] [CrossRef] [Green Version]
  21. Tuan, N.H.; Mohammadi, H.; Rezapour, S. A mathematical model for COVID-19 transmission by using the Caputo fractional derivative. Chaos Solitons Fractals 2020, 140, 110107. [Google Scholar] [CrossRef] [PubMed]
  22. Yi, M.; Huang, J. Wavelet operational matrix method for solving fractional differential equations with variable coefficients. Appl. Math. Comput. 2014, 230, 383–394. [Google Scholar] [CrossRef]
  23. Yi, M.; Wang, L.; Huang, J. Legendre wavelets method for the numerical solution of fractional integro-differential equations with weakly singular kernel. Appl. Math. Model. 2016, 40, 3422–3437. [Google Scholar] [CrossRef]
  24. Razzaghi, M.; Yousefi, S. Sine-cosine wavelets operational matrix of integration and its applications in the calculus of variations. Int. J. Syst. Sci. 2002, 33, 805–810. [Google Scholar] [CrossRef]
  25. Saeed, A.; Saeed, U. Sine-cosine wavelet method for fractional oscillator equations. Math. Methods Appl. Sci. 2019, 42, 6960–6971. [Google Scholar] [CrossRef]
  26. Wang, Y.; Yin, T.; Zhu, L. Sine-cosine wavelet operational matrix of fractional order integration and its applications in solving the fractional order Riccati differential equations. Adv. Differ. Equ. 2017, 2017, 222. [Google Scholar] [CrossRef] [Green Version]
  27. Kajani, M.T.; Ghasemi, M.; Babolian, E. Numerical solution of linear integro-differential equation by using sine–cosine wavelets. Appl. Math. Comput. 2006, 180, 569–574. [Google Scholar] [CrossRef]
  28. Kilicman, A.; Al Zhour, Z.A.A. Kronecker operational matrices for fractional calculus and some applications. Appl. Math. Comput. 2007, 187, 250–265. [Google Scholar] [CrossRef]
  29. Keshavarz, E.; Ordokhani, Y. Bernoulli wavelets method for solution of fractional differential equations in a large interval. Math. Res. 2016, 2, 17–32. [Google Scholar] [CrossRef]
  30. Rahimkhani, P.; Ordokhani, Y.; Babolian, E. Numerical solution of fractional pantograph differential equations by using generalized fractional-order Bernoulli wavelet. J. Comput. Appl. Math. 2017, 309, 493–510. [Google Scholar] [CrossRef]
  31. Tohidi, E.; Bhrawy, A.; Erfani, K. A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation. Appl. Math. Model. 2013, 37, 4283–4294. [Google Scholar] [CrossRef]
  32. Soltanpour Moghadam, A.; Arabameri, M.; Baleanu, D.; Barfeie, M. Numerical solution of variable fractional order advection-dispersion equation using Bernoulli wavelet method and new operational matrix of fractional order derivative. Math. Methods Appl. Sci. 2020, 43, 3936–3953. [Google Scholar] [CrossRef]
Figure 1. Comparison of the numerical solutions of Sp(t) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Figure 1. Comparison of the numerical solutions of Sp(t) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Axioms 10 00122 g001
Figure 2. Comparison of the numerical solutions of e p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Figure 2. Comparison of the numerical solutions of e p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Axioms 10 00122 g002
Figure 3. Comparison of the numerical solutions of i p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Figure 3. Comparison of the numerical solutions of i p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Axioms 10 00122 g003
Figure 4. Comparison of the numerical solutions of a p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Figure 4. Comparison of the numerical solutions of a p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Axioms 10 00122 g004
Figure 5. Comparison of the numerical solutions of r p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Figure 5. Comparison of the numerical solutions of r p ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Axioms 10 00122 g005
Figure 6. Comparison of the numerical solutions of w ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Figure 6. Comparison of the numerical solutions of w ( t ) by using the sine–cosine wavelets method, Bernoulli wavelets method, ALPs and RK4 methods for υ = 1 .
Axioms 10 00122 g006
Figure 7. The effect of υ on the distribution of s p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Figure 7. The effect of υ on the distribution of s p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Axioms 10 00122 g007
Figure 8. The effect of υ on the distribution of e p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Figure 8. The effect of υ on the distribution of e p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Axioms 10 00122 g008
Figure 9. The effect of υ on the distribution of i p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Figure 9. The effect of υ on the distribution of i p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Axioms 10 00122 g009
Figure 10. The effect of υ on the distribution of a p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Figure 10. The effect of υ on the distribution of a p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Axioms 10 00122 g010
Figure 11. The effect of υ on the distribution of r p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Figure 11. The effect of υ on the distribution of r p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Axioms 10 00122 g011
Figure 12. The effect of υ on the distribution of w p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Figure 12. The effect of υ on the distribution of w p t . Black: υ = 0.2 , Blue: υ = 0.4 , Red: υ = 0.6 , Green: υ = 0.8 .
Axioms 10 00122 g012
Table 1. Variables and Parameters.
Table 1. Variables and Parameters.
Variables and ParametersDefinitionVariables and ParametersDefinition
npThe birth rate of people.µpThe shedding coefficients from Ap to W.
mpThe death rate of people.δpThe proportion of asymptomatic infection rate of people.
1 ω p The incubation period of people.βpThe transmission rate from Ip to Sp.
1 ω p The latent period of people.βWThe transmission rate from W to Sp.
1 γ p The infectious period of symptomatic infection in people.kThe multiple of the transmissibility of Ap to that of Ip.
1 γ p The infectious period of asymptomatic infection in people. 1 ε The lifetime of the virus in W.
µpThe shedding coefficients from Ip to W.cThe relative shedding rate of Ap compared to Ip.
Table 2. The values of 1010 Sp and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
Table 2. The values of 1010 Sp and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
t1010 Sp
M = 4
1010 Sp
M = 5
1010 Sp
M = 6
1010 Sp
M = 7
ROC
M = 7
1010 Sp
M = 8
ROC
M = 8
1010 Sp
M = 9
ROC
M = 9
1010 Sp
M = 10
ROC
M = 10
0.21.081.772.572.62−18.952.24−0.732.03−0.292.131.16
0.42.280.8042.091.424.691.473.831.81−0.731.350.16
0.61.301.291.291.29−0.011.29−2341.290.021.29−20
0.81.121.121.121.123.301.121.101.12−0.861.121.30
Table 3. The values of ep and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
Table 3. The values of ep and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
tep
M = 4
ep
M = 5
ep
M = 6
ep
M = 7
ROC
M = 7
ep
M = 8
ROC
M = 8
ep
M = 9
ROC
M = 9
ep
M = 10
ROC
M = 10
0.25.45705.45605.45655.45643.175.4562−0.775.4560−0.175.45592.09
0.45.25605.25195.25355.25290.945.25293.795.2529−0.255.25280.94
0.65.10565.10545.10535.10530.725.10520.795.10520.855.10520.88
0.84.98614.98604.98594.98591.134.98590.794.98590.874.98590.88
Table 4. The values of ip and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
Table 4. The values of ip and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
tip
M = 4
ip
M = 5
ip
M = 6
ip
M = 7
ROC
M = 7
ip
M = 8
ROC
M = 8
ip
M = 9
ROC
M = 9
ip
M = 10
ROC
M = 10
0.23.00583.00583.00583.0057−0.883.00580.393.00580.453.00585.47
0.42.99742.99762.99752.99751.792.99751.592.9975−0.0252.9975−0.8
0.62.98782.98782.98782.98780.952.98770.902.98770.882.98770.89
0.82.97752.97752.97752.97750.742.97750.882.97750.872.97750.88
Table 5. The values of ap and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
Table 5. The values of ap and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
tap
M = 4
ap
M = 5
ap
M = 6
ap
M = 7
ROC
M = 7
ap
M = 8
ROC
M = 8
ap
M = 9
ROC
M = 9
ap
M = 10
ROC
M = 10
0.23.00583.00583.00583.0057−0.883.00580.393.00580.453.00585.47
0.42.99742.99762.99752.99751.792.99751.592.9975−0.0252.9975−0.8
0.62.98782.98782.98782.98780.952.98770.902.98770.882.98770.88
0.82.97752.97752.97752.97750.742.97750.882.97750.872.97750.89
Table 6. The values of rp and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
Table 6. The values of rp and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
trp
M = 4
rp
M = 5
rp
M = 6
rp
M = 7
ROC
M = 7
rp
M = 8
ROC
M = 8
rp
M = 9
ROC
M = 9
rp
M = 10
ROC
M = 10
0.23.14033.14223.14123.14126.153.1417−0.803.1421−0.113.14232.20
0.43.54763.55633.55293.55421.003.55434.183.5542−0.383.55440.47
0.63.84943.84973.84993.85000.723.85010.783.85010.853.85010.88
0.84.08694.08704.08714.08711.214.08720.784.0872−0.874.08720.89
Table 7. The values of wp and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
Table 7. The values of wp and rate of convergence (ROC) at different points for k = 2 and various values of M by Bernoulli wavelets.
twp
M = 4
wp
M = 5
wp
M = 6
wp
M = 7
ROC
M = 7
wp
M = 8
ROC
M = 8
wp
M = 9
ROC
M = 9
wp
M = 10
ROC
M = 10
0.23.56083.56093.56093.56099.393.5609−0.773.5609−0.103.56092.31
0.43.57833.57893.57863.57871.063.57875.13.5787−0.533.57870.21
0.63.58983.58983.58983.58980.713.58980.783.58980.853.58980.88
0.83.59753.59763.59763.59761.363.59760.763.5976−0.873.59760.89
Table 8. The values of 1010 Sp and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
Table 8. The values of 1010 Sp and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
t1010 Sp
L = 4
1010 Sp
L = 5
1010 Sp
L = 6
1010 Sp
L = 7
ROC
L = 7
1010 Sp
L = 8
ROC
L = 8
1010 Sp
L = 9
ROC
L = 9
1010 Sp
L = 10
ROC
L = 10
0.20.3681.263.063.62−1.682.44−0.631.170.101.37−24.1
0.41.161.362.200.8740.312.11−0.131.398.581.407.31
0.61.531.520.7530.920.080.821−0.141.478.871.468.38
0.81.961.910.8400.246−0.190.865−0.070.789.241.75−9.38
Table 9. The values of ep and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
Table 9. The values of ep and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
tep
L = 4
ep
L= 5
ep
L = 6
ep
L = 7
ROC
L = 7
ep
L = 8
ROC
L = 8
ep
L = 9
ROC
L = 9
ep
L = 10
ROC
L = 10
0.25.38335.39515.47665.5134−0.415.4717155.41891.915.4222−11.6
0.45.23575.23855.28575.21750.125.2783−0.315.24414.965.24496.59
0.65.12275.11985.07295.14030.135.0800−0.305.11395.185.11316.50
0.85.05485.04434.96664.9308−0.384.9707−0.135.02192.358.0187−11.0
Table 10. The values of ip and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
Table 10. The values of ip and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
tip
L = 4
ip
L = 5
ip
L = 6
ip
L = 7
ROC
L = 7
ip
L = 8
ROC
L = 8
ip
L = 9
ROC
L = 9
ip
L = 10
ROC
L = 10
0.23.00323.00363.00653.0077−0.403.0063−0.143.00452.093.0046−11.5
0.42.99692.99702.99872.99620.132.9984−0.342.99724.692.99726.39
0.62.98842.98832.98652.98900.132.9868−0.352.98104.782.98806.22
0.82.98032.97982.97682.9754−0.452.9769−0.142.97881.812.9787−12.1
Table 11. The values of ap and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
Table 11. The values of ap and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
tap
L = 4
ap
L = 5
ap
L = 6
ap
L = 7
ROC
L = 7
ap
L = 8
ROC
L = 8
ap
L = 9
ROC
L = 9
ap
L = 10
ROC
L = 10
0.23.0803.32033.40463.4426−0.413.3996−0.153.34491.913.3483−11.6
0.43.15953.16243.21113.14080.123.2035−0.313.16825.003.16906.60
0.63.05273.04973.00153.07090.133.0088−0.303.04375.203.04296.53
0.82.99423.98352.90402.8672−0.382.9082−0.132.96092.392.9577−11.0
Table 12. The values of rp and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
Table 12. The values of rp and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
trp
L = 4
rp
L = 5
rp
L = 6
rp
L = 7
ROC
L = 7
rp
L = 8
ROC
L = 8
rp
L = 9
ROC
L = 9
rp
L = 10
ROC
L = 10
0.23.29073.26663.10003.0248−0.413.1099−0.153.21781.913.2111−11.6
0.43.58933.58363.48723.62630.123.5023−0.313.57214.973.57066.59
0.63.81443.82043.91603.77850.133.9016−0.303.83245.183.83406.51
0.84.94694.96834.12654.1995−0.384.1182−0.134.01372.364.0201−11.0
Table 13. The values of wp and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
Table 13. The values of wp and rate of convergence (ROC) at different points for k = 0 and various values of L by sine–cosine wavelets.
twp
L = 4
wp
L = 5
wp
L = 6
wp
L = 7
ROC
L = 7
wp
L = 8
ROC
L = 8
wp
L = 9
ROC
L = 9
wp
L = 10
ROC
L = 10
0.23.56753.56643.55913.5557−0.413.5595−0.153.56431.883.5640−11.6
0.43.58023.58003.57583.58190.123.5764−0.303.57955.013.57946.64
0.63.58833.58853.59273.58670.133.5921−0.293.58905.283.58916.55
0.83.59153.59243.59933.9024−0.373.5989−0.133.59432.4483.5946−10.8
Table 14. Comparison of CPU running time in seconds, used in Bernoulli for k = 2 and sine–cosine wavelets for k = 0.
Table 14. Comparison of CPU running time in seconds, used in Bernoulli for k = 2 and sine–cosine wavelets for k = 0.
Bernoulli WaveletsM = 4
m ^ = 8  
M = 5
m ^ = 10  
M = 6
m ^ = 12  
L = 7
m ^ = 14  
M = 8
m ^ = 16  
M = 9
m ^ = 18  
M = 10
m ^ = 20  
CPU running time0.280.671.712.553.919.0067.9
Sine-cosine waveletsL = 4
m ^ = 9  
L = 5
m ^ = 11  
L = 6
m ^ = 13  
L = 7
m ^ = 15  
L = 8
m ^ = 17  
L = 9
m ^ = 19  
L = 10
m ^ = 21
CPU running time0.821.733.085.358.869.6713.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hedayati, M.; Ezzati, R.; Noeiaghdam, S. New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method. Axioms 2021, 10, 122. https://doi.org/10.3390/axioms10020122

AMA Style

Hedayati M, Ezzati R, Noeiaghdam S. New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method. Axioms. 2021; 10(2):122. https://doi.org/10.3390/axioms10020122

Chicago/Turabian Style

Hedayati, Maryamsadat, Reza Ezzati, and Samad Noeiaghdam. 2021. "New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method" Axioms 10, no. 2: 122. https://doi.org/10.3390/axioms10020122

APA Style

Hedayati, M., Ezzati, R., & Noeiaghdam, S. (2021). New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method. Axioms, 10(2), 122. https://doi.org/10.3390/axioms10020122

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