An integral operational matrix based on Jacobi polynomials for solving fractional-order differential equations

This paper aims to construct a general formulation for the Jacobi operational matrix of fractional integral operator. Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. Therefore, a reliable and efficient technique for the solution of them is too important. For the concept of fractional derivative we will adopt Caputo’s definition by using Riemann–Liouville fractional integral operator. Our main aim is to generalize the Jacobi integral operational matrix to the fractional calculus. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

Chebyshev operational matrix method for solving multi-order fractional ordinary differential equations

The aim of this article is to present an analytical approximation solution for linear and nonlinear multi-order fractional differential equations (FDEs) by extending the application of the shifted Chebyshev operational matrix. For this purpose, we convert FDE into a counterpart system and then using proposed method to solve the resultant system. Our results in solving four different linear and nonlinear FDE, confirm the accuracy of proposed method.     

The convergence and stability analysis of the Jacobi collocation method for solving nonlinear fractional differential equations with integral boundary conditions

In this paper, we apply the Jacobi collocation method for solving nonlinear fractional differential equations with integral boundary conditions. Due to existence of integral boundary conditions, after reformulation of this equation in the integral form, the method is proposed for solving the obtained integral equation. Also, the convergence and stability analysis of the proposed method are studied in two main theorems. Furthermore, the optimum degree of convergence in the L₂ norm is obtained for this method. Furthermore, some numerical examples are presented in order to illustrate the performance of the presented method. Finally, an application of the model in control theory is introduced. Copyright © 2015 John Wiley & Sons, Ltd.     

A computational matrix method for solving systems of high order fractional differential equations

In this paper, we introduced an accurate computational matrix method for solving systems of high order fractional differential equations. The proposed method is based on the derived relation between the Chebyshev coefficient matrix A of the truncated Chebyshev solution u(t)$\mathbf{u}(\mathbf{t})$ and the Chebyshev coefficient matrix A^(ν)${\mathbf{A}}^{(\nu )}$ of the fractional derivative u^(ν)${\mathbf{u}}^{(\nu )}$. The fractional derivatives are presented in terms of Caputo sense. The matrix method for the approximate solution for the systems of high order fractional differential equations (FDEs) in terms of Chebyshev collocation points is presented. The systems of FDEs and their conditions (initial or boundary) are transformed to matrix equations, which corresponds to system of algebraic equations with unknown Chebyshev coefficients. The remaining set of algebraic equations is solved numerically to yield the Chebyshev coefficients. Several numerical examples for real problems are provided to confirm the accuracy and effectiveness of the present method.     

The general Jacobi matrix method for solving some nonlinear ordinary differential equations

In this paper, we obtain the approximate solutions for some nonlinear ordinary differential equations by using the general Jacobi matrix method. Explicit formulae which express the Jacobi expansion coefficients for the powers of derivatives and moments of any differentiable function in terms of the original expansion coefficients of the function itself are given in the matrix form. Three test problems are discussed to illustrate the efficiency of the proposed method.     

Bernoulli wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations

In this paper, a new numerical method for solving fractional differential equations is presented. The fractional derivative is described in the Caputo sense. The method is based upon Bernoulli wavelet approximations. The Bernoulli wavelet is first presented. An operational matrix of fractional order integration is derived and is utilized to reduce the initial and boundary value problems to system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique.     

Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations

In this paper, we state and prove a new formula expressing explicitly the derivatives of shifted Chebyshev polynomials of any degree and for any fractional-order in terms of shifted Chebyshev polynomials themselves. We develop also a direct solution technique for solving the linear multi-order fractional differential equations (FDEs) with constant coefficients using a spectral tau method. The spatial approximation with its fractional-order derivatives (described in the Caputo sense) are based on shifted Chebyshev polynomials T_L,n(x) with x ∈ (0, L), L > 0 and n is the polynomial degree. We presented a shifted Chebyshev collocation method with shifted Chebyshev–Gauss points used as collocation nodes for solving nonlinear multi-order fractional initial value problems. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed techniques and to compare with the existing results.     

Bernstein operational matrix of fractional derivatives and its applications

In this paper, Bernstein operational matrix of fractional derivative of order α in the Caputo sense is derived. We also apply this matrix to the collocation method for solving multi-order fractional differential equations. The numerical results obtained by the present method compares favorably with those obtained by various collocation methods earlier in the literature.     

Convergence analysis of the Jacobi spectral-collocation method for fractional integro-differential equations

We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L²-norm. The numerical examples are given to illustrate the theoretical results.     

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.     

A reliable algorithm of homotopy analysis method for solving nonlinear fractional differential equations

In this paper, based on the homotopy analysis method (HAM), a powerful algorithm is developed for the solution of nonlinear ordinary differential equations of fractional order. The proposed algorithm presents the procedure of constructing the set of base functions and gives the high-order deformation equation in a simple form. Different from all other analytic methods, it provides us with a simple way to adjust and control the convergence region of solution series by introducing an auxiliary parameter ?$?$. The analysis is accompanied by numerical examples. The algorithm described in this paper is expected to be further employed to solve similar nonlinear problems in fractional calculus.     

