计算物理 homework11
李明达 PB18020616
1
*
摘要
这是计算物理第 11 次作业,作业题目是数值研究 d(d=1, 2, 3) 维空间中随机行走返回原点的几率
P
d
,讨论它随步数 N 的变化关系 P
d
(N),能否定义相关的指数值?
关键词
随机行走自回归
1
中国科学技术大学物理学院 * 作者: dslmd@mail.ustc.edu.cn
1. 模型和算法
1.1 模型
在本实验中,我们考虑 d 维空间中的粒子在 N
次随机游走恰好能回到原点的几率。由讲义的知识,
一维随机行走之后,其分布趋向于高斯分布,我们
有
x
2
(t) = 2Dt
f(x) ∼ N(0, σ
2
)
可以看到,随着时间(步数)的增加,回到原点的几
率越来越小。我们可以将其扩展到二维网格、三维
网格上,在不考虑自规避的前提下,可以看出,显然
几率要更小一些。
理论上的分析告诉我们,P
d
(N) 应该可以定义
一个指数值,定义如下:
P
d
(N) = A × N
−ν
这个指数值可以通过如下变形看的更清楚:
lnP
d
(N) = lnA − νlnN
对一维 d=1,理论值为 0.5;对二维 d=2,理论
值为 1.0;对三维 d=3,理论值为 1.5。
下面是我的算法:
1.2 一维随机行走
单个粒子的位置用数组 x[i] 来表述,我们采用
一下 MC 模拟方法,先产生 [0, 1] 随机数:
如果随机数 ∈ [0, 0.5],粒子往 x 正方向走,x+1
如果随机数 ∈ [0.5, 1],粒子往 x 正方向走,x-1
1.3 二维随机行走
单个粒子的位置用数组 x[i], y[i] 来表述,我们
采用一下 MC 模拟方法,先产生 [0, 1] 随机数:
如果随机数 ∈ [0, 0.25],粒子往 x 正方向走,x+1
如果随机数 ∈ [0.25, 0.5],粒子往 x 正方向走,
x-1
如果随机数 ∈ [0.5, 0.75],粒子往 x 正方向走,
y+1
如果随机数 ∈ [0.75, 1],粒子往 x 正方向走,y-1
1.4 三维随机行走
单个粒子的位置用数组 x[i], y[i], z[i] 来表述,我
们采用一下 MC 模拟方法,先产生 [0, 1] 随机数:
如果随机数 ∈ [0, 1/6],粒子往 x 正方向走,x+1
如果随机数 ∈ [1/6, 1/3],粒子往 x 正方向走,
x-1
如果随机数 ∈ [1/3, 1/2],粒子往 x 正方向走,
y+1
如果随机数 ∈ [1/2, 2/3],粒子往 x 正方向走,
y-1
如果随机数 ∈ [2/3, 5/6],粒子往 x 正方向走,
z+1
如果随机数 ∈ [5/6, 1],粒子往 x 正方向走,z-1
1.5 程序
在程序中,我通过 d
i
(intseed, intN) 函数来模
拟 i 维的随机行走自回归概率,函数的参数 N 是指
N 次随机行走,每一个 N 次随机行走进行 n 次(n
是函数内的的变量,可以进行调节)的重复模拟来
计算随机行走自回归概率。
2. 实验结果讨论
评论0