Numerical Schemes for Fractional Ordinary Differential Equations
Weihua Deng,Can Li
DOI: https://doi.org/10.5772/34965
2012-01-01
Abstract:Will-be-set-by-IN-TECHto make sure the discussed equation has a unique solution.These additional conditions, in many situations, describe some properties of the solution at the initial time (Heymans & Podlubny, 2005), the fractional derivative does not have convenient used physical meaning (there are already some progress in the geometric and the physical interpretation of the fractional calculus (Podlubny, 2002) and the physical interpretation of the initial conditions in terms of the Riemann-Liouville fractional derivatives of the unknown function has also been discussed in (Podlubny, 2002)).Just like the classic calculus and differential equations, the theories of fractional differentials, integrals and differential equations have been developing.With the development of the theories of fractional calculus, many research monographs are published, e.g., (Oldham & Spanier, 1974;Podlubny, 1999;Samko et al., 1993).In the literatures, several analytical methodologies, such as, Laplace transform, Mellin transform, Fourier transform, are restored to obtain the analytical solutions of the fractional equations by many authors (Metzler & Klafter, 2000;Podlubny, 1999;Samko et al., 1993;Zaslavsky, 2002;Mainrdi, 2008), however, similar to treating classical differential equations, they can mainly deal with linear fractional differential equation with constant coefficients.Usually, for nonlinear systems these techniques do not work.So in many cases the more reasonable option is to find its numerical solution.As is well known, the difficulty of solving fractional differential equations is essentially because fractional calculus are non-local operators.This non-local property means that the next state of a system not only depends on its current state but also on its historical states starting from the initial time.This property is closer to reality and is the main reason why fractional calculus has become more and more useful and popular.In other words, this non-local property is good for modeling reality, but a challenge for numerical computations.Much effort has been devoted during the recent years to the numerical investigations of fractional calculus and fractional dynamics of (1) (Lubich, 1985;1986;Podlubny, 1999).More recently, Diethelm et al successfully presented the numerical approximation of (1) using Adams-type predictor-corrector approach (Diethelm & Ford, 2002a) and the detailed error analysis of this method was given in (Diethelm et al., 2004).The convergent order of Diethelm's predictor-corrector approach was proved to be min(2, 1 + α) (Diethelm et al., 2004).Because of the non-local property of the fractional derivatives, the arithmetic complexity of their algorithm with step size h is O(h -2 ), whereas a comparable algorithm for a classical initial value problem only gives rise to O(h -1 ).To improve the accuracy and reduce the arithmetic complexity, some techniques such as the Richardson extrapolation, short memory principle and corresponding mixed numerical schemes are developed.In (Deng, 2007a), we present an improved version of the predictor-corrector algorithm with the accuracy increased to min(2, 1 + 2α) and half of the computational cost is reduced comparing to the original one in (Diethelm et al., 2004).Furthermore, we apprehend the short memory principle from a new viewpoint (Deng, 2007b); after using the nested meshes presented in (Ford & Simpson, 2001) and combining the short memory principle and the predictor-corrector approach, we minimize the computational complexity to O(h -1 log(h -1 )) at preserving the order of accuracy 2. This chapter briefly reviews the recent development of the predictor-corrector approach for fractional dynamic systems.The plan of this chapter is as follows.In Section 2, we briefly discuss the short memory principle and the nested meshes.In Section 3, the predictor-corrector schemes and its improved versions are presented, meanwhile the convergent order and arithmetic complexity are also proposed.In Section 4, we provide two 356 Numerical Modelling www.intechopen.com