Numerical solution of multiterm variable‐order fractional differential equations via shifted Legendre polynomials

In this paper, shifted Legendre polynomials will be used for constructing the numerical solution for a class of multiterm variable‐order fractional differential equations. In the proposed method, the shifted Legendre operational matrix of the fractional variable‐order derivatives will be investigated. The fundamental problem is reduced to an algebraic system of equations using the constructed matrix and the collocation technique, which can be solved numerically. The error estimate of the proposed method is investigated. Some numerical examples are presented to prove the applicability, generality, and accuracy of the suggested method.     

A fractional‐order Jacobi Tau method for a class of time‐fractional PDEs with variable coefficients

This paper presents a shifted fractional‐order Jacobi orthogonal function (SFJF) based on the definition of the classical Jacobi polynomial. A new fractional integral operational matrix of the SFJF is presented and derived. We propose the spectral Tau method, in conjunction with the operational matrices of the Riemann–Liouville fractional integral for SFJF and derivative for Jacobi polynomial, to solve a class of time‐fractional partial differential equations with variable coefficients. In this algorithm, the approximate solution is expanded by means of both SFJFs for temporal discretization and Jacobi polynomials for spatial discretization. The proposed tau scheme, both in temporal and spatial discretizations, successfully reduced such problem into a system of algebraic equations, which is far easier to be solved. Numerical results are provided to demonstrate the high accuracy and superiority of the proposed algorithm over existing ones. Copyright © 2015 John Wiley & Sons, Ltd.     

The construction of operational matrix of fractional derivatives using B-spline functions

Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. For that reason we need a reliable and efficient technique for the solution of fractional differential equations. Here we construct the operational matrix of fractional derivative of order α in the Caputo sense using the linear B-spline functions. The main characteristic behind the approach using this technique is that it reduces such problems to those of solving a system of algebraic equations thus we can solve directly the problem. The method is applied to solve two types of fractional differential equations, linear and nonlinear. Illustrative examples are included to demonstrate the validity and applicability of the new technique presented in the current paper.     

Numerical solutions to initial and boundary value problems for linear fractional partial differential equations

In this article, Haar wavelets have been employed to obtain solutions of boundary value problems for linear fractional partial differential equations. The differential equations are reduced to Sylvester matrix equations. The algorithm is novel in the sense that it effectively incorporates the aperiodic boundary conditions. Several examples with numerical simulations are provided to illustrate the simplicity and effectiveness of the method.     

Haar wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations

Haar wavelet operational matrix has been widely applied in system analysis, system identification, optimal control and numerical solution of integral and differential equations. In the present paper we derive the Haar wavelet operational matrix of the fractional order integration, and use it to solve the fractional order differential equations including the Bagley-Torvik, Ricatti and composite fractional oscillation equations. The results obtained are in good agreement with the existing ones in open literatures and it is shown that the technique introduced here is robust and easy to apply.     

Application of generalized differential transform method to multi-order fractional differential equations

In a recent paper [Odibat Z, Momani S, Erturk VS. Generalized differential transform method: application to differential equations of fractional order, Appl Math Comput. submitted for publication] the authors presented a new generalization of the differential transform method that would extended the application of the method to differential equations of fractional order. In this paper, an application of the new technique is applied to solve fractional differential equations of the form y^(μ)(t)=f(t,y(t),y^(β₁)(t),y^(β₂)(t),…,y^(β_n)(t)) with μ>β_n>β_n-1>…>β₁>0, combined with suitable initial conditions. The fractional derivatives are understood in the Caputo sense. The method provides the solution in the form of a rapidly convergent series. Numerical examples are used to illustrate the preciseness and effectiveness of the new generalization.     

Analysis of differential equations of fractional order

This paper provides a robust convergence checking method for nonlinear differential equations of fractional order with consideration of homotopy perturbation technique. The differential operators are taken in the Caputo sense. Some theorems to prove the existence and uniqueness of the series solutions are presented. Results show that the proposed theoretical analysis is accurate.     

Fractional‐order orthogonal Bernstein polynomials for numerical solution of nonlinear fractional partial Volterra integro‐differential equations

In this paper, a new two‐dimensional fractional polynomials based on the orthonormal Bernstein polynomials has been introduced to provide an approximate solution of nonlinear fractional partial Volterra integro‐differential equations. For this aim, the fractional‐order orthogonal Bernstein polynomials (FOBPs) are constructed, and its operational matrices of integration, fractional‐order integration, and derivative in the Caputo sense and product operational matrix are derived. These operational matrices are utilized to reduce the under study problem to a nonlinear system of algebraic equations. Using the approximation of FOBPs, the convergence analysis and error estimate associated to the proposed problem have been investigated. Finally, several examples are included to clarify the validity, efficiency, and applicability of the proposed technique via FOBPs approximation.