离散大系统非线性比较方程的稳定性

用矢量李雅普诺夫函数解决大系统的稳定性问题必须要判断矢量比较方程的稳定性.对离散系统,过去只研究过线性驻定比较方程的稳定性.本文全面建立了离散非线性驻定比较方程的各种稳定性判别准则,其中渐近稳定的准则既是充分也是必要的,并由此推得了一个用于C₁类函数的准则,两者均可用来判断离散非线性(驻定或非驻定)系统的非指数稳定以至全局非指数稳定.所有准则均具有简单的代数形式,便于应用.     

变时滞退化Lurie控制系统的绝对稳定性

研究了变时滞退化Lurie控制系统的绝对稳定性问题.基于Lyapunov稳定性理论,利用线性矩阵不等式方法给出了系统绝对稳定的判别准则.讨论了变时滞退化Lurie直接控制系统和间接控制系统的绝对稳定性,得到绝对稳定性的充分条件仅依赖于时滞导数的大小,且时滞可以是无界函数：最后给出了实例说明本文结果的有效性.     

线性随机变时滞微分方程指数Euler方法的收敛性和稳定性

文应用指数Euler方法研究了线性随机变时滞微分方程的收敛性和稳定性;首先,证明了指数Euler方法是$\frac{1}{2}$阶均方收敛的;其次,在解析解均方稳定的前提下,通过跟Euler-Maruyama方法比较发现指数Euler方法在大步长下依然保持解析解的均方稳定性;最后,用数值试验验证了收敛和稳定的结果.     

预估-校正方法的绝对稳定性讨论

预估-校正方法,即PECE方法,常被用于求解常微分方程的初值问题.而一般文献中常只讨论了单个线性多步法公式的稳定性问题,很少涉及由一个显式公式和一个隐式公式组合而成的PECE方法的稳定性.本文应用根轨迹法和对分法讨论了常用的PECE方法的稳定性,求出了一些常用PECE方法的组合公式的绝对稳定区间和绝对稳定区域,并用数值...     

关于一类非线性控制系统绝对稳定性的研究

通过构造Liapunov函数讨论了一类非线性非定常控制系统并得出其零解绝对稳定性的条件,推广了相关结论.     

求解一类非线性偏微分方程的异步并行算法

本文针对一类非线性偏微分方程,把区域分裂法与异步混合算法结合起来,产生了一种异步混合单调算法,证明了收敛性定理。我们已经用区域分裂技术建立了一类求解偏微分方程的异步并行算法。这类方法成功地用于带有线性或非线性边界条件的线性或拟线性方程的定解问题,那么怎样处理较一般的非线性问题呢?     

中立型延迟微分方程一般线性方法的非线性稳定性

本文研究一类非线性中立型延迟微分方程一般线性方法的数值稳定性.证明了一般线性方法为(k,p,O)-代数稳定时,在一定的约束条件下,其数值解保持微分方程理论解的稳定性质,特别是证明了在约束网格情形代数靛的-般线性方法能无条件保持解析解的稳定性.     

大系统渐近稳定的一般判别定理

利用分解集结和向量 V 函数判别大系统的稳定性是一个广泛采用的有效方法.但过去一般限于集结成常系数线性比较方程,判定的只是指数稳定.本文提出一个非线性比较方程的构造定理,推广了 Bailey 方法,并且进一步推广了作者的前期工作,判定了大系统的非指数稳定.     

多时滞退化Lurie控制系统的绝对稳定性

本文研究退化Lurie直接控制时滞系统和退化Lurie间接控制时滞系统的绝对稳定性,基于Lyapunov稳定性理论和线性矩阵不等式得到绝对稳定的判据,即绝对稳定性仅仅依赖于时滞导数的大小.特别地,时滞可以是无界的函数,仿真示例同时说明了此方法的有效性.     

非自治刚性随机微分方程正则EM分裂方法的收敛性和稳定性

本文研究了数值求解非自治随机微分方程的正则Euler-Maruyama分裂（CEMS）方法,该方程的漂移项系数带有刚性且允许超线性增长,扩散项系数满足全局Lipschitz条件.首先,证明了CEMS方法的强收敛性及收敛速度.其次,证明了在适当条件下CEMS方法是均方稳定的.进一步,利用离散半鞅收敛定理,研究了CEMS方法的几乎必然指数稳定性.结果表明,CEMS方法在漂移系数的刚性部分满足单边Lipschitz条件下可保持几乎必然指数稳定性.最后通过数值实验,检验了CEMS方法的有效性并证实了我们的理论结果.     

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…     

Numerical Stability and Oscillations of Runge-Kutta Methods for Differential Equations with Piecewise Constant Arguments of Advanced Type

