In this paper, we propose a high-precision numerical algorithm for a fractional integro-differential equation based on the shifted Legendre polynomials and the idea of Gauss-Legendre quadrature rule and spectral collocation method. The error analysis of this method is also given in detail. Some numerical examples are give to illustrate the exponential convergence of our method.
An efficient spectral method is proposed for a new Steklov eigenvalue problem in inverse scattering. Firstly, we establish the weak form and the associated discrete scheme by introducing an appropriate Sobolev space and a corresponding approximation space. Then, according to the Fredholm Alternative, the corresponding operator forms of weak formulation and discrete formulation are derived. After that, the error estimates of approximated eigenvalues and eigenfunctions are proved by using the spectral approximation results of completely continuous operators and the approximation properties of orthogonal projection operators. We also construct an appropriate set of basis functions in the approximation space and derive the matrix form of the discrete scheme based on the tensor product. In addition, we extend the algorithm to the circular domain. Finally, we present plenty of numerical experiments and compare them with some existing numerical methods, which validate that our algorithm is effective and high accuracy.
In this paper, the time fractional diffusion equations optimal control problem is solved by order with uniform accuracy scheme in time and finite element method (FEM) in space. For the state and adjoint state equation, the piecewise linear polynomials are used to make the space variables discrete, and obtain the semidiscrete scheme of the state and adjoint state. The priori error estimates for the semidiscrete scheme for state and adjoint state equation are established. Furthermore, the order uniform accuracy scheme is used to make the time variable discrete in the semidiscrete scheme and construct the full discrete scheme for the control problems based on the first optimal condition and ‘first optimize, then discretize’ approach. The fully discrete scheme’s stability and truncation error are analyzed. Finally, two numerical examples are denoted to show that the theoretical analysis are correct.
In this paper, we present a high-order approximate solution with uniform accuracy for nonlinear 3D Volterra integral equations. This numerical scheme is constructed based on the three-dimensional block cubic Lagrangian interpolation method. At the same time, we give the local truncation error analysis of the numerical scheme based on Taylor’s theorem. Through theoretical analysis, we reach the conclusion that the optimal convergence order of this high-order numerical scheme is 4. Finally, we verify the effectiveness and applicability of the method through four numerical examples.
1 June 2020 in Journal of Applied Analysis & Computation