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.     

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

Linear stability of general linear methods for systems of neutral delay differential equations

This paper is concerned with the numerical solution of delay differential equations (DDEs). We focus on the stability of general linear methods for systems of neutral DDEs with multiple delays. A type of interpolation procedure is considered for general linear methods. Linear stability properties of general linear methods with this interpolation procedure are investigated. Many extant results are unified.     

On the simultaneous diagonal stability of a pair of positive linear systems

In this paper we derive a necessary and sufficient condition for the existence of a diagonal common quadratic Lyapunov function (CQLF) for a pair of positive linear time-invariant (LTI) systems.     

Numerical modelling of linear and nonlinear diffusion equations by compact finite difference method

In this work, accurate solutions to linear and nonlinear diffusion equations were introduced. A combination of a sixth-order compact finite difference scheme in space and a low-storage third-order total variation diminishing Runge-Kutta scheme in time have been used for treatment of these equations. The computed results with the use of this technique have been compared with the exact solution to show the accuracy of it. Here, the approximate solution to the diffusion equations has been obtained easily and elegantly with neither transforming nor linearizing the equation. The present method is seen to be a very good alternative method to some existing techniques for realistic problems.     

Continuous approximation of linear impulsive systems and a new form of robust stability

The time-scale tolerance for linear ordinary impulsive differential equations is introduced. How large the time-scale tolerance is directly reflects the degree to which the qualitative dynamics of the linear impulsive system can be affected by replacing the impulse effect with a continuous (as opposed to discontinuous, impulsive) perturbation, producing what is known as an impulse extension equation. Theoretical properties related to the existence of the time-scale tolerance are given for periodic systems, as are algorithms to compute them. Some methods are presented for general, aperiodic systems. Additionally, sufficient conditions for the convergence of solutions of impulse extension equations to the solutions of their associated impulsive differential equation are proven. Counterexamples are provided.     

A new augmented Lyapunov-Krasovskii functional approach for stability of linear systems with time-varying delays

This paper proposes improved delay-dependent conditions for asymptotic stability of linear systems with time-varying delays. The proposed method employs a suitable Lyapunov-Krasovskii’s functional for new augmented system. Based on Lyapunov method, delay-dependent stability criteria for the systems are established in terms of linear matrix inequalities (LMIs) which can be easily solved by various optimization algorithms. Three numerical examples are included to show that the proposed method is effective and can provide less conservative results.