Symmetric and symplectic exponential integrators for nonlinear Hamiltonian systems

This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show that these conditions are extensions of the symmetry and symplecticity conditions of Runge–Kutta methods. Based on these conditions, some symmetric and symplectic exponential integrators up to order four are derived. Two numerical experiments are carried out and the results demonstrate the remarkable numerical behaviour of the new exponential integrators in comparison with some symmetric and symplectic Runge–Kutta methods in the literature.     

Symplectic rigidity,symplectic fixed points,and global perturbations of Hamiltonian systems

In this paper we study a generalized symplectic fixed‐point problem, first considered by J. Moser in [20], from the point of view of some relatively recently discovered symplectic rigidity phenomena. This problem has interesting applications concerning global perturbations of Hamiltonian systems. © 2007 Wiley Periodicals, Inc.     

Hamilton系统MSRK算法的数值色散关系

Numerical dispersion relation of the multi-symplectic Runge-Kutta (MSRK) method for linear Hamiltonian PDEs is derived in the present paper, which is shown to be a discrete counterpart to that possessed by the differential equation. This provides further understanding of MSRK methods. However, much still remains to be investigated further.     

Lyapunov inequalities and stability for discrete linear Hamiltonian systems

In this paper, we establish several new Lyapunov type inequalities for discrete linear Hamiltonian systems when the end-points are not necessarily usual zeros, but rather, generalized zeros, which generalize and improve almost all related existing ones. Applying these inequalities, an optimal stability criterion is obtained.     

Preservation of equilibria for symplectic methods applied to Hamiltonian systems

In this paper, the linear stability of symplectic methods for Hamiltonian systems is studied. In par- ticular, three classes of symplectic methods are considered： symplectic Runge-Kutta （SRK） methods, symplectic partitioned Runge-Kutta （SPRK） methods and the composition methods based on SRK or SPRK methods. It is shown that the SRK methods and their compositions preserve the ellipticity of equilibrium points uncondi- tionally, whereas the SPRK methods and their compositions have some restrictions on the time-step.     

A remark on focal points of recessive solutions of discrete symplectic systems

We investigate nonoscillatory and controllable symplectic difference systems. We show that the recessive solution of such a system at +∞ has the same number of focal points (counting multiplicities) as the recessive solution at −∞.     

线性矩阵哈密顿系统的振动性

研究了线性矩阵 Hamilton系统X′=A( t) X + B( t) YY′=C( t) X -A*( t) Y　　 t≥ 0的振动性 .其中 A( t) ,B( t) ,C( t) ,X,Y为实 n× n矩阵值函数 ,B,C为对称矩阵 ,B正定 .借助于正线性泛函 ,采用加权平均法 ,得到了该系统的非平凡预备解的振动性 .这些结果推广、改进了许多已知的结果     

Regular reduction of controlled Hamiltonian system with symplectic structure and symmetry

In this paper, our goal is to study the regular reduction theory of regular controlled Hamiltonian (RCH) systems with symplectic structure and symmetry, and this reduction is an extension of regular symplectic reduction theory of Hamiltonian systems under regular controlled Hamiltonian equivalence conditions. Thus, in order to describe uniformly RCH systems defined on a cotangent bundle and on the regular reduced spaces, we first define a kind of RCH systems on a symplectic fiber bundle. Then we introduce regular point and regular orbit reducible RCH systems with symmetry by using momentum map and the associated reduced symplectic forms. Moreover, we give regular point and regular orbit reduction theorems for RCH systems to explain the relationships between RpCH-equivalence, RoCH-equivalence for reducible RCH systems with symmetry and RCH-equivalence for associated reduced RCH systems. Finally, as an application we regard rigid body and heavy top as well as them with internal rotors as the regular point reducible RCH systems on the rotation group SO(3) and on the Euclidean group SE(3), as well as on their generalizations, respectively, and discuss their RCH-equivalence. We also describe the RCH system and RCH-equivalence from the viewpoint of port Hamiltonian system with a symplectic structure.     

Conservation properties of numerical integrators for highly oscillatory Hamiltonian systems

^** Email: David.Cohen{at}math.unige.ch. Present address: Mathematisches Institut, Universität Tübingen, D-72076 Tübingen, Germany (cohen{at}na.uni-tuebingen.de) Modulated Fourier expansion is used to show long-time near-conservationof the total and oscillatory energies of numerical methods forHamiltonian systems with highly oscillatory solutions. The numericalmethods considered are an extension of the trigonometric methods.A brief discussion of conservation properties in the continuousproblem and in the multi-frequency case is also given.     

Focal points and principal solutions of linear Hamiltonian systems revisited

In this paper we present a novel view on the principal (and antiprincipal) solutions of linear Hamiltonian systems, as well as on the focal points of their conjoined bases. We present a new and unified theory of principal (and antiprincipal) solutions at a finite point and at infinity, and apply it to obtain new representation of the multiplicities of right and left proper focal points of conjoined bases. We show that these multiplicities can be characterized by the abnormality of the system in a neighborhood of the given point and by the rank of the associated T-matrix from the theory of principal (and antiprincipal) solutions. We also derive some additional important results concerning the representation of T-matrices and associated normalized conjoined bases. The results in this paper are new even for completely controllable linear Hamiltonian systems. We also discuss other potential applications of our main results, in particular in the singular Sturmian theory.     

Spectral properties of singular discrete linear Hamiltonian systems

