分形读书报告

读书报告近期，我阅读了一些有关于分形方面的文献，大致的了解了一下分形以及分形的量度一—分数维，下面就是一些我对于分形的理解。在经典的欧几里德几何中，可以用直线、圆、球等这一类规则的形状去描述诸如墙、车轮、卫星等人造物体，因为这些物体本来就是根据欧氏几何的规则图形生成的。然而自然界中，却存在着许许多多极其复杂的形状，女口：山不是锥体，云不是球体，闪电不是折线，雪花边缘也不是圆等等，再如宇宙中点点繁星所构成的集合更非经典几何所能描述的，它们不再具有人们早已熟知的数学分析中的连续、光滑（可导）这一基本性质了，而是非线性的。为了描述这些问题，哈佛大学数学系教授曼德布罗特（BenoitB.Mandelbort）在1975年首次提出分形（Fractal）概念。1982年Mandelbrot著作《TheFractalGeometryofNature》的出版，标志着分形理论的产生。分形理论的建立为研究无序结构和探索复杂事物提供了一极有力的工具。分形理论与耗散结构理论、协同学、混沌理论都是同一时期在非线性科学研究中取得的重要成果。所谓分形就是事物组成部分以某种方式与整体相似的形其整体具有自相似性。分形研究的对象是具有自相似性的无序系统其维数的变化是连续的，而非欧氏几何中的整数维，如空间的欧氏维数是3。（琚正挺,2006）从Mandelbrot在《英国的海岸线有多长》一文中提出分数维概念以后，分形几何学逐渐发展成为专门研究复杂、非规则现象的新理论，并已被证实在研究过去常被认为的无规律体,如地质体的内在规律方面行之有效，分形能够对自然世界和表面的复杂性作出更精确的表达。分形具有自相似性和无标度性。（1）自相似性：一个分形的某种结构或过程从不同的空间或时间尺度来看都是自相似的。事实上，在标度区内具有对称性，即表征自相似系统或结构的定量性质如分维数，并不会因为放大或缩小等而变化，所改变的只是系统的外部形式，即系统的部分和整体之间存在自相似性。虽然这种定义不完备，但抓住了分形的本质特征一一自相似性。通常所说的自相似可以分为两类:一类是完全相似，由数学模型生成，如图1所示科赫曲线的构造:设E0为一单位线段，将其三等分，中间的1/3用边长为1/3的等边三角形向上指的另两条边代替,得到的集记为E1，它包含四条线段。对气的每条线段重复这一过程得到E2。归

