Research Article  Open Access
ShuLi Mei, DeHai Zhu, "Interval Shannon Wavelet Collocation Method for Fractional FokkerPlanck Equation", Advances in Mathematical Physics, vol. 2013, Article ID 821820, 12 pages, 2013. https://doi.org/10.1155/2013/821820
Interval Shannon Wavelet Collocation Method for Fractional FokkerPlanck Equation
Abstract
Metzler et al. introduced a fractional FokkerPlanck equation (FFPE) describing a subdiffusive behavior of a particle under the combined influence of external nonlinear force field and a Boltzmann thermal heat bath. In this paper, we present an interval Shannon wavelet numerical method for the FFPE. In this method, a new concept named “dynamic interval wavelet” is proposed to solve the problem that the numerical solution of the fractional PDE is usually sensitive to boundary conditions. Comparing with the traditional wavelet defined in the interval, the Newton interpolator is employed instead of the Lagrange interpolation operator, so, the extrapolation points in the interval wavelet can be chosen dynamically to restrict the boundary effect without increase of the calculation amount. In order to avoid unlimited increasing of the extrapolation points, both the error tolerance and the condition number are taken as indicators for the dynamic choice of the extrapolation points. Then, combining with the finite difference technology, a new numerical method for the time fractional partial differential equation is constructed. A simple FokkerPlanck equation is taken as an example to illustrate the effectiveness by comparing with the GrunwaldLetnikov central difference approximation (GLCDA).
1. Introduction
Due to the fact that 1/ signal gains the increasing interests in the field of biomedical signal processing and engineering systems [1], the differential equations of fractional order appear more and more frequently in various research areas and engineering applications [2, 3]. As a matter of fact, the applications of fractional differential equations and their corresponding time series have been developed in various fields of sciences and technologies [4, 5] in recent years, ranging from computer science to physics [6, 7]. An effective and easytouse method for solving such equations is needed. However, known methods have certain disadvantages. Methods, described in detail in [3] for fractional differential equations of rational order, do not work in the case of arbitrary real order. On the other hand, there is an iteration method described in [8], which allows solution of fractional differential equations of arbitrary real order but it works effectively only for relatively simple equations, in addition to the series method. Up to now, most studies on the numerical methods for the fractional PDEs concentrate on the finite difference methods. Li [9] proposed an analytical method taking the fractal time series as the solution to a differential equation of fractional order or a response of a fractional system or a fractional filter driven with a white noise in the domain of stochastic processes and gave the exact solution of impulse response to a class of fractional oscillators [10]. According to this idea, Li and his coresearchers solved many problems in science and technology [11–14]. In addition, Wavelet numerical method is another way to get the solution of the fractional PDEs. In fact, the wavelet transform theory has been widely used in numerical analysis such as PDEsbased image processing [15–17], option pricing model [18], integrodifferential operators [19–23], and other nonlinear PDEs [24–28]. The wavelet functions possess many excellent numerical properties, such as orthogonality, interpolation, smoothness, and compact support, which are helpful in improving numerical accuracy and efficiency. In recent decades, many wavelets which have compact support, smoothness, and other properties have been constructed. Among these wavelets, Shannon wavelet is paid little attention in applications as it does not possess compact support property although it possesses orthogonality, smooth continuity, and analytical expression. Cattani studied the properties of the Shannon wavelet function, which possesses many advantages such as orthogonality, continuous and differentiable. It also has the advantage over the Hermite DAF in that it is an interpolating function, producing matrix equations that have the potential to be relatively sparse. In addition, the second order approximation of a C^{2}function, based on Shannon wavelet functions, is given [29]. The approximation is compared with the wavelet reconstruction formula and the error of approximation is explicitly computed [30].
A perceived disadvantage of the Shannon scaling function is that it tends to zero quite slowly as. A direct consequence of this is that when calculating the derivatives, a large number of the nodal values will contribute significantly. It is for this reason that Hoffman et al. [31] have suggested using the ShannonGabor wavelet, which introduces the Gaussian window function to improve the compact support property of Shannon wavelet function in required precision range. However, the presence of the Gaussian window destroys the orthogonal properties possessed by the Shannon wavelet, effectively worsening the approximation to a Dirac delta function.
Comparing with the common PDEs, the solutions of the fractional PDEs are more sensitive to the boundary condition. Using the wavelet transform defined in infinite domain to solve the engineering problems in finite interval, the wavelet transform coefficients at the boundary are usually very large. It will bring server boundary effect which affects the calculation accuracy and efficiency. Vasilyev and Paolucci [32] construct an interval wavelet using external wavelets, which can decrease the boundary effect to some extent. Based on the same principle, a more general construction method for the interval interpolation wavelet [33, 34] was given in the framework of generalized variational principle and has been widely used in many areas [35–37]. But the choice of parameter(that is the amount of the external collocation points) was not discussed in detail. It just points out that the value ofshould be taken between 1 and 3 based on experience. In fact, the value ofdepends on the smoothness and derivative of the approximated function at boundary points. That is, if the approximated function is the solution of the diffusion PDEs with respect to the time parameter, the value ofshould be taken dynamically. In addition, we should take into account that the impact of the external collocation points to the condition number of the system of the discretized algebraic equations. So, it is necessary to construct a dynamic interval wavelet in solving the PDEs with dynamic boundary conditions such as the fractional PDEs.
In this paper, a dynamic interval Shannon wavelet collocation method for the fractional FPDs is proposed. In this method, the relation between the parameterand the wavelet approximation error was discussed based on the interpolation error theory, and an adaptive choice procedure onwas constructed. Therefore, the socalled dynamic interval Shannon wavelet is constructed. Next, based on the GrünwaldLetnikov definition of the fractional order derivative, we construct a Shannon wavelet numerical method for the fraction FokkerPlanck equation.
2. Fractional FokkerPlanck Equation
The fractional FokkerPlanck equation has been used in many physical transport problems which take place under the influence of an external force field [2, 38].
In the presence of an external force field, the evolution of a test particle is usually described in terms of the FokkerPlanck equation (FPE)
which defines the probabilityof finding the particle at a certain position at a given time.denotes the mass of the diffusing particle,denotes the generalized diffusion coefficient with dimension, andis the generalized friction coefficient with dimension. The corresponding initial condition is and the boundary conditions are Equation (1) uses the RiemannLiouville fractional derivative of order, defined by whereis the gamma function.
According to the properties of the RiemannLiouville fractional derivative, it is easy to know that, if, (1) can be rewritten as follows: Metzler et al. [2] proposed three implicit approximations for solving (5) as follows.
The GrünwaldLetnikovexpansion and the backward Euler implicit approximation (GLBDIA) where,,, and andare positive integers. The local truncation error is.
approximation and the central difference implicit approximation (CDIA) The local truncation error is.
approximation and the backward difference implicit approximation (BDIA) The local truncation error is.
In fact, (6)–(8) are not perfect approximation as the boundary effect is not taken into account. So, it will introduce boundary effect in solving the PDEs with the Nuemann boundary conditions. It is well known that the finite difference method is equivalent to the FaberSchauder wavelet collocation method, so the construction method of the dynamic interval wavelet introduced in this paper can also be used to deal with the boundary problem in the finite difference method.
According to the Shannon sample theory, it can improve the calculation precision that combining the GrünwaldLetnikov expansion orapproximation of the fractional derivative in (5) with the Shannon scaling function as the weight function instead of the various difference operators as follows:
3. Construction of the Interval Interpolation Wavelet
3.1. Shannon Wavelet and ShannonGabor Wavelet
The representation of Shannon wavelet is based upon approximating the Dirac delta function as a bandlimited function and is given by and the ShannonGabor scaling function is defined as [17] where is the window size.
Consider a onedimensional function,. A discrete point sequence of the variableis defined as and the corresponding discrete point sequence of the scaling function andcan be defined, respectively, as where.
The first and second order derivatives ofat the discrete point are
And the first and second order derivatives ofat the discrete point are In fact, there is no difference between the construction method of the Interval Shannon wavelet and the interval ShannonGabor wavelet. So, we just take one uniform symbolas the representation of the Shannon wavelet and the ShannonGabor wavelet in the following.
3.2. Interval Interpolation Wavelet
According to the definition of the interval wavelet, the interval interpolation basis functions can be expressed as:
where, whereis the amount of the external collocation points, the amount of discrete points in the definition domain isis the definition domain of the approximated function.
Equations (16) and (17) show that the interval wavelet is derived from the domain extension. The supplementary discrete points in the extended domain are called external points. The value of the approximated function at the external points can be obtained by Lagrange extrapolation method. Using the interval wavelet to approximate a function, the boundary effect can be left in the supplementary domain; that is, the boundary effect is eliminated in the definition domain.
According to (16) and (17), the interval wavelet approximant of the functioncan be expressed as is the given value at the discrete point. At the external points, can be obtained by extrapolation; that isSo the interval wavelet approximant ofcan be rewritten as Let then andcorrespond to the left and the right external points, respectively. They are obtained by Lagrange extrapolation using the internal collocation points near the boundary. So, the interval wavelet’s influence on the boundary effect can be attributed to Lagrange extrapolation. It should be pointed out that we did not care about the reliability of the extrapolation. The only function of the extrapolation is enlarging the definition domain of the given function which can avoid the boundary effect occurred in the domain. Therefore, we can discuss the choice ofby means of Lagrange innerand extrapolation error polynomial as follows: Equation (23) indicates that the approximation error is both related to the smoothness and the gradient of the original function near the boundary. Setting differentcan satisfy the error tolerance.
3.3. Adaptive Interval Interpolation Wavelet
The interval interpolation wavelet is often used to solve the diffusion PDEs with Neumann boundary conditions. The smoothness and gradient of the PDE’s solution usually vary with the time parameter. If the parameteris a constant, we have to take a bigger value in order to obtain results with higher calculation precision. But the biggerusually introduces the famous Gibbs phenomena into the numerical solution, which usually results in the algorithm becoming invalid. In addition, the biggerwill bring much more calculation. To keep higher numerical precision and save calculation, the best way is to design a procedure thatcan vary with the curve’s smoothness and gradient dynamically.
In this dynamic procedure, the error estimation equation (23) can be taken as the criterion about. But in most cases, we cannot know the smoothness and the derivative’s order of the original function. This can be solved by substituting the difference coefficient for the derivative. This is coincident with the Newton interpolation equation which is equivalent with Lagrange interpolation equation. In addition, the Lagrange interpolation algorithm has no inheritance which is the key feature of Newton interpolation. So, the basis function has to be calculated repeatedly as interpolation points are added into the calculation, which increases the computation complexity greatly. In contracst to the Lagrange method, the advantage of Newton interpolation method is that the Newton divided difference form is employed, which can produce a mathematically equivalent result by using recurrence relations, which reduces the number of compute operation, especially the multiplication. So it is convenient using the Newton interpolation method to construct the dynamic procedure.
3.3.1. Newton Interpolation
The expression of Newton interpolation can be written as Substituting the Newton interpolation instead of the Lagrange interpolation into (22) can be rewritten as where
3.3.2. Relation between the Newton Interpolation Error and the Choice of
It is well known that the Newton interpolation is equivalent to the Lagrange interpolation. The corresponding error estimation can be expressed as And the simplest criterion to terminate the dynamic choice onis (is the absolute error tolerance). Obviously, it is difficult to definewhich should meet with the precision requirement of all approximated curves. In fact, the difference coefficient ) can be used directly as the criterion; that is As mentioned above, once the curves with lower order smoothness are approximated by higher order polynomial expression, the errors will become bigger on the contrary. In fact, even if theis infinite, the computational precision cannot be satisfied except by increasing computational complexity. To avoid this, we design the termination procedure of dynamic choice aboutas follows: If , then else if or , then else if or , then …
3.3.3. and the Condition Number of the System of Algebraic Equations
In the field of numerical analysis, the condition number of a function with respect to an argument measures how much the output value of the function can change for a small change in the input argument. This is used to measure how sensitive a function is to changes or errors in the input and how much error in the output results from an error in the input. It is no doubt that the choice ofcan change the condition number of the system of algebraic equations discretized by the wavelet interpolation operator or the finite difference method. Therefore, the choice ofshould take the condition number into account. In fact, if the condition number, then you may lose up todigits of accuracy on top of what would be lost to the numerical method due to loss of precision from arithmetic methods [24]. According to the general rule of thumb, the choice should follow the rule as follows:
3.3.4. Relation between and Computation Complexity
The computational complexity of interpolation calculation is not proportional to the increasing points. The former is mainly up to the computation amount of and the derivative operations. Obviously, according to (5), the increase in computational complexity iswhen the number of extension pointsincreases by 1. But the computational complexity of adaptively increasing collocation points is related to the different wavelet functions. For the wavelet with compact support property such as Daubechies wavelet and Shannon wavelet, the value ofis impossible to be infinite. For Haar wavelet which has no smoothness property,can be taken as 0 at most since it need not to be extended. For FaberSchauder wavelet,can be taken as 1 at most. For Daubechies wavelet,can be taken as different values according to the order of vanishing moments, but it must be finite. For the wavelets without compact support property,can take value dynamically, such as Shannon wavelet. The computational complexity of increasing points is mainly up to the wavelet function of itself.
4. Numerical Results and Discussion
Fractional FokkerPlanck equation is a typical fractional PDE, which is often used to describe a subdiffusive behavior of a particle under the combined influence of external nonlinear force field, and a Boltzmann thermal heat bath. This section considers the accuracy and efficiency of the proposed method for a fractional FokkerPlanck equation. Comparisons are made with results obtained with Chen’s finite difference approximations and the exact analytic solution.
It has been pointed out that the finite difference approximation formats proposed in [2] are not perfect as they do not take the boundary problems into account. In this section, we take the GrünwaldLetnikovexpansion and the central difference implicit approximation (GLCDIA) to solve the example. That is, According to the wavelet collocation method, the fractional FokkerPlanck equation can be approximately represented as where . Let Then the system of (31) can be expressed as the matrix format: Next, we will discuss the precision of the method proposed in this paper with numerical experience. Consider the FokkerPlanck equation as follows: with the initial condition and the boundary conditions The exact analytic solution is
All the comparisons in this section are made qualitatively by comparing the calculation precision in the same time step and space mesh grid size. The first measure of error is given by which provides a measure of the accuracy of the solution near the boundary. The second measure of erroris given by which provides a general measure of the accuracy of the solution over the main body of the distribution and was often used to investigate the accuracy of the FEM.
The comparisons between the static interval ShannonGabor wavelets withandare showen in Figure 1. The boundary effect of the interval wavelet with(Figure 1(a)) is almost eliminated compared to (Figure 1(b)). FFPE is a 2order PDEs with respect to, sois the necessary condition for the interval wavelet satisfying the requirement of FFPE. We also noticed that the condition number of FFPE from the Table 1 that the condition number ofincreases more rapid thanwith the increase ofand the decrease of . It has been mentioned in Section 2 that the larger condition number can decrease the calculation precision greatly. This also can be illustrated in Figure 2. The condition number in Figure 2(a) is greatly larger than in Figure 2(b), although the approximation ofis better than. The former has failed to solve FFPE obviously. In fact, this explained the reason why we construct the dynamic interval wavelet.

