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.     

Legendre polynomial solutions of high-order linear Fredholm integro-differential equations

In this study, a Legendre collocation matrix method is presented to solve high-order Linear Fredholm integro-differential equations under the mixed conditions in terms of Legendre polynomials. The proposed method converts the equation and conditions to matrix equations, by means of collocation points on the interval [−1, 1], which corresponding to systems of linear algebraic equations with Legendre coefficients. Thus, by solving the matrix equation, Legendre coefficients and polynomial approach are obtained. Also examples that illustrate the pertinent features of the method are presented and by using the error analysis, the results are discussed.     

Legendre wavelets method for constrained optimal control problems

A numerical method for solving non‐linear optimal control problems with inequality constraints is presented in this paper. The method is based upon Legendre wavelet approximations. The properties of Legendre wavelets are first presented. The operational matrix of integration and the Gauss method are then utilized to reduce the optimal control problem to the solution of algebraic equations. The inequality constraints are converted to a system of algebraic equalities; these equalities are then collocated at the Gauss nodes. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2002 John Wiley & Sons, Ltd.     

The L<Superscript>2</Superscript>-convergence of the Legendre spectral Tau matrix formulation for nonlinear fractional integro differential equations

The operational Tau method, a well known method for solving functional equations is employed to approximate the solution of nonlinear fractional integro-differential equations. The fractional derivatives are described in the Caputo sense. The unique solvability of the linear Tau algebraic system is discussed. In addition, we provide a rigorous convergence analysis for the Legendre Tau method which indicate that the proposed method converges exponentially provided that the data in the given FIDE are smooth. To do so, Sobolev inequality with some properties of Banach algebras are considered. Some numerical results are given to clarify the efficiency of the 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.     

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 technique for solving fractional variational problems by Müntz–Legendre polynomials

This paper presents an efficient approximate method for solving a class of fractional variational problems (FVPs). The fractional derivatives are described in the Caputo sense. In the our method, we use the Müntz–Legendre orthonormal basis and the properties of Rayleigh–Ritz method for fractional calculus to reduce FVPs to solve a system of algebraic equations which solved using Newton’s iterative method. Also, an estimation of the error is given in the sense of Sobolev norms. The illustrative examples are provided to demonstrate the applicability and simplicity of the new technique.     

Fredholm Integro——differential型方程的Legendre小波方法

研究Legendre小波方法求解具有一阶导和二阶导类型的线性Fredholm integro-differential型方程。应用Legendre小波逼近法把这两类方程分别化为代数方程求解．实例说明。Legendre小波在解决这两类方程时的可行性和有效性．     

The modulus-based matrix splitting iteration methods for second-order cone linear complementarity problems

This paper presents a new computational technique for solving fractional pantograph differential equations. The fractional derivative is described in the Caputo sense. The main idea is to use Müntz-Legendre wavelet and its operational matrix of fractional-order integration. First, the Müntz-Legendre wavelet is presented. Then a family of piecewise functions is proposed, based on which the fractional order integration of the Müntz-Legendre wavelets are easy to calculate. The proposed approach is used this operational matrix with the collocation points to reduce the under study problem to a system of algebraic equations. An estimation of the error is given in the sense of Sobolev norms. The efficiency and accuracy of the proposed method are illustrated by several numerical examples.     

Legendre wavelets approach for numerical solutions of distributed order fractional differential equations

In this study, a new numerical method for the solution of the linear and nonlinear distributed fractional differential equations is introduced. The fractional derivative is described in the Caputo sense. The suggested framework is based upon Legendre wavelets approximations. For the first time, an exact formula for the Riemann–Liouville fractional integral operator for the Legendre wavelets is derived. We then use this formula and the properties of Legendre wavelets to reduce the given problem into a system of algebraic equations. Several illustrative examples are included to observe the validity, effectiveness and accuracy of the present numerical method.     

A computational method for solving optimal control of a system of parallel beams using Legendre wavelets

The optimal control of transverse vibration of two Euler–Bernoulli beams coupled in parallel by discrete springs is considered. An index of performance is formulated which consists of a modified energy functional of two coupled structures at a specified time and penalty functions involving the point control forces. The minimization of the performance index over these forces is subject to the equation of motion governing the structural vibrations, the imposed initial condition as well as the boundary conditions. By use of the modal space technique, the optimal control of distributed parameter systems is simplified into the optimal control of a linear time-invariant lumped-parameter systems. A computationally attractive method based on Legendre wavelets in time domain for solving the optimal control of the lumped parameter systems for any finite interval is proposed. Legendre wavelet integral operational matrix and the properties of a Kronecker product are used to find the approximated optimal trajectory and optimal law of the linear systems with respect to a quadratic cost function by only solving a linear system of algebraic equations. This method provides a straightforward and convenient approach for digital computation. A numerical example is provided to demonstrate the applicability and effectiveness of the proposed method.     

Shifted Legendre direct method for variational problems

The shifted Legendre polynomial series is employed to solve variational problems. The solution is carried out by using an operational matrix for integrating the shifted Legendre polynomial vector. Variational problems are reduced to solving algebraic equations. Two illustrative examples are given, and the computational results obtained by Legendre series direct method are compared with the exact solutions.     

A quadrature tau method for fractional differential equations with variable coefficients

