蒙特卡罗法对布朗运动的模拟.pdfVIP

经济论丛 蒙特卡罗法对布朗运动的模拟 ◇ 谭 　波 　霍海波 　郑州航空工业管理学院数理系 　　【摘 　要 】布朗运动是分子物理的重要内容 , 它对人们认识物质的结构, 热运动的本质及其规律都有重要意义 , 教材中一般只做定性的讨论 。 本 文通过采用随机游走的模型 , 对布朗运动进行了二维空间的模拟 , 所得结论与爱因斯坦给出的理论公式趋向一致 , 并与自由粒子的匀速 直线运送进 行了比较 。 可以使同学们比较形象直观的理解分子的无规则运动 。 也为同学们初步体会学习用计算机来模拟物理问题 。 【关键词 】蒙特卡洛法 　布朗运动 　模拟 位移为 我们取 布朗粒子运动的方向为 　　一 、引言 ■r=k, k=0.6, θ=2πξ2, ξ∈ (0, 近年来 , 由于计算方法的深入发展和计算机的 高速发展及普 及, 随 1), 则可以推出下一步粒子所在位置为 着物理学基础理论的进一步突 破 , 物理 学家们 逐步可 以应用 一些更 严 xi+1 =xi+■rcosθ 格和更全面的复杂模型 , 来定 量研究 实际的 复杂体 系的 物理性 质 。 基 yi+1 =yi+■rcosθ 于物理学基本原理的数值计算和模拟已经成为将 理论物理和实验 物理 　　设 t=0 时粒子的位置取为 x0 =0, y0 =0。 每隔 1000步 取一次矢径 紧密联系在一起的一座重要的 桥梁 , 它 不仅能 够弥补 简单的 解析理 论 平方平均值 , 总共经历了 100000步 , 则相对出发点的矢径平方的平均值 模型难以完全描述复杂物理现 象的不 足 , 而且 可以克 服实验 物理中 遇 可以表示为 1000 1000 到的许多困难 , 例如直接模拟实验上不能实现或技术条件要 求很高 、实 2 1 2 1 2 2 r1 = ∑ ri= ∑ (xi +yi) 1000i=1 1000i=1 验代价昂贵的物理系统等