`International Journal of Differential EquationsVolume 2010 (2010), Article ID 968186, 13 pageshttp://dx.doi.org/10.1155/2010/968186`
Research Article

## The Use of Fractional B-Splines Wavelets in Multiterms Fractional Ordinary Differential Equations

School of Mathematical and Computer Sciences, Fuzhou University, Fuzhou 350002, China

Received 31 July 2009; Revised 2 November 2009; Accepted 4 November 2009

Copyright © 2010 X. Huang and X. Lu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We discuss the existence and uniqueness of the solutions of the nonhomogeneous linear differential equations of arbitrary positive real order by using the fractional B-Splines wavelets and the Mittag-Leffler function. The differential operators are taken in the Riemann-Liouville sense and the initial values are zeros. The scheme of solving the fractional differential equations and the explicit expression of the solution is given in this paper. At last, we show the asymptotic solution of the differential equations of fractional order and corresponding truncated error in theory.

#### 1. Introduction

Recently,there have been several schemes devoted to the solution of fractional differential equations. These schemes can be broadly classified into two classes, numerical and analytical ([1]). As we know, with the help of some special functions, such as Mittag-Leffler function and Green function, Miller and Ross have obtained the explicit representations of solutions of some classes of homogeneous linear fractional differential equations (FDEs) in [2]; through the use of the technique of Laplace and Fourier Transform, the analytical solutions have been given by Podlubny in [3]. The numerical scheme we have encountered can be divided into two groups. In the first group, the solution is approximated over the entire domain using approximating functions such as polynomials and orthogonal functions. In the second group, the entire domain is divided into several small domains like in a finite element technique, and the solution is obtained for variables at the node points ([1]).As Edwards et al. declaimed in [4], several forms of fractional differential equations have been proposed in standard models, and there has been significant interest in developing numerical schemes for their solution. Thus, several papers have been presented in dealing with approximate numerical techniques for FDEs. Among these, the papers of Diethelm, Edwards, Ford, Freed, and Simpson are noteworthy (see, e.g., [1, 46]). One of these schemes which should be mentioned is the use of a Predictor-Corrector (more precisely, PECE) method in [6]. In order to obtain higher precision, they have replaced the PECE method by a method with in [7]. In particular, the PECE method is an important numerical scheme which has been applied in many fields; for example, Yang and Liu have applied the PECE method for simulating fraction order dynamical control system in [8]. In addition, several numerical schemes have also been proposed by other authors (see [9, 10]). In this paper, we present a new scheme which contains the features of both two groups by using the fractional B-splines wavelet. In paper [11], Unser and Blu have proved that the fractional B-splines generates a valid multiresolution analyses of for , which means that the orthogonal fractional B-splines could be obtained by applying the standard technique described in [12]. And they also obtained that the fractional B-splines deriving the asymptotic development of the have a fractional order of approximation.

In fact, as the theories of wavelets analyses improve day by day, the wavelet has become a powerful mathematical tool which widely used in signal processing, image compression and enhancement, pattern recognition, control systems, and other fields in the past two decades. But almost no papers or books have applied the theories of wavelets to solve the fractional differential equations. And our fundamental purpose of this paper is applying the fractional B-Splines wavelets to prove the existence and uniqueness of the solution of the nonhomogeneous linear fractional differential equations (also so-called linear multiterm fractional differential equations) with its initial conditions. Let us begin to discuss the solution of multiterms fractional ordinary differential equations with the following form: For convenience, we consider the initial values: where are taken in the Riemann-Liouville sense, , and . The function belongs to the space ; without loss of generality, in this paper, we consider the interval as ,.

The plan of this paper is as follows. In Section 2, we recall the definitions of fractional derivative and integral and related properties which will be used in this paper, give the representation of Mittag-Leffler function and generalized Mittag-Leffler function, and then introduce the fractional B-splines and some related properties of wavelets. In Section 3, by applying the technique of the Laplace Transforms, and considering the proprieties of the generalized Mittag-Leffler function, we prove the lemma of the differential equations of arbitrary positive real order, which make sure the solution belongs to the space ; by virtue of the solution that can be expressed as the form of wavelet series and the basis function that is the orthogonal fractional B-splines wavelet which yields the Riesz basis for the space ,we can prove the uniqueness of the coefficients of the representation of the solution, which gain the uniqueness of the solution of the fractional differential equations and validate the representation of solution. Thus, we have finished the proof of the theorem. In Section 4, the asymptotic solution of the differential equations of fractional order and corresponding truncated error will be discussed.

