作业十二Cellular Automata

作业

FinalExam Sel-organise: Sandpile model
2013301020084 许晗

说明

主程序:sandpile.py;

任意曲线方程拟合数据:data_analyse.py;

有一次N=10000的计算结果存入了:data.txt中;

sandavalanch01.png: collapes caused by one drop

100_avalanch_10000.png: the amount of collapes happend in 10,000 times randomadd and evolve distributed in the lattice

1ei_iej_distribution.png: the fitcurve of the data from a N=10,000 calculation, selecting the $N \in [10^i,10^j)$ to fit.

不好储存的数据

N:总坍塌次数, D:每次投点所引起的平均坍塌次数的倒数, iN:在$100\times 100$的点阵里有$i$个$10\times 10$的方形边界, L:总点阵大小, l:内部方形边界大小$l\times l$

总点阵大小

correspond to "size_distribution"

特殊边界条件

correspond to "avaboundary_10_400_i", which represent boundary-iN said as above and are averaged for 400 times random add and evolve


and

correspond to "avaboundary_i_400", which represent boundary-il*il said as above and are averaged for 400 times random add and evolve