python三角型递归图案_Python数据结构与算法29：递归可视化：谢尔宾斯基三角

注：本文如涉及到代码，均经过Python 3.7实际运行检验，保证其严谨性。

本文阅读时间约为6分钟。

谢尔宾斯基三角

谢尔宾斯基(Sierpinski)是另外一种分形构造，如果是平面就称谢尔宾斯基三角形，如果是立体则称谢尔宾斯基金字塔。

真正的谢尔宾斯基三角形是完全不可见的，面积为0，周长无穷，是介于一维和二维之间的分数维(约1.585维)构造，一个很奇特的几何形状。

根据自相似特性，谢尔宾斯基三角形是有3个尺寸减半的谢尔宾斯基三角形，按照品字形拼叠而成。

由于无法真正作出谢尔宾斯基三角形(真正的谢尔宾斯基的degree趋向∞)，只能作近似谢尔宾斯基三角形(degree是有限的)。

谢尔宾斯基三角形的作图思路

在degree有限的情况下，degree=n的三角形，是由3个degree=n-1的三角形按照品字拼叠而成。同时，这3个degree=n-1的三角形边长均为degree=n的三角形的一半(规模减小)。

当degree=0，则就是一个等边三角形，这也是递归基本结束条件。

Pic-405-1 谢尔宾斯基三角形作图思路：degree=0, 1, 2, 3时的不同情况

经过上述分析，得出谢尔宾斯基三角形的绘制代码如下：

# 谢尔宾斯基三角形的作图。

import turtle

def sierpinski(degree, points):

colormap = ['blue', 'red', 'green', 'white', 'yellow', 'orange']

drawTriangle(points, colormap[degree]) # 绘制等边三角形。

if degree > 0: # 最小规模为degree=0，那时就直接退出。

# 以下三个情况都是调用自身，按照左、上、右的次序使得三角形的边长减半。

sierpinski(degree - 1,

{'left': points['left'],

'top': getMid(points['left'], points['top']),

'right': getMid(points['left'], points['right'])})

sierpinski(degree - 1,

{'left': getMid(points['left'], points['top']),

'top': points['top'],

'right': getMid(points['top'], points['right'])})

sierpinski(degree - 1,

{'left': getMid(points['left'], points['right']),

'top': getMid(points['top'], points['right']),

'right': points['right']})

def drawTriangle(points, color):

t.fillcolor(color)

t.penup()

t.goto(points['top'])

t.pendown()

t.begin_fill()

t.goto(points['left'])

t.goto(points['right'])

t.goto(points['top'])

t.end_fill()

def getMid(p1, p2): # 减小规模，边长减半。

return((p1[0] + p2[0]) / 2, (p1[1] + p2[1]) / 2)

t = turtle.Turtle()

# 外轮廓的三个顶点。

points = {'left': (-200, -100),

'top': (0, 200),

'right': (200, -100)}

sierpinski(5, points) # 绘制degree=5的谢尔宾斯基三角形。

# degree越大，海龟绘图时间越长；degree每增加1，绘图时间可能是原来的好几倍。

turtle.done()

degree=5的谢尔宾斯基三角形绘图成功时，如下图所示：

Pic-405-2 degree=5的谢尔宾斯基三角形

To be continued.