计算物理 ›› 2005, Vol. 22 ›› Issue (1): 1-6.

• 研究论文 •    下一篇

一般阻尼随机微分方程的拟局部振荡算法

宋艳丽1, 包景东1,2   

  1. 1. 北京师范大学物理系, 北京 100875;
    2. 北京师范大学低能核物理研究所, 北京 100875
  • 收稿日期:2004-07-12 修回日期:2004-09-25 出版日期:2005-01-25 发布日期:2005-01-25
  • 作者简介:宋艳丽(1977-),女,河北邯郸,博士生,从事涨落与耗散动力学研究.
  • 基金资助:
    国家自然科学基金重点(10235020);教育部跨世纪优秀人才基金资助项目

A Quasi-local Vibration Algorithm for Damped Stochastic Differential Equations

SONG Yan-li1, BAO Jing-dong1,2   

  1. 1. Department of Physics, Beijing Normal University, Beijing 100875, China;
    2. Institute of Low-Energy Nuclear Physics, Beijing Normal University, Beijing 100875, China
  • Received:2004-07-12 Revised:2004-09-25 Online:2005-01-25 Published:2005-01-25

摘要: 提出了一种模拟随机微分方程的拟局部振荡算法,即利用算符劈裂方法和势函数的泰勒展开,对噪声作用下耗散粒子的时间演化算符进行分解,得到了对应涨落行为的扩散算符和对应确定轨迹的漂移算符.其中局部简谐势场的涨落过程可获得解析解,而剩余的确定项则利用简单的Euler算法积分.应用到几个算例并与常用的两种算法相比较,结果表明:本算法随时间步长最稳定,可使用较大的时间步长.

关键词: 算符劈裂, 朗之万方程, 局域振荡, 稳定性

Abstract: A quasi-local vibration algorithm is proposed to simulate the damped stochastic differential equations.Using the potential Taylor expansion and the operator splitting technique,the time evolution operator of a dissipative particle is factorized into two parts: an operator which describes the deterministic path of the system and a diffusion operator which incorporates fluctuation away from path.The differential equation corresponding to the fluctuation process is solved analytically, and the deterministic trajectory equation is calculated with the Euler algorithm.Compared with other algorithms, it is shown that the calculated results by the present algorithm are more stable and converge to the correct value when the time step is large.

Key words: operator splitting, Langevin equation, quasi-local vibration, stableness

中图分类号: