布朗运动的理论分析与蒙特卡罗模拟.docVIP

布朗运动的理论分析与蒙特卡罗模拟 徐建强 卢氏县第一高级中学 472200 摘要：本文简要给出布朗运动的郎之万理论，利用蒙特卡罗方法中的随机游走理论，通过Matlab编程.直观、形象地模拟了三维空间布郎运动，同时计算了布朗粒子位移平方的平均值与时间的关系，其结果与理论分析一致。 关键词：布朗运动；蒙特卡罗方法；Matlab 中图分类号：O552.1 1 引言 1827 年,苏格兰植物学家布朗在显微镜下观察到花粉在水中作无规则运动,并发现其它悬浮微粒也不停息地作无规则运动。微粒在流体介质中这种随机的连续运动,称为布朗运动;该微粒可以叫做布朗粒子。布朗运动发现之后50 年内,其实质一直不为人知。1877 年,德索耳提出,微粒受到周围媒质分子不平衡的碰撞是产生布朗运动的原因。1905 年爱因斯坦发表了关于布朗运动的论文,并被普遍接受。随后斯莫陆绰斯基(1906 年) 、朗之万(1908年) 都发表了各自关于布朗运动的理论,证明布朗粒子位移平方的平均值正比于时间。 [1]现在对布朗运动二维可视化研究已很多，本文着重从三维立体空间实现布朗运动的可视化。 2 理论分析 由郎之万方程： （1） F(t)为涨落力，为斯托克斯公式，a为颗粒（看作小球）半径，η为粘滞系数。对（1）式两边乘x得 （2） 将（2）式对大量颗粒求平均，即把各布郎粒子的运动方程相加然后用粒子数去除。因F(t)的平均值为零，并对（2）式第二项应用能量均分定理，经整理可得 （3） 式（3）是的二阶常系数线性非齐次微分方程，通解为 （4） 经过近似处理，并假设所有粒子在t=0时都处在x=0处，则C2=0.因此得 （5） （5）式表明位移平方的平均值与时间成正比，正是理论结论。 3 蒙特卡罗模拟3.1蒙特卡罗方法 蒙特卡罗（Monte-Carlo）方法是利用连续产生的无规数求解的一种计算方法。蒙特卡罗方法在物理工程中有着广泛的应用，例如对于分子运动学，粒子输运现象，布朗运动，放射性衰变等现象，由于问题本身有一定的无规性，同时又具有一定的统计规律性，用这个方法很合适。, 同时发展成为随机过程研究的一个重要分支.蒙特卡罗中随机行走就是采用伪随机数发生器所产生的伪随机数序列对布朗运动这一类物理过程作计算机仿真的方法. 为简化起见，我们研究具有代表性的二维随机行走。设一醉汉从原点出发行走N步, 每一步所行的方向与其前一步的方向完全无关具有随机性,每步的步长则在指定的区间里涨落。这样，各步在x方向和y方向的位移分别为 ， ，……， （6） 在行走N步之后, 行走者所在位置与原点的距离R就有如下的关系: （7） 只要行走是真正随机的, 则行走者向前后左右行走的概率是一样的, 走了相当多的步数之后, 上式中的交叉项就相互抵消, 于是得到: （8） 于是有: （9） 其中r为步长。从此式可以看出当步长一定的时候，R2与步数n成正比。对于该式成立的条件是n足够大。[3] 3.2计算机模拟 为了直观在此模拟三维空间中布朗运动，图1为空间中任意一点，坐标图示， 图1空间点坐标表示图示 用极坐标表示则 ，，， ，，， 为了让这些动起来，我们需要给出运动的方向和步长。这里我们选择步长ss=，theta1, a=theta2并在极坐标下通过随机数产生函数给出分子随机的运动方向，之后再回到笛卡儿坐标系中换算出每个分子位置的x、y分量。这里我们将通过计算机模拟00个分子在维内的随机运动情况，在Matlab中定义一个00×3的矩阵来存放这些的位置。其中每一行代表一个分子，第一列存放位置的x分量，类似的第二列存放其y分量相应的Matlab程序段如下：rand(state,sum(100*clock)); %产生随机数 ns=200; %设定时间 np=200; %分