博客
关于我
04.沙堆模型
阅读量:461 次
发布时间:2019-03-06

本文共 1293 字,大约阅读时间需要 4 分钟。

"""沙堆模型:1. 在一个二维正方格子上,每个格点上都可以放上一定的沙子;2. 当某个格子 (i,j) 里的沙子数目大于 4 时,就会发生崩塌(或者也可以设定为按照一定的概率崩塌);3. 格子 (i,j) 里面的 4 颗沙子则转移到 (i,j-1),(i,j+1),(i-1,j),(i+1,j) 四个格子中;4. 随着这四个格子中沙粒数目的变化,上述的演化规则可以进一步逐层扩展到其它相邻的格点,最终蔓延整个沙堆。"""import numpy as npimport matplotlib.pyplot as plt  def iterate(grid):    changed = False    for ii, arr in enumerate(grid):  # ii:矩阵行号, jj:矩阵列号, 便于定位        for jj, val in enumerate(arr):  # 沙子数目大于4发生崩塌            if val >= 4:                grid[ii, jj] -= 4                # 满足条件的前提下移动沙子                if ii >= 1:                    grid[ii - 1, jj] += 1                if ii <= len(arr) - 1:                    grid[ii + 1, jj] += 1                if jj >= 1:                    grid[ii, jj - 1] += 1                if jj <= len(arr) - 1:                    grid[ii, jj + 1] += 1                # 矩阵是否还在变化的标志                changed = True    return grid, changed def simulate(grid):    # 循环递归,直到矩阵不再变化为止    while True:        grid, changed = iterate(grid)        if not changed:            return grid  if __name__ == "__main__":    start_grid = np.zeros((10, 10))    start_grid[3:4, 6:7] = 64    final_grid = simulate(start_grid.copy())    plt.figure()    plt.gray()    plt.imshow(start_grid)    plt.show()    plt.figure()    plt.gray()    plt.imshow(final_grid)    plt.show()

 

转载地址:http://rekbz.baihongyu.com/

你可能感兴趣的文章
MySQL 常见的开放性问题
查看>>
Mysql 常见错误
查看>>
mysql 常见问题
查看>>
MYSQL 幻读(Phantom Problem)不可重复读
查看>>
mysql 往字段后面加字符串
查看>>
mysql 快照读 幻读_innodb当前读 与 快照读 and rr级别是否真正避免了幻读
查看>>
MySQL 快速创建千万级测试数据
查看>>
mysql 快速自增假数据, 新增假数据,mysql自增假数据
查看>>
MySql 手动执行主从备份
查看>>
Mysql 批量修改四种方式效率对比(一)
查看>>
Mysql 报错 Field 'id' doesn't have a default value
查看>>
MySQL 报错:Duplicate entry 'xxx' for key 'UNIQ_XXXX'
查看>>
Mysql 拼接多个字段作为查询条件查询方法
查看>>
mysql 排序id_mysql如何按特定id排序
查看>>
Mysql 提示:Communication link failure
查看>>
mysql 插入是否成功_PDO mysql:如何知道插入是否成功
查看>>
Mysql 数据库InnoDB存储引擎中主要组件的刷新清理条件:脏页、RedoLog重做日志、Insert Buffer或ChangeBuffer、Undo Log
查看>>
mysql 数据库中 count(*),count(1),count(列名)区别和效率问题
查看>>
mysql 数据库备份及ibdata1的瘦身
查看>>
MySQL 数据库备份种类以及常用备份工具汇总
查看>>