Application of Haar wavelet method to eigenvalue problems of high order differential equations

In this work, we present a computational method for solving eigenvalue problems of high-order ordinary differential equations which based on the use of Haar wavelets. The variable and their derivatives in the governing equations are represented by Haar function and their integral. The first transform the spectral coefficients into the nodal variable values. The second, solve the obtained system of algebraic equation. The efficiency of the method is demonstrated by four numerical examples.     

Extreme eigenvalues of real symmetric Toeplitz matrices

We exploit the even and odd spectrum of real symmetric Toeplitz matrices for the computation of their extreme eigenvalues, which are obtained as the solutions of spectral, or secular, equations. We also present a concise convergence analysis for a method to solve these spectral equations, along with an efficient stopping rule, an error analysis, and extensive numerical results.

     

一类二阶变系数线性微分方程的解法

介绍如何通过变换把二阶变系数线性微分方程转化为一阶非线性微分方程,进而利用待定系数法对其求解,并对二阶变系数线性微分方程与一阶常系数非线性微分方程的内在的关系进行讨论.     

A new Jacobi rational–Gauss collocation method for numerical solution of generalized pantograph equations

This paper is concerned with a generalization of a functional differential equation known as the pantograph equation which contains a linear functional argument. In this article, a new spectral collocation method is applied to solve the generalized pantograph equation with variable coefficients on a semi-infinite domain. This method is based on Jacobi rational functions and Gauss quadrature integration. The Jacobi rational-Gauss method reduces solving the generalized pantograph equation to a system of algebraic equations. Reasonable numerical results are obtained by selecting few Jacobi rational–Gauss collocation points. The proposed Jacobi rational–Gauss method is favorably compared with other methods. Numerical results demonstrate its accuracy, efficiency, and versatility on the half-line.     

Controllability method for acoustic scattering with spectral elements

We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improvements in efficiency due to the higher order spectral elements. For a given accuracy, the controllability technique with spectral element method requires fewer computational operations than with conventional finite element method. In addition, by using higher order polynomial basis the influence of the pollution effect is reduced.     

Solving 2D and 3D Poisson equations and biharmonic equations by the Haar wavelet method

In this paper, we present a computational method for solving 2D and 3D Poisson equations and biharmonic equations which based on the use of Haar wavelets. The highest derivative appearing in the differential equation is expanded into the Haar series, this approximation is integrated while the boundary conditions are incorporated by using integration constants. In 2D the first transform the spectral coefficients into the nodal variable values and then use Kronecker products to construct the approximations for derivatives over a tensor product grid of the horizontal and vertical blocks. Finally, solutions to four test problems are investigated.     

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 modification of the invariant imbedding method for a singular boundary value problem

We consider a dynamically-consistent analytical model of a 3D topographic vortex. The model is governed by equations derived from the classical problem of the axisymmetric Taylor–Couette flow. Using linear expansions, these equations can be reduced to a differential sixth-order equation with variable coefficients. For this differential equation, we formulate a boundary value problem, which has a number of issues for numerical solving. To avoid these issues and find the eigenvalues and eigenfunctions of the boundary value problem, we suggest a modification of the invariant imbedding method (the Riccati equation method). In this paper, we show that such a modification is necessary since the boundary conditions possess singular matrices, which sufficiently complicate the derivation of the Riccati equation. We suggest algebraic manipulations, which permit the initial problem to be reduced to a problem with regular boundary conditions. Also, we propose a method for obtaining a numerical solution of the matrix Riccati equation by means of recurrence relations, which allow us to obtain a matrizer converging to the required eigenfunction. The suggested method is tested by calculating the corresponding eigenvalues and eigenfunctions, and then, by constructing fluid particle trajectories on the basis of the eigenfunctions.     

A new version of the Smith method for solving Sylvester equation and discrete-time Sylvester equation

Recently, Xue etc. \cite{28} discussed the Smith method for solving Sylvester equation $AX+XB=C$, where one of the matrices $A$ and $B$ is at least a nonsingular $M$-matrix and the other is an (singular or nonsingular) $M$-matrix. Furthermore, in order to find the minimal non-negative solution of a certain class of non-symmetric algebraic Riccati equations, Gao and Bai \cite{gao-2010} considered a doubling iteration scheme to inexactly solve the Sylvester equations. This paper discusses the iterative error of the standard Smith method used in \cite{gao-2010} and presents the prior estimations of the accurate solution $X$ for the Sylvester equation. Furthermore, we give a new version of the Smith method for solving discrete-time Sylvester equation or Stein equation $AXB+X=C$, while the new version of the Smith method can also be used to solve Sylvester equation $AX+XB=C$, where both $A$ and $B$ are positive definite. % matrices. We also study the convergence rate of the new Smith method. At last, numerical examples are given to illustrate the effectiveness of our methods     

Numerical solution of a linear elliptic partial differential equation with variable coefficients: A complex variable boundary element approach

This article presents a complex variable boundary element method for the numerical solution of a second order elliptic partial differential equation with variable coefficients. To assess the validity and accuracy of the method, it is applied to solve some specific problems with known solutions. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Efficient spectral ultraspherical-Galerkin algorithms for the direct solution of 2nth-order linear differential equations