In this article, we develop a direct solution technique for solving multi-order fractional differential equations (FDEs) with variable coefficients using a quadrature shifted Legendre tau (Q-SLT) method. The spatial approximation is based on shifted Legendre polynomials. A new formula expressing explicitly any fractional-order derivatives of shifted Legendre polynomials of any degree in terms of shifted Legendre polynomials themselves is proved. Extension of the tau method for FDEs with variable coefficients is treated using the shifted Legendre–Gauss–Lobatto quadrature. Numerical results are given to confirm the reliability of the proposed method for some FDEs with variable coefficients.     

A New Approach for Solving a Class of Delay Fractional Partial Differential Equations

In this article, a new numerical approach has been proposed for solving a class of delay time-fractional partial differential equations. The approximate solutions of these equations are considered as linear combinations of Müntz–Legendre polynomials with unknown coefficients. Operational matrix of fractional differentiation is provided to accelerate computations of the proposed method. Using Padé approximation and two-sided Laplace transformations, the mentioned delay fractional partial differential equations will be transformed to a sequence of fractional partial differential equations without delay. The localization process is based on the space-time collocation in some appropriate points to reduce the fractional partial differential equations into the associated system of algebraic equations which can be solved by some robust iterative solvers. Some numerical examples are also given to confirm the accuracy of the presented numerical scheme. Our results approved decisive preference of the Müntz–Legendre polynomials with respect to the Legendre polynomials.     

Numerical solution of singular IVPs of Lane–Emden type using a modified Legendre-spectral method

In this paper, a numerical method for singular initial value problems of the Lane–Emden type in the second-order ordinary differential equations is proposed. The method changes solving the equation to solving a Volterra integral equation. We have applied the improved Legendre-spectral method to solve Lane–Emden type equations. The Legendre Gauss points are used as collocation nodes and Lagrange interpolation is employed in Volterra term. The results reveal that the method is effective, simple and accurate.     

An optimal adaptive wavelet method without coarsening of the iterands

In this paper, an adaptive wavelet method for solving linear operator equations is constructed that is a modification of the method from [Math. Comp, 70 (2001), pp. 27-75] by Cohen, Dahmen and DeVore, in the sense that there is no recurrent coarsening of the iterands. Despite this, it will be shown that the method has optimal computational complexity. Numerical results for a simple model problem indicate that the new method is more efficient than an existing alternative adaptive wavelet method.

     

A space-time Legendre spectral tau method for the two-sided space-time Caputo fractional diffusion-wave equation

The space-time fractional diffusion-wave equation (FDWE) is a generalization of classical diffusion and wave equations which is used in modeling practical phenomena of diffusion and wave in fluid flow, oil strata and others. This paper reports an accurate spectral tau method for solving the two-sided space and time Caputo FDWE with various types of nonhomogeneous boundary conditions. The proposed method is based on shifted Legendre tau (SLT) procedure in conjunction with the shifted Legendre operational matrices of Riemann-Liouville fractional integral, left-sided and right-sided fractional derivatives. We focus primarily on implementing this algorithm in both temporal and spatial discretizations. In addition, convergence analysis is provided theoretically for the Dirichlet boundary conditions, along with graphical analysis for several special cases using other conditions. These suggest that the Legendre Tau method converges exponentially provided that the data in the given FDWE are smooth. Finally, several numerical examples are given to demonstrate the high accuracy of the proposed method.     

A Legendre–Galerkin Chebyshev collocation method for the Burgers equation with a random perturbation on boundary condition

In the paper, we apply the generalized polynomial chaos expansion and spectral methods to the Burgers equation with a random perturbation on its left boundary condition. Firstly, the stochastic Galerkin method combined with the Legendre–Galerkin Chebyshev collocation scheme is adopted, which means that the original equation is transformed to the deterministic nonlinear equations by the stochastic Galerkin method and the Legendre–Galerkin Chebyshev collocation scheme is used to deal with the resulting nonlinear equations. Secondly, the stochastic Legendre–Galerkin Chebyshev collocation scheme is developed for solving the stochastic Burgers equation; that is, the stochastic Legendre–Galerkin method is used to discrete the random variable meanwhile the nonlinear term is interpolated through the Chebyshev–Gauss points. Then a set of deterministic linear equations can be obtained, which is in contrast to the other existing methods for the stochastic Burgers equation. The mean square convergence of the former method is analyzed. Numerical experiments are performed to show the effectiveness of our two methods. Both methods provide alternative approaches to deal with the stochastic differential equations with nonlinear terms.     

The comparison of two reliable methods for accurate solution of time‐fractional Kaup–Kupershmidt equation arising in capillary gravity waves

In this paper, we compared two different methods, one numerical technique, viz Legendre multiwavelet method, and the other analytical technique, viz optimal homotopy asymptotic method (OHAM), for solving fractional‐order Kaup–Kupershmidt (KK) equation. Two‐dimensional Legendre multiwavelet expansion together with operational matrices of fractional integration and derivative of wavelet functions is used to compute the numerical solution of nonlinear time‐fractional KK equation. The approximate solutions of time fractional Kaup–Kupershmidt equation thus obtained by Legendre multiwavelet method are compared with the exact solutions as well as with OHAM. The present numerical scheme is quite simple, effective, and expedient for obtaining numerical solution of fractional KK equation in comparison to analytical approach of OHAM. Copyright © 2016 John Wiley & Sons, Ltd.     

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.