Spline collocation methods for linear multi-term fractional differential equations

Some regularity properties of the solution of linear multi-term fractional differential equations are derived. Based on these properties, the numerical solution of such equations by piecewise polynomial collocation methods is discussed. The results obtained in this paper extend the results of Pedas and Tamme (2011) [15] where we have assumed that in the fractional differential equation the order of the highest derivative of the unknown function is an integer. In the present paper, we study the attainable order of convergence of spline collocation methods for solving general linear fractional differential equations using Caputo form of the fractional derivatives and show how the convergence rate depends on the choice of the grid and collocation points. Theoretical results are verified by some numerical examples.     

Piecewise polynomial collocation for linear boundary value problems of fractional differential equations

We consider a class of boundary value problems for linear multi-term fractional differential equations which involve Caputo-type fractional derivatives. Using an integral equation reformulation of the boundary value problem, some regularity properties of the exact solution are derived. Based on these properties, the numerical solution of boundary value problems by piecewise polynomial collocation methods is discussed. In particular, we study the attainable order of convergence of proposed algorithms and show how the convergence rate depends on the choice of the grid and collocation points. Theoretical results are verified by two numerical examples.     

On the convergence of spline collocation methods for solving fractional differential equations

In the first part of this paper we study the regularity properties of solutions of initial value problems of linear multi-term fractional differential equations. We then use these results in the convergence analysis of a polynomial spline collocation method for solving such problems numerically. Using an integral equation reformulation and special non-uniform grids, global convergence estimates are derived. From these estimates it follows that the method has a rapid convergence if we use suitable nonuniform grids and the nodes of the composite Gaussian quadrature formulas as collocation points. Theoretical results are verified by some numerical examples.     

Numerical solution of multi-order fractional differential equations with multiple delays via spectral collocation methods

This paper discusses a general framework for the numerical solution of multi-order fractional delay differential equations (FDDEs) in noncanonical forms with irrational/rational multiple delays by the use of a spectral collocation method. In contrast to the current numerical methods for solving fractional differential equations, the proposed framework can solve multi-order FDDEs in a noncanonical form with incommensurate orders. The framework can also solve multi-order FDDEs with irrational multiple delays. Next, the framework is enhanced by the fractional Chebyshev collocation method in which a Chebyshev operation matrix is constructed for the fractional differentiation. Spectral convergence and small computational time are two other advantages of the proposed framework enhanced by the fractional Chebyshev collocation method. In addition, the convergence, error estimates, and numerical stability of the proposed framework for solving FDDEs are studied. The advantages and computational implications of the proposed framework are discussed and verified in several numerical examples.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

A new Jacobi operational matrix: An application for solving fractional differential equations

In this paper, we derived the shifted Jacobi operational matrix (JOM) of fractional derivatives which is applied together with spectral tau method for numerical solution of general linear multi-term fractional differential equations (FDEs). A new approach implementing shifted Jacobi operational matrix in combination with the shifted Jacobi collocation technique is introduced for the numerical solution of nonlinear multi-term FDEs. The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations which greatly simplifying the problem. The proposed methods are applied for solving linear and nonlinear multi-term FDEs subject to initial or boundary conditions, and the exact solutions are obtained for some tested problems. Special attention is given to the comparison of the numerical results obtained by the new algorithm with those found by other known methods.     

Collocation methods for fractional integro-differential equations with weakly singular kernels

In this paper, the piecewise polynomial collocation methods are used for solving the fractional integro-differential equations with weakly singular kernels. We present that a suitable transformation can convert fractional integro-differential equations to one type of second kind Volterra integral equations (VIEs) with weakly singular kernels. Then we solve the VIEs by standard piecewise polynomial collocation methods. It is shown that such kinds of methods are able to yield optimal convergence rate. Finally, some numerical experiments are given to show that the numerical results are consistent with the theoretical results.     

一类分数阶多项延迟微分方程的Jacobi谱配置方法

本文利用Jacobi谱配置方法数值求解了一类分数阶多项延迟微分方程,并证明了该方法是收敛的,通过若干数值算例验证了相应的理论结果,结果表明Jacobi谱配置方法求解这类方程是非常高效的,同时也为这类分数阶延迟微分方程的数值求解提供了新的选择,对分数阶泛函方程的数值方法的研究有一定的指导意义.     

A fast numerical algorithm based on the Taylor wavelets for solving the fractional integro‐differential equations with weakly singular kernels

In this paper, a fast numerical algorithm based on the Taylor wavelets is proposed for finding the numerical solutions of the fractional integro‐differential equations with weakly singular kernels. The properties of Taylor wavelets are given, and the operational matrix of fractional integration is constructed. These wavelets are utilized to reduce the solution of the given fractional integro‐differential equation to the solution of a linear system of algebraic equations. Also, convergence of the proposed method is studied. Illustrative examples are included to demonstrate the validity and applicability of the technique.     

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.     

Stability of piecewise polynomial collocation for computing periodic solutions of delay differential equations