For differential equations with piecewise constant arguments of advanced type, numerical stability and oscillations of Runge-Kutta methods are investigated. The necessary and sufficient conditions under which the numerical stability region contains the analytic stability region are given. The conditions of oscillations for the Runge-Kutta methods are obtained also. We prove that the Runge-Kutta methods preserve the oscillations of the analytic solution. Moreover, the relationship between stability and oscillations is discussed. Several numerical examples which confirm the results of our analysis are presented.     

Fourth-Order Splitting Methods for Time-Dependant Differential Equations

This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.     

New open modified Newton Cotes type formulae as multilayer symplectic integrators

New modified open Newton Cotes integrators are introduced in this paper. For the new proposed integrators the connection between these new algorithms, differential methods and symplectic integrators is studied. Much research has been done on one step symplectic integrators and several of them have obtained based on symplectic geometry. However, the research on multistep symplectic integrators is very poor. Zhu et al. [1] studied the well known open Newton Cotes differential methods and they presented them as multilayer symplectic integrators. Chiou and Wu [2] studied the development of multistep symplectic integrators based on the open Newton Cotes integration methods. In this paper we introduce a new open modified numerical method of Newton Cotes type and we present it as symplectic multilayer structure. The new obtained symplectic schemes are applied for the solution of Hamilton’s equations of motion which are linear in position and momentum. An important remark is that the Hamiltonian energy of the system remains almost constant as integration proceeds. We have applied also efficiently the new proposed method to a nonlinear orbital problem and an almost periodic orbital problem.     

On the Nonlinear Stability of Symplectic Integrators

The modified Hamiltonian is used to study the nonlinear stability of symplectic integrators, especially for nonlinear oscillators. We give conditions under which an initial condition on a compact energy surface will remain bounded for exponentially long times for sufficiently small time steps. While this is easy to achieve for non-critical energy surfaces, in some cases it can also be achieved for critical energy surfaces (those containing critical points of the Hamiltonian). For example, the implicit midpoint rule achieves this for the critical energy surface of the Hénon–Heiles system, while the leapfrog method does not. We construct explicit methods which are nonlinearly stable for all simple mechanical systems for exponentially long times. We also address questions of topological stability, finding conditions under which the original and modified energy surfaces are topologically equivalent. This revised version was published online in July 2006 with corrections to the Cover Date.     

Geometric integration of non-autonomous linear Hamiltonian problems

Symplectic integration of autonomous Hamiltonian systems is a well-known field of study in geometric numerical integration, but for non-autonomous systems the situation is less clear, since symplectic structure requires an even number of dimensions. We show that one possible extension of symplectic methods in the autonomous setting to the non-autonomous setting is obtained by using canonical transformations. Many existing methods fit into this framework. We also perform experiments which indicate that for exponential integrators, the canonical and symmetric properties are important for good long time behaviour. In particular, the theoretical and numerical results support the well documented fact from the literature that exponential integrators for non-autonomous linear problems have superior accuracy compared to general ODE schemes.     

Volume preserving RK methods for linear systems

1. IntroductionIn recent yearss there has been a great interest in constructing numerical integrationschemes for ODEs in such a way that some qualitative geometrical properties of the solutionof the ODEs are exactly preserved. R.th[ll and Feng Kang[2'31 has proposed symplectic algorithms for Hamiltollian systems, and since then st ruct ure s- preserving me t ho ds fordynamical systems have been systematically developed[4--7]. The symplectic algorithms forHamiltonian systems, the volume-pre…     

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.     

Conserved quantities of some Hamiltonian wave equations after full discretization

Hamiltonian PDEs have some invariant quantities, which would be good to conserve with the numerical integration. In this paper, we concentrate on the nonlinear wave and Schrödinger equations. Under hypotheses of regularity and periodicity, we study how a symmetric space discretization makes that the space discretized system also has some invariants or `nearly' invariants which well approximate the continuous ones. We conjecture some facts which would explain the good numerical approximation of them after time integration when using symplectic Runge-Kutta methods or symmetric linear multistep methods for second-order systems.     

Symplectic integration of Hamiltonian wave equations

Summary The numerical integration of a wide class of Hamiltonian partial differential equations by standard symplectic schemes is discussed, with a consistent, Hamiltonian approach. We discretize the Hamiltonian and the Poisson structure separately, then form the the resulting ODE's. The stability, accuracy, and dispersion of different explicit splitting methods are analyzed, and we give the circumstances under which the best results can be obtained; in particular, when the Hamiltonian can be split into linear and nonlinear terms. Many different treatments and examples are compared.