摘要：提出一种朗之万动力学方法获取处于热平衡态耦合系统内部振子坐标，数值模拟了单端固定简谐振子链的时间演化行为，并将其平衡性质与解析解进行了比较．结果表明了朗之万动力学方法的有效性．推广应用于非简谐四次方型耦合系统，模拟得到振子的四次方均坐标，与理论值验证；以模拟结果作为样本点计算哈密顿量，其能量分布与Boltzmann分布相符．

关键词：耦合振子链 langevin方程 稳态分布 

单位：北京师范大学物理学系 北京100875