In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermit...In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.展开更多
An Euler wavelets method is proposed to solve a class of nonlinear variable order fractional differential equations in this paper.The properties of Euler wavelets and their operational matrix together with a family of...An Euler wavelets method is proposed to solve a class of nonlinear variable order fractional differential equations in this paper.The properties of Euler wavelets and their operational matrix together with a family of piecewise functions are first presented.Then they are utilized to reduce the problem to the solution of a nonlinear system of algebraic equations.And the convergence of the Euler wavelets basis is given.The method is computationally attractive and some numerical examples are provided to illustrate its high accuracy.展开更多
Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds...Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds in C-n. The Plemelj formula and composite formula of higher order singular integral are obtained. Differential integral equations on smooth closed orientable manifolds are treated by using the composite formula.展开更多
By the generalized Riccati transformation and the integral averaging technique, some sufficient conditions of oscillation of the solutions for second order nonlinear differential equations with damping were discussed....By the generalized Riccati transformation and the integral averaging technique, some sufficient conditions of oscillation of the solutions for second order nonlinear differential equations with damping were discussed.Some sufficient oscillation criteria for previous equations were built up.Some oscillation criteria have been expanded and strengthened in some other known results.展开更多
This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether me...This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether method and the Poisson method. Then the solution of the higher-order equation can be obtained by integrating the solution of the second-order equation.展开更多
This paper studies the existence and uniqueness of solutions for a class of boundary value problems of nonlinear fractional order differential equations involving the Caputo fractional derivative by employing the Ban...This paper studies the existence and uniqueness of solutions for a class of boundary value problems of nonlinear fractional order differential equations involving the Caputo fractional derivative by employing the Banach’s contraction principle and the Schauder’s fixed point theorem. In addition, an example is given to demonstrate the application of our main results.展开更多
In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubi...In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.展开更多
An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
In this paper, a new numerical method for solving fractional differential equations(FDEs) is presented. The method is based upon the fractional Taylor basis approximations. The operational matrix of the fractional int...In this paper, a new numerical method for solving fractional differential equations(FDEs) is presented. The method is based upon the fractional Taylor basis approximations. The operational matrix of the fractional integration for the fractional Taylor basis is introduced. This matrix is then utilized to reduce the solution of the fractional differential equations to a system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of this technique.展开更多
The aim of this paper deals with the study of the Horn matrix function of two complex variables. The convergent properties, an integral representation of H2(A,A′,B,B′;C;z,w) is obtained and recurrence matrix relatio...The aim of this paper deals with the study of the Horn matrix function of two complex variables. The convergent properties, an integral representation of H2(A,A′,B,B′;C;z,w) is obtained and recurrence matrix relations are given. Some result when operating on Horn matrix function with the differential operator D and a solution of certain partial differential equations are established. The Hadamard product of two Horn’s matrix functions is studied, certain results as, the domain of regularity, contiguous functional relations and operating with the differential operator D and D2 are established.展开更多
Based on the linear theories of thin cylindrical shells and viscoelastic materials, a governing equation describing vibration of a sandwich circular cylindrical shell with a viscoelastic core under harmonic excitation...Based on the linear theories of thin cylindrical shells and viscoelastic materials, a governing equation describing vibration of a sandwich circular cylindrical shell with a viscoelastic core under harmonic excitation is derived. The equation can be written as a matrix differential equation of the first order, and is obtained by considering the energy dissipation due to the shear deformation of the viscoelastic core layer and the interaction between all layers. A new matrix method for solving the governing equation is then presented With an extended homogeneous capacity precision integration approach. Having obtained these, vibration characteristics and damping effect of the sandwich cylindrical shell can be studied. The method differs from a recently published work as the state vector in the governing equation is composed of displacements and internal forces of the sandwich shell rather than displacements and their derivatives. So the present method can be applied to solve dynamic problems of the kind of sandwich shells with various boundary conditions and partially constrained layer damping. Numerical examples show that the proposed approach is effective and reliable compared with the existing methods.展开更多
In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation method...In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation meth- ods. We discuss in detail the connections of EFCMs with trigonometric Fourier colloca- tion methods （TFCMs）, the well-known Hamiltonian Boundary Value Methods （HBVMs）, Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an es- sential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first- order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM（2,2） as an illustrative example. The numerical experiments using EFCM（2,2） are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM（2,2）.展开更多
The asymptotic behavior and oscillation of the solutions of second order integro-differential equations with deviating argumentis studied. Our technique depends on an integral inequality containing a deviating argumen...The asymptotic behavior and oscillation of the solutions of second order integro-differential equations with deviating argumentis studied. Our technique depends on an integral inequality containing a deviating argument. From this we obtain some sufficient conditions under which all solutions of Eq.(1.4) have some asymptotic behavior and oscillation.展开更多
By using the integral averaging technique, oscillating criteria for certain even order neutral differential equation with deviating arguments are established. These results extend some known oscillation criteria due t...By using the integral averaging technique, oscillating criteria for certain even order neutral differential equation with deviating arguments are established. These results extend some known oscillation criteria due to Xu and Xia.展开更多
Using the generalized Riccati technique and the averaging technique, we establish several oscillation criteria for certain even order neutral differential equation. The results obtained in this paper extend and improv...Using the generalized Riccati technique and the averaging technique, we establish several oscillation criteria for certain even order neutral differential equation. The results obtained in this paper extend and improve some known results in the previous literatures.展开更多
文摘In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.
基金The authors are grateful to the editor,the associate editor and the anonymous reviewers for their constructive and helpful comments.This work was supported by the Zhejiang Provincial Natural Science Foundation of China(No.LY18A010026),the National Natural Science Foundation of China(No.11701304,11526117)Zhejiang Provincial Natural Science Foundation of China(No.LQ16A010006)+1 种基金the Natural Science Foundation of Ningbo City,China(No.2017A610143)the Natural Science Foundation of Ningbo City,China(2018A610195).
文摘An Euler wavelets method is proposed to solve a class of nonlinear variable order fractional differential equations in this paper.The properties of Euler wavelets and their operational matrix together with a family of piecewise functions are first presented.Then they are utilized to reduce the problem to the solution of a nonlinear system of algebraic equations.And the convergence of the Euler wavelets basis is given.The method is computationally attractive and some numerical examples are provided to illustrate its high accuracy.
基金the Bilateral Science and Technology Collaboration Program of Australia 1998 the Natural Science Foundation of China (No. 1
文摘Using integration by parts and Stokes' formula, the authors give a new definition of Hadamard principal value of higher order singular integrals with Bochner-Martinelli kernel on smooth closed orientable manifolds in C-n. The Plemelj formula and composite formula of higher order singular integral are obtained. Differential integral equations on smooth closed orientable manifolds are treated by using the composite formula.
文摘By the generalized Riccati transformation and the integral averaging technique, some sufficient conditions of oscillation of the solutions for second order nonlinear differential equations with damping were discussed.Some sufficient oscillation criteria for previous equations were built up.Some oscillation criteria have been expanded and strengthened in some other known results.
基金Project supported by the National Natural Science Foundation of China(Grant No10572021)Doctoral Programme Foundation of Institution of Higher Education of China(Grant No20040007022)
文摘This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether method and the Poisson method. Then the solution of the higher-order equation can be obtained by integrating the solution of the second-order equation.
文摘This paper studies the existence and uniqueness of solutions for a class of boundary value problems of nonlinear fractional order differential equations involving the Caputo fractional derivative by employing the Banach’s contraction principle and the Schauder’s fixed point theorem. In addition, an example is given to demonstrate the application of our main results.
文摘In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
文摘In this paper, a new numerical method for solving fractional differential equations(FDEs) is presented. The method is based upon the fractional Taylor basis approximations. The operational matrix of the fractional integration for the fractional Taylor basis is introduced. This matrix is then utilized to reduce the solution of the fractional differential equations to a system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of this technique.
文摘The aim of this paper deals with the study of the Horn matrix function of two complex variables. The convergent properties, an integral representation of H2(A,A′,B,B′;C;z,w) is obtained and recurrence matrix relations are given. Some result when operating on Horn matrix function with the differential operator D and a solution of certain partial differential equations are established. The Hadamard product of two Horn’s matrix functions is studied, certain results as, the domain of regularity, contiguous functional relations and operating with the differential operator D and D2 are established.
基金supported by the National Natural Science Foundation of China (No. 10662003)the Doctoral Fund of Ministry of Education of China (No. 20040787013)
文摘Based on the linear theories of thin cylindrical shells and viscoelastic materials, a governing equation describing vibration of a sandwich circular cylindrical shell with a viscoelastic core under harmonic excitation is derived. The equation can be written as a matrix differential equation of the first order, and is obtained by considering the energy dissipation due to the shear deformation of the viscoelastic core layer and the interaction between all layers. A new matrix method for solving the governing equation is then presented With an extended homogeneous capacity precision integration approach. Having obtained these, vibration characteristics and damping effect of the sandwich cylindrical shell can be studied. The method differs from a recently published work as the state vector in the governing equation is composed of displacements and internal forces of the sandwich shell rather than displacements and their derivatives. So the present method can be applied to solve dynamic problems of the kind of sandwich shells with various boundary conditions and partially constrained layer damping. Numerical examples show that the proposed approach is effective and reliable compared with the existing methods.
文摘In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation meth- ods. We discuss in detail the connections of EFCMs with trigonometric Fourier colloca- tion methods （TFCMs）, the well-known Hamiltonian Boundary Value Methods （HBVMs）, Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an es- sential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first- order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM（2,2） as an illustrative example. The numerical experiments using EFCM（2,2） are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM（2,2）.
文摘The asymptotic behavior and oscillation of the solutions of second order integro-differential equations with deviating argumentis studied. Our technique depends on an integral inequality containing a deviating argument. From this we obtain some sufficient conditions under which all solutions of Eq.(1.4) have some asymptotic behavior and oscillation.
文摘By using the integral averaging technique, oscillating criteria for certain even order neutral differential equation with deviating arguments are established. These results extend some known oscillation criteria due to Xu and Xia.
文摘Using the generalized Riccati technique and the averaging technique, we establish several oscillation criteria for certain even order neutral differential equation. The results obtained in this paper extend and improve some known results in the previous literatures.
