Legendre wavelets method for solving fractional integro-differential equations

In this paper, based on the constructed Legendre wavelets operational matrix of integration of fractional order, a numerical method for solving linear and nonlinear fractional integro-differential equations is proposed. By using the operational matrix, the linear and nonlinear fractional integro-differential equations are reduced to a system of algebraic equations which are solved through known numerical algorithms. The upper bound of the error of the Legendre wavelets expansion is investigated in Theorem 5.1. Finally, four numerical examples are shown to illustrate the efficiency and accuracy of the approach.     

Using reproducing kernel for solving a class of singular weakly nonlinear boundary value problems

Cui et al. [M. Cui and F. Geng, Solving singular two point boundary value problems in reproducing kernel space, J. Comput. Appl. Math. 205 (2007), pp. 6–15; H. Yao and M. Cui, A new algorithm for a class of singular boundary value problems, Appl. Math. Comput. 186 (2007), pp. 1183–1191] presents an algorithm to solve a class of singular linear boundary value problems in the reproducing kernel space. In this paper, we will present three new algorithms to solve a class of singular weakly nonlinear boundary value problems in reproducing kernel space. The algorithms are efficiently applied to solving some model problems. It is demonstrated by the numerical examples that those algorithms are highly accurate.     

Quasi-wavelet based numerical method for fourth-order partial integro-differential equations with a weakly singular kernel

In this paper, we study the numerical solution of initial boundary-value problem for the fourth-order partial integro-differential equations with a weakly singular kernel. We use the forward Euler scheme for time discretization and the quasi-wavelet based numerical method for space discretization. Detailed discrete formulations are given to the treatment of three different boundary conditions, including clamped-type condition, simply supported-type condition and a transversely supported-type condition. Some numerical experiments are included to demonstrate the validity and applicability of the discrete technique. The comparisons of present results with analytical solutions show that the quasi-wavelet based numerical method has a distinctive local property. Especially, the method is easy to implement and produce very accurate results.     

An efficient fractional-order wavelet method for fractional Volterra integro-differential equations

In this paper, an efficient and robust numerical technique is suggested to solve fractional Volterra integro-differential equations (FVIDEs). The proposed method is mainly based on the generalized fractional-order Legendre wavelets (GFLWs), their operational matrices and the Collocation method. The main advantage of the proposed method is that, by using the GFLWs basis, it can provide more efficient and accurate solution for FVIDEs in compare to integer-order wavelet basis. A comparison between the achieved results confirms accuracy and superiority of the proposed GFLWs method for solving FVIDEs. Error analysis and convergence of the GFLWs basis is provided.     

Numerical approach for solving fractional Fredholm integro-differential equation

In this article, we present a new method which is based on the Taylor Matrix Method to give approximate solution of the linear fractional Fredholm integro-differential equations. This method is based on first taking the truncated Taylor expansions of the functions in the linear fractional differential part and Fredholm integral part then, substituting their matrix forms into the equation. We solve this matrix equation with the assistance of Maple 13. In addition, illustrative examples are presented to demonstrate the effectiveness of the proposed method.     

A wavelet approach for the multi-term time fractional diffusion-wave equation

In this paper, a Galerkin method based on the second kind Chebyshev wavelets (SKCWs) is established for solving the multi-term time fractional diffusion-wave equation. To do this, a new operational matrix of fractional integration for the SKCWs must be derived and in order to improve the computational efficiency, the hat functions are proposed to create a general procedure for constructing this matrix. Implementation of these wavelet basis functions and their operational matrix of fractional integration simplifies the problem under consideration to a system of linear algebraic equations, which greatly decreases the computational cost for finding an approximate solution. The main privilege of the proposed method is adjusting the initial and boundary conditions in the final system automatically. Theoretical error and convergence analysis of the SKCWs expansion approve the reliability of the approach. Also, numerical investigation reveals the applicability and accuracy of the presented method.     

应用Laguerre小波算子矩阵求二重数值积分

基于Laguerre小波函数及其对应的积分算子矩阵给出了一个求二重积分的数值方法。该方法通过对被积函数进行恰当的离散化,将二重数值积分问题转化为矩阵运算,从而易于求解、方便计算。该方法不仅适用于积分区域是矩形区域,也适用于二重变限积分的情况。数值算例验证了该方法的可行性及有效性。     

Extrapolation for solving a system of weakly singular nonlinear Volterra integral equations of the second kind

This article discusses an extrapolation method for solving a system of weakly singular nonlinear Volterra integral equations of the second kind. Based on a generalization of the discrete Gronwall inequality and Navot's quadrature rule, the modified trapeziform quadrature algorithm is presented. The iterative algorithm for solving the discrete system possesses a high accuracy order O(h ^2+α). After the asymptotic expansion of errors is proved, we can obtain an approximation with a higher accuracy order using extrapolation. An a posteriori error estimation is provided. Some numerical results are presented to illustrate the efficiency of our methods.     

A direct method to solve integral and integro-differential equations of convolution type by using improved operational matrix

In this article, we use improved operational matrix of block pulse functions on interval [0,?1) to solve Volterra integral and integro-differential equations of convolution type without solving any system and projection method. We first obtain Laplace transform of the problem and then we find numerical inversion of Laplace transform by improved operational matrix of integration. Numerical examples show that the approximate solutions have a good degree of accuracy.     