The present paper is essentially based on the works of the Unser and Blu in [11] and Podlubny in [3], to which more general classes of the fractional differential equations we shall refer in the following research. For more related review of fractional B-spline wavelet, see the papers of the Unser and Blu ([11, 13, 14]), which have discussed some more important wavelet properties, such as Riesz bounds and two-scale relation.

#### 2. Preliminary and Definitions

##### 2.1. Definitions

We may recall the definition of the Left Riemann-Liouville differential operators of arbitrary order , which take the form where is the integer defined by (see [3]), and is gamma function. Similarly, the left Riemann-Liouville integral operators of order is defined as And then we should give the following expression for the Laplace transform of the Riemann-Liouville differential operators of the order , which is

Let us now introduce the case of the Caputo differential operators of arbitrary order , which is defined as And the formula of its Laplace transform can be expressed as

##### 2.2. Mittag-Leffler Functions and Generalized

The Mittag-Leffler functions and its generalized forms have played a special role in solving the fractional differential equations. In this section, we just give the definition of the following series of representation of the Mittag-Leffler function with , which validates in the whole complex plane: And for the generalized Mittag-Leffler function, we use the following definition: Let us now consider the Laplace transforms of the function (see [3]), which is defined by

##### 2.3. Fractional B-Splines Wavelet

Splines have had a significant impact on the early development of the theory of the wavelet transform (see [13]). And Unser and Blu have first mentioned the fractional B-Splines in [13], who extended Schoenberg's family of polynomial splines to all fractional degrees and defined the fractional causal B-splines by taking the fractional difference of the one-sided power function: where the one-side power function is defined as follows: Then we introduce the fractional B-splines autocorrelation sequence as follows: From [13], we know the explicit form of the fractional B-splines wavelet: which yields a Riesz basis for . Thus, using the standard orthogonalize technique described in [12], we have where is the Fourier Transform of the function of . Thus we obtain the orthogonal fractional B-splines which also yields a Riesz basis for . According to the theories of wavelet analyses, the function can be expressed as where the coefficients are constants.

#### 3. The Existence and Uniqueness of the Solution ofMultiterms Fractional Ordinary Differential Equations

In this section, we will prove the existence and uniqueness properties of the solutions of the nonhomogeneous linear differential equations of arbitrary real order .

##### 3.1. An Important Lemma

Under the hypothesis of the existence for the solutions of (1.1), we have the following lemma.

Lemma 3.1. Let and , and , then the solution of the initial value problem (1.1) and (1.2) is also in .

Proof. In order to prove this lemma, we divide the proof into two steps.
Firstly, we consider the case of , then (1.1) can be rewritten as Taking the transform of (3.1) to both sides, we obtain By virtue of the functions and that is a bounded linear operator (see [15]), we can easily derive the function .
And then, let us consider the case of .
To (1.1), using the formula (2.2) and taking the Laplace transform to both sides, we conclude that and then To prove the functions , we change the functions into the following form: Then replacing the factor with in (3.5), and expanding the second factor into a formal of series, we have Using the method described in [3], we obtain the expression of , which is Substituting in (3.4), considering the formula (2.8), and taking the inverse Laplace transform term-by-term, we obtain that the analytical solution of initial values problem (1.1) and (1.2) in the following form: with where is the multinomial coefficient, the representation of is the convolution of functions and , and is the derivative of the Mittag-Leffler function with parameters and (see [16]).
From the representation (3.9), the function is bounded in , and the index satisfies ; so we can easily gain that the function is bounded in . And then it yields that the functions , and so we have finished the proof of the lemma.

##### 3.2. The Proof of Existence and Uniqueness

Considering the solution of the fractional differential equation and the fractional B-spline wavelet which generates a Riesz basis for , we can prove the following theorem.

Theorem 3.2. Let and , and , then the initial value problems (1.1) and (1.2) have a unique solution. Further more, the solution has explicit representation of fractional B-splines wavelets series.