This paper is concerned with spectral properties of singular discrete linear Hamiltonian systems. It is shown that properties of the essential spectrum of each self-adjoint subspace extension (SSE) of the corresponding minimal subspace are independent of the values of the coefficients of the system on any finite subinterval. The analyticity of the Weyl function is studied by employing the Schwarz reflection principle for the system in the limit point case. Based on the above work, several sufficient conditions are obtained for each SSE to have no essential spectrum points in an interval of the real line in the strong limit point case, and then a sufficient condition for the essential spectrum to be bounded from below (above) and a sufficient condition for the pure discrete spectrum are presented, respectively. As a direct consequence, the related spectral properties of singular higher order symmetric vector difference expressions are given.     

Stability of linear neutral differential equations with delays and impulses established by the fixed points method

In this paper, we consider the impulsive effects on the stability of the zero solution of the linear neutral differential equations with variable delays. By transforming the equations into ones without impulses and using fixed point theory, some sufficient conditions for asymptotic stability and exponential stability of the zero solution are obtained. The paper extends and improves results on sufficient conditions obtained by Jin and Luo (2008) [17], and Ardjouni and Djoudi (2011) [18], which is shown clearly in Example 1. This paper also shows that the impulse intensity and the impulse time both influence the decay rate of the convergence to zero of the solutions. Finally, two examples are given to show applications of some results obtained.     

Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models

In this paper we are concerned with the fractional-order predator-prey model and the fractional-order rabies model. Existence and uniqueness of solutions are proved. The stability of equilibrium points are studied. Numerical solutions of these models are given. An example is given where the equilibrium point is a centre for the integer order system but locally asymptotically stable for its fractional-order counterpart.     

On the stability of iterative approximations to fixed points of nonexpansive mappings

We study iterative retraction approximations to fixed points of the nonexpansive self-mapping given on the closed convex set G in a Banach space B. The conditions which guarantee weak and strong convergence and stability of these approximations with respect to perturbations of both operator A and constraint set G are considered. The results of this paper are new even in a Hilbert space for the iterative projection approximations.     

The asymptotic stability of theoretical and numerical solutions for systems of neutral multidelay-differential equations

The asymptotic stability of theoretical and numerical solutions for neutral multidelay-differential equations (NMDEs) is dealt with. A sufficient condition on the asymptotic stability of theoretical solutions for NMDEs is obtained. On the basis of this condition, it is proved that A-stability of the multistep Runge-Kutta methods for ODEs is equivalent to NGP_k-stability of the induced methods for NMDEs. Project supported by the National Natural Science Foundation of China (Grant No. 19771034).     

On data-dependence of exponential stability and stability radii for linear time-varying differential-algebraic systems

This paper is addressed to some questions concerning the exponential stability and its robustness measure for linear time-varying differential-algebraic systems of index 1. First, the Bohl exponent theory that is well known for ordinary differential equations is extended to differential-algebraic equations. Then, it is investigated that how the Bohl exponent and the stability radii with respect to dynamic perturbations for a differential-algebraic system depend on the system data. The paper can be considered as a continued and complementary part to a recent paper on stability radii for time-varying differential-algebraic equations [N.H. Du, V.H. Linh, Stability radii for linear time-varying differential-algebraic equations with respect to dynamic perturbations, J. Differential Equations 230 (2006) 579-599].     

On the uniform exponential stability of a wide class of linear time-delay systems

This paper deals with the global uniform exponential stability independent of delay of time-delay linear and time-invariant systems subject to point and distributed delays for the initial conditions being continuous real functions except possibly on a set of zero measure of bounded discontinuities. It is assumed that the delay-free system as well as an auxiliary one are globally uniformly exponentially stable and globally uniform exponential stability independent of delay, respectively. The auxiliary system is typically a part of the overall dynamics of the delayed system but not necessarily the isolated undelayed dynamics as usually assumed in the literature. Since there is a great freedom in setting such an auxiliary system, the obtained stability conditions are very useful in a wide class of practical applications.     

On stability of equilibrium points in nonlinear fractional differential equations and fractional Hamiltonian systems

In this article, a brief stability analysis of equilibrium points in nonlinear fractional order dynamical systems is given. Then, based on the first integral concept, a definition of planar Hamiltonian systems with fractional order introduced. Some interesting properties of these fractional Hamiltonian systems are also presented. Finally, we illustrate two examples to see the differences between fractional Hamiltonian systems with their classical order counterparts.© 2014 Wiley Periodicals, Inc. Complexity 21: 93–99, 2015     

On the stability of linear systems with an exact constraint set

This paper deals with the stability of the intersection of a given set with the solution, , of a given linear system whose coefficients can be arbitrarily perturbed. In the optimization context, the fixed constraint set X can be the solution set of the (possibly nonlinear) system formed by all the exact constraints (e.g., the sign constraints), a discrete subset of (as or { 0,1} ⁿ , as it happens in integer or Boolean programming) as well as the intersection of both kind of sets. Conditions are given for the intersection to remain nonempty (or empty) under sufficiently small perturbations of the data. Research supported by Fondecyt Grant 1020(7020)-646. Research supported by DGES and FEDER, Grant BFM2002-04114-C02-01     

Stochastic stability of positive Markov jump linear systems with fixed dwell time

This paper studies the stochastic stability of positive Markov jump linear systems with a fixed dwell time. By constructing an auxiliary system that originated from the initial system with state jumps, sufficient and necessary conditions of stochastic stability for positive Markov jump linear systems are obtained with both exactly known and partially known transition rates. The main idea in the latter case is applying a convex combination to convert bilinear programming into linear programming problems. On this basis, multiple piecewise linear co-positive Lyapunov functions are provided to achieve less conservative results. Then state feedback controller is designed to stabilize the positive Markov jump linear systems by solving linear programming problems. Numerical examples are presented to illustrate the viability of our conclusions.