Solving fractional integral equations by the Haar wavelet method
Introduction
Although the conception of the fractional derivatives was introduced already in the middle of the 19th century by Riemann and Liouville, the first work, devoted exclusively to the subject of fractional calculus, is the book by Oldham and Spanier [1] published in 1974. After that the number of publications about the fractional calculus has rapidly increased. The reason for this is that some physical processes as anomalous diffusion, complex viscoelasticity, behaviour of mechatronic and biological systems, rheology etc. cannot be described adequately by the classical models.
At the present time we possess several excellent monographs about fractional calculus for example the book [2] by Kilbas et al., to which is also included a rather large and up-to-date Bibliography (928 items). Because of the enormous number of the papers about this topic we shall cite here only some papers which are more close to subject of this paper.
In a number of papers fractional differential equations are discussed; mostly these equations are transformed to fractional Volterra integral equations. For solution different techniques, as Fourier and Laplace transforms, power spectral density, Adomian decomposition method, path integration etc., are applied.
One-dimensional fractional harmonic oscillator is analysed in [3], [4], [5], [6]. In [3], [4] the solution is obtained in terms of Mittag–Leffler functions using Laplace transforms; several cases of the forcing function equation are considered. In [5] the fractional equation of motion is solved by the path integral method. In [6] the case, where the fractional derivatives only slightly differ from the ordinary derivatives, is analysed. Fractional Hamilton’s equations are discussed in [7]. In [8] multiorder fractional differential equations are solved by using the Adomian decomposition. In several papers fractional chaotic systems are discussed. In [8] a three-dimensional fractional chaotic oscillator model is proposed. Chaotic dynamics of the fractionally damped Duffing equation is investigated in [9]. Two chaotic models for third-order chaotic nonlinear systems are analysed in [10].
It is somewhat surprising that among different solution techniques the wavelet method has not attained much attention. We found only one paper [11] in which the wavelet method is applied for solving fractional differential equations; for this purpose the Daubechies wavelet functions are used.
Among the different wavelet families mathematically most simple are the Haar wavelets. Due to the simplicity the Haar wavelets are very effective for solving differential and integral equations (see e.g. [12], [13], [14], [15]). Therefore the idea, to apply Haar wavelet technique also for solving problems of fractional calculus, arises. This is the main aim of the present paper.
The paper is organized as follows. Sections 2 About the fractional calculus, 3 Haar wavelets are preparative: in Section 2, basic equations of the fractional calculus are briefly reviewed, in Section 3, the Haar wavelet method is described. In Section 4, two error estimates for the results, obtained by the Haar wavelet method, are introduced. In Section 5, the wavelet solution for fractional Volterra integral equations is presented. The one-dimensional fractional harmonic vibrations are discussed in Section 6. The solution for the fractional Fredholm integral equation is presented in Section 7.
Section snippets
About the fractional calculus
Let us briefly consider some basic formulae about the fractional calculus.
The Riemann–Liouville fractional integrals of order are defined byandHere is the gamma function and the integer part of . The integrals (1), (2) are called left-sided and right-sided fractional integrals.
As to the fractional derivatives, then in this paper we shall use the Caputo derivatives defined as
Haar wavelets
Usually the Haar wavelets are defined for the interval . In this paper the more general case , is considered. Let us define the quantity , where J is the maximal level of resolution. We shall divide the interval into subintervals of equal length; each subinterval has the length . Next two parameters are introduced: the dilatation parameter j for which and the translation parameter (here the notation is introduced). The wavelet
Error estimates
It is essential to estimate the exactness of the obtained solutions, for this purpose in the following two error estimates are defined. Here we have to distinguish the two following situations.
- (i)
If the exact solution of the problem is known we shall calculate the differences and define the error estimates as (local estimate) or (global estimate).
- (ii)
Mostly the exact solution is unknown. For this case the following procedure is
Fractional Volterra integral equation
The fractional Volterra integral equation has the form [2]The kernel and the right-side function are given, is a real number. The value corresponds to the ordinary (nonfractional) Volterra equation.
According to the Haar wavelet method, the solution of (15) is sought in the form (8). Replacing (8), (15) and satisfying this equation in the collocation points we obtainHere the symbol
Fractional harmonic vibrations
Consider the equationwhere are prescribed constants, is the forcing term. To (23) belong initial conditions . If and we get the usual differential equation of the harmonic oscillator.
The symbols denote left-sided Caputo derivatives, which are defined by (3). Since in the present case the equation (23) gets the formThis is a
Fractional Fredholm integral equation
The usual Fredholm integral equation has the formIf we want to get its fractional analogue we must take into account the fact that in the case of non-integer the expression has sense only for . Therefore it seems reasonable to define the fractional Fredholm integral equation in the formIn (36) stand the left-sided and right-sided Riemann–Liouville integrals, respectively. As in Section 4
Conclusion
In the present paper three nonfractional equations are solved by the Haar wavelet method, these problems must be considered as examples for the recommended method of solution, since this approach can be easily carried over to solving some other problems, as e.g.
- (i)
first kind linear integral equations [12]:
- (ii)
integro–differential equations [12], [14];
- (iii)
weakly singular integral equations [12].
In the present paper only linear equations are considered, but the method is applicable also for nonlinear
Acknowledgement
Financial support from the Estonian Science Foundation under Grant ETF-6697 is gratefully acknowledged.
References (18)
- et al.
Dynamics of the fractional oscillator
Physica A
(2001) - et al.
Path integral representation of fractional harmonic oscillator
Physica A
(2006) On a almost free damping vibration equation using N-fractional calculus
J. Comp. Appl. Math.
(2002)- et al.
The Hamilton formalism with fractional derivatives
J. Math. Anal. Appl.
(2007) - et al.
Solving a multi-order fractional differential equation using Adomian decomposition
Appl. Math. Comp.
(2007) - et al.
Chaotic dynamics of the fractionally damped Duffing equation
Chaos Soliton Fract.
(2007) - et al.
Chaos in fractional-order autonomous nonlinear systems
Chaos Soliton Fract.
(2003) - et al.
Identification of fractional chaotic system parameters
Chaos Soliton Fract.
(2004) Haar wavelet method for nonlinear integro–differential equations
Appl. Math. Comp.
(2006)
Cited by (123)
A novel L1-Predictor-Corrector method for the numerical solution of the generalized-Caputo type fractional differential equations
2024, Mathematics and Computers in SimulationExtension of Haar wavelet technique for numerical solution of three-dimensional linear and nonlinear telegraph equations
2024, Partial Differential Equations in Applied MathematicsOn the numerical solution of fractional differential equations under white noise processes
2023, Probabilistic Engineering MechanicsSolving an inverse problem for a time-fractional advection-diffusion equation with variable coefficients by rationalized Haar wavelet method
2022, Journal of Computational ScienceFractional modelling and numerical simulations of variable-section viscoelastic arches
2021, Applied Mathematics and ComputationDynamics of flight of the fragments with higher order Haar wavelet method
2024, Proceedings of the Estonian Academy of Sciences