Symplectic Difference Schemes for Linear Hamiltonian Canonical Systems

In this paper, we present some results of a study, specifically within the framework of symplectic geometry, of difference schemes for numerical solution of the linear Hamiltonian systems. We generalize the Cayley transform with which we can get different types of symplectic schemes. These schemes are various generalizations of the Euler centered scheme. They preserve all the invariant first integrals of the linear Hamiltonian systems.     

Explicit Symplectic Difference Schemes for Separable Hamiltonian Systems

This paper is to develop explicit fourth order symplectic difference schemes for separable Hamiltonian systems.     

线性Hamilton系数的几个辛格式

本文利用二次Pad啨逼近的几个结果构造了线性Ｈａｍｉｌｔｏｎ系数的几个辛格式     

Symplectic Collocation Schemes for Hamiltonian Systems

The symplectic collocation schemes, which are based on the framework established by Feng Kang [1], are proposed for numerical solution of Hamiltonian systems. The sufficient and necessary conditions for various collocation schemes to be symplectic are obtained. Some examples of symplectic collocation schemes are also given.     

A Note on Conservation Laws of Symplectic Difference Schemes for Hamiltonian Systems

In this paper we consider the necessary conditions of conservation laws of symplectic difference schemes for Hamiltonian systems and give an example which shows that there does not exist any centered symplectic difference scheme which preserves all Hamiltonian energy.     

A Note on Convergence of Symplectic Schemes for Hamiltonian Systems

In this note we prove that all canonical (or symplectic) schemes for Hamiltonian systems constructed in [1-3] are convergent.     

Symplectic Structure of Discrete Hamiltonian Systems

This paper is concerned with the symplectic structure of discrete nonlinear Hamiltonian systems. The results are related to an open problem that was first proposed by C. D. Ahlbrandt [J. Math. Anal. Appl.180 (1993), 498-517] discussed elsewhere in the literature. But we give a different statement and different proof. Under a solvable condition, we show that the solution operator of a discrete nonlinear Halmiltonian system is symplectic. Then its phase flow is a discrete one-parameter family of symplectic transformations and preserves the phase volume.     

Poisson Difference Schemes for Hamiltonian Systems on Poisson Manifolds

In this paper a systematical method for the construction of Poisson difference schemes with arbitrary order of accuracy for Hamiltonian systems on Poisson manifolds is considered. The transition of such difference schemes from one time-step to the next is a Poisson map. In addition, these schemes preserve all Casimir functions and, under certain conditions, quadratic first integrals of the original Hamiltonian systems. Especially, the arbitrary order centered schemes preserve all Casimir functions and all quadratic first integrals of the original Hamiltonian systems.     

Symplectic Multistep Methods for Linear Hamiltonian Systems