Summary. We prove numerical stability of a class of piecewise polynomial collocation methods on nonuniform meshes for computing asymptotically stable and unstable periodic solutions of the linear delay differential equation by a (periodic) boundary value approach. This equation arises, e.g., in the study of the numerical stability of collocation methods for computing periodic solutions of nonlinear delay equations. We obtain convergence results for the standard collocation algorithm and for two variants. In particular, estimates of the difference between the collocation solution and the true solution are derived. For the standard collocation scheme the convergence results are “unconditional”, that is, they do not require mesh-ratio restrictions. Numerical results that support the theoretical findings are also given. Received June 9, 2000 / Revised version received December 14, 2000 / Published online October 17, 2001     

On spline collocation methods for boundary integral equations in the plane

Nowadays boundary elemen; methods belong to the most popular numerical methods for solving elliptic boundary value problems. They consist in the reduction of the problem to equivalent integral equations (or certain generalizations) on the boundary Γ of the given domain and the approximate solution of these boundary equations. For the numerical treatment the boundary surface is decomposed into a finite number of segments and the unknown functions are approximated by corresponding finite elements and usually determined by collocation and Galerkin procedures. One finds the least difficulties in the theoretical foundation of the convergence of Galerkin methods for certain classes of equations, whereas the convergence of collocation methods, which are mostly used in numerical computations, has yet been proved only for special equations and methods. In the present paper we analyse spline collocation methods on uniform meshes with variable collocation points for one-dimensional pseudodifferential equations on a closed curve with convolutional principal parts, which encompass many classes of boundary integral equations in the plane. We give necessary and sufficient conditions for convergence and prove asymptotic error estimates. In particular we generalize some results on nodal and midpoint collocation obtained in [2], [7] and [8]. The paper is organized as follows. In Section 1 we formulate the problems and the results, Section 2 deals with spline interpolation in periodic Sobolev spaces, and in Section 3 we prove the convergence theorems for the considered collocation methods.     

Gegenbauer spectral method for time‐fractional convection–diffusion equations with variable coefficients

In this paper, we study the numerical solution to time‐fractional partial differential equations with variable coefficients that involve temporal Caputo derivative. A spectral method based on Gegenbauer polynomials is taken for approximating the solution of the given time‐fractional partial differential equation in time and a collocation method in space. The suggested method reduces this type of equation to the solution of a linear algebraic system. Finally, some numerical examples are presented to illustrate the efficiency and accuracy of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

Analysis of a moving collocation method for one-dimensional partial differential equations

A moving collocation method has been shown to be very efficient for the adaptive solution of second- and fourth-order time-dependent partial differential equations and forms the basis for the two robust codes MOVCOL and MOVCOL4.In this paper,the relations between the method and the traditional collocation and finite volume methods are investigated.It is shown that the moving collocation method inherits desirable properties of both methods: the ease of implementation and high-order convergence of the traditional collocation method and the mass conservation of the finite volume method.Convergence of the method in the maximum norm is proven for general linear two-point boundary value problems.Numerical results are given to demonstrate the convergence order of the method.     

Numerical comparison of methods for solving linear differential equations of fractional order

In this article, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving linear differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper will present a numerical comparison between the two methods and a conventional method such as the fractional difference method for solving linear differential equations of fractional order. The numerical results demonstrates that the new methods are quite accurate and readily implemented.     

Efficient Mittag-Leffler Collocation Method for Solving Linear and Nonlinear Fractional Differential Equations

In this paper, a new approximation method for fractional differential equations based on Mittag-Leffler function is developed. Finite Mittag-Leffler function and its fractional-order derivatives are investigated. An efficient technique for solving linear and nonlinear fractional order differential equations is developed. The proposed method combines Mittag-Leffler collocation method and optimization technique. Error estimation of the approximation is stated and proved. We present numerical results and comparisons of previous treatments to demonstrate the efficiency and applicability of the proposed method. Making use of small number of unknowns, the resulting solution converges to the exact one in the linear case and it has a very small error in the nonlinear case.     

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.     

变系数线性微分方程初值问题数值解的小波方法

通过利用小波尺度函数的正交性并结合配点法 ,本文给出了一种求解变系数线性微分方程初值问题数值解的小波算法 .在一定的假设条件下 ,对算法的收敛性进行了理论分析 .最后 ,我们还给出了一个具体的数值计算例子 .     

A highly accurate Jacobi collocation algorithm for systems of high‐order linear differential–difference equations with mixed initial conditions

In this paper, a shifted Jacobi–Gauss collocation spectral algorithm is developed for solving numerically systems of high‐order linear retarded and advanced differential–difference equations with variable coefficients subject to mixed initial conditions. The spatial collocation approximation is based upon the use of shifted Jacobi–Gauss interpolation nodes as collocation nodes. The system of differential–difference equations is reduced to a system of algebraic equations in the unknown expansion coefficients of the sought‐for spectral approximations. The convergence is discussed graphically. The proposed method has an exponential convergence rate. The validity and effectiveness of the method are demonstrated by solving several 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 more easier. Copyright © 2015 John Wiley & Sons, Ltd.     

Legendre spectral method for the fractional Bratu problem

In this paper, the Legendre spectral collocation method (LSCM) is applied for the solution of the fractional Bratu's equation. It shows the high accuracy and low computational cost of the LSCM compared with some other numerical methods. The fractional Bratu differential equation is transformed into a nonlinear system of algebraic equations for the unknown Legendre coefficients and solved with some spectral collocation methods. Some illustrative examples are also given to show the validity and applicability of this method, and the obtained results are compared with the existing studies to highlight its high efficiency and neglectable error.