Proof. From Lemma 3.1, we have obtained the functions , which can be expressed as where the index satisfies .
Then substituting (2.14) into (1.1), and taking the Fourier transform to both sides, we obtain where , and (3.12) is equivalent to the following form: And then, by taking the inverse Fourier transform to (3.13), we have Because the function is the orthogonal fractional B-splines, hence, the representation of the coefficient is defined by where is the inner product. Then substituting (3.15) in (2.14), we obtain the solution of the initial value problems (1.1) and (1.2).
Let us suppose that the initial value problems (1.1) and (1.2) have another solution, which can be expressed as Utilizing the similar scheme to , we acquire the representation of , which is Taking (3.17) into (3.16), we obtain the representation of . Obviously, we have it derives that , which means that the solution of the initial value problems (1.1) and (1.2) is unique. Then substituting with its coefficient in (1.1), we can easily check that the equation is correct. It indicates that the function is a solution of the differential equation (1.1), which yields the existence of solution. Finally, we have completed the proof of the solution of existence and uniqueness of the nonhomogeneous linear differential equations of arbitrary order.

#### 4. The Asymptotic Solution of the Fractional Differential Equations and Error Estimation

##### 4.1. The Asymptotic Solution of the Fractional Differential Equations

The purpose of this section is to discuss the case of in (1.1), which appeared in most practical applications we have encountered, and show the asymptotic solution of the fractional differential equations and error estimation. In fact, the case of in Theorem 3.2 is equivalent to the following corollary.

Corollary 4.1. Let , and , then the nonhomogeneous linear fractional differential equations on the initial values where,  ,  ,  , have a unique solution.

In Section 3, we have proved the existence and uniqueness of the nonhomogeneous linear fractional differential equations of arbitrary positive real order. Obviously, it also satisfies the case of positive rational order. To (1.1), let be the least common multiple of the denominators of , so , where . And to (4.1), let and , then considering the conditions of , where . Hence, the initial value problems (4.1) and (4.2) have been changed into of the problems (1.1) and (1.2) with the order , respectively, which means that Theorem 3.2 in the case of is equivalent to Corollary 4.1.

Similarity to the scheme of Theorem 3.2, we can easily prove Corollary 4.1 and give the representation of explicit solution of the initial value problems (4.1) and (4.2), which can be defined by where , and is the inverse of the Laplace Transform of the function , where To obtain the asymptotic solution of (4.1), we will give the explicit formulation of the function . Let and ; thus the is a polynomial of the degree . Moreover, we suppose that are distinct zeros of with the order of , respectively; by applying the theories of polynomials, the can be rewritten in the following form where , Then, substituting the function in (4.2) and expanding it to the sum of partial fractions, we have where the coefficients are constants.

Thus, taking the inverse Laplace Transform to (4.6) and using the formulation (2.8), we obtain explicit formulation of the function , defined by According to Unser and Blu in [13], the wavelet base generated by fractional B-splines wavelet is denseness of the representation in . Let be the inverse of Fourier Transform of , where With the help of the theories of wavelet analyses, we know that the functions are an orthogonal fractional B-splines wavelet and dense in with . Thus the solutions of initial value problems (4.1) and (4.2) have the following form: Then substituting in (4.1), we obtain the representation of the coefficients , defined by Finally, by combining (4.9) and (4.10), the obtained function is the asymptotic solution of initial value problems (4.1) and (4.2) that we are looking for.

##### 4.2. Order of the Error Estimation

To estimate the error of asymptotic solution of initial value problems (4.1) and (4.2), we introduce the following properties of the fractional B-spines , which have been proved by Unser and Blu (see [13, Theorem ]).

Proposition 4.2. For all , we have when tends to .

To (4.8), the function is -periodic and symmetric, and so we can restrict its study to . In particular, one has since is strictly decreasing over (see [13]).

Combining (4.12) and (4.8) and taking the inverse of Fourier Transform, we obtain Then according to (4.10), (4.11), and (4.13), the asymptotic solution can be defined by To calculate the truncated error of the asymptotic solution, let be the truncated sum corresponding to asymptotic solution (4.8), where Thus the truncated error will be obtained as follows: where the function is bounded in and ; with the help of (4.12), the inequality of (4.16) can be amplified, which means that Then substituting (4.10) in (4.16), we have To the right side of (4.17), we divided the representation into three parts for discussion. Firstly, because the function is bounded in and , there exists a constant which holds Secondly, we consider the series function which is convergent in ; thus there exists a constant which satisfies At last, we denote a constant which defined by Then combining (4.19), (4.21), and (4.22), the inequality of (4.18) can be amplified to the following form: when tends to .

