200字范文,内容丰富有趣,生活中的好帮手!
200字范文 > Python递归绘制谢尔宾斯基三角形

Python递归绘制谢尔宾斯基三角形

时间:2019-04-20 04:20:37

相关推荐

Python递归绘制谢尔宾斯基三角形

首先,回顾递归的概念:

一个函数不停地调用它本身,我们就认为这个函数使用了递归。

先来看一个谢尔宾斯基三角形的示例:

所以首先,我们需要一个方法,它能帮我们画三角形,在这里画一个三角形,需要引入Turtle模块,画三角形需要三个点的坐标,不做过多的叙述。注;规定左下角顶点为p[0],上方顶点和右顶点分别为p[1]p[2]

from turtle import *def drawTriangle(points, color, myTurtle):myTurtle.fillcolor(color)myTurtle.up()myTurtle.goto(points[0])myTurtle.down()myTurtle.begin_fill()myTurtle.goto(points[1])myTurtle.goto(points[2])myTurtle.goto(points[0])myTurtle.end_fill()

其次,要绘制一个谢尔宾斯基三角形,我们需要从一个大三角形开始,找到每个边的中点,连接起来形成一个三角形(粉色部分)

在把中点连接起来之前,我们要找到中点,找到中点的方法大家高中都学过:

其中p1,p2是坐标,在这里我们用元组的形式表达。

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

然后我们要开始绘制谢尔宾斯基三角形了,下面考虑它的具体画法:如果我们只简单地分一层,则:

不难注意到,层数为0,就是一个三角形。

画一层,最主要的是找到它的三个中点,然后去绘制一个如上图所示的三角形。

如果层数为2,我们需要将周围的三个三角形再分别进行一次如上的处理,即,找到三个蓝色三角形的中点,然后把这三个中点连上线,忽略中间的三角形。我们不妨把上图看作一个单位,然后进一步地发现,高层数的谢尔宾斯基三角形就是用很多个这个单位构成的:如degree == 2:

所以考虑到它在不断地画相同的三角形,我们在这里当然要使用递归的方法来进行处理。

首先,我们要先画外轮廓的大三角形,代码如下:

def sierpinski(points, degree, myTurtle):colormap = ['blue', 'red', 'green', 'white', 'yellow','violet', 'orange']drawTriangle(points, colormap[degree], myTurtle)

解释一下,colormap是颜色的一个列表,不同的层次用不同的颜色,points接受一个由元组构成的列表,是整个外轮廓大三角形的坐标,然后调用上面的绘制三角形的函数进行绘制外轮廓。

前面提到,如果层数为0,就是一个三角形,此时我们当然没有进行绘制的必要,因此我们需要一个判断语句,再进入到递归的部分:

if degree>0:

接下来进入到了绘制谢尔宾斯基三角形的环节,我们可以先让Python绘制左下角的三角形,像分形树先走右子树一样,我们先去画左下角的三角形直到degree==0,则对左下角三角的递归完成,绘制左下角的三角形,需要的是大三角形左顶点的坐标和两个中点的坐标,这里可以使用getMid()完成

sierpinski([points[0],getMid(points[0], points[1]),getMid(points[0], points[2])],degree - 1, myTurtle)

像分形树一样,我们已经绘制出了左下角的三角

然后去绘制p[1]处的三角,我们需要其和p[0]p[2]的中点

sierpinski([points[1],getMid(points[0], points[1]),getMid(points[1], points[2])],degree - 1, myTurtle)

最后:

sierpinski([points[2],getMid(points[2], points[1]),getMid(points[0], points[2])],degree - 1, myTurtle)

绘制完成:

from turtle import *def drawTriangle(points, color, myTurtle):myTurtle.fillcolor(color)myTurtle.up()myTurtle.goto(points[0])myTurtle.down()myTurtle.begin_fill()myTurtle.goto(points[1])myTurtle.goto(points[2])myTurtle.goto(points[0])myTurtle.end_fill()def getMid(p1, p2):return ( (p1[0]+p2[0]) /2, (p1[1] + p2[1]) / 2)def sierpinski(points, degree, myTurtle):colormap = ['blue', 'red', 'green', 'white', 'yellow','violet', 'orange']drawTriangle(points, colormap[degree], myTurtle)if degree > 0:sierpinski([points[0],getMid(points[0], points[1]),getMid(points[0], points[2])],degree - 1, myTurtle)sierpinski([points[1],getMid(points[0], points[1]),getMid(points[1], points[2])],degree - 1, myTurtle)sierpinski([points[2],getMid(points[2], points[1]),getMid(points[0], points[2])],degree - 1, myTurtle)myTurtle = Turtle()myWin = myTurtle.getscreen()myTurtle.speed(1)myPoints = [(-200, -100), (0, 200), (200, -100)]sierpinski(myPoints, 3, myTurtle)myWin.exitonclick()

整体代码如上面所示,可以考虑把递归中只留下一行 代码,看看小乌龟 的行动轨迹,比如为什么单独运行第一个递归,代码就一直在画左下角的三角形,而剩下的两个不行呢?(考虑绘制三角形函数的定义)

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。