The operational matrix of fractional integration for shifted Chebyshev polynomials

A new shifted Chebyshev operational matrix (SCOM) of fractional integration of arbitrary order is introduced and applied together with spectral tau method for solving linear fractional differential equations (FDEs). The fractional integration is described in the Riemann–Liouville sense. The numerical approach is based on the shifted Chebyshev tau method. The main characteristic behind the approach using this technique is that only a small number of shifted Chebyshev polynomials is needed to obtain a satisfactory result. Illustrative examples reveal that the present method is very effective and convenient for linear multi-term FDEs.     

The alternative Legendre Tau method for solving nonlinear multi-order fractional differential equations

In this paper, the alternative Legendre polynomials (ALPs) are used to approximate the solution of a class of nonlinear multi-order fractional differential equations (FDEs). First, the operational matrix of fractional integration of an arbitrary order and the product operational matrix are derived for ALPs. These matrices together with the spectral Tau method are then utilized to reduce the solution of the mentioned equations into the one of solving a system of nonlinear algebraic equations with unknown ALP coefficients of the exact solution. The fractional derivatives are considered in the Caputo sense and the fractional integration is described in the Riemann-Liouville sense. Numerical examples illustrate that the present method is very effective for linear and nonlinear multi-order FDEs and high accuracy solutions can be obtained only using a small number of ALPs.     

A numerical-analytical solution of multi-term fractional-order differential equations

In this paper, a solution to initial value problems for fractional-order linear commensurate multi-term differential equations with Caputo derivatives is presented. The solution is obtained in the form of a finite sum of the Mittag-Leffler–type functions and the meta-trigonometric cosine function by using a numerical-analytical method. The results of presented numerical experiments show that for high accuracy calculations of these functions, the multi-precision arithmetic must be applied. The approach for solving of the initial value problems for generalized Basset equation, generalized Bagley-Torvik equation, and multi-term fractional equation is demonstrated.     

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.     

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.     

一类变系数分数阶微分方程的数值解法

本文研究了一类变系数分数阶微分方程的数值解法问题. 利用Cheyshev小波推导出的分数阶微分方程的算子矩阵把分数阶微分方程转换为代数方程组. 同时给出了Cheyshev小波基的收敛性和误差估计表达式, 并给出数值算例说明所提方法的精确性和有效性     

On the applicability of Genocchi wavelet method for different kinds of fractional‐order differential equations with delay

A novel collocation method based on Genocchi wavelet is presented for the numerical solution of fractional differential equations and time‐fractional partial differential equations with delay. In this work, to achieve the approximate solution with height accuracy, we employed the operational matrix of integer derivative and the pseudo‐operational matrix of fractional derivative in Caputo sense. Also, based on Genocchi function properties, we presented delay and pantograph operational matrices of Genocchi wavelet functions (GWFs). Due to operational and pseudo‐operational matrices, the equations under this study can be turned into nonlinear algebraic equations with the unknown GWF coefficients. For illustrating the upper bound of error for the proposed method, we estimate the error in the sense of Sobolev space. In addition, to demonstrate the efficacy of the pseudo‐operational matrix of fractional derivative, we investigate the upper bound of error for the mentioned matrix. Finally, the algorithm based on the proposed approach is implemented for some numerical experiments to confirm accuracy and applicability.     

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.     

A numerical scheme based on Bernoulli wavelets and collocation method for solving fractional partial differential equations with Dirichlet boundary conditions

In this paper, an efficient and accurate numerical method is presented for solving two types of fractional partial differential equations. The fractional derivative is described in the Caputo sense. Our approach is based on Bernoulli wavelets collocation techniques together with the fractional integral operator, described in the Riemann‐Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved and the solution of fractional partial differential equations is achieved. Some results concerning the error analysis are obtained. The validity and applicability of the method are demonstrated by solving four numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions much easier.     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method together with the fractional integral operator, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved, and the solution of fractional nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

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 of distributed order fractional differential equations in the time domain using the Müntz–Legendre wavelets approach

In this paper, a numerical method is presented to obtain and analyze the behavior of numerical solutions of distributed order fractional differential equations of the general form in the time domain with the Caputo fractional derivative. The suggested method is based on the Müntz–Legendre wavelet approximation. We derive a new operational vector for the Riemann–Liouville fractional integral of the Müntz–Legendre wavelets by using the Laplace transform method. Applying this operational vector and collocation method in our approach, the problem can be reduced to a system of linear and nonlinear algebraic equations. The arising system can be solved by the Newton method. Discussion on the error bound and convergence analysis for the proposed method is presented. Finally, seven test problems are considered to compare our results with other well‐known methods used for solving these problems. The results in the tabulated tables highlighted that the proposed method is an efficient mathematical tool for analyzing distributed order fractional differential equations of the general form.     

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.     

A new operational vector approach for time-fractional subdiffusion equations of distributed order based on hybrid functions

This research study deals with the numerical solutions of linear and nonlinear time-fractional subdiffusion equations of distributed order. The main aim of our approach is based on the hybrid of block-pulse functions and shifted Legendre polynomials. We produce a novel and exact operational vector for the fractional Riemann–Liouville integral and use it via the Gauss–Legendre quadrature formula and collocation method. Consequently, we reduce the proposed equations to systems of equations. The convergence and error bounds for the new method are investigated. Six problems are tested to confirm the accuracy of the proposed approach. Comparisons between the obtained numerical results and other existing methods are provided. Numerical experiments illustrate the reliability, applicability, and efficiency of the proposed method.     

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.     

Investigating a Class of Nonlinear Fractional Differential Equations and Its Hyers-Ulam Stability by Means of Topological Degree Theory

The aim of this article is to seek some adequate conditions via a prior estimate method (topological degree method) to derive the existence of solution to a nonlinear boundary value problem of fractional differential equations (FDEs). With the help of topological degree method which has been applied in many articles, we establish the required results for existence and uniqueness of solution to a class of FDEs. Moreover, we also formulate sufficient conditions for Hyers-Ulam stability to the solution of the considered problem. Finally, an appropriate example is provided to justify the relevant results.     

Simulation of nonlinear fractional dynamics arising in the modeling of cognitive decision making using a new fractional neural network

By the rapid growth of available data, providing data-driven solutions for nonlinear (fractional) dynamical systems becomes more important than before. In this paper, a new fractional neural network model that uses fractional order of Jacobi functions as its activation functions for one of the hidden layers is proposed to approximate the solution of fractional differential equations and fractional partial differential equations arising from mathematical modeling of cognitive-decision-making processes and several other scientific subjects. This neural network uses roots of Jacobi polynomials as the training dataset, and the Levenberg-Marquardt algorithm is chosen as the optimizer. The linear and nonlinear fractional dynamics are considered as test examples showing the effectiveness and applicability of the proposed neural network. The numerical results are compared with the obtained results of some other networks and numerical approaches such as meshless methods. Numerical experiments are presented confirming that the proposed model is accurate, fast, and feasible.