详解Python 实现元胞自动机中的生命游戏(Game of life)


Posted in Python onJanuary 27, 2018

简介

细胞自动机(又称元胞自动机),名字虽然很深奥,但是它的行为却是非常美妙的。所有这些怎样实现的呢?我们可以把计算机中的宇宙想象成是一堆方格子构成的封闭空间,尺寸为N的空间就有NN个格子。而每一个格子都可以看成是一个生命体,每个生命都有生和死两种状态,如果该格子生就显示蓝色,死则显示白色。每一个格子旁边都有邻居格子存在,如果我们把33的9个格子构成的正方形看成一个基本单位的话,那么这个正方形中心的格子的邻居就是它旁边的8个格子。

每个格子的生死遵循下面的原则:

1. 如果一个细胞周围有3个细胞为生(一个细胞周围共有8个细胞),则该细胞为生(即该细胞若原先为死,则转为生,若原先为生,则保持不变) 。

2. 如果一个细胞周围有2个细胞为生,则该细胞的生死状态保持不变;

3. 在其它情况下,该细胞为死(即该细胞若原先为生,则转为死,若原先为死,则保持不变)

设定图像中每个像素的初始状态后依据上述的游戏规则演绎生命的变化,由于初始状态和迭代次数不同,将会得到令人叹服的优美图案。

代码

"""
元胞自动机 Python 实现
"""
import numpy as np
import matplotlib.pyplot as plt


class GameOfLife(object):

  def __init__(self, cells_shape):
    """
    Parameters
    ----------
    cells_shape : 一个元组,表示画布的大小。

    Examples
    --------
    建立一个高20,宽30的画布
    game = GameOfLife((20, 30))
    
    """

    # 矩阵的四周不参与运算
    self.cells = np.zeros(cells_shape)

    real_width = cells_shape[0] - 2
    real_height = cells_shape[1] - 2
    
    self.cells[1:-1, 1:-1] = np.random.randint(2, size=(real_width, real_height))
    self.timer = 0
    self.mask = np.ones(9)
    self.mask[4] = 0
  
  def update_state(self):
    """更新一次状态"""
    buf = np.zeros(self.cells.shape)
    cells = self.cells
    for i in range(1, cells.shape[0] - 1):
      for j in range(1, cells.shape[0] - 1):
        # 计算该细胞周围的存活细胞数
        neighbor = cells[i-1:i+2, j-1:j+2].reshape((-1, ))
        neighbor_num = np.convolve(self.mask, neighbor, 'valid')[0]
        if neighbor_num == 3:
          buf[i, j] = 1
        elif neighbor_num == 2:
          buf[i, j] = cells[i, j]
        else:
          buf[i, j] = 0
    self.cells = buf
    self.timer += 1
  
  def plot_state(self):
    """画出当前的状态"""
    plt.title('Iter :{}'.format(self.timer))
    plt.imshow(self.cells)
    plt.show()

  def update_and_plot(self, n_iter):
    """更新状态并画图
    Parameters
    ----------
    n_iter : 更新的轮数
    """
    plt.ion()
    for _ in range(n_iter):
      plt.title('Iter :{}'.format(self.timer))
      plt.imshow(self.cells)
      self.update_state()
      plt.pause(0.2)
    plt.ioff()
          

if __name__ == '__main__':
  game = GameOfLife(cells_shape=(60, 60))
  game.update_and_plot(200)

效果图

详解Python 实现元胞自动机中的生命游戏(Game of life)

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Python3中条件控制、循环与函数的简易教程
Nov 21 Python
浅析Python数据处理
May 02 Python
Python下使用Scrapy爬取网页内容的实例
May 21 Python
解决python3 Pycharm上连接数据库时报错的问题
Dec 03 Python
python3.x+pyqt5实现主窗口状态栏里(嵌入)显示进度条功能
Jul 04 Python
python 五子棋如何获得鼠标点击坐标
Nov 04 Python
python打印n位数“水仙花数”(实例代码)
Dec 25 Python
关于Pytorch的MNIST数据集的预处理详解
Jan 10 Python
python字符串,元组,列表,字典互转代码实例详解
Feb 14 Python
解决python便携版无法直接运行py文件的问题
Sep 01 Python
python_tkinter弹出对话框创建
Mar 20 Python
python Tkinter模块使用方法详解
Apr 07 Python
Python实现的栈(Stack)
Jan 26 #Python
使用python实现链表操作
Jan 26 #Python
Python中optparser库用法实例详解
Jan 26 #Python
python利用socketserver实现并发套接字功能
Jan 26 #Python
Django的HttpRequest和HttpResponse对象详解
Jan 26 #Python
Python编程实现的简单神经网络算法示例
Jan 26 #Python
Django使用httpresponse返回用户头像实例代码
Jan 26 #Python
You might like
php中3种方法删除字符串中间的空格
2014/03/10 PHP
PHP中使用array函数新建一个数组
2015/11/19 PHP
PHP切割汉字的常用方法实例总结
2019/04/27 PHP
用js实现的抽象CSS圆角效果!!
2007/05/03 Javascript
很多人都是用下面的js刷新站IP和PV
2008/09/05 Javascript
js 单引号 传递方法
2009/06/22 Javascript
iframe自适应宽度、高度 ie6 7 8,firefox 3.86下测试通过
2010/07/29 Javascript
js循环改变div颜色具体方法
2013/06/25 Javascript
js窗口关闭提示信息(兼容IE和firefox)
2015/10/23 Javascript
图解js图片轮播效果
2015/12/20 Javascript
详解webpack require.ensure与require AMD的区别
2017/12/13 Javascript
浅谈redux以及react-redux简单实现
2018/08/28 Javascript
详解React 的几种条件渲染以及选择
2018/10/23 Javascript
详解小程序毫秒级倒计时(适用于拼团秒杀功能)
2019/05/05 Javascript
详解微信小程序回到顶部的两种方式
2019/05/09 Javascript
Vue实现点击显示不同图片的效果
2019/08/10 Javascript
解决layer.prompt无效的问题
2019/09/24 Javascript
[00:03]DOTA2新版本PA至宝展示
2014/11/19 DOTA
Python实现信用卡系统(支持购物、转账、存取钱)
2016/06/24 Python
Python实现全角半角字符互转的方法
2016/11/28 Python
pip安装python库的方法总结
2019/08/02 Python
Python中typing模块与类型注解的使用方法
2019/08/05 Python
python redis 批量设置过期key过程解析
2019/11/26 Python
将keras的h5模型转换为tensorflow的pb模型操作
2020/05/25 Python
Pytorch实现将模型的所有参数的梯度清0
2020/06/24 Python
Python实现Canny及Hough算法代码实例解析
2020/08/06 Python
Python如何测试stdout输出
2020/08/10 Python
C++如何引用一个已经定义过的全局变量
2014/08/25 面试题
廉洁自律承诺书2015
2015/01/22 职场文书
感谢信范文大全
2015/01/23 职场文书
出生公证书
2015/01/23 职场文书
趣味运动会开幕词
2015/01/28 职场文书
店长岗位职责
2015/02/11 职场文书
司机岗位职责范本
2015/04/10 职场文书
公积金贷款承诺书
2015/04/30 职场文书
vue递归实现树形组件
2022/07/15 Vue.js