哈密尔顿系统的有限元法

利用常微分方程的连续有限元法,结合函数的M-型展开,对非线性哈密尔顿系统证明了连续一、二次有限元分在3阶量、5阶量意义下近似保辛,且保持能量守恒.在数值实验中结合庞加莱截面,哈密尔顿混沌数值试验结果与理论相吻合.     

A new approximate analytical solution of the ideal potential flow around a circular cylinder between two parallel flat plates

In ref. [1], Lin obtained an approximate analytical solution of the ideal potential flow around a circular cylinder between two parallel flat flates.In this paper, the author shows that one may obtain the result coinciding with that obtained in ref. [1] by making use of the Shvez's method^[2]. Morever, we can obtain a more accurate result than that obtained in ref. [1], if we make use of the improved Shvez's method^[2]. Some calculating examples are presented.     

SINGULAR PERTURBATION OF NONLINEAR BOUNDARY VALUE PROBLEMS

In this paper we consider the boundary value problem where ε.μ are two positive parameters. Under f_y≤-k<0 and other suitable restrictions, there exists a solution and it satisfied where y_(0,0)(x) is solution of reduced problem while y_i-j,j(x)(j=0,1,...,i;i=1,2,...,m) can be obtained successively from certain linear equations.     

两类辛算法的稳定性比较

主要对一阶隐式Euler法、二阶隐式Euler中点公式分别嵌入到一阶显辛Euler法和二阶leapfrog显辛积分器所得算法M1和M2进行了稳定性比较。对于坐标与速度耦合的线性哈密顿系统,理论分析寻找稳定区比较困难,依靠数值方法定性地比较了他们的稳定性。结果表明M1算法一般优于M2算法。     

双色激光场中一维氢分子的经典动力学性质研究

应用经典理论方法,采用辛算法,数值求解双色激光场作用下的一维共线氢分子的哈密顿正则方程,得到氢分子在激光场作用下的经典轨迹,并计算H₂⁺,H₂²⁺,H,H⁺等产物的产生概率随时间的演化,分析了双色场的光强、倍频、相位的变化对氢分子动力学性质的影响并给出相关的物理解释.     

一类带不平坦界面声波导中的共轭特征算子构造

无界且带有不平坦界面的声波导经坐标变换和完美匹配层截断,波的传播计算问题近似转化为在有界且带有平坦界面的声波导中数值求解改进Helmholtz方程,由于用步进算法解此改进复Helmholtz方程所涉及到的算子的特征函数一般不具有正交性,故在数值步进求解复方程时存在着局部基下坐标计算的困难.本文一方面推导出此复偏微分方程的共轭特征函数所满足的方程,并论证方程的特征函数与共轭特征函数正交之性质;另一方面,给出局部基下坐标计算的简便公式,它可使步进计算保持高效率.数值模拟结果表明所提方法切实可行、有效.     

李级数算法和显式辛算法的相位分析

以线性可分Hamilton动力学系统为例,研究了李级数算法和显式辛算法的相位精度,研究了李级数算法的保辛精度及其保辛精度的提高方法;指出了显式辛算法相位精度与算法阶次的不协调性,印辛算法的阶次高并不意味着其相位精度也高,李级数算法不存在这种问题,指出了一个算法的相位可能超前也可能滞后.分析结果表明三阶显式辛算法具有比较高的相位精度.     

一种光学层析图像的多准则重建方法

光学层析成像是一个病态重建问题,为克服重建过程的病态性,提出将多准则优化理论引入到图像重建中。利用了三个用于光学层析图像重建的准则:平方误差函数、图像熵函数和局部平滑函数。采用向量优化方法将多准则优化问题转化为单准则优化问题求解。为了确定各个目标函数间的权重系数,提出一种动态权重系数求解方法。重建过程目标函数关于光学参量的梯度计算是关键,因此提出一种基于梯度树的计算方法。实验过程中对多准则重建结果和基于平方误差函数的单准则重建结果做了比较,证明该方法能够克服传统的偏重单一目标的单准则重建的不足,有效地重建光学层析图像,提高图像重建质量。     

Perfect plane-wave source for a high-order symplectic finite-difference time-domain scheme

The method of splitting a plane-wave finite-difference time-domain (SP-FDTD) algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field (TF/SF) formulation of high-order symplectic finite-difference time-domain (SFDTD) scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain (FDTD), the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.     

Completeness of eigenfunction systems for the product of two symmetric operator matrices and its application in elasticity

The completeness theorem of the eigenfunction systems for the product of two 2×2 symmetric operator matrices is proved. The result is applied to 4×4 infinite-dimensional Hamiltonian operators. A modified method of separation of variables is proposed for a separable Hamiltonian system. As an application of the theorem, the general solutions for the plate bending equation and the free vibration of rectangular thin plates are obtained. Finally, a numerical test is analysed to show the correctness of the results.