Thus, from the above discussion, it is evident to derive that the truncated error converges as , where . Note that we should choose a suitable which should be much more greater than , where, . Finally, we have obtained the error order of the asymptotic truncated sum in theory.

Remark 4.3. Noting that the process of the proof of the existence and uniqueness of the solution of the initial value problems (1.1) and (1.2) in the case of Riemann-Liouville fractional differential operator, and considering the formula of the Laplace transform of Caputo differential operator, we can be easily replaced the case of Riemann-Liouville fractional differential operator by Caputo sense with its initial values. It means the following corollary is correct.

Corollary 4.4. Let and , and , then the multiterms fractional ordinary differential equations with its initial values where , have a unique solution.

#### 5. Conclusion

In this paper, we have proved the existence and uniqueness of the solution of the differential equations of arbitrary positive real order. And the representation of the solution of (1.1) has been given in the process of proof. We have obtained the asymptotic solution of the differential equations of fractional order and corresponding error estimation. The most notable feature is the order of the asymptotic truncated error, namely, , which is effective to calculate the numerical solution of (4.1). In particular, the case of Riemann-Liouville differential operator is replaced by Caputo sense with its initial values in Corollary 4.4. Similarly to the scheme of the proof as showen in the paper, by considering the relationship between Riemann-Liouville and Caputo differential operator, we can easily complete the proof of the corollary in the case of the Caputo differential operator, which shows that the method we have discussed can be applied more widely.

#### References

1. P. Kumar and O. P. Agrawal, “Numerical scheme for the solution of fractional differential equations of order greater than one,” Journal of Computational and Nonlinear Dynamics, vol. 1, no. 2, 8 pages, 2006.
2. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 1993.
3. I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
4. J. T. Edwards, N. J. Ford, and A. C. Simpson, “The numerical solution of linear multi-term fractional differential equations; systems of equations,” Journal of Computational and Applied Mathematics, vol. 148, no. 2, pp. 401–418, 2002.
5. K. Diethelm and N. J. Ford, “Analysis of fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 265, no. 2, pp. 229–248, 2002.
6. K. Diethelm, N. J. Ford, and A. D. Freed, “A predictor-corrector approach for the numerical solution of fractional differential equations,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 3–22, 2002.
7. K. Diethelm, “Efficient solution of multi-term fractional differential equations using $P{\left(EC\right)}^{\text{m}}E$ methods,” Computing, vol. 71, no. 4, pp. 305–319, 2003.
8. C. Yang and F. Liu, “A computationally effective predictor-corrector method for simulating fractional order dynamical control system,” The ANZIAM Journal, vol. 47, pp. C168–C184, 2005.
9. K. Diethelm and N. J. Ford, “Multi-order fractional differential equations and their numerical solution,” Applied Mathematics and Computation, vol. 154, no. 3, pp. 621–640, 2004.
10. V. Daftardar-Gejji and A. Babakhani, “Analysis of a system of fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 293, no. 2, pp. 511–522, 2004.
11. M. Unser and T. Blu, “Construction of fractional spline wavelet bases,” in Wavelets Applications in Signal and Image Processing VII, vol. 3813 of Proceedings of SPIE, pp. 422–431, Denver, Colo, USA, July 1999.
12. S. Mallat, A Wavelet Tour of Signal Processing, Academic Press, San Diego, Calif, USA, 2nd edition, 1998.
13. M. Unser and T. Blu, “Fractional splines and wavelets,” SIAM Review, vol. 42, no. 1, pp. 43–67, 2000.
14. M. Unser, A. Aldroubi, and M. Eden, “A family of polynomial spline wavelet transforms,” Signal Processing, vol. 30, pp. 141–162, 1993.
15. V. J. Ervin and J. P. Roop, “Variational formulation for the stationary fractional advection dispersion equation,” Numerical Methods for Partial Differential Equations, vol. 22, no. 3, pp. 558–576, 2005.
16. F. Mainardi and R. Gorenflo, “On Mittag-Leffler-type functions in fractional evolution processes,” Journal of Computational and Applied Mathematics, vol. 118, no. 1-2, pp. 283–299, 2000.