An improved reproducing kernel method for Fredholm integro-differential type two-point boundary value problems

In this work, an improved reproducing kernel method to find the numerical solution of Fredholm integro-differential equation type boundary value problems has been developed. Based on the good properties of reproducing kernel function and the conjugate operator, the solution representation is obtained. Meanwhile, we prove that the approximation converges to the exact solution uniformly. After that the convergence estimates are derived.     

Wavelets Galerkin method for solving stochastic heat equation

In this paper, a new computational method based on the second kind Chebyshev wavelets (SKCWs) together with the Galerkin method is proposed for solving a class of stochastic heat equation. For this purpose, a new stochastic operational matrix for the SKCWs is derived. A collocation method based on block pulse functions is employed to derive a general procedure for forming this matrix. The SKCWs and their operational matrices of integration and stochastic Itô-integration are used to transform the under consideration problem into the corresponding linear system of algebraic equations which can be simply solved to achieve the solution of the problem. The proposed method is very convenient for solving such problems, since the initial and boundary conditions are taken into account automatically. Moreover, the efficiency of the proposed method is shown for some concrete examples. The results reveal that the proposed method is very accurate and efficient.     

Legendre wavelet method for solving variable-order nonlinear fractional optimal control problems with variable-order fractional Bolza cost

This paper deals with a numerical method for solving variable-order fractional optimal control problem with a fractional Bolza cost composed as the aggregate of a standard Mayer cost and a fractional Lagrange cost given by a variable-order Riemann–Liouville fractional integral. Using the integration by part formula and the calculus of variations, the necessary optimality conditions are derived in terms of two-point variable-order boundary value problem. Operational matrices of variable-order right and left Riemann–Liouville integration are derived, and by using them, the two-point boundary value problem is reduced into the system of algebraic equations. Additionally, the convergence analysis of the proposed method has been considered. Moreover, illustrative examples are given to demonstrate the applicability of the proposed method.     

Numerical solution of integro-differential equations of high-order Fredholm by the simplified reproducing kernel method

The key of the reproducing kernel method (RKM) to solve the initial boundary value problem is to construct the reproducing kernel meeting the homogenous initial boundary conditions of the considered problems. The usual method is that the initial boundary conditions must be homogeneous and put them into space. Another common method is to put homogeneous or non-homogeneous conditions directly into the operator. In addition, we give a new numerical method of RKM for dealing with initial boundary value problems, homogeneous conditions are put into space, and for nonhomogeneous conditions, we put them into operators. The focus of this paper is to further verify the reliability and accuracy of the latter two methods. Through solving three numerical examples of integral–differential equations and comparing with other methods, we find that the two methods are useful.     

Chebyshev collocation method for solving singular integral equation with cosecant kernel

A collocation method based on Chebyshev polynomials is proposed for solving cosecant-type singular integral equations (SIE). For solving SIE, difficulties lie in its singular term. In order to remove singular term, we introduce Gauss–Legendre integration and integral properties of the cosecant kernel. An advantage of this method is to approximate the best uniform approximation by the best square approximation to obtain the unknown coefficients in the method. On the other hand, the convergence is fast and the accuracy is high, which is verified by the final numerical experiments compared with the existing references.     

Numerical solution for a class of space-time fractional equation by the piecewise reproducing kernel method

ABSTRACT

Due to the non-locality of fractional derivative, the analytical solution and good approximate solution of fractional partial differential equations are usually difficult to get. Reproducing kernel space is a perfect space in studying this type of equations, however the numerical results of equations by using the traditional reproducing kernel method (RKM) isn't very good. Based on this problem, we present the piecewise technique in the reproducing kernel space to solve this type of equations. The focus of this paper is to verify the stability and high accuracy of the present method by comparing the absolute error with traditional RKM and study the effect on absolute error for different values of α. Furthermore, we can study the distribution of entire space at a particular time period. Three numerical experiments are provided to verify the efficiency and stability of the proposed method. Meanwhile, it is tested by experiments that the change of the value of α has little effect on its accuracy.     

Numerical solution of an integral equations system of the first kind by using an operational matrix with block pulse functions

This article proposes a simple efficient method for solving a Volterra integral equations system of the first kind. By using block pulse functions and their operational matrix of integration, a first kind integral equations system can be reduced to a linear system of algebraic equations. The coefficient matrix of this system is a block matrix with lower triangular blocks. Numerical examples show that the approximate solutions have a good degree of accuracy.     

An efficient method based on the second kind Chebyshev wavelets for solving variable-order fractional convection diffusion equations

In this paper, a class of variable-order fractional convection diffusion equations have been solved with assistance of the second kind Chebyshev wavelets operational matrix. The operational matrix of variable-order fractional derivative is derived for the second kind Chebyshev wavelets. By implementing the second kind Chebyshev wavelets functions and also the associated operational matrix, the considered equations will be reduced to the corresponding Sylvester equation, which can be solved by some appropriate iterative solvers. Also, the convergence analysis of the proposed numerical method to the exact solutions and error estimation are given. A variety of numerical examples are considered to show the efficiency and accuracy of the presented technique.