图2谢尔宾斯基垫片(李伯奎等,2004)纳得到Ek和Ek+1。当k充分大时,Ek+1与Ek只在精细的细节上不同。当kTa时，极限曲线称为科契曲线，被人们用来做典型的海岸线模型，它可以刻画出真实海岸线的复杂性和粗糙程度。又如谢尔宾斯基垫的构造，数学过程简单描述为:在每步构造中都将前次的正三角形等分成4个小正三角形并去掉中间的一个，这一构造过程的极限图形是一曲线,称为谢尔宾斯基垫,如图2所示。另外一类就是自然界中的分形,如蜿蜒曲折的海岸线、云彩的形状等，其相似性并不是严格的图2谢尔宾斯基垫片(李伯奎等,2004) Eti J £1凸 i£图1Koch曲线(李伯奎等,2004)(2)无标度性：在具有分形性质的物体上任选一局部区域由于其自身具有自相似性，对它进行放大后,得到的放大图形会显示出原图的形态特性，即它的形态、内在的复杂程度、不规则形等各种特性，与原图相比均不会发生变化，如上面讨论的科契曲线的性质这种特性称为无标度性，又称为伸缩对称性。(李伯奎等,2004)分形维作为分形的标准量度，必然有它作为标准的经典模型维数一一自相似维数(Sel-fsimilarDimension),豪斯道夫维数(HausdorffDimension),盒计数维数(Box-countingDimension),功率谱维数(PowerspectrumDimension),结构函数法维数(StructurefunctionDimen-sion)。自相似维数(Sel-fsimilarDimension)自相似维数的引入受到规则形体如线段、正方形、立方体的启发。如果把线段、正方形和立方体的边分成两等份，这时线段是原来一半长度的两个线段,正方形被分成四个全等的小正方形，立方体则被分成八个全等的小立方体。也就是说，线段、正方形和立方体可被看成是由2、4、8个与整体相似的图形组成。2、4、8个这些数字可以改写成21、22、23,这里出现的指数分别与图形的欧氏维数与拓扑维数一致。一般地若把某个图形的长度(或标度)缩小1/r时得到N个和原图形相似的图形，有N=r-D，这里的指数D就具有维数的意义，称为自相似维数，用数学语言描述如下：如果一个集F由m个相等的且与F相似的部分组成，则称F为自相似集。若部分与F的相似比为r，则定义自相似维数为：D=-logm/logr自相似维数只对严格自相似的均匀一致的线性分形集有意义为了刻画更广泛的集类，需要引入更一般的维数Hausdorff维数。豪斯道夫维数(HausdorffDimension)其计算的基本原理为:分形集都遵循一定的标度律000，即测度M(5)随测量尺度5按照一种幕指数规律而变化，即M(5)正比5k将M(5)和5在双对数坐标中作图，并进行最小二乘拟合得一直线,其斜率K与分形维数D之间有如下关系：D=f(K)采用不同的测度，对应的函数也不同，在后面将分别进行说明。盒计数维数(Box-countingDimension)对于分形集F,N(5)是覆盖F的直径至多为5的集的个数,N(5)和5之间有幕律关系：N(5)正比5k在双对数坐标中拟合的logN(5)-log5直线的斜率K与分形维数D的关系为D=1-K盒维数是最简单也最明了的分形维数。在不同的标度下,用盒计数法来分析实际分形集的方法适用范围广，无论分形集是不连通的点集，还是曲线、曲面或立体都可用这一方法，除了自仿射分形。功率谱维数(PowerspectrumDimension)分形曲线若以功率谱Pg)为测度，以频率®为尺度，则有P(①)正比®k那么，所拟合的logP(®)-log®直线的斜率K与D的关系为D=(5+K)/2功率谱法适合于自仿射分形曲线,但在用于工程表面轮廓曲线的分维计算时，其幕律关系不很明显，误差较大,使用场合受到很大的限制。结构函数法维数(StructurefunctionDimen-sion)自相关函数(ACF)已经成为描述空间变量的最流行的方法。它毫无疑问地包含了有用的空间信息,然而当我们用ACF来研究已知形貌在磨损、变形或者某些类似过程中的变化时,许多变化由于集合平均而被掩盖着。如图3表示某一表面在磨损前后的轮廓形貌。图4表示该表面未磨损轮廓和已磨损轮廓的自相关函数曲线，图中磨损几乎在整波长上出现，因而在自相关函数图中二者变化甚微是显然的。［汕未磨损轮啣程样芒矍//rm(b)已磨损轮邮图3磨损前后的表面形貌比较(李伯奎等,2004)另外当从轮廓扩展到表面时，尤其是要求三维描述各向异性结构时，自相关函数(ACF)也会出现问题了。标准化的自相关函数必须利用轮廓方差，但是对各个角度上的轮廓来说轮廓方差是不同的，而且在自协方差函数上出现奇异点。如果用表面方差标准化ACF,那么在坐标原点会有奇异点，只有在相对于同一个中位面测量各轮廓时才会消除。以结构函数或者方差函数(SF)的方式来描述时就会消除这些问题，Sayles和Thomas(1977)对于轮廓结构函数定义为S(t)=E[Z(x+1)-Z(x)]2=f+TOS(w)(ejwt-1)dw=c/4~2D)—g这个函数是描述任意空间距离S上高度差均方的期望值。对平稳结构来说它含有与ACF相同的信息，该函数的两个主要优点是:它的意义不局限于平稳的情况，其次它不依赖于中位面。因此与中线有关的任何轮廓结构函数SF是表面结构函数SF的一个部分，这一性质ACF是不具备的。图5表示图3中的磨损和未磨损轮廓的结构函数SF曲线，从图中可以清晰地看到磨损和未磨损表面变化明显,表明了这种描述方法优于图4的自相关函数ACF方法。同时，拟合的logS(i)-logI直线,其斜率K与D有关系为D=(4-K)/2很多的研究表明,对于自仿射轮廓曲线，最具有意义的量就是结构函数,因此结构函数法特别适用于具有统计自仿射分形特征的工程表面轮廓曲线的分形维数计算，而机加工表面大多都具有自仿射分形的特征。(李伯奎等,2004)图4磨损(I)和未磨损(II)轮廓的自相关函数(李伯奎等,2004)图5磨损和未磨损轮廓的结构函数(李伯奎等,2004)然而，理论只有与实践联系在一起才具有真实的意义，则在观察自然界中的具有统计自相似性的景象时，我们需要一些方法去计算它们的分数维，经过许多科学家研究与努力,大致总结出了以下的几种方法：量规法量规法的思路是使用不同长度的尺子去度量同一段曲线(以海岸线为例,海岸线的长度L(r)由尺子长度r和尺子测量的次数N(r)来决定：L(r)=N(r)xr当海岸线的弯曲程度、复杂程度不同，且尺子长度r也出现变化,那么被测海岸线的长度也必然出现相应的变化，尺子长度越小,则所测得的海岸线长度值越接近被测海岸线长度的真实值。然而，随着尺子尺度的缩小，海岸线边上的各种小岛屿也就被纳入了测量的范围，所以,分形理论告诉人们，与传统的认知不一样的是，海岸线的长度更确切地来说是一个变量，它并不是描绘海岸线的一个完好的量度，而必须找到一个表征海岸线性质的客观量度，这就是分维数。根据Mandelbrot的研究，有下式成立：L(r)=Mxr1-D式中,L(r)为被测海岸线的长度;r为标度;M为待定常数;D为被测海岸线的分维数。

