不连续曲线 highcharts_【分形几何】02.Koch曲线及推广

本文主要是对经典分形中Koch曲线和Koch雪花作简单的介绍，给出分形的生成规则，并对其几何性质进行分析。另外，讨论Koch曲线的推广形式。

1. 原始Koch曲线

Koch曲线是早期被描述的一种分形曲线。它由瑞典数学家Helge von Koch在1904年发表的一篇题为“从初等几何构造的一条没有切线的连续曲线”的论文中提出。虽然这条曲线处处连续，但处处不光滑，处处不可微。构造这类曲线是从初始元出发，如取一条线段或一个多边形。原则上是按一定的规则变线段为折线，形成生成元，分形的结构也就确定了。

设初始元是[0,1]，根据下列规则构造一个Koch曲线：

① 三等分一条线段；

② 用一个等边三角形两条边替代第①步划分三等分的中间部分；

③ 在每一条直线上，重复第②步。

Koch曲线是以上步骤地无限重复的极限结果，这就是Koch本人的原始构造（图1、2）。

图1 原始Koch曲线，第1到4级的构造

图2 原始Koch曲线

为讨论方便，设Koch曲线初始元的边长为

，边数为

，长度为

，依次所得第

级Koch曲线构造的边长为

，边数为

，长度为

①边长

的计算：

，得

②边数

的计算：

，得

③长度

的计算：

，得

由上面的计算，Koch曲线的长度是无穷大。另外，可证明在曲线的任一点都没有切线，即处处不可微。显然Koch曲线是自相似的且无自重叠，那么其分形维数可用相似维数来表示，我们先叙述一下相似维数的定义。

图3 相似维数的定义

由上述定义，因为Koch曲线可由4个非重叠且全等的部分组成，即

，相似比

，Koch曲线的相似维数

另外，Koch曲线可由很简单的L系统代码生成，如下：

1. 角度：δ=60°
2. 初始元：F
3. 生成元：F→F+F--F+F

有兴趣的读者，请参考Lindenmayer系统(简称为L系统)的相关资料。

IFS算法如下：

图4 Koch曲线的IFS算法

通过化简后，有

图5 化简的IFS算法

2. Koch曲线的Cesàro构造

Koch曲线的另一种构造是由意大利数学家蔡查罗(Ernesto Cesàro，1859—1906)创造，图6中是Cesàro构造的第1级到第4级。生成元是夹角为120°等长线段，其中红色细线是上一级的构造，黑色粗线是本级的构造。每一条红色线和两条黑色线组成一个等腰三角形，这样最终的Koch曲线总是包含在全部等腰三角组成的区域之中。

原始Koch曲线的构造只是Cesàro构造的子序列，相当于取n=2,4,6,8,…，时的Cesàro构造。

图6 Koch曲线的Cesàro构造

设初始元的边长为

，边数为

，长度为

，依次所得第

级Koch曲线构造的边长为

，边数为

，长度为

①边长

的计算：

，得

②边数

的计算：

，得

③长度

的计算：

，得

Cesàro构造的L系统代码生成，如下：

1. 角度：δ=30°
2. 初始元：L
3. 左生成元：L→+R--R+
4. 右生成元：R→-L++L-

3. Koch雪花

Koch雪花（也称为Koch岛）是由三个全等的Koch曲线组成的封闭曲线（图7），同样是处处连续，处处不可微的自回避曲线。

图7 Koch雪花第1到3级的构造

设初始三角形

的边长为

，边数为

，周长为

，面积为

；依次所得第

级Koch雪花

的边长为

，边数为

，周长为

，面积为

①边长

的计算：

，得

②边数

的计算：

，得

③周长

的计算：

，得

④面积

的计算：

因

是在

的每条边上再生成一个小三角形，

；同理，

是在

的每条边上再生成一个小三角形，

；综合①②可推出

，通过化简得：

Koch雪花的面积：

Koch雪花是周长无穷大，而面积有限的图形。并且可平面密铺（至少需要两种不同比例），图8中由6小1大Koch雪花组成一个更大的Koch雪花，长度比例是

等价于生成元由7部分组成，相似比

，维数方程

，解得

，分形是充满平面的。

这里不做过多介绍，后面会在Peano曲线的文中会有详细叙述。

图8 Koch雪花

Koch雪花的L系统代码生成，如下：

1. 角度：δ=60°
2. 初始元：F++F++F
3. 生成元：F→F+F--F+F

4. 改变生成元线段的角度

通过改变Koch曲线的Cesàro构造二条线段的角度

，可得到一条自回避的分形曲线（图9、10），在前面的一篇文章有介绍：

木木：【分形几何】01.等长二线段生成元​zhuanlan.zhihu.com

，相似比

，分形维数

特别的，当

时，曲线充满一个等腰直角三角形，分维

。

图9 角δ=43°的前6级构造

图10 角δ=43°的分形

L系统代码生成，如下：

1. 角度：δ=43°
2. 初始元：L
3. 左生成元：L→+R--R+
4. 右生成元：R→-L++L-

IFS算法如下：

图11 两线段的IFS算法

5. 改变相似比

改变原始Koch曲线的相似比（图12），变成有两种相似比的四线段生成元，第1相似比

， 第2相似比

，有

和

图12 四线段生成元

1. 当
   时，分形是自回避的（图13中
   ）；
2. 当
   时，有
   ，分形是恰接触的（图14、15）；
3. 当
   时，分形是自回避的（图16中
   ）。

图14的维数方程

，令

，方程变换为

，解得

分形维数

图13 第1相似比等于0.3

图14 第1相似比等于0.25

图15 分形第1级到第n级的叠加

图16 第1相似比等于0.2

另外，分形和0轴围的面积公式：

，其中

这里简单说一下，推导过程。定义一个序列

，实际就是第n级新增的面积系数 ：

面积系数

面积

当

时，系数

，面积

（图15）

当

时，系数

，面积

（图17）

图17 Koch曲线的面积

IFS算法如下：

图17 两种相似比的IFS算法

图14中四线段分形的外边界同样是分形（图18），生成元是由7条线段组成，相似比

，成立

维数方程

，令

，方程变换为

，解得

分形维数

图18 分形的外边界