python分形几何_分形造型（一）

一、引言

欧氏几何方法

以规则几何图形为研究对象，物体形状由方程来描述

欧氏几何主要描述工具

点，直线，平滑的曲线、平面

欧氏几何不能模拟自然景物

山脉、云、树、草、火、雾、波浪等

为解决复杂图形生成问题，分形造型应运而生

二、分形的历史

分形(Fractal)一词，是曼德布罗特创造出来的，其原意具有不规则、支离破碎等意义

分形几何学是一门以非规则几何形态为研究对象的几何学

分形理论创始人是美籍法国数学家B.B.曼德布罗特。他1924年11月20日生于波兰华沙。1985年获巴纳德(Barnard)奖章。1986年获富兰克林(Franklin)奖章。1967年发表于美国《科学》杂志上的“英国的海岸线有多长”的划时代论文，是他的分形思想萌芽的重要标志。1975年冬他创造了fractal一词。著作《分形：形状、机遇和维数》法文版于同年出版。

1967年B.B.Mandelbrot提出了“英国的海岸线有多长？”的问题，以1km为单位测量海岸线，得到的近似长度将短于1km的迂回曲折都忽略掉了，若以1m为单位测量，则能测出被忽略掉的迂回曲折，长度将变大，测量单位进一步变小，测得的长度将愈来愈大，这些愈来愈大的长度将趋近于一个确定值，这个极限值就是海岸线的长度。海岸线的长度是不确定的，长度依赖于测量单位。

三、分形几何与传统几何相比

从整体上看，分形几何图形是处处不规则的。例如，海岸线和山川形状，从远距离观察，其形状是极不规则的。

在不同尺度上，图形的规则性又是相同的。上述的海岸线和山川形状，从近距离观察，其局部形状又和整体形态相似，它们从整体到局部，都是自相似的。

四、分数维D定义

N相当于把原图复制出若干份

如果某图形是由把原图缩小为\(\frac 1 S\)的相似的N个图形所组成，有：

\[D=\frac {lnN}{lnS}

\]

的关系成立，则指数D称为相似性维数，D可以是整数，也可以是分数。

五、经典分形

Koch曲线

import turtle as t

def koch(size, n):

if n == 0:

t.fd(size)

else:

for angle in [0, 60, -120, 60]:

t.left(angle)

koch(size / 3, n - 1)

if __name__ == '__main__':

t.setup(1000, 1000)

t.pen(speed=0, pendown=False, pencolor='blue')

a, n = 400, 4

t.goto(-a / 2, a / 2 / pow(3, 0.5))

t.pd()

for i in range(3):

koch(a, n)

t.right(120)

t.ht()

t.done()

原来的一条边变成4条边，所以N=4

现在的一个小三角的边长是原来的\(\frac 1 3\)，所以S=\(\frac 1 3\)

所以\(D=\frac {ln4}{ln3}=1.2618\)，D介于1和2之间

Cantor集

Cantor集是由不断去掉线段中间的三分之一开区间而得出。

#!/apps/bin/python

from tkinter import *

class Cantor(Frame):

limit = 1

def __init__(self, master=None):

Frame.__init__(self, master)

self.draw = Canvas(self, width=800, height=500)

self.grid()

self.createWidgets()

def createWidgets(self):

# self.quitButton = Button(self, text="QUIT", command=self.quit)

# self.quitButton.pack(side=BOTTOM, fill=BOTH)

self.draw.pack(side=LEFT)

self.drawCanvas(100, 100, 700, 100)

def drawCanvas(self, ax, ay, bx, by):

self.draw.create_line(ax, ay, bx, by)

if (bx - ax) > self.limit:

cx = ax + (bx - ax) / 3

cy = ay + 50

dx = bx - (bx - ax) / 3

dy = by + 50

ay = ay + 50

by = by + 50

self.drawCanvas(ax, ay, cx, cy)

self.drawCanvas(dx, dy, bx, by)

app = Cantor()

app.master.title("Cantor (recursive)")

app.mainloop()

一条线变成2条线，所以D=2，长度是原来的\(\frac 1 3\)，所以S=\(\frac 1 3\)

所以

\[D=\frac {ln2}{ln3}=0.6309

\]

Sierpinski垫片

将一个等边三角形四等分，得到四个小等边三角形，去掉中间一个，保留它的三条边，将剩下的三个小等边三角形在分别进行四等分，并分别去掉中间一个，保留它们的边。

"""

功能：绘制谢尔宾斯基三角形

环境：python3.7

日期：2019/1/14 21:49

作者：指尖魔法师

版本：1.0

"""