Python实现的matplotlib动画演示之细胞自动机
更新时间：2022年04月21日 09:40:07 作者：小小明-代码实体
这篇文章主要介绍了Python实现的matplotlib动画演示之细胞自动机,用python来模拟，首先尝试表示Beacon，本文通过实例代码给大家介绍的非常详细，需要的朋友可以参考下
维基百科上有个有意思的话题叫细胞自动机：https://en.wikipedia.org/wiki/Cellular_automaton
在20世纪70年代，一种名为生命游戏的二维细胞自动机变得广为人知，特别是在早期的计算机界。由约翰 · 康威发明，马丁 · 加德纳在《科学美国人》的一篇文章中推广，其规则如下：
- Any live cell with fewer than two live neighbours dies, as if caused by underpopulation.
- Any live cell with two or three live neighbours lives on to the next generation.
- Any live cell with more than three live neighbours dies, as if by overpopulation.
- Any dead cell with exactly three live neighbours becomes a live cell, as if by reproduction.
总结就是：任何活细胞在有两到三个活邻居时能活到下一代，否则死亡。任何有三个活邻居的死细胞会变成活细胞，表示繁殖。
在Conway’s Game of Life中，展示了几种初始状态：
下面我们用python来模拟，首先尝试表示Beacon：
import numpy as np import matplotlib.pyplot as plt universe = np.zeros((6, 6), "byte") # Beacon universe[1:3, 1:3] = 1 universe[3:5, 3:5] = 1 print(universe) im = plt.imshow(universe, cmap="binary")
[[0 0 0 0 0 0] [0 1 1 0 0 0] [0 1 1 0 0 0] [0 0 0 1 1 0] [0 0 0 1 1 0] [0 0 0 0 0 0]]
可以看到已经成功的打印出了Beacon的形状，下面我们继续编写细胞自动机的演化规则：
def cellular_auto(universe):
universe_new = universe.copy()
h, w = universe.shape
for y in range(h):
for x in range(w):
neighbor_num = universe[x-1:x+2, y-1:y+2].sum()-universe[x, y]
# 任何有三个活邻居的死细胞都变成了活细胞，繁殖一样。
if universe[x, y] == 0 and neighbor_num == 3:
universe_new[x, y] = 1
# 任何有两到三个活邻居的活细胞都能活到下一代，否则就会死亡。
if universe[x, y] == 1 and neighbor_num not in (2, 3):
universe_new[x, y] = 0
return universe_new
universe = cellular_auto(universe)
print(universe)
plt.axis("off")
im = plt.imshow(universe, cmap="binary")
[[0 0 0 0 0 0] [0 1 1 0 0 0] [0 1 0 0 0 0] [0 0 0 0 1 0] [0 0 0 1 1 0] [0 0 0 0 0 0]]
ArtistAnimation动画
基于此我们可以制作matplotlib的动画，下面直接将Blinker、Toad、Beacon都放上去：
from matplotlib import animation
import numpy as np
import matplotlib.pyplot as plt
%matplotlib notebook
def cellular_auto(universe):
universe_new = universe.copy()
h, w = universe.shape
for y in range(h):
for x in range(w):
neighbor_num = universe[x-1:x+2, y-1:y+2].sum()-universe[x, y]
# 任何有三个活邻居的死细胞都变成了活细胞，繁殖一样。
if universe[x, y] == 0 and neighbor_num == 3:
universe_new[x, y] = 1
# 任何有两到三个活邻居的活细胞都能活到下一代，否则就会死亡。
if universe[x, y] == 1 and neighbor_num not in (2, 3):
universe_new[x, y] = 0
return universe_new
universe = np.zeros((12, 12), "byte")
# Blinker
universe[2, 1:4] = 1
# Beacon
universe[4:6, 5:7] = 1
universe[6:8, 7:9] = 1
# Toad
universe[8, 2:5] = 1
universe[9, 1:4] = 1
fig = plt.figure()
plt.axis("off")
im = plt.imshow(universe, cmap="binary")
frame = []
for _ in range(2):
frame.append((plt.imshow(universe, cmap="binary"),))
universe = cellular_auto(universe)
animation.ArtistAnimation(fig, frame, interval=500, blit=True)
然后我们画一下Pulsar：
# Pulsar
universe = np.zeros((17, 17), "byte")
universe[[2, 7, 9, 14], 4:7] = 1
universe[[2, 7, 9, 14], 10:13] = 1
universe[4:7, [2, 7, 9, 14]] = 1
universe[10:13, [2, 7, 9, 14]] = 1
fig = plt.figure()
plt.axis("off")
im = plt.imshow(universe, cmap="binary")
frame = []
for _ in range(3):
frame.append((plt.imshow(universe, cmap="binary"),))
universe = cellular_auto(universe)
animation.ArtistAnimation(fig, frame, interval=500, blit=True)
FuncAnimation动画
另一种创建matplotlib动画的方法是使用FuncAnimation，完整代码：
from matplotlib import animation
import numpy as np
import matplotlib.pyplot as plt
from IPython.display import HTML
# %matplotlib notebook
def cellular_auto(universe):
universe_new = universe.copy()
h, w = universe.shape
for y in range(h):
for x in range(w):
neighbor_num = universe[x-1:x+2, y-1:y+2].sum()-universe[x, y]
# 任何有三个活邻居的死细胞都变成了活细胞，繁殖一样。
if universe[x, y] == 0 and neighbor_num == 3:
universe_new[x, y] = 1
# 任何有两到三个活邻居的活细胞都能活到下一代，否则就会死亡。
if universe[x, y] == 1 and neighbor_num not in (2, 3):
universe_new[x, y] = 0
return universe_new
def update(i=0):
global universe
im.set_data(universe)
universe = cellular_auto(universe)
return im,
# Pulsar
universe = np.zeros((17, 17), "byte")
universe[[2, 7, 9, 14], 4:7] = 1
universe[[2, 7, 9, 14], 10:13] = 1
universe[4:7, [2, 7, 9, 14]] = 1
universe[10:13, [2, 7, 9, 14]] = 1
fig = plt.figure()
plt.axis("off")
im = plt.imshow(universe, cmap="binary")
plt.show()
anim = animation.FuncAnimation(
fig, update, frames=3, interval=500, blit=True)
HTML(anim.to_jshtml())
这种动画生成速度较慢，好处是可以导出html文件：
with open("out.html", "w") as f:
f.write(anim.to_jshtml())
还可以保存MP4视频：
anim.save("out.mp4")
或gif动画：
anim.save("out.gif")
注意：保存MP4视频或GIF动画，需要事先将ffmpeg配置到环境变量中
ffmpeg下载地址：
链接: https://pan.baidu.com/s/1aioB_BwpKb6LxJs26HbbiQ?pwd=ciui
提取码: ciui
随机生命游戏
接下来，我们创建一个50*50的二维生命棋盘，并选取其中1500个位置作为初始活细胞点，我们看看最终生成的动画如何。
完整代码如下：
from matplotlib import animation
import numpy as np
import matplotlib.pyplot as plt
%matplotlib notebook
def cellular_auto(universe):
universe_new = universe.copy()
h, w = universe.shape
for y in range(1, h-1):
for x in range(1, w-1):
neighbor_num = universe[x-1:x+2, y-1:y+2].sum()-universe[x, y]
# 任何有三个活邻居的死细胞都变成了活细胞，繁殖一样。
if universe[x, y] == 0 and neighbor_num == 3:
universe_new[x, y] = 1
# 任何有两到三个活邻居的活细胞都能活到下一代，否则就会死亡。
if universe[x, y] == 1 and neighbor_num not in (2, 3):
universe_new[x, y] = 0
# 边缘置零
universe[[0, -1]] = 0
universe[:, [0, -1]] = 0
return universe_new
boardsize, pad = 50, 2
universe = np.zeros((boardsize+pad, boardsize+pad), "byte")
# 随机选取1500个点作为初始活细胞
for i in range(1500):
x, y = np.random.randint(1, boardsize+1, 2)
universe[y, x] = 1
fig = plt.figure()
plt.axis("off")
im = plt.imshow(universe, cmap="binary")
frame = []
for _ in range(200):
frame.append((plt.imshow(universe, cmap="binary"),))
universe = cellular_auto(universe)
animation.ArtistAnimation(fig, frame, interval=50, blit=True)
到此这篇关于Python实现的matplotlib动画演示之细胞自动机的文章就介绍到这了,更多相关python matplotlib动画内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
相关文章
simple-pytest框架主要参考了httprunner的yaml数据驱动部分设计思路,是基于Pytest + Pytest-html+ Log + Yaml + Mysql 实现的简易版接口自动化框架，这篇文章主要介绍了simple-pytest 框架使用指南,需要的朋友可以参考下2024-02-02
这篇文章主要介绍了Python中的CURL PycURL使用例子,需要的朋友可以参考下2014-06-06
这篇文章主要为大家详细介绍了python运用sklearn实现KNN分类算法，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2019-10-10
这篇文章主要介绍了python异常和文件处理机制,详细分析了Python异常处理的常用语句、使用方法及相关注意事项,需要的朋友可以参考下2016-07-07
大家好，本篇文章主要讲的是Python中引用传参四种方式介绍，感兴趣的同学赶快来看一看吧，对你有帮助的话记得收藏一下2021-12-12
今天带大家学习怎么用python实现一个无界面的小型图书管理系统,文中有非常详细的图文解说及代码示例,对正在学习python的小伙伴们有很好地帮助,需要的朋友可以参考下2021-05-05
这篇文章主要介绍了python 如何实现跳过异常继续执行，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教2022-07-07
今天小编就为大家分享一篇对dataframe进行列相加,行相加的实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2018-06-06
本篇文章将会使大家了解Python 语句、表达式以及它们之间的区别。还包含几个示例来更清楚地解释这个概念。接下来，我们将解释如何在 Python 编程中使用多行语句和缩进,需要的朋友可以参考一下2021-09-09
常见的数据增强操作有：按比例放大或缩小图片、旋转、平移、水平翻转、改变图像通道等。本文将通过Python OpenCV实现这些操作，需要的可以参考一下2022-02-02
最新评论