分形理论及技术4.docVIP

一、关于分形的生长模型 客观世界大量存在的实际物体的几何结构大都是具有无规则自相似的分形，然而这些分形结构是如何形成的呢？也就是说，物质粒子是如何一步一步堆砌起来而生长成分形结构？这些年来有不少人提出了一些无规则分形结构生长的模型，使得该问题得到了比较满意的解答。关于无规则分形生长有基本物理机制、DLA模型、DBM模型等。 分形生长的基本物理机制： DLA模型： DBM模型： 二、分形的计算机模拟 由于计算机技术的发展，特别是计算机图形学的建立和发展，使得所有规则分形图案也都可以用构造相应的算法并通过计算机绘制出来，图形为分形理论的研究提供了最直观的几何图形，也可以从多姿多彩的分形图案中验证理论并总结出规律，不断丰富和完善分形理论。用计算机绘制或模拟分形图也是研究分形的一个重要手段。 这里以规则分形为研究对象，介绍几种用计算机绘制分形的算法等。 （1）迭代函数系统的迭代算法 迭代函数系统（Iteration Function System，IFS）是分形几何学的重要分支之一，它是由Barnsley于1981年提出的一个分形结构系统，现已成为分形几何中的重要研究内容之一。IFS是以仿射变换为框架，根据几何对象的整体与局部具有自相似结构，经过迭代而产生的。 仿射变换是对图形所作的绕原点旋转、比例放大（缩小）及平移等操作（方向改变、大小改变、位置改变）的几何变换： （在数学中，线性映射(也叫做线性变换或线性算子)是在两个向量空间之间的函数，它保持向量加法和标量乘法的运算。几何上，两个向量空间之间的一个仿射变换由一个线性变换加上一个平移组成） ，也可以记为 其中，为实数，此即为一个二维仿射变换。有下述四种平面仿射变换具有明显的几何意义： （a）缩放变换：；(b)伸长变换：； (c)剪切变换：；(d)旋转变换： 对应的仿射变换如下图示： 另外，式中的表示坐标系中的原点被平移到。 仿射变换的几何特性： （1）仿射变换的逆变换仍为仿射变换； （2）为线性变换，直线段经仿射变换后仍为直线段，并且保持线段上点的比例关系不变； （3）二平行直线变换后仍保持平行性； （4）任意平面图形经仿射变换后，其面积发生变化，为变换前的倍， 几何图形的全貌与局部，在上述仿射变换下，具有自相似结构，即若将图形局部放大会发现与原图是相似的。也就是说，图形的局部是整体的一个小复制品。因此，如果找到了描述图形过程的仿射变换，那末就可以把握住相应图形的构成要素，从而可以得到分形图案各个层次上的细节。 基于这一思想，Barnsley得出了将图形简化为一系列仿射变换的一般过程，提出了IFS。这一算法被广泛应用在计算机生成图形及图像数据压缩等方面。 例如，有下列迭代函数： 于是，若从单位间隔开始，则迭代函数将单位间隔变换成间隔，而将单位间隔变换成间隔。 再变换下去，，； ， ， (不断下去，变换的吸引子就是康托尔集合。 正像上式一样，对于给定的一个初始图像，通过个迭代函数，将变换成许多小的拷贝，那末这些变换的和集 就产生了新的结构。 在一个变换过程中，不断重复应用若干变换公式的变换系统成为迭代函数系统。 IFS迭代算法有两种：确定性迭代算法与随机迭代算法。 确定性迭代算法思路：（1）首先任意在平面上选择一个集合，例如是矩形或单位圆等。然后顺序使用每个仿射变换得到。令；再对类似地顺序使用每个变换，得到。如此反复，将得到一个集合序列，且该序列一定收敛于该IFS的吸引子。 具体步骤： （1）初始化，在计算机上任意给定一些点，设定最大迭代步数； （2）对这些点顺序计算他们经过各个变换后的点集并保存。 （3）擦去屏幕上原来的点，显示由（2）计算出的点。 （4）返回到第（2）步，直到达到最大的迭代步数为止。 以希尔宾斯基三角形为例，进行算例说明，对应的IFS映射见下表，所产生的迭代函数系统包含三个仿射变换。 a b c d e f 1 0.5 0 0 0.5 0 0 2 0.5 0 0 0.5 0.5 0 3 0.5 0 0 0.5 0 0.5 初始集合选取以为顶点的直角三角形，算法过程如下图所示。 第一次迭代过程与结果： 第二次迭代过程（部分）与结果： 继续迭代，可知，该迭代函数系统的吸引子就是希尔宾斯基三角形。 随机迭代系统：是一种随机地选取IFS中仿射变换的方法。 确定性算法的原理简单明了。但由于方法不考虑的区别，而在迭代中对他们的使用处于平等的地位（没有考虑不同仿射变换的面积压缩比的不同），执行起来比较浪费机时（如果某一个仿射变换的面积压缩比特别小，那末该变换对初始集合作用后，面积以几何级数缩小，由于计算机分辨率所限，很快会缩成一个点）。所以通常IFS的实现不采用确定性算法，而较多地使用随机迭代算法。即随机迭代算法的思想为：为了使压缩能同步进行，