(a) , condition number is 76.7074, , and
(b) , condition number is 31.2950, , and
(a) , and
(b) , and
The numerical errors comparisons among the dynamic, static interval wavelet method and the interval finite difference method are showen in Figure 3. The result also can be illustrated in Table 2.

(a) Adapt interval wavelet collocation method
(b) Intervalfinite difference method
(c) Interval wavelet collocation method with constant
The robustness of the dynamic interval wavelet collocation method (DIWCM) is the best compared to the interval FDM and the static interval WCM, as it avoids both of the larger condition number and the error of the approximation simultaneity. The varied process of is showen in Table 3. It shows that the value of is fixed atafter a short time of vibration. This reflects the properties of the FFPE to some extent.

In addition, it also has to be noticed that we can get the higher precision solution with the interval finite difference method (FDM) as the amount of the collocation points decreases (Figure 4). It is well known that increasing the collocation points can impove the approximation although it can increase the condition number in FFPE. In fact, it profits from the smoothness of the solution, which would not work in solving the nonlinear problems.
(a) , , , and condition number is 39.1651
(b) , , , and condition number is 1.2202
All above numerical experiments are done with the ShannonGabor wavelet. It is well known that the presence of the Gaussian window destroys the orthogonal properties possessed by the Shannon wavelet, effectively worsening the approximation efficiency to a Dirac delta function. Comparing with the Shannon wavelet collocation method (Figure 5), the ShannonGabor wavelet numerical method has higher precision and more complicated calculation amount. But it is showen in Figure 5 that dynamic interpolation wavelet construction scheme can be applied to both of the ShannonGabor wavelet and the Shannon wavelet. As a matter of fact, the dynamic scheme is designed for the interpolation wavelet, which has no connection with certain concrete wavelet function.
(a) Static interval Shannon wavelet collocation method, , and ,
(b) Dynamic interval Shannon wavelet collocation method, , , and
5. Conclusions
In solving the fractional FokkerPlanck equations, there are two factors related to the choice of. The first factor is the condition number, which relates to the parameters , and the time step. The largercan decrease the calculation precision greatly. Another factor is the approximation of the function and its derivatives, especially near the boundary. Using the interval wavelet with constant to solve the fraction FokkerPlanck equations cannot eliminate the boundary effect completely as the condition number is sensitive to the parameter. With regard to the accuracy and time complexity of the solution in comparison with those obtained with other algorithms, the dynamic interval wavelet onconstructed in this paper is more reasonable. The numerical experiments illustrate that it is necessary to construct the dynamic interval wavelet collocation method for the fractional PDEs.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (no. 41171337) and the National High Technology Research and Development Program of China (no. 2006AA10Z235).
References
 M. Li and W. Zhao, “On $1/f$ noise,” Mathematical Problems in Engineering, vol. 2012, Article ID 673648, 23 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Metzler, E. Barkai, and J. Klafter, “Anomalous diffusion and relaxation close to thermal equilibrium: a fractional fokkerplanck equation approach,” Physical Review Letters, vol. 82, no. 18, pp. 3563–3567, 1999. View at: Publisher Site  Google Scholar
 S. Chen, F. Liu, P. Zhuang, and V. Anh, “Finite difference approximations for the fractional fokkerplanck equation,” Applied Mathematical Modelling, vol. 33, no. 1, pp. 256–273, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. G. Bakhoum and C. Toma, “Modeling transitions in complex systems by multiplicative effect of temporal patterns extracted from signal flows,” Mathematical Problems in Engineering, vol. 2012, Article ID 409856, 11 pages, 2012. View at: Publisher Site  Google Scholar
 C. Toma, “Advanced signal processing and command synthesis for memorylimited complex systems,” Mathematical Problems in Engineering, vol. 2012, Article ID 927821, 13 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. G. Bakhoum and C. Toma, “Transient aspects of wave propagation connected with spatial coherence,” Mathematical Problems in Engineering, vol. 2013, Article ID 691257, 5 pages, 2013. View at: Publisher Site  Google Scholar
 S. X. Hu, Z. W. Liao, and W. F. Chen, “Sinogram restoration for lowdosed xray computed tomography using fractionalorder peronamalik diffusion,” Mathematical Problems in Engineering, vol. 2012, Article ID 391050, 13 pages, 2012. View at: Publisher Site  Google Scholar
 J. H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Li, “Fractal time series—a tutorial review,” Mathematical Problems in Engineering, vol. 2010, Article ID 157264, 26 pages, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Li, S. C. Lim, and S. Y. Chen, “Exact solution of impulse response to a class of fractional oscillators and its stability,” Mathematical Problems in Engineering, vol. 2011, Article ID 657839, 9 pages, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Li, “Approximating ideal filters by systems of fractional order,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 365054, 6 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Li and W. Zhao, “Solving Abel's type integral equation with Mikusinski's operator of fractional order,” Advances in Mathematical Physics, vol. 2013, Article ID 806984, 4 pages, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Li, C. Cattani, and S. Y. Chen, “Viewing sea level by a onedimensional random function with long memory,” Mathematical Problems in Engineering, vol. 2011, Article ID 654284, 13 pages, 2011. View at: Publisher Site  Google Scholar
 M. Li and W. Zhao, “Visiting power laws in cyberphysical networking systems,” Mathematical Problems in Engineering, vol. 2012, Article ID 302786, 13 pages, 2012. View at: Publisher Site  Google Scholar
 I. Daubechies and G. Teschke, “Variational image restoration by means of wavelets: simultaneous decomposition, deblurring, and denoising,” Applied and Computational Harmonic Analysis, vol. 19, no. 1, pp. 1–16, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. L. Mei, “Construction of target controllable image segmentation model based on homotopy perturbation technology,” Abstract and Applied Analysis, vol. 2013, Article ID 131207, 8 pages, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S.L. Mei, “HPMbased dynamic wavelet transform and its application in image denoising,” Mathematical Problem in Engineering, vol. 2013, Article ID 309418, 10 pages, 2013. View at: Publisher Site  Google Scholar
 H. Yan, “Adaptive wavelet precise integration method for nonlinear blackscholes model based on variational iteration method,” Abstract and Applied Analysis, vol. 2013, Article ID 735919, 6 pages, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. W. Liu, “Interval wavelet numerical method on fokkerplanck equations for nonlinear random system,” Advances in Mathematical Physics, vol. 2013, Article ID 651357, 7 pages, 2013. View at: Publisher Site  Google Scholar
 C. Cattani and A. Ciancio, “Separable transition density in the hybrid model for tumorimmune system competition,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 610124, 6 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Cattani, A. Ciancio, and B. Lods, “On a mathematical model of immune competition,” Applied Mathematics Letters, vol. 19, no. 7, pp. 678–683, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Cattani, “Fractional calculus and shannon wavelet,” Mathematical Problems in Engineering, vol. 2012, Article ID 502812, 26 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Cattani, “Shannon wavelets for the solution of integrodifferential equations,” Mathematical Problems in Engineering, vol. 2010, Article ID 408418, 22 pages, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Y. Xing, “Waveletbased homotopy analysis method for nonlinear matrix system and its application in Burgers equation,” Mathematical Problems in Engineering, vol. 2013, Article ID 982810, 7 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 S.L. Pang, “Wavelet numerical method for nonlinear random system,” Transactions of the Chinese Society of Agricultural Machinery, vol. 38, no. 3, pp. 168–170, 2007. View at: Google Scholar
 S. M. Bibic, “Harmonic wavelet analysis—connection coefficients for nonlinear PDE,” UPB Scientific Bulletin, Series A, vol. 73, no. 1, pp. 27–36, 2011. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 M. Bieri, “A sparse composite collocation finite element method for elliptic SPDEs,” SIAM Journal on Numerical Analysis, vol. 49, no. 6, pp. 2277–2301, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. A. Libre, A. Emdadi, E. J. Kansa, M. Shekarchi, and M. Rahimian, “A fast adaptive wavelet scheme in RBF collocation for nearly singular potential PDEs,” Computer Modeling in Engineering and Sciences, vol. 38, no. 3, pp. 263–284, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 C. Cattani, “Second order shannon wavelet approximation of C^{2}functions,” UPB Scientific Bulletin, Series A, vol. 73, no. 3, pp. 73–84, 2011. View at: Google Scholar  MathSciNet
 C. Cattani and L. M. S. Ruiz, “Discrete differential operators in multidimensional haar wavelet spaces,” International Journal of Mathematics and Mathematical Sciences, vol. 2004, no. 44, pp. 2347–2355, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. K. Hoffman, G. W. Wei, D. S. Zhang, and D. J. Kouri, “ShannonGabor wavelet distributed approximating functional,” Chemical Physics Letters, vol. 287, no. 12, pp. 119–124, 1998. View at: Google Scholar
 O. V. Vasilyev and S. Paolucci, “A dynamically adaptive multilevel wavelet collocation method for solving partial differential equations in a finite domain,” Journal of Computational Physics, vol. 125, no. 2, pp. 498–512, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. L. Mei, H. L. Lv, and Q. Ma, “Construction of interval wavelet based on restricted variational principle and its application for solving differential equations,” Mathematical Problems in Engineering, vol. 2008, Article ID 629253, 14 pages, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. L. Mei, C. J. Du, and S. W. Zhang, “Asymptotic numerical method for multidegreeoffreedom nonlinear dynamic systems,” Chaos, Solitons and Fractals, vol. 35, no. 3, pp. 536–542, 2008. View at: Publisher Site  Google Scholar
 S. L. Mei, Q. S. Lu, S. W. Zhang, and L. Jin, “Adaptive interval wavelet precise integration method for partial differential equations,” Applied Mathematics and Mechanics, vol. 26, no. 3, pp. 364–371, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. L. Mei, C. J. Du, and S. W. Zhang, “Linearized perturbation method for stochastic analysis of a rill erosion model,” Applied Mathematics and Computation, vol. 200, no. 1, pp. 289–296, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S.L. Mei, “HAMbased adaptive multiscale meshless method for Burgers equation,” Journal of Applied Mathematics, vol. 2013, Article ID 248246, 10 pages. View at: Publisher Site  Google Scholar
 T. D. Frank, “Autocorrelation functions of nonlinear fokkerplanck equations,” European Physical Journal B, vol. 37, no. 2, pp. 139–142, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 ShuLi Mei and DeHai Zhu. 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.