学习Python中turtle模块的基本用法（4：绘制科赫曲线和谢尔宾斯基三角形）

  科赫曲线和谢尔宾斯基三角形是常见的分形图形（详细介绍见参考文献1），本文使用turtle库绘制这两类图形。

科赫曲线

  科赫曲线的详细介绍见参考文献2，其中的绘图思路是“画正三角形，并把每一边三等分,取三等分后的一边中间一段为边向外作正三角形，并把这“中间一段”擦掉,重复上述两步，画出更小的三角形。”考虑到绘制完三角形再处理每条边的话，使用turtle的绘制函数不便于计算坐标，于是调整为绘制每条边时都将其三等分，只处理要显示的前后两个子边，重复该过程至指定的迭代次数（看下面录制的绘图动画比较直观）。最终的绘图代码及绘图效果如下所示。

from turtle import *

def Kohn(length, n):
    if n==0:
        forward(length)
    else:
        thirdLen=length/3
        Kohn(thirdLen, n-1)
        right(60)
        Kohn(thirdLen, n-1)        
        left(120)
        Kohn(thirdLen, n-1)  
        right(60)
        Kohn(thirdLen, n-1)  

screensize(600,400)
tracer(0, 0)
penup()
goto(-180,-90)
pendown()
Kohn(360,8)
left(120)
Kohn(360,8)
left(120)
Kohn(360,8)
left(120)
hideturtle()
        
exitonclick()

谢尔宾斯基三角形

  参考文献3中介绍的谢尔宾斯基三角形的画法有去掉中心法、Chaos Game法和L系统，本文中采用Chaos Game法绘制。最初采用参考文献3中的构造方法绘图，绘图代码及绘出的图形如下所示（判断点是否在三角形内的代码来自参考文献4）。：

def Sierpinski_1(x1,y1,x2,y2,x3,y3):
    maxX=max((x1,x2,x3))
    minX=min((x1,x2,x3))
    maxY=max((y1,y2,y3))
    minY=min((y1,y2,y3))        
    rx=random.randint(minX+1,maxX-1)
    ry=random.randint(minY+1,maxY-1)

    while(not IsInside(x1,y1,x2,y2,x3,y3,rx,ry)):
          rx=random.randint(minX+1,maxX-1)
          ry=random.randint(minY+1,maxY-1)
       
    tmp=((x1,y1),(x2,y2),(x3,y3))[random.randint(0,2)]
    mx=(tmp[0]+rx)/2
    my=(tmp[1]+ry)/2
    penup() 
    goto(mx-1,my)
    pendown()
    right(90)
    circle(1)
    left(90)

screensize(600,400)
tracer(0, 0)
x1=-250
y1=-150
x2=150
y2=-150
x3=0
y3=100

for i in range(0,5000):    
    Sierpinski_1(x1,y1,x2,y2,x3,y3)

hideturtle()

  估计应该是绘图思路不对，百度谢尔宾斯基三角形的混沌游戏画法，发现步骤都大致相同，但重复的步骤不一致，最终基于参考文献5-6中的介绍，调整绘图思路为：
  1、随机生成三个点A,B,C，组成三角形；
  2、随机生成三角形的一点P；
  3、绘制点P，同时计算P与A,B,C中任意一点的中点，并将中点坐标赋予P；
  4、重复步骤3。
基于上述步骤，重新调整代码，绘图代码及绘图效果如下：

def Sierpinski_2(x1,y1,x2,y2,x3,y3,n):
    maxX=max((x1,x2,x3))
    minX=min((x1,x2,x3))
    maxY=max((y1,y2,y3))
    minY=min((y1,y2,y3))        
    px=random.randint(minX+1,maxX-1)
    py=random.randint(minY+1,maxY-1)

    while(not IsInside(x1,y1,x2,y2,x3,y3,px,py)):
          px=random.randint(minX+1,maxX-1)
          py=random.randint(minY+1,maxY-1)
    penup()
    goto(px-1,py)
    pendown()
    right(90)
    circle(1)
    left(90)

    points=((x1,y1),(x2,y2),(x3,y3))
    for index in range(0,n):
        tmp=points[random.randint(0,2)]
        px=(tmp[0]+px)/2
        py=(tmp[1]+py)/2
        penup() 
        goto(px-1,py)
        pendown()
        right(90)
        circle(1)
        left(90)

screensize(600,400)
tracer(0, 0)
x1=-250
y1=-150
x2=150
y2=-150
x3=0
y3=100

Sierpinski_2(x1,y1,x2,y2,x3,y3,5000)

hideturtle()

  最后是一个完成的随机版谢尔宾斯基三角形绘制代码

from turtle import *
import random

def IsTrangleOrArea(x1,y1,x2,y2,x3,y3):
    return abs((x1 * (y2 - y3) + x2 * (y3 - y1) + x3 * (y1 - y2)) / 2.0)

def IsInside(x1,y1,x2,y2,x3,y3,x,y): 
    ABC = IsTrangleOrArea(x1,y1,x2,y2,x3,y3)
    PBC = IsTrangleOrArea(x,y,x2,y2,x3,y3)
    PAC = IsTrangleOrArea(x1,y1,x,y,x3,y3)
    PAB = IsTrangleOrArea(x1,y1,x2,y2,x,y) 
    return (ABC == PBC + PAC + PAB)

    maxX=max((x1,x2,x3))
    minX=min((x1,x2,x3))
    maxY=max((y1,y2,y3))
    minY=min((y1,y2,y3))        
    px=random.randint(minX+1,maxX-1)
    py=random.randint(minY+1,maxY-1)

    while(not IsInside(x1,y1,x2,y2,x3,y3,px,py)):
          px=random.randint(minX+1,maxX-1)
          py=random.randint(minY+1,maxY-1)
    penup()
    goto(px-1,py)
    pendown()
    right(90)
    circle(1)
    left(90)

    points=((x1,y1),(x2,y2),(x3,y3))
    for index in range(0,n):
        tmp=points[random.randint(0,2)]
        px=(tmp[0]+px)/2
        py=(tmp[1]+py)/2
        penup() 
        goto(px-1,py)
        pendown()
        right(90)
        circle(1)
        left(90)

width=600
height=400

screensize(width,400)
tracer(0, 0)

x1=random.randint(1-width/2,width/2-1)
y1=random.randint(1-height/2,height/2-1)
x2=random.randint(1-width/2,width/2-1)
y2=random.randint(1-height/2,height/2-1)
x3=random.randint(1-width/2,width/2-1)
y3=random.randint(1-height/2,height/2-1)
Sierpinski_2(x1,y1,x2,y2,x3,y3,5000)
hideturtle()
exitonclick()

参考文献
[1]https://baike.baidu.com/item/%E5%88%86%E5%BD%A2/85449?fr=aladdin
[2]https://baike.baidu.com/item/%E7%A7%91%E8%B5%AB%E6%9B%B2%E7%BA%BF/7090673?fromModule=lemma_inlink
[3]https://baike.baidu.com/item/%E8%B0%A2%E5%B0%94%E5%AE%BE%E6%96%AF%E5%9F%BA%E4%B8%89%E8%A7%92%E5%BD%A2?fromModule=lemma_inlink
[4]https://blog.csdn.net/mxw322/article/details/56026607
[5]https://zhuanlan.zhihu.com/p/103610498?utm_source=wechat_session&ivk_sa=1024320u
[6]https://wiki.swarma.org/index.php?title=%E8%B0%A2%E5%B0%94%E5%AE%BE%E6%96%AF%E5%9F%BA%E4%B8%89%E8%A7%92%E5%BD%A2