Wavelet approximation scheme for distributed order fractional differential equations
Introduction
In mathematical science, the topic of fractional differential equations (FDEs) is ancient. FDEs are generalization of ordinary differential equations (ODEs) and partial differential equations (PDEs) to an arbitrary non-integer order. Fractional derivative is a great tool, as it is used to explain the consciousness and biological characteristics of various techniques and materials. This is the main reason why FDEs are used to model electrical and mechanical properties of actual materials and archaeological properties of rocks and several other fields. For the development history of fractional differential operator (FDO) one can see [1]. The FDEs are widely used in different fields of science and engineering such as viscoelastic damped structure model [2], hysteretic damping model [3], transient analysis of viscoelastic structure [4], electrochemistry [5] and ultraslow kinetics [6]. To the best of our fundamental understanding, two separate approaches based on spectral methods were carried out for FODEs and FPDEs. The first approach contains the classical orthogonal functions as a trail function to approximate the corresponding FDEs. In the second approach, Jacobi polyfractonomials were used as a trail function to find the approximate solutions of FODEs and FPDEs (see for example [7] and references therein). The discussion of fractional-order control and system has attracted many researchers to work on as an important application field of fractional calculus.
In sciences and engineering there are number of existing problems which are governed by the distributed order FDEs [8], [9]. If the fractional calculus operators operate on , and we integrate with respect to the order , then distributed-order differential/integral equations can be obtained. In 1969, Caputo was the first who presented the concept of distributed-order equation and solved it later in 1995. The distributed order equation is designed to model the input–output relationship of a linear time-invariant system based on observation of the frequency domain response. In other words, the distributed-order equation is the time-domain representation of the input–output relationship perceived and constructed in the frequency domain. The distributed-order differential equations and their applications in engineering fields have recently gained much attention. For instance, Bagley and Torvik [10] discussed the general solution of linear distributed-order differential equation and also examined the existence of solution of DOFDEs. In [11], distributed-order equations were introduced in the constitutive dielectric media. In 2016, the fractional kinetics of distributed-order equations was discussed by Sokolov et al. [12]. In 2006, Umarov and Steinberg [13] proposed multi-dimensional random walk models governed by distributed fractional order differential equations. In particular, the distributed-order operator is a more reliable method for understanding and defining certain actual physical processes such as the complexities of nonlinear systems [14], [15], networked structures [16] and nonhomogeneous phenomena [9].
In modeling many fields of mathematical physics, chemistry, biology and electronics, the study of distributed order fractional differential equations (DFDEs) plays major role; such as in the distributed order oscillator [17], complex system [18], viscoelastic model [19], PI-control systems [20], distributed order system identification [21], frequency control system [22], and many more (see [23] and references therein). It has already been shown that both integer and fractional order differential equations are the special cases of DFDEs. One can see [10], [24] regarding the uniqueness and existence of the solution of DFDEs. Since most of the DFDEs do not have precise analytical solutions therefore an efficient numerical methods are needed to solve them. For this purpose, computational methods such as Petrov–Galerkin and spectral collocation methods [25], Legendre collocation method [26], pseudo-spectral method [27], improved composite collocation method [28], spectral method [29], an improved meshless method [30], finite element method [31], Legendre spectral element method [32], implicit difference scheme [33], Chebyshev collocation method [34], Legendre–Gauss collocation method [35], block-pulse wavelet [36] are widely used.
The operational matrix method has been proved simple and efficient computational technique for solving fractional differential equations. In 2010, Saadatmandi and Dehghan [37] proposed and derived the operating matrix of fractional derivative (OMFD) for the shifted Legendre polynomials (SLPs). Li et al. [38] developed the haar wavelet OMFD in 2010. Keshavarz et al. [39] also launched the Bernoulli wavelets OMFD in 2014. The OMFD for B-spline functions was also constructed in 2012 by Lakestani et al. [40]. In 2013, Bhrawy et al. [41] obtained the OMFD for the generalize Laguerre polynomials. In 2014 Kayedi-Bardeh et al. [42] developed the OMFD using the fractional jacobi function. The interested readers can also see [43], [44] for more details of operational matrices.
Recently, M. Pourbabaee and Saadatmandi [45] provided a useful technique based on the operational matrices of Legendre polynomials for finding the approximate solution of DFDE along with initial conditions. Based on the facts and to the best of our knowledge, there is no numerical technique available based on wavelets operational matrix for solving DFDEs. Motivated by the work of M. Pourbabaee and Saadatmandi [45] and usefulness of operational matrices, the aim of this paper is to create the distributed order fractional derivative operational matrix for Legendre and Chebyshev wavelets. The purpose of using Legendre and Chebyshev wavelets in this paper are the properties of their basis functions which are constructed by orthogonal polynomials. Thus, Legendre and Chebyshev wavelets basis have advantages of both infinitely differentiable functions and small compact supports. Furthermore, the operational matrices of Legendre and Chebyshev wavelets have the sparsity and thus contain many zeros, which make computation simple and fast. These operational matrices are utilized along with Gauss–Legendre quadrature rule and standard Tau approach to solve linear DFDEs and distributed order time fractional diffusion equations (DTFDEs). The motivation of this method is the provide good experimental results with less computational complexity.
The remainder of the paper is presented as follows. Section 2 briefly presents some basic definitions which include fractional derivative, the distributed differential operator, Legendre wavelet, Chebyshev wavelet and their approximation properties. Also, we define Gauss Legendre quadrature formula in this section. Section 3 is devoted to the operational matrices. In this section, operational matrices of derivative of Legendre and Chebyshev wavelets for integer order and distributed order are derived. The implementation of the method to one and two-variables distributed order equations is described in Section 4. Convergence analysis and error bounds with the help of Legendre as well as Chebyshev wavelets for the mentioned scheme are given in Section 5. Error estimation of the presented method is provided in Section 6. Finally, numerical experiments are performed alongwith discussion of numerical stability in Section 7.
Section snippets
Definition of fractional derivative in the Caputo sense
The Caputo fractional derivative with order is defined as [1]
Here is the Gamma function. This Caputo operator has some following properties:
- •
, where, C is any arbitrary constant.
- •
The Caputo derivative of is given as: where, is the ceiling function.
- •
satisfies the linearity property, i.e. where, are arbitrary constants for
Legendre wavelets integer order derivative matrix
Theorem 3.1 Let be the Legendre wavelets vector defined in Section 2.3. The derivative of can be expressed as [48] where is the operational matrix defined as follows: in which F is matrix and its th element is defined as follows:
The above matrix can be generalized for higher orders as
Chebyshev wavelets integer order derivative matrix
Theorem 3.2 Let be the Chebyshev wavelets
Numerical method: Application to the distributed order linear equations
In this section we discuss numerical method for solving distributed order linear fractional differential equations in one and two variables. We provide two schemes namely scheme (I) and scheme (II) for dealing with distributed order fractional differential equations in one and two-variable, respectively.
Convergence analysis and error bound by using Legendre wavelets
Theorem 5.1 Let be the approximation of continuous functions , respectively defined on with bounded mixed second derivative, say for some positive constant , where, and . Then
can be expanded as an infinite sum of Legendre wavelets and the series converges uniformly to the function , that is
Error estimation for one variable case
In this part, we discuss an error estimation for one variable distributed order fractional differential equation. We have one-variable distributed fractional differential equation as under initial conditions First, we assume be the error function, where, is exact and is approximate solutions of Eq. (6.1). Hence satisfies the following equations under initial
Numerical examples
In this part, we consider five test functions in order to verify the usefulness and accuracy of the method. The numerical results obtained via the presented method are provided in the form of Tables and Figures for different values of and . Numerical results of -errors, -errors and pointwise errors along with used CPU time are provided in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11 and 18, 19. Comparison of numerical results produced by
Conclusion
In this paper, a new approach based on operational matrices of wavelets and standard Tau method is developed for solving distributed order fractional differential equations. Operational matrices of Legendre and Chebyshev wavelets together with standard Tau method and Legendre Gauss quadrature formula convert the original problem into a system of linear algebraic equations. The testing of method is done on several test problems. To verify the efficiency and accuracy of the method -errors,
CRediT authorship contribution statement
Yashveer Kumar: Conceptualization, Software, Writing. Somveer Singh: Methodology, Editing. Nikhil Srivastava: Visualization, Investigation. Aman Singh: Writing - review. Vineet Kumar Singh: Supervision.
Acknowledgments
The first & fourth authors acknowledge the financial support from Council of Scientific & Industrial Research (CSIR), India, under Junior Research Fellow (JRF) scheme. The second author acknowledges the financial support from National Board for Higher Mathematics, Department of Atomic Energy, India, with sanction order no. 0204/17/2019/R & D-II/9814. The third author acknowledges the financial support from Ministry of Human Resource and Development (MHRD) , under JRF scheme.
References (53)
Fractional differential equations in electrochemistry
Adv. Eng. Softw.
(2010)- et al.
Numerical solution of distributed order fractional differential equations by hybrid functions
J. Comput. Phys.
(2016) - et al.
On a nonlinear distributed order fractional differential equation
J. Math. Anal. Appl.
(2007) - et al.
Semilinear ordinary differential equation coupled with distributed order fractional differential equation
Nonlinear Anal. TMA
(2010) - et al.
Design, implementation and application of distributed order pi control
ISA Trans.
(2013) - et al.
Fractional-order system identification based on continuous order-distributions
Signal processing
(2003) - et al.
Numerical analysis for distributed-order differential equations
J. Comput. Appl. Math.
(2009) - et al.
An improved composite collocation method for distributed-order fractional differential equations based on fractional chelyshkov wavelets
Appl. Numer. Math.
(2019) - et al.
Numerical solution for diffusion equations with distributed order in time using a Chebyshev collocation method
Appl. Numer. Math.
(2017) - et al.
Error analysis of the legendre-gauss collocation methods for the nonlinear distributed-order fractional differential equation
Appl. Numer. Math.
(2019)
A new operational matrix for solving fractional-order differential equations
Comput. Math. Appl.
Haar wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations
Appl. Math. Comput.
Bernoulli wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations
Appl. Math. Model.
The construction of operational matrix of fractional derivatives using b-spline functions
Commun. Nonlinear Sci. Numer. Simul.
Operational matrix approach for the solution of partial integro-differential equation
Appl. Math. Comput.
Numerical solution of nonlinear weakly singular partial integro-differential equation via operational matrices
Appl. Math. Comput.
A novel legendre operational matrix for distributed order fractional differential equations
Appl. Math. Comput.
Application of wavelet collocation method for hyperbolic partial differential equations via matrices
Appl. Math. Comput.
A new legendre wavelet operational matrix of derivative and its applications in solving the singular ordinary differential equations
J. Franklin Inst. B
A stable numerical inversion of abel’s integral equation using almost bernstein operational matrix
J. Quant. Spectrosc. Radiat. Transfer
Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications
Fractional calculus-a different approach to the analysis of viscoelastically damped structures
AIAA J.
Damping description using fractional derivatives
Mech. Syst. Signal Process.
Fractional calculus in the transient analysis of viscoelastically damped structures
AIAA J.
Fractional fokker-planck equation for ultraslow kinetics
Europhys. Lett.
Cauchy and nonlocal multi-point problems for distributed order pseudo-differential equations: Part one
J. Anal. Appl
Cited by (21)
Weakly singular Volterra integral equation with combined logarithmic-power-law kernel: Analytical and computational consideration
2024, Applied Numerical MathematicsA novel method based on fractional order Gegenbauer wavelet operational matrix for the solutions of the multi-term time-fractional telegraph equation of distributed order
2024, Mathematics and Computers in SimulationNumerical solution of nonlinear stochastic differential equations with fractional Brownian motion using fractional-order Genocchi deep neural networks
2023, Communications in Nonlinear Science and Numerical SimulationFractional shifted Morgan–Voyce neural networks for solving fractal-fractional pantograph differential equations
2023, Chaos, Solitons and FractalsNumerical simulation of variable-order fractal-fractional delay differential equations with nonsingular derivative
2023, Engineering Science and Technology, an International JournalDistributed order Gauss-Quadrature scheme for distributed order fractional sub-diffusion model
2023, Chaos, Solitons and Fractals