A wavelet operational method for solving fractional partial differential equations numerically

Fractional calculus is an extension of derivatives and integrals to non-integer orders, and a partial differential equation involving the fractional calculus operators is called the fractional PDE. They have many applications in science and engineering. However not only the analytical solution existed for a limited number of cases, but also the numerical methods are very complicated and difficult. In this paper, we newly establish the simulation method based on the operational matrices of the orthogonal functions. We formulate the operational matrix of integration in a unified framework. By using the operational matrix of integration, we propose a new numerical method for linear fractional partial differential equation solving. In the method, we (1) use the Haar wavelet; (2) establish a Lyapunov-type matrix equation; and (3) obtain the algebraic equations suitable for computer programming. Two examples are given to demonstrate the simplicity, clarity and powerfulness of the new method.

[1]  Chih-Fan Chen,et al.  Numerical inversion of Laplace transform using Haar wavelet operational matrices , 2001 .

[2]  Mihály Kovács,et al.  Numerical solutions for fractional reaction-diffusion equations , 2008, Comput. Math. Appl..

[3]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[4]  P. N. Paraskevopoulos,et al.  Chebyshev series approach to system identification, analysis and optimal control , 1983 .

[5]  Chih-Fan Chen,et al.  Design of piecewise constant gains for optimal control via Walsh functions , 1975 .

[6]  Jelena Kovacevic,et al.  Wavelets and Subband Coding , 2013, Prentice Hall Signal Processing Series.

[7]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[8]  G. Jumarie,et al.  Modified Riemann-Liouville derivative and fractional Taylor series of nondifferentiable functions further results , 2006, Comput. Math. Appl..

[9]  I. Podlubny.  Fractional-Order Systems and -Controllers , 1999 .

[10]  G. Jumarie.  FOURIER'S TRANSFORM OF FRACTIONAL ORDER VIA MITTAG-LEFFLER FUNCTION AND MODIFIED RIEMANN-LIOUVILLE DERIVATIVE , 2008 .

[11]  W. Press,et al.  Numerical Recipes in C++: The Art of Scientific Computing (2nd edn)1 Numerical Recipes Example Book (C++) (2nd edn)2 Numerical Recipes Multi-Language Code CD ROM with LINUX or UNIX Single-Screen License Revised Version3 , 2003 .

[12]  Guy Jumarie,et al.  Fractional partial differential equations and modified Riemann-Liouville derivative new methods for solution , 2007 .

[13]  Spyridon G. Mouroutsos,et al.  Taylor series approach to system identification, analysis and optimal control , 1985 .

[14]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[15]  Spyridon G. Mouroutsos,et al.  The Fourier series operational matrix of integration , 1985 .

[16]  C. S. Lu.  Solution of the matrix equation AX+XB = C , 1971 .

[17]  I. Podlubny.  Fractional differential equations , 1998 .

[18]  Fan-Chu Kung,et al.  Laguerre series approach to the analysis of a linear control system incorporating observers , 1986 .

[19]  Rong-Yeu Chang,et al.  Legendre polynomials approximation to dynamic linear state equations with initial or boundary value conditions , 1984 .

[20]  Igor Podlubny,et al.  The Laplace Transform Method for Linear Differential Equations of the Fractional Order , 1997, funct-an/9710005.

[21]  M. Meerschaert,et al.  Finite difference approximations for two-sided space-fractional partial differential equations , 2006 .

[22]  C. F. Chen,et al.  Haar wavelet method for solving lumped and distributed-parameter systems , 1997 .

[23]  R. Haddad,et al.  Multiresolution Signal Decomposition: Transforms, Subbands, and Wavelets , 1992 .

[24]  Y. T. Tsay,et al.  Walsh operational matrices for fractional calculus and their application to distributed systems , 1977 .