1.IntroductionProhaorFengKangadvancedtheprincipleforconstructiollofsymplecticalgrvrithm8forHarniloniansystemsI11andpointedout.thatsymplecticalgoritlunscanre-ffedmainhauresofHtalltonianSystems,thereforetheyaremoreavailable.Plentyoft~talandrnunericalresultshaveprovedthesepoints.PrO~FengKangalsodiscussedtheaPproalmationproblemsbyalgebraicfull-tfonS.Theconclusionsaxestatedasfolfows[2l:1.wenoteop(f)=p(()/a(f).AmulistepmethodM(p,a)issymplecticforlinearIhailonhaSy8tems(wecallitlinearsymPlectic…     

辛映射与Hamilton控制系统

给定等价辛流形 ,即辛同态或形变等价的辛流形 ,研究了建立在这些辛流形上的Hamilton控制系统之间的一些性质的联系 ,诸如 (局部 )能观测性 ,强可接近性 ,(拟 )极小性等 .而且 ,利用Cort啨ｓ介绍的 (弱 )外等价系统的概念 ,给出使得两个Hamilton控制系统是辛同态的一个充分条件     

Essentially Symplectic Boundary Value Methods for Linear Hamiltonian Systems

1.IlltroductiollInmanyareasofphysics,mechanics,etc.,HamiltoniansystemsofODEsplayaveryimportantrole.Suchsystemshavethefollowinggeneralform:where,bydenotingwithOfandimthenullmatrixandtheidentitymatrixofordermarespectively,SimplepropertiesofthematrixJZmarethefollowingones:Inequation(1)AH(~,t)isthegradientofascalarfunctionH(y,t),usuallycalledHamiltonian.InthecasewhereH(y,t)=H(y),thenthevalueofthisfunctionremainsconstantalongt.hesollltion7/(t),t,hatis'*ReceivedFebruaryI3,1995.l)Worksupporte…     

线性Hamilton系统边值问题的保辛数值方法

以Hamilton系统的正则变换和生成函数为基础研究线性时变Hamilton系统边值问题的保辛数值求解算法.根据第二类生成函数系数矩阵与状态传递矩阵的关系,构造了生成函数系数矩阵的区段合并递推算法,并进一步将递推算法推广到线性非齐次边值问题中;然后利用生成函数的性质将边值问题转化为初值问题,最后采用初值问题的保辛算法求解以达到整个Hamilton系统保辛的目的.数值算例表明该方法能够有效地求解线性齐次与非齐次问题,并能很好地保持Hamilton系统的固有特性.     

Construction of Canonical Difference Schemes for Hamiltonian Formalism via Generating Functions

This paper discusses the relationship between canonical maps and generating functions and gives the general Hamilton-Jacobi theory for time-independent Hamiltonian systems. Based on this theory, the general method — the generating function method — of the construction of difference schemes for Hamiltonian systems is considered. The transition of such difference schemes from one time-step to the next is canonical. So they are called the canonical difference schemes. The well known Euler centered scheme is a canonical difference scheme. Its higher order canonical generalisations and other families of canonical difference schemes are given. The construction method proposed in the paper is also applicable to time-dependent Hamiltonian systems.     

SYMPLECTIC SCHEMES FOR NONAUTONOMOUS HAMILTONIAN SYSTEM

ＳＹＭＰＬＥＣＴＩＣＳＣＨＥＭＥＳＦＯＲＮＯＮＡＵＴＯＮＯＭＯＵＳ　ＨＡＭＩＬＴＯＮＩＡＮ　ＳＹＳＴＥＭＱＩＮＭＥＮＧＺＨＡＯ（秦孟兆）（ＩｎｓｔｉｔａｔｅｏｆＣｏｍｐｕｔａｔｉｏｎａｌＭａｔｈｅｍａｔｉｃｓａｎｄＳｃｉｅｎｔｉｆｉｃ－Ｅｎｇｉｎｅｅ...     

High-Order Symplectic and Symmetric Composition Methods for Multi-Frequency and Multi-Dimensional Oscillatory Hamiltonian Systems

The multi-frequency and multi-dimensional adapted Runge-Kutta-Nyström (ARKN) integrators, and multi-frequency and multi-dimensional extended Runge-Kutta-Nyström(ERKN) integrators have been developed to efficiently solve multi-frequency oscillatory Hamiltonian systems. The aim of this paper is to analyze and derive high-order symplectic and symmetric composition methods based on the ARKN integrators and ERKN integrators. We first consider the symplecticity conditions for the multi-frequency and multi-dimensional ARKN integrators. We then analyze the symplecticity of the adjoint integrators of the multi-frequency and multi-dimensional symplectic ARKN integrators and ERKN integrators, respectively. On the basis of the theoretical analysis and by using the idea of composition methods, we derive and propose four new high-order symplectic and symmetric methods for the multi-frequency oscillatory Hamiltonian systems. The numerical results accompanied in this paper quantitatively show the advantage and efficiency of the proposed high-order symplectic and symmetric methods.     

Symplectic Topology of Integrable Hamiltonian Systems, II: Topological Classification

The main purpose of this paper is to give a topological and symplectic classification of completely integrable Hamiltonian systems in terms of characteristic classes and other local and global invariants.     

斜对角无穷维Hamilton算子的点谱和特征函数系辛结构的非退化性

本文研究斜对角无穷维Hamilton算子$H=\begin{pmatrix}0&B\\C&0\end{pmatrix}$的点谱和特征函数系辛结构的非退化性, 给出斜对角无穷维Hamilton算子$H$的特征函数系具有非退化辛结构的充分必要条件. 基于此, 进一步刻画了斜对角无穷维Hamilton算子$H$的点谱分别包含于实轴、虚轴以及其它区域的充分必要条件. 最后, 以板弯曲问题和弦振动问题中导出的斜对角无穷维Hamilton算子为例, 验证了所得结论的正确性.     

Difference Method of General Schemes with Intrinsic Parallelism for One-dimensional Quasilinear Parabolic Systems with Bounded Measurable Coefficients

The difference method of the general finite difference schemes with intrinsic parallelism for the boundary value problem of the quasilinear parabolic system is studied without assuming heuristically that the original boundary value problem has the unique smooth vector solution.     

Hamilton系统下基于相位误差的精细辛算法

Hamilton系统是一类重要的动力系统，辛算法（如生成函数法、SRK法、SPRK法、多步法等）是针对Hamilton系统所设计的具有保持相空间辛结构不变或保Hamilton函数不变的算法.但是，时域上，同阶的辛算法与Runge-Kutta法具有相同的数值精度，即辛算法在计算过程中也存在相位误差，导致时域上解的数值精度不高.经过长时间计算后，计算结果在时域上也会变得“面目全非”.为了提高辛算法在时域上解的精度，将精细算法引入到辛差分格式中，提出了基于相位误差的精细辛算法（HPD-symplectic method），这种算法满足辛格式的要求，因此在离散过程中具有保Hamilton系统辛结构的优良特性.同时，由于精细化时间步长，极大地减小了辛算法的相位误差，大幅度提高了时域上解的数值精度，几乎可以达到计算机的精度，误差为O（10-13）.对于高低混频系统和刚性系统，常规的辛算法很难在较大的步长下同时实现对高低频精确仿真，精细辛算法通过精细计算时间步长，在大步长情况下，没有额外增加计算量，实现了高低混频的精确仿真.数值结果验证了此方法的有效性和可靠性.