兰顿蚂蚁python实现
由来:
兰顿蚂蚁(英语:Langton’s ant)是细胞自动机的例子。它由克里斯托夫·兰顿在1986年提出,它由黑白格子和一只“蚂蚁”构成,是一个二维图灵机。兰顿蚂蚁拥有非常简单的逻辑和复杂的表现。在2000年兰顿蚂蚁的图灵完备性被证明。兰顿蚂蚁的想法后来被推广,比如使用多种颜色。
释义:
在平面上的正方形格被填上黑色或白色。在其中一格正方形有一只“蚂蚁”。它的头部朝向上下左右其中一方。
若蚂蚁在黑格,右转90度,将该格改为白格,向前移一步;
若蚂蚁在白格,左转90度,将该格改为黑格,向前移一步。
很多时,蚂蚁刚刚开始时留下的路线都会有接近对称、像是会重复。但不论起始状态如何,蚂蚁的路线必然是无限长的。
简单来说,你可以想象一个无限延伸的二维平面.平面被无数个等大黑色或白色的小方格分割,在其中的一个方格中存在一只蚂蚁.这只蚂蚁的初始方向任意,然后,这只蚂蚁按着如下的方式爬行:
若蚂蚁在白格,右转90度,将该格改为黑格,向前移一步;
若蚂蚁在黑格,左转90度,将该格改为白格,向前移一步。
按照这种逻辑进行下去,蚂蚁的初始路线会出现许多对称或重复的形状,然后是混沌的假随机,最后进行到一万步时会出现一条104周期的’高速公路’(在已知的实验中,都出现了这种现象,但尚无法证明这是必然现象).
知道了其中的道理,我们就可以通过编程实现了.
兰顿蚂蚁的实现可以分为两个方面.一是绘图,这方面比较容易,只需要让turtle向前向左向右就行了;二是黑白的判断与变化,虽然无法直观的表现,但我们可以使用numpy提供的矩阵来模拟,用0代表黑色,用1代表白色,接着用0或1填满这个二维矩阵并实现相应的变化.ok,思路有了,我们开始实现.
代码如下:
import numpy
i = int(input('输入屏幕的长\n'))
j = int(input('输入屏幕的宽\n'))
arr = numpy.ones((i,j))
建立一个用户指定大小的白矩阵(黑矩阵也行)
代码如下:
def draw(i, j):
turtle.setup(0.5, 0.5)
turtle.speed(0)
i = i//2
j = j//2
sign = ['d', 'w', 'a', 's'] # 用于判断方向的列表
a = 0
count = 0 # 记录步数
ward = sign[a] # 设定初始方向为向右
while True:
if arr[i][j]: # 如果为白块
arr[i][j] = 0 # 改为黑块
turtle.color('white', 'black')
turtle.begin_fill()
fill_squ(ward) # 调用自定义的填充函数
turtle.end_fill()
turtle.left(90) # 向左转90°并前进
turtle.forward(10)
if a == 3: # 周期性边界条件判断
a = 0
else:
a += 1
ward = sign[a]
if ward == 'd': # 对矩阵上的位置移动
i += 1
if ward == 'w':
j -= 1
if ward == 'a':
i -= 1
if ward == 's':
j += 1
count += 1
print(count) # 打印步数
else: # 如果为黑块(与上注释类似)
arr[i][j] = 1
turtle.color('white', 'white')
turtle.begin_fill()
fill_squ(ward)
turtle.end_fill()
turtle.right(90)
turtle.forward(10)
if a == 0:
a = 3
else:
a -= 1
ward = sign[a]
if ward == 'd':
i += 1
if ward == 'w':
j -= 1
if ward == 'a':
i -= 1
if ward == 's':
j += 1
count += 1
print(count)
else:
turtle.mainloop()
在这里我们遇到了一个问题,那就是在蚂蚁移动的过程当中,无法判断它的方向,就不知道它在矩阵上的移动方式.我的解决方案是,用一个元素为四个方向的列表追踪蚂蚁的方向,如果是白块向右,就是顺时针,列表元素递增,如果是黑块向左,就是逆时针,列表元素递减;并设置一个边界条件,让它循环.
def fill_squ(ward):
if ward == 'd':
for i in range(4):
turtle.forward(10)
turtle.left(90)
if ward == 'w':
for i in range(4):
turtle.forward(10)
turtle.right(90)
if ward == 'a':
turtle.right(180)
for i in range(4):
turtle.forward(10)
turtle.left(90)
turtle.left(180)
if ward == 's':
turtle.left(90)
for i in range(4):
turtle.forward(10)
turtle.left(90)
turtle.right(90)
方块可不会填充它自己,我们手动给他加上,由于每次蚂蚁的方向不同,故填充方式也需要经过判断后执行.
import numpy
import random
import turtle
def set_screen(i, j):
arr = numpy.empty((i, j), dtype=int)
for a in range(i):
for b in range(j):
arr[a][b] = random.randint(0, 1)
print(arr)
return arr
def fill_squ(ward):
if ward == 'd':
for i in range(4):
turtle.forward(10)
turtle.left(90)
if ward == 'w':
for i in range(4):
turtle.forward(10)
turtle.right(90)
if ward == 'a':
turtle.right(180)
for i in range(4):
turtle.forward(10)
turtle.left(90)
turtle.left(180)
if ward == 's':
turtle.left(90)
for i in range(4):
turtle.forward(10)
turtle.left(90)
turtle.right(90)
def draw(i, j):
turtle.setup(0.5, 0.5)
turtle.speed(0)
i = i//2
j = j//2
sign = ['d', 'w', 'a', 's']
a = 0
count = 0
ward = sign[a]
while True:
if arr[i][j]:
arr[i][j] = 0
turtle.color('white', 'black')
turtle.begin_fill()
fill_squ(ward)
turtle.end_fill()
turtle.left(90)
turtle.forward(10)
if a == 3:
a = 0
else:
a += 1
ward = sign[a]
if ward == 'd':
i += 1
if ward == 'w':
j -= 1
if ward == 'a':
i -= 1
if ward == 's':
j += 1
count += 1
print(count)
else:
arr[i][j] = 1
turtle.color('white', 'white')
turtle.begin_fill()
fill_squ(ward)
turtle.end_fill()
turtle.right(90)
turtle.forward(10)
if a == 0:
a = 3
else:
a -= 1
ward = sign[a]
if ward == 'd':
i += 1
if ward == 'w':
j -= 1
if ward == 'a':
i -= 1
if ward == 's':
j += 1
count += 1
print(count)
else:
turtle.mainloop()
if __name__ == '__main__':
i = int(input('输入屏幕的长\n'))
j = int(input('输入屏幕的宽\n'))
#arr = set_screen(i,j)
arr = numpy.ones((i,j))
draw(i, j)
这里我加入了一个set_screen(i,j),用于模拟完全随机的黑白块组合的效果,如果想用,可以在主函数那里把注释改掉.
这就是最后的结果,可以看到有一条长形的直线,只要计算机运算能力足够,它可以一直延长.(这个结果的计算过程比较久,我算了半个小时左右,也有可能因计算机不同而异).
兰顿蚂蚁的实现过程十分有趣,作为一名新人程序员我捣鼓了老长时间才出成果.关于兰顿蚂蚁,这种混沌中诞生的秩序令人深思,这究竟是否为必然结果?还是一种偶然?在给定的规则下,虽然一开始表现出随机与混沌,但最终还是表现出趋同性,这背后是否暗示了某种原理与宿命?我觉得它具有一定的深邃性,有待探讨.
本人为新人博主,文章中有何不妥与错误还望大佬们斧正.