问将一个群体随机分成6个

EN

我在研究三角网格上元胞自动机粒子扩散的算法。这意味着每个细胞有6个邻居。

每个细胞都有一定数量的粒子。

在每一次迭代中，每个细胞都将其粒子分散到相邻的细胞中。

我很难有效地做到这一点，因为有数十万个细胞(有时是几百万个)，每个细胞中都有大量的粒子n (n >> 100)。

，我在寻找一种算法，把一个数字随机分成6部分，

一种有效但效率低下的办法：

从区间(0,6)上的均匀分布中产生随机数，就像单元中有粒子一样。

• 如果数字在(0，1 )中:将粒子扩展到相邻的1。
• 如果数字在(1，2 )中:将粒子分散到相邻的2。
• 如果数字在(2，3 )中:将粒子分散到相邻的3。
• 等等。

这适用于“小”粒子数(n < 50)，但得到非常计算密集。

我的理论方法：

调用要分配的粒子数n。

从均值为0和方差为1的正态分布(高斯)中生成5个随机数。将这些数字称为r0、r1、r2、r3、r4

r0 = n/2 + r0*(n/4) // this transforms r0 to a random number drawn from a normal distribution with mean n/2 and variance n/2

r0有效地将n个粒子的粒子数分成两组，每个组分配给三个邻域。大小r0，n- r0大小的

r1 = r0/3 + r0*(r0/9) // this transforms r1 to a random number drawn from a normal distribution with mean r0/3 and variance r0/3

r1有效地将r0粒子的粒子数分成两组，一组分配给单个邻居，另一组分配给两个邻居。前者为size r1，后者为size r0 - r1

r2 = (r0 - r1)/2 + r2*((r0 - r1)/4) // this transforms r2 to a random number drawn from a normal distribution with mean (r0 - r1)/2 and variance (r0 - r1)/2

r2有效地将(r0 - r1)粒子的粒子数分成两组，每一组分配给一个单独的邻居。

数字r0、r1和r2现在应该从n个粒子的总体中分离出3个随机部分，每个部分按正态分布。

以同样的方式继续，将人口(n - r0)分成三部分。

这种方法对我来说似乎是有意义的，但我相信我的差异可能离这里很远，所以我得到了奇怪的结果，有太多的粒子被“分离”给一组邻居，而没有一个被留给其他的邻居。这就引入了奇怪的各向异性效应。

背景:许多均匀分布的组合被高斯分布很好地近似。该算法是试图修改Bastien在“元胞自动机物理系统建模”第5.7章(第213页)中描述的一种算法。

如果能在我的方法中看到一个错误或一个不同的、同样有效的错误，将非常感谢。

我还没有指定im编码语言，因为我只是在寻找一般的算法。我正在使用java (Process3.5)，但是如果您可以提供任何语言，这对我来说是可以的。