含分布时滞的时滞微分系统多步龙格-库塔方法的时滞相关稳定性

本文研究了一类含分布时滞的时滞微分系统的多步龙格-库塔方法的稳定性.基于辐角原理，本文给出了多步龙格-库塔方法弱时滞相关稳定性的充分条件，并通过数值算例验证了理论结果的有效性.     

A Runge-Kutta–based WENO reconstruction on moving mesh for compressible fluids

In this paper, we construct a high-order moving mesh method based on total variation diminishing Runge-Kutta and weighted essential nonoscillatory reconstruction for compressible fluid system. Beginning with the integral form of fluid system, we get the semidiscrete system with an arbitrary mesh velocity. We use weighted essential nonoscillatory reconstruction to get the space accuracy on moving meshes, and the time accuracy is obtained by modified Runge-Kutta method; the mesh velocity is determined by moving mesh method. One- and two-dimensional numerical examples are presented to demonstrate the efficient and accurate performance of the scheme.     

On a certain class of locally conformal symplectic structures of the second kind

We study locally conformal symplectic (LCS) structures of the second kind on a Lie algebra. We show a method to construct new examples of Lie algebras admitting LCS structures of the second kind starting with a lower dimensional Lie algebra endowed with a LCS structure and a suitable extension. Moreover, we characterize all LCS Lie algebras obtained with our construction. Finally, we study the existence of lattices in the associated simply connected Lie groups in order to obtain compact examples of manifolds admitting this kind of structure.     

传输线方程的高精度龙格−库塔数值求解方法

提出了一种求解传输线方程的高精度龙格-库塔(RK)方法。此方法在空间上采取高阶泰勒展开,提高了对空间微分的近似精度,减少了数值色散所带来的误差。与传统的时域有限差分法(FDTD)方法相比,在每波长采样数相同时,RK方法的计算精度更高。同时,根据Taylor模型,对外界平面波激励源进行离散,成功利用RK方法对外部场激励传输线进行求解,扩大了龙格-库塔方法在求解传输线方程时的应用范围。通过编程对平面波辐照下无限大地平面上的单导体与双导体的算例分别应用FDTD方法与RK方法进行了计算,验证了RK方法的正确性。结果表明同等计算条件下RK方法的计算精度更高。     

两类五阶解非线性方程组的迭代算法

本文首先根据Runge-Kutta方法的思想,结合Newton迭代法,提出了一类带参数的解非线性方程组F(x)=0的迭代算法,然后基于解非线性方程f(x)=0的King算法,给出第二类解非线性方程组的迭代算法,收敛性分析表明这两类算法都是五阶收敛的.其次给出了本文两类算法的效率指数,以及一些已知算法的效率指数,并且将本文算法的效率指数与其它方法进行详细的比较,通过效率比率R_(i,j)可知本文算法具有较高的计算效率.最后给出了四个数值实例,将本文两类算法与现有的几种算法进行比较,实验结果说明本文算法收敛速度快,迭代次数少,有明显的优势.     

Low-Storage Runge-Kutta Method for Simulating Time-Dependent Quantum Dynamics

A wide range of quantum systems are time-invariant and the corresponding dynamics is dictated by linear differential equations with constant coefficients.Although simple in mathematical concept,the integration of these equations is usually complicated in practice for complex systems,where both the computational time and the memory storage become limiting factors.For this reason,low-storage Runge-Kutta methods become increasingly popular for the time integration.This work suggests a series of s-stage sth-order explicit RungeKutta methods specific for autonomous linear equations,which only requires two times of the memory storage for the state vector.We also introduce a 13-stage eighth-order scheme for autonomous linear equations,which has optimized stability region and is reduced to a fifth-order method for general equations.These methods exhibit significant performance improvements over the previous general-purpose low-stage schemes.As an example,we apply the integrator to simulate the non-Markovian exciton dynamics in a 15-site linear chain consisting of perylene-bisimide derivatives.     

Dominant and recessive solutions at infinity and genera of conjoined bases for discrete symplectic systems

In this paper we introduce the theory of dominant solutions at infinity for nonoscillatory discrete symplectic systems without any controllability assumption. Such solutions represent an opposite concept to recessive solutions at infinity, which were recently developed for such systems by the authors. Our main results include: (i) the existence of dominant solutions at infinity for all ranks in a given range depending on the order of abnormality of the system, (ii) construction of dominant solutions at infinity with eventually the same image, (iii) classification of dominant and recessive solutions at infinity with eventually the same image, (iv) limit characterization of recessive solutions at infinity in terms of dominant solutions at infinity and vice versa, and (v) Reid’s construction of the minimal recessive solution at infinity. These results are based on a new theory of genera of conjoined bases for symplectic systems developed for this purpose in this paper.     

陀螺系统与反对称矩阵辛本征解的计算

用状态向量法,引出陀螺线性系统的广义本征问题,证明了本征向量之间的加权共轭辛正交关系,以及用本征向量对任意状态向量的展开定理。运用反对称矩阵胞块组成的LDL~T分解,将本征方程导向辛本征问题的标准型。这套方法适用于陀螺系统K阵不正定的情形。对于辛本征问题用SH变换将矩阵化为半边三对角线胞块阵或三对角线胞块阵,然后再求解其全部本征解。为陀螺系统的模态分析打下了基础。     

粘弹性厚壁筒问题的辛本征解方法

将哈密顿体系引进到粘弹性力学厚壁筒问题中,在辛体系下重新描述了基本问题,即建立了正则方程组。借助于积分变换,得到了拉伸、扭转和弯曲等问题的解以及有边界局部效应的解。将原问题归结为辛几何空间中的零本征值本征解和非零本征值本征解问题,从而建立了一种有效的分析问题方法和数值方法。为解决同类问题提供了一条可行的路径。     

Numerical solutions of linear quadratic control for time-varying systems via symplectic par conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented.It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation.The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.