HERMITE MATRIX POLYNOMIALS AND SECOND ORDER MATRIX DIFFERENTIAL EQUATIONS

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.

Solving the Nonlinear Variable Order Fractional Differential Equations by Using Euler Wavelets

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 DIFFERENTIAL INTEGRAL EQUATIONS ONSMOOTH 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.

OSCILLATION CRITERIA FOR SECOND ORDER NONLINEAR DIFFERENTIAL EQUATION WITH DAMPING

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.

Noether's and Poisson's methods for solving differential equation x_s^((m))=F_s(t,x_k^((m-2)) ,x_k^((m-1)))

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.

Existence and Uniqueness for the Boundary Value Problems of Nonlinear Fractional Differential 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.

AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

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.

Improved precise integration method for differential Riccati equation

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

Numerical Solutions of Fractional Differential Equations by Using Fractional Taylor Basis

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.

On Horn Matrix Function <i>H₂(A,A′,B,B′;C;z,w)</i>of Two Complex Variables under Differential Operator

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.

精细积分方法的发展与扩展应用

钟万勰院士于1991年首先提出计算矩阵指数的精细积分方法,其要点是2N类算法和增量存储。精细积分方法可给出矩阵指数在计算机意义上的精确解,为常微分方程的数值计算提供了高精度、高稳定性的算法,现已成功应用于结构动力响应、随机振动、热传导以及最优控制等众多领域。本文首先介绍矩阵指数精细积分方法的提出、基本思想和发展;然后依次介绍在时不变/时变线性微分方程、非线性微分方程以及大规模问题求解中发展起来的各种精细积分方法,分析了其优缺点和适用范围;最后介绍了精细积分方法的基本思想在两点边值问题、椭圆函数和病态代数方程等问题的扩展应用,进一步展示了该思想的特色。

结构动力问题的高阶精确时步群积分方法

高阶精确时间积分方法可为与时间相关的高频复杂动力行为提供高精度的预测结果,但既有高阶精确时间积分方法普遍存在计算工作量偏大的问题,难以满足实际结构线性和非线性动力分析日益增长的计算需求。该文提出了一种基于时步群的高阶精确时间积分方法,将p(p≥2)个相邻的未知时步组成待求解的时步群,以结构动力方程积分解为基础构建逐时步群求解结构动态响应的时间积分方案。在对每个时步群进行积分的过程中,无需联立求解方程,仅通过矩阵乘法运算即可一次性地计算得到时步群内全部p个时步的动态响应。数值特性分析以及线性与非线性算例试验均表明,该文算法精度高、稳定性好、数值耗散可控,在选择较大的时间步距情形下依然能够稳定地获得高精度的计算结果。相较传统二阶精度时间积分方法,该文算法的计算效率亦有较大幅度地提高。

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

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

Fredholm Integro-Differential型方程的Legendre小波方法

研究Legendre小波方法求解具有一阶导和二阶导类型的线性Fredhol mintegro-differential型方程,应用Leg-endre小波逼近法将这2类方程分别化为代数方程求解.实例说明,Legendre小波在解决这2类方程时具可行性和有效性.

非线性二阶变系数微分方程的三点边值问题

研究了非线性二阶变系数微分方程的三点边值问题.首先,对非线性二阶变系数微分方程多次积分得到与之等价的Fredholm-Hammerstein积分方程;其次,利用分段泰勒级数得到Fredholm-Hammerstein积分方程的数值解;最后,通过具体算例验证此方法的可行性与有效性,并给出相应的误差估计.

New matrix method for analyzing vibration and damping effect of sandwich circular cylindrical shell with viscoelastic core

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.

EXPONENTIAL FOURIER COLLOCATION METHODS FOR SOLVING FIRST-ORDER DIFFERENTIAL EQUATIONS

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）.

基于高阶动力格式的基底隔震结构地震响应分析

隔震支座在强烈地震作用下将产生复杂的动力响应,精确掌握地震响应行为及其规律是研发高质量隔震支座和进行隔震结构设计的基础。在进行结构地震响应分析,特别是非线性分析时,采用传统逐步积分方法以及低阶格式建立起的动力时程分析方法都难以达到高精度,而采用传统高阶时间积分算法则会出现计算过程繁琐,工作量大等情况。本文采用一种基于高阶动力格式的显式时间积分新式算法,应用于基底隔震结构的地震响应分析中。首先通过自由度凝聚缩减原理,建立设置隔震支座钢筋混凝土框架结构简化计算模型;同时采用隔震支座恢复力双切线本构模型构建基底隔震体系运动方程,并转化为状态方程;进而基于微分求积(DQ)原理编制高阶动力格式算法求解状态方程;最后分析一高层基底隔震结构,在高、中、低不同频率地震波下地震响应,验证该方法适用性。结果表明,①利用静态凝聚技术建立隔震结构模型,从而构建地震响应状态方程,是高阶动力分析方法实施的基础。②结构高阶动力分析方法仅通过矩阵乘法运算即可获得结构动力响应,可较好且高效地用于解决非线性地震响应分析问题。③将结构高阶动力分析方法应用于基底隔震结构,能获得隔震结构较为精确的分析结果。

常数变易法在微分方程求解中的应用探究

常数变易法是解常微分方程行之有效的一种方法,是拉格朗日历经十一年研究的一种特殊的变量代换法。为探究常数变易法的教学拓展,将常数变易法应用于求解线性微分方程组和高阶线性微分方程,通过常数变易过程,给出简洁推演,建立通解公式,并以典型示例,阐明了公式的实际操作过程。

带有p-Laplacian算子的分数阶积分-微分方程边值问题正解的存在性

研究了带有p-Laplacian算子以及变Riemann-Liouville分数阶积分的分数阶积分-微分方程的边值问题,利用锥上的不动点定理,得到了该边值问题正解的存在性结果.