Some efficient and accurate algorithms based on the ultraspherical-Galerkin method are developed and implemented for solving 2nth-order linear differential equations in one variable subject to homogeneous and nonhomogeneous boundary conditions using a spectral discretization. We extend the proposed algorithms to solve the two-dimensional 2nth-order differential equations. The key to the efficiency of these algorithms is to construct appropriate base functions, which lead to linear systems with specially structured matrices that can be efficiently inverted, hence greatly reducing the cost and roundoff errors.     

A numerical study for the KdV and the good Boussinesq equations using Fourier Chebyshev tau meshless method

The current paper presents a scheme, which combines Fourier spectral method and Chebyshev tau meshless method based on the highest derivative (CTMMHD) to solve the nonlinear KdV equation and the good Boussinesq equation. Fourier spectral method is used to approximate the spatial variable, and the problem is converted to a series of equations with Fourier coefficients as unknowns. Then, CTMMHD is applied blockwise in time direction. For the long time computing of solitons, we introduce the computational area moving technique. The numerical results show that the accuracy of Fourier-CTMMHD is good and the computational area moving technique makes the long-time numerical behavior well for the problems with solitons moving towards the same direction.     

Boundary‐value problems for the third‐order loaded equation with noncharacteristic type‐change boundaries

The purpose of this paper is to establish unique solvability for a certain generalized boundary‐value problem for a loaded third‐order integro‐differential equation with variable coefficients. Moreover, the method of integral equations is applied to obtain an equation related to the Riemann‐Liouville operators.     

Operational matrices of Chebyshev cardinal functions and their application for solving delay differential equations arising in electrodynamics with error estimation

In this paper, a new and effective direct method to determine the numerical solution of pantograph equation, pantograph equation with neutral term and Multiple-delay Volterra integral equation with large domain is proposed. The pantograph equation is a delay differential equation which arises in quite different fields of pure and applied mathematics, such as number theory, dynamical systems, probability, mechanics and electrodynamics. The method consists of expanding the required approximate solution as the elements of Chebyshev cardinal functions. The operational matrices for the integration, product and delay of the Chebyshev cardinal functions are presented. A general procedure for forming these matrices is given. These matrices play an important role in modelling of problems. By using these operational matrices together, a pantograph equation can be transformed to a system of algebraic equations. An efficient error estimation for the Chebyshev cardinal method is also introduced. Some examples are given to demonstrate the validity and applicability of the method and a comparison is made with existing results.     

一阶线性双曲方程的耗散谱元法

讨论了二维一阶线性变系数双曲方程的耗散谱元法,得到拟最优估计.数值结果表明,耗散谱元法对于具有较复杂边界条件的问题同样有效,对于有限光滑问题,耗散谱元法能够得到比传统的谱元法更好的结果.     

Iterative solution of some nonlinear differential equations

We propose an iterative method to solve some non-linear ordinary differential equations. Comparing on the Mathieu, van der Pol and Hill equation of fourth order, we see that this method is much more efficient than the well known methods by Lyapunov or Picard.     

On an inverse problem of reconstructing a subdiffusion process from nonlocal data

We consider a problem of modeling the thermal diffusion process in a closed metal wire wrapped around a thin sheet of insulation material. The layer of insulation is assumed to be slightly permeable. Therefore, the temperature value from one side affects the diffusion process on the other side. For this reason, the standard heat equation is modified, and a third term with an involution is added. Modeling of this process leads to the consideration of an inverse problem for a one‐dimensional fractional evolution equation with involution and with periodic boundary conditions with respect to a space variable. This equation interpolates heat equation. Such equations are also called nonlocal subdiffusion equations or nonlocal heat equations. The inverse problem consists in the restoration (simultaneously with the solution) of the unknown right‐hand side of the equation, which depends only on the spatial variable. The conditions for overdefinition are initial and final states. Existence and uniqueness results for the given problem are obtained via the method of separation of variables.     

A lattice Boltzmann model for the Fokker-Planck equation

In this paper, a lattice Boltzmann model is presented for solving one and two-dimensional Fokker-Planck equations with variable coefficients. In particular, it is efficient to simulate one-dimensional stochastic processes governed by the Fokker-Planck equation. Numerical results agree well with the exact solutions, which indicates that the proposed model is suitable for solving the Fokker-Planck equation.     

二阶非线性中立型时滞差分方程的振动性

研究了一类具有多个变滞量的变系数的二阶非线性中立型时滞差分方程的振动性,得到了该类方程振动及其解的一阶差分振动的充分条件,推广了现有文献中的某些结果.     

A Hermite collocation method for the approximate solutions of high‐order linear Fredholm integro‐differential equations

In this study, a Hermite matrix method is presented to solve high‐order linear Fredholm integro‐differential equations with variable coefficients under the mixed conditions in terms of the Hermite polynomials. The proposed method converts the equation and its conditions to matrix equations, which correspond to a system of linear algebraic equations with unknown Hermite coefficients, by means of collocation points on a finite interval. Then, by solving the matrix equation, the Hermite coefficients and the polynomial approach are obtained. Also, examples that illustrate the pertinent features of the method are presented; the accuracy of the solutions and the error analysis are performed. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1707–1721, 2011