对式(2)两边同取双对数，可得:LgL(r)=(l-D)lgr+C式中,C为待定常数;该式斜率值等于1-D,即分维数D=l-K(该式的斜率值)。(朱晓华,2002)盒维数法将曲线用一边长等于1的方盒子覆盖，将此方盒分割成含有2n个小方盒的网格集,小方盒的边长为2-n,用这个网格集覆曲线，统计出与相交的小盒子数量M(n)，则曲线的分形维数为(王东升等,1995)D=limD=limnslogM(n)nlog2含量面积法如果把地球化学元素的数据记为x,y,z,其中x,y代表地理位置，z代表元素TOC\o"1-5"\h\zii i i i i含量，贝收,y,z构成的曲面称为含量曲面.先将含量曲面的投影平面用矩形网络分割i i i为边长为5x^5y的矩形，第k个矩形记为abcd,这4个点的投影高度为hhhak bk ckhdk(即4个点处元素的含量),选取含量尺度r,当hk,hbk,hk,hdk均大于或等于r时，计dk akbk ck dk算该投影网格对应的小曲面面积，近似面积公式为Sk心1[J6Sk心1[J6x2+(h-h)2x丿 akdk-h}+.6x2+(h-h^2ck' bkck-hbk则整个投影网络对应在曲面上的覆盖总面积可近似为S(r)=TNS(r)kk=1式中,N是小矩形数目。当取不同的r值,将得到不同的S(r)。为了求出分维数D,将观测数据S(r),S(r),,,S(r)和r,r,,,r绘在双对数1 2 n 12 n坐标图上，用最小乘法进行分段拟合，求出斜率D的估计量，即为分维数。(陈聆等，2004)为使含量-面积分形法所计算的结果更加精确,在估计的拐点两侧可适当加密分类的面积数。4）三角棱柱表面积法三角棱柱表面积法是通过比较在采用不同尺度观测图像“表面”时表面积的大小的变化情况来计算曲面的分形维数的一种方法。此方法利用栅格的4个角点（A、B、C、D）像素值来计算，而中心点的像素值为这4个角点像素平均值。中心点将栅格正方形分成4个三角形，分别计算这4个三角形的面积,4个三角形的面积之和即为栅格的表面积（图6）。改变栅格尺寸，重复上述计算，从而得到图像表面积与栅格尺度之间的关系，即可计算出图像的分维。（陈文凯等，2010）图6三角棱柱表面积法示意图（陈文凯等，2010）5）结构函数法将表面轮廓曲线视为一个时间序列z（x），则具有分形特征的时间序列能使其采样数据的结构函数满足E[z（x+九）—Z（x）1二C|九|4-2D（1〈D〈2）式中,E[Z（x+入）-Z（x）]2表示方差的算术平均值，t是数据间隔的任意选择值，即尺度标准。针对若干尺度t对轮廓曲线的离散信号计算出相应的差方的算术平均值，然后在对数坐标中得logE[Z（x+入）-Z（x）]2-logt直线的斜率k,则分形维数D与斜率k的转换关系为D=2-k/2也就是说得到k值，就可以得到分形维数D。（王建军等，2006）具有分形特征的是复杂系统，复杂程度可以用非整数维一一分数维来描述。各种不同的分形维数是集合划分不同层次的层次标号，它们从不同的角度对集合进行层次的划分。分形维数D度量了系统填充空间的能力，它从测度论和对称理论方面刻画了系统的无序性,是描述复杂对象的最基本特征。比如豪斯道夫维数是描述点集规则与不规则的几何尺度，同时其整数部分反映出图形的空间状态。对于动力学系统,豪斯道夫维数大体上表示了独立变量的数目。在工程表面评价中已经证明分形维数是不受仪器精度和测量基准影响的重要指标;在材料科学中已经发现分形维数与材料的某些性质参数有关;在化学领域，分形维数同催化剂的催化作用和选择性相关。在信号处理、地震的预报、石油的开采、生物生态学、经济因素的分析等方面，分形维数都有其独特的含义。（李伯奎等，2004）以上就是对于分形的一些了解，经过探讨,我们发现对于本次我们的课题长石环带中的微量元素的分析可以运用到分形中的盒维数法，通过盒维数法我们可以计算出北京周口店地区的房山岩体中的分数维,并推测出各种元素的大致分布,从而为后来者对于房山岩体的研究提供帮助。主要通过两种的盒维数法来计算分数维。1）剖面曲线的盒维数计算方法：通过我们已经获得的房山岩体的长石环带薄片上的各种元素的含量，我们将其制作成元素含量剖面曲线（如图7），并将元素含量剖面曲线嵌入平面空间，再将平面空间划分成边长为a的盒子，数出曲线所占据的盒子数Na,则曲线的长度L可以近似表示为：L=N0然后改变盒子的边长为b数出曲线所占据的盒子数Nb，获a aaA b得曲线的长度L，如此改变盒子的边长计算下去。若将盒子的边长记为r则获得的结果b满足如下关系：L=Nr正比口-drr式中D即为剖面曲线的分形维数，若将上式两边取对数，则有：lgN=-Dlg（r）+C（常数），r

则上述剖面曲线所获得的盒子数与所采用的盒子尺度在双对数坐标中将呈线性分布，从拟合直线的斜率（slope）即可获得剖面曲线的分维值（D=-slope）。11瓷至k11瓷至k图7湖南某地坞含量剖面曲线（龚庆杰等，2002）2）空间曲面的盒维数计算方法:将我们已知的各种元素的含量制作成元素含量空间曲面（如图8），再将曲面嵌入立体空间，并将立体空间划分为边长为a的立方体盒子，数出曲面所占据的盒子数N,则曲面的面积A可以近似表示为：A=Nxa2a a aa然后改变盒子的边长为b,数出曲面所占据的盒子数Nb,获得曲面的面积Ab；如此类推计b b算下去。若将盒子的边长记为r,则获得的结果满足如下关系：A=Nxr2正比（r2）2-drr式中D即为曲面的分形维数。上式两边取对数有：LgA=lg（Nrx"）=（2-D）lg3）+C（常数）lg（Nr）=（1-D）lg（r2）+C（常数）则由上述曲面所获得的盒子数或面积与所采用的盒子面