曼德布洛特集
若一个复二次多项式 [math]\displaystyle{ f_c(z)=z^2+c }[/math],从[math]\displaystyle{ z=0 }[/math]的前提下,开始递归计算,存在复数[math]\displaystyle{ c }[/math]使得该方程无限迭代后的结果能保持有界(即不发散),将满足上述条件的复数[math]\displaystyle{ c }[/math]的集合视为一种特殊集。
阿德里安 · 杜阿迪 Adrien Douady为纪念数学家伯努·瓦曼德布洛特 Benoît B. Mandelbrot而将该特殊集命名为曼德布洛特集 Mandelbrot set 。[1]曼德布洛特集中的复数[math]\displaystyle{ c }[/math]使得数列[math]\displaystyle{ f_c(0) , f_c(f_c(0)) , \dotsc }[/math]等在取绝对值后仍保持有界。该集合与朱利亚集 Julia set有着很深的内在联系,由于他们都使用相同的复二次多项式来进行迭代。朱利亚集也会产生与曼德布洛特集相类似的分形图案。
曼德布洛特集图像可以通过选取不同的复数[math]\displaystyle{ c }[/math],并观察其数列 [math]\displaystyle{ f_c(0) , f_c(f_c(0)) , \dotsc }[/math]是否达到无穷大(实践中通常指其是否在预设迭代次数后离开预设的某个含0在内的有界邻域)来实现:将c的实部作为复平面图像 Complex Plane的横坐标,虚部作为复平面图像的纵坐标,然后根据数列[math]\displaystyle{ |f_c(0)|,|f_c(f_c(0))|,\dotsc }[/math]超过主观设定阈值的速度来给每个点染色,若复数[math]\displaystyle{ c }[/math]的值在经过预设迭代次数之后没有超过阈值(注意:这里是将曼德布洛特集图像与其补集图像区分开来的标志)则染成特殊的颜色(通常是黑色)。如果复数[math]\displaystyle{ c }[/math]保持不变,而[math]\displaystyle{ z }[/math]的初始值——通常由[math]\displaystyle{ z_0 }[/math]来表示——却是一个变量,那么就能得到在简单函数的参数空间 Parameter Space中复数[math]\displaystyle{ c }[/math]对应的朱利亚集。
曼德布洛特集的图像展示了一条无比精致又无限复杂的分界线,将其不断的放大,就可看到更加精密的、基于递归的细节部分。也就是说,曼德布洛特集的分界线是一条分形曲线。重复出现的细部“样式”由所观测的集合区域所决定。曼德布洛特集的分界线上也包含了集合整体形状的较小版本,因此自相似性的分形特性不仅适用于该集合的局部,还适用于整个集合。
因为曼德布洛特集的图像在美学上拥有独特的吸引力,并且还是一个根据简单规则产生复杂结构的代表例子,这使得曼德布洛特集在数学之外的其他领域中也十分流行。它也是数学可视化和表现数学之美的最著名例子之一。
历史
曼德尔洛特集起源于20世纪初由法国数学家皮埃尔费托 Pierre Fatou 和加斯顿茱莉亚 Gaston Julia 首先研究的复杂动力学。首次确切定义分形,并绘制出可视化的分形图案得益于 罗伯特·W·布鲁克斯 Robert W. Brooks和彼得·马特尔斯基 Peter Matelski在1978年对克莱尼群 Kleinian Groups 的部分研究工作。[2] 在此基础上,1980年3月1日,在位于纽约的约克敦海茨 Yorktown Heights的 IBM的汤玛士·J·华生研究中心 Thomas J. Watson Research Center,伯努·瓦曼德布洛特 Benoît B. Mandelbrot首次绘制出曼德布洛特集的可视化图形。[3]且Benoît B. Mandelbrot在1980年发表了一篇关于二次多项式的参数空间 Parameter Space的研究论文。[4]
对曼德布洛特集的数学研究真正始于数学家阿德里安· 杜阿迪 Adrien Douady和约翰·H·哈伯德 John H. Hubbard [1] 的一系列研究工作。他们探明了曼德布洛特集的许多基本性质,Adrien Douady 为纪念Benoît B. Mandelbrot在分形几何中做出的杰出贡献,将该集合命名为曼德布洛特集。
数学家海因茨-奥托·佩特根 Heinz-Otto Peitgen 和彼得·里希特 Peter Richter通过照片、书籍[5],在德国歌德学院 Goethe-Institut举办国际巡回展览等宣传方式,让曼德布洛特集进入大众的视野中,受到广泛的关注。 [6][7]该封面展示了一幅以(-0.909,-0.275)为坐标的曼德布洛特集图形,该图形由Peitgen创作。[8][9] 1985年8月《科学美国人 Scientific American 》的封面文章向广大读者介绍了计算曼德布洛特集的算法。20世纪80年代中期,当个人计算机的功能变得足够强大,可以绘制图形并以高分辨率显示这些图形时,曼德布洛特集被运用到计算机图形学的图像演示中,并日益凸显了它的重要性。[10]
Adrien Douady和 John H. Hubbard 的研究工作不断取得新的成果,同时感兴趣于复杂动力学和抽象数学 Abstract Mathematics领域的队伍快速壮大,自此,深入了解曼德布洛特集一直是这些领域的核心研究。包括米哈伊尔·柳比奇 Mikhail Lyubich,[11][12] 科特·麦克马伦 Curt McMullen, 约翰·米尔诺 John Milnor,石仓光博 Mitsuhiro Shishikur和'让-克里斯托夫·约科兹 Jean-Christophe Yoccoz在内的许多人在曼德布洛特集的研究工作中,都作出了大大小小的贡献。
定义
曼德布洛特集是满足使得复二次多项式[math]\displaystyle{ z_{n+1}=z_{n}^{2}+c }[/math]中,从0开始进行迭代将方程进行无限迭代,并保持在有界范围内的复平面上的[math]\displaystyle{ c }[/math]的集合。[13]
因此,该复数[math]\displaystyle{ c }[/math]满足以下条件:如果从[math]\displaystyle{ z_0=0 }[/math] 开始,进行重复迭代,所取的复数[math]\displaystyle{ c }[/math]使得[math]\displaystyle{ z_n }[/math] 的绝对值都保持有界(不发散)。(对任意[math]\displaystyle{ n\gt 0 }[/math])
例如,令[math]\displaystyle{ c=1 }[/math] ,函数值数列为0,1,2,5,26,... ,发散,所以复数1 ∉ 曼德布洛特集。
令 [math]\displaystyle{ c=-1 }[/math],函数值数列是0,-1,0,-1,0,... ,收敛,所以复数-1 ∈ 曼德布洛特集。
曼德布洛特集也可以定义为一族多项式的连通轨迹 Connectedness Locus。
基本性质
由于曼德布洛特集是一个封闭图形,且包含在以原点为中心,以2为半径的封闭圆盘中。故曼德布洛特集是一个紧集 Compact Set,更具体地说,对于 [math]\displaystyle{ n\geq 0. }[/math] 而言,当且仅当[math]\displaystyle{ |P_c^n(0)|\leq 2 }[/math]时,复数[math]\displaystyle{ c }[/math]属于曼德布洛特集。也就是说,若[math]\displaystyle{ |P_c^n(0)| }[/math]大于2,则说明其函数值数列发散,所对应的复数[math]\displaystyle{ c }[/math]不属于曼德布洛特集。
曼德布洛特集与复平面的实轴交点所构成的区间为[math]\displaystyle{ [−2, 1/4] }[/math]。属于该区间的参数c满足 逻辑斯蒂映射 Logistic Map的双射条件。 其中,逻辑斯蒂映射是一种二次的多项式映射(递归关系式),是一个由简单非线性方程序产生混沌现象的经典范例。
- [math]\displaystyle{ x_{n+1} = r x_n(1-x_n),\quad r\in[1,4]. }[/math]
其对应关系为 [math]\displaystyle{ z = r\left(\frac12 - x\right),\quad c = \frac{r}{2}\left(1-\frac{r}{2}\right). }[/math]
事实上,这反映了逻辑斯蒂映射和曼德布洛特集的参数空间 Parameter Space之间的对应关系。
关于曼德布洛特集是否连通,许多数学家都有过研究:
Adrien Douady和 John H. Hubbard在曼德布洛特集的补集与闭合单位圆盘(以原点为中心,以1为半径做圆)的补集之间构造了一个显式的共形同构 Conformal Isomorphism ,由此证明了曼德布洛特集是连通的。 Benoît B. Mandelbrot 最初猜测曼德布洛特集是不连通的。这是由于当时计算机程序的局限性,导致程序无法检测到所生成的曼德布洛特集图形中连接不同细部的微小连线。通过进一步的实验后,Benoît B. Mandelbrot 纠正了之前的看法,认为曼德布洛特集是连通的。
此外,杰里米 · 卡恩 Jeremy Kahn在2001年利用严格的拓扑证明,论证了曼德布洛特集的连通性。[14]
由Adrien Douady和 John H. Hubbard证明曼德布洛特集连通时,所用到的曼德布洛特集的补集均匀化的动力学公式,引出了曼德布洛特集的外部尾迹射线。可将这些射线进行组合来研究曼德布洛特集,形成了 Yoccoz拼图的组合技术。[15]
曼德布洛特集的分界线是二次族 bifurcation locus的分岔轨迹 bifurcation locus,即参数[math]\displaystyle{ c }[/math]在极微小的变化下会产生很突然的动力学变化。该分界线是一类多项式双曲线,可被构造成一系列平面代数曲线的极限集。 通过设置 p0 = z, pn+1 = pn2 + z,然后将复平面上的点集pn(z)| = 2 转化为实笛卡尔平面 Cartesian Plane上度为2n+1的在 x 和 y 上的曲线,从而对该曲线进行定义。使用下面所提及的“逃逸时间算法”计算并绘制出来的曼德布洛特集图像中也显示了该分界线。
其他性质
心脏形结构和圆盘形“芽苞”Main cardioid and period bulbs
观察曼德布洛特集的图像,会首先看到中间的心脏形结构。(即图11中的区域1)该心形曲线内部的参数[math]\displaystyle{ c }[/math]若满足[math]\displaystyle{ c = \frac\mu2\left(1-\frac\mu2\right) }[/math], 其中 [math]\displaystyle{ \mu }[/math] 属于单位开圆盘(即在单位圆内的开集区域,开集定义:设A是度量空间X的一个子集。如果A中的每一个点都有一个以该点为中心的邻域包含于A,则称A是度量空间X中的一个开集。),则称这样的点[math]\displaystyle{ c }[/math]为复二次映射的周期点,记作[math]\displaystyle{ P_c }[/math]。
类似这样的“芽苞”突起有无穷多个,它们都与主心脏形结构曲线相切。
且满足以下定义:
对于任意有理数 [math]\displaystyle{ \tfrac{p}{q} }[/math] ,[math]\displaystyle{ p }[/math]、[math]\displaystyle{ q }[/math]互素,若参数[math]\displaystyle{ c }[/math]满足:
- [math]\displaystyle{ c_{\frac{p}{q}} = \frac{e^{2\pi i\frac pq}}2\left(1-\frac{e^{2\pi i\frac pq}}2\right). }[/math]
则该点处存在一个与主心脏形结构曲线相切的"芽苞",且记作“ [math]\displaystyle{ \tfrac{p}{q} }[/math]-芽苞”。
公式中的参数其由吸性循环周期的周期值q和组合旋转数 [math]\displaystyle{ \tfrac{p}{q} }[/math] 组成。包含周期为q的吸性周期循环的Factou 集 Fatou components 都在吸性不动点相交。若记分量[math]\displaystyle{ U_0,\dots,U_{q-1} }[/math]为逆时针方向,[math]\displaystyle{ P_c }[/math]将分量[math]\displaystyle{ U_j }[/math] 映射到分量[math]\displaystyle{ U_{j+p\,(\operatorname{mod} q)} }[/math]
[math]\displaystyle{ c_{\frac{p}{q}} }[/math] 为分枝 Bifurcation点,在该点吸性的不动点与斥性的周期循环发生相反的变化。当经过分枝点进入[math]\displaystyle{ \tfrac{p}{q} }[/math]-“芽苞”时,即[math]\displaystyle{ c_{\frac{p}{q}} }[/math] 从心脏形结构曲线内部变为开圆内部时,不动点[math]\displaystyle{ \alpha }[/math]由吸性变为斥性,周期循环由斥性变为吸性。
双曲分量 Hyperbolic components
上述所提及的圆盘形“芽苞”都是曼德布洛特集的内部分量,其中映射[math]\displaystyle{ P_{c} }[/math]具有吸性周期循环,这样的分量称为双曲分量。
推测这些双曲分量是否仅存在于曼德布洛特集内部区域,称为双曲密度问题。这或许也是复动力学领域中最值得关注的公开问题。曼德布洛特集中假设的非双曲分量通被称为“怪胎”或幽灵分量。[16][17]对于实二次多项式,在1990年,双曲密度问题得到解决。柳比奇 Lyubich 、格拉奇克 Graczyk、斯维亚特克 Świątek分别独立证明了双曲分量仅存在于曼德布洛特集的内部区域。
并不是每个双曲分量都可由曼德布洛特集的主心脏形结构经过一系列的直接分叉即可得到。但像图(15)的双曲分量可由小的曼德布洛特集副本的主心脏结构曲线经过一系列的直接分叉得到。
每个双曲分量都有一个中心,该中心记作点[math]\displaystyle{ c }[/math],使得[math]\displaystyle{ P_{c}(z) }[/math]的内部Fatou区域具有一个超吸性周期循环,即其吸引力是无穷的。这意味着该循环包含临界点0,因此经过数次迭代后,0会回到本身。
也就是说,将0带入[math]\displaystyle{ P_{c}(z) }[/math],存在[math]\displaystyle{ n }[/math],使得:
- [math]\displaystyle{ P_c }[/math]n[math]\displaystyle{ (0) = 0 }[/math]。
如果将公式的变量固定为[math]\displaystyle{ c }[/math]而不是[math]\displaystyle{ z }[/math],得到:
- [math]\displaystyle{ Q^{n}(c) }[/math]。
由[math]\displaystyle{ Q^{n+1}(c) = Q^{n}(c)^{2} + c }[/math],且令[math]\displaystyle{ Q^{n}(c)的纬度为2^{n-1} }[/math]。
因此,我们可以通过依次求解方程[math]\displaystyle{ Q^{n}(c) = 0, n = 1, 2, 3, ... }[/math]来构造双曲分量的中心。每一步产生的新中心的个数由斯隆的OEIS: A00074给出。OEIS的全称为:The On-Line Encyclopedia of Integer Sequences™ (OEIS™),它是一个关于整数序列(数列)的专业型网站,是一个关于组合数学研究的重要的网站,里面包含了众多数列的研究成果。
连通性 Local connectivity
在Adrien Douady 和 John H. Hubbard的研究工作中,提出了著名的MLC猜想:他们推测曼德布洛特集是局部连通的。这一猜想将曼德布洛特集简单抽象成一个“压缩圆盘”。 其中,它还体现了上述所提及的双曲分量的思想。
让-克里斯托夫·约科兹 Jean-Christophe Yoccoz证明了在所有有限可重整化参数下的曼德布洛特集的局部连通性; 也就是说,该种局部连通性只体现在有限多个小曼德布洛特集副本中。[18] 从那时起,在曼德布洛特集的其他参数点也体现了局部连通性。
自相似 Self-similarity
在 Misiurewicz 点附近,对曼德布洛特集进行放大,能够观察到自相似性。将其收敛于一个极限集后,我们还推测在广义 Feigenbaum 点(例如-1.401155或-0.1528 + 1.0397 i)周围能够观察到自相似的特征。[19][20]
总的来说,曼德布洛特集合并不是严格意义上的具有自相似特征的集合,但它具有准自相似性,因为在任意小的空间尺度上,都可以找到与自身略有不同的小副本。小副本之间的细微差别体现在它们与整体集合之间的连接细线上。
进一步的计算结果
Mitsuhiro Shishikura计算出曼德布洛特集分界线的豪斯多夫维数 Hausdorff dimension 为2。[21] 现在还不知道曼德布洛特集分界线是否具有正平面上的勒贝格测度 Lebesgue Measure。
利用可在实数域内计算的 Blum-Shub-Smale 模型 Blum-Shub-Smale Model (简称BBS模型),可得出以下结论:虽然无法对于曼德布洛特集进行计算,但可计算其补集。BBS模型也不可计算包括求幂图在内的许多比较简单的图像。目前,在基于可计算性分析 Computable Analysis(该种分析方式更接近于直观概念上的“由计算机绘制集合”)的实际计算模型中,曼德布洛特集是否可算尚未有定论。由该种模型计算集合的方式更接近于直观概念上的“由计算机绘制集合”。 赫特林 Hertling 证明了若双曲性猜想是正确的,则基于可计算性分析的实数域模型中的曼德布洛特集也是可计算的。
与朱利亚集之间的关系
朱利亚集 Julia set可由[math]\displaystyle{ f_c(z)=z^2+c }[/math]反复迭代得到。与曼德布洛特集不同,其对于复数[math]\displaystyle{ c }[/math]进行固定,取某一个[math]\displaystyle{ z }[/math]值,(如[math]\displaystyle{ z=z_0 }[/math]),可以得到数列[math]\displaystyle{ z_0,f_c(z_0),f_c(f_c(z_0)),f_c(f_c(f_c(z_0))).... }[/math]。该数列可能趋于无穷大或者始终处于某一范围之内并收敛于某一值。将使该数列不发散的[math]\displaystyle{ z }[/math]值集合称为朱利亚集。
由曼德布洛特集的定义,在给定点处曼德布洛特集的几何结构与相应的朱利亚集的几何结构有着密切的对应关系。 例如,当一个点在曼德布洛特集聚集时,相应的朱利亚集是连通的。
事实上,该原理被运用到曼德布洛特集的所有深层结果中。 例如,Shishikura 证明了对于在曼德布洛特集分界线上的一组稠密的参数,相对应的朱利亚集的豪斯多夫维数为2,然后将这些信息传递到参数平面上。[21] 同样,Yoccoz 首先证明了朱利亚集的局部连通性,然后在相应的参数点处进一步证明曼德布洛特集的局部连通性。[18]Adrien Douady将这一原则表述为:在动力平面上耕耘,在参数空间中收获。
几何结构
对于任意有理数[math]\displaystyle{ \tfrac{p}{q} }[/math],其中p 和q互素,周期 q 的一个双曲分量从主心脏形结构分支出来。 在该分枝点上与主心脏形结构相连的曼德布洛特集部分称为 p / q-limb。 计算机实验表明,分枝体半径像[math]\displaystyle{ \tfrac{1}{q^2} }[/math]一样趋于0,目前最佳的拟合形式是约科兹不等式 Yoccoz-inequality,其尺寸像[math]\displaystyle{ \tfrac{1}{q} }[/math]一样趋于0。 周期为[math]\displaystyle{ q }[/math]的分枝体顶端有[math]\displaystyle{ q-1 }[/math]个天线。我们可以通过计数这些触角反过来推算该特定分枝体的周期为多少。
曼德布洛特集中的π
为了证明 p / q-limb 的厚度为零,大卫·鲍尔 David Boll 在1991年利用计算机计算了使得 [math]\displaystyle{ z =-\tfrac{3}{4} + i\epsilon }[/math] ([math]\displaystyle{ -\tfrac{3}{4} }[/math]的级数发散所需的迭代次数。([math]\displaystyle{ -\tfrac{3}{4} }[/math]是所处的位置)。由于z = [math]\displaystyle{ -\tfrac{3}{4} }[/math]为确切值,并不发散。则需要进一步缩小 ε、增加迭代次数,更加接近目标值。譬如:当ε = 0.0000001时,需迭代31415928次,此时输出值的π值为3.1415928。[22]
曼德布洛特集中的斐波那契数列
可以看出,斐波那契数列 Fibonacci Sequence位于曼德布洛特集内,且与主心脏形结构曲线和法雷数列图存在关联。
将主心脏结构曲线映射到一个圆盘上时,会观察到由下一个最大的双曲分量延伸出来的天线,其位于之前选择的两个分量之间。天线的数目变化是一个斐波那契数列。
天线数目也与法雷数列图有关,由于n级法雷数列是指对任意给定的一个自然数n,将分母小于等于n的不可约的真分数按升序排列,并且在第一个分数之前加上0/1,在最后一个分数之后加上1/1这个序列称为n级法雷数列,以[math]\displaystyle{ F_n }[/math]表示。如[math]\displaystyle{ F_5 }[/math]为:0/1, 1/5, 1/4, 1/3, 2/5, 1/2, 3/5, 2/3, 3/4, 4/5, 1/1。
可得到:在曼德布洛特集中的法雷数列,其分母数量在相应的分数值内,分数值与圆盘周围的距离有关。这两个部分的分数值本身可以在[math]\displaystyle{ {\frac{1}{3}} }[/math]之后相加,以计算出序列中下一个双曲分量的位置。因此,在曼德尔勃特集合中可以找到1、2、3、5、8、13和21的斐波那契数列。
不同缩放比率下的图像库
当一个人看得越近或放大图像时,达到的放大效果可以让曼德布洛特集显示出更复杂的细节。当将下图不断的放大,缩放到选定的[math]\displaystyle{ c }[/math]值形成一个反映其变化的图集时,会给人一种不同几何结构中蕴藏着无限想象力的感觉。下也对于它们的一些典型规则进行说明。
zoom 03相对于zoom 00的放大倍率之比约为[math]\displaystyle{ 10^{10}:1 }[/math]。以普通的显示器显示出该图像的完整图形,其表示直径为400万公里的的曼德布洛特集合的一部分,其分界线将显示出数量庞大的不同分形结构。
海马的“身体”由25个“辐条”组成。25个“辐条”被分为三组,其中两组中各含有12个“辐条”,另一组单独由一个连接到主心脏形结构的“辐条”组成。这各含有12个“辐条”的两组可以通过某种变形变集合的“上臂”为曼德布洛特集的两根“手指”;中心是通常说的 Misiurewicz 点。在海马的“身体的上半部分”和“尾巴”之间,可发现一个扭曲的曼德布洛特集的小副本。该集合称为“卫星集”,也就是附属集。
上面的岛屿似乎是由无限多的部分组成,就像康托集 Cantor Sets一样。但实际上这是对应的朱利亚集[math]\displaystyle{ J_c }[/math]的情况。然而,它们通过微小的结构进行连接,其整体代表了一个单连通集。
这些微小的结构在中间的卫星集处相遇,在这样的放大倍率下,该卫星集太小以至于无法被识别。相应的[math]\displaystyle{ J_c }[/math] 的[math]\displaystyle{ c }[/math]值不是该图像中心所对应的[math]\displaystyle{ c }[/math]值,而是第六个缩放步骤中,相对于曼德布洛特集的主体,其中心所示的相同位置上的卫星集所对应的[math]\displaystyle{ c }[/math]值。
曼德布洛特集和朱利亚集的3D图像
除了能绘制出曼德布洛特集的二维图像外,还可以利用各种技术绘制出曼德布洛特集和朱利亚集的三维图像。具体做法为:将其二维图像中的每个像素赋予一个高度值,即可生成三维图像。
- 最简单的三维图像渲染方法是将每个像素的迭代值作为高度值,从而生成具有不同高度值的“台阶”的图像。
- 如果使用分数迭代值(也称为势函数)来计算每个点的高度值,则可以避免在生成的图像中执行步骤。 然而,使用分数迭代数据渲染的三维图像看起来比较粗糙且不太美观。
- 另一种方法是使用每个点的距离估计[23] (DE)数据来计算高度值。利用指数函数对距离估计值进行非线性映射,可以得到感观较好的图像。使用 DE 数据绘制的图像往往在视觉上更引人关注。更重要的是,三维图像使得连接图像中各点的细“卷须”更易于观察。 在“分形图像的科学”第121页的图29、30显示了使用外部距离估计绘制的二维和三维图像。
- 以下是前面“不断缩放的曼德布洛特集图集”的3D版本,在高度值图中利用距离估计将其进行渲染,并使用类似的裁剪和着色。
下图类似于上图中的“zoom 05”,它来自于《分形之美》这本书的第85页图44的3D图像版本。[24]
扩展
多重曼德布洛特集 Multibrot sets
在数学中,Multibrot集是复平面上的有限集,集中的元素取自一般一元多项式递归族中。如:
- [math]\displaystyle{ z \mapsto z^d + c.\ }[/math]
这个名字是多重和曼德布洛特集的混合词。同样的也适用于朱利亚集合,称为多重朱利亚集 Multi-Julia set。
对于整数[math]\displaystyle{ d }[/math],这些集合是由同一公式构造的朱利亚的连通轨迹。 对全三次连通轨迹进行研究,这里考虑双参数递归 [math]\displaystyle{ z \mapsto z^3 + 3kz + c }[/math],两个临界点是参数k的复数平方根。若两个临界点都固定,则说明参数在三次连通轨迹中。 [25] 对于一般的全纯函数族,将曼德布洛特集的分界线推广到分支轨迹,即使在连通轨迹不存在的情况下,分支轨迹也是一个自然的研究对象。
多重曼德布洛特集是通过改变指数[math]\displaystyle{ d }[/math] 的值得到的。其中[math]\displaystyle{ d≥2 }[/math]。指数[math]\displaystyle{ d }[/math]可以进一步推广到负数和小数。 右侧有一视频,其显示了从 d=0到7的发展过程,在这个过程中,周边有6个(即 d-1)凸起部分。类似的负指数发展也会导致(1-d)条裂缝出现在环的内侧。
点击此处观看d从0到7的多重曼德布洛特集视频
更高维下的曼德布洛特集
由于曼德布洛特集无法将复数扩展到三维来进行迭代,故曼德布洛特集不能完美的扩展到三维图形。但四元数 Quaternions的方法可将复数扩展到四维。[26] 其能够将曼德布洛特集和朱利亚集成功扩展到四维,再利用投影或横切成三维。
其他非解析性质的映射
特别有趣的是三角分形 Tricorn fractal ,其反全纯族的连通轨迹为:[math]\displaystyle{ z \mapsto \bar{z}^2 + c. }[/math]
其定义与曼德布洛特集一样,区别在于映射方式。约翰·米尔诺 John Milnor在研究实三次多项式的参数切片时发现了Tricorn Fractal(有时也称为曼德尔巴分形 Mandelbar Fractal))。 它不是局部连接的。 这一性质是由实三次多项式的连通轨迹继承过来的。
另一种非解析泛化是燃烧船分形,它是通过迭代以下得到的:[math]\displaystyle{ z \mapsto (|\Re \left(z\right)|+i|\Im \left(z\right)|)^2 + c. }[/math]
绘制曼德布洛特集
主要文章:绘制曼德尔勃特集合的算法
可以利用多种算法通过计算设备来绘制曼德布洛特集。在这里,展示应用最广泛且最简单的算法:即朴素的“逃逸时间算法”。在逃逸算法中,对于绘图区域中每一个(x,y)点重复计算,根据计算后的不同表现为某一像素区域选择不同的颜色。每个点(x,y)用作重复计算或迭代计算中的起始值(下面将详细描述)。 每次迭代的结果用作下一次迭代的起始值。在每次迭代期间都会检测这些值,以查看它们是否已经达到了关键的“逃逸”条件或“跳出”条件。 如果达到该条件,则停止计算,选择像素的颜色并绘制,并检测下一个(x,y)点。每个点的颜色值反应了逃逸点处的逃逸速度。通常,黑色用于显示在迭代限制之前未能逃逸的值,而逐渐逃逸的点使用较亮的颜色表示。这很直观的表示了在达到逃逸条件之前需要多少个循环。
为了渲染这样的图像,考虑将复平面区域分割成一定数量的点区域。要为任意点区域着色,设置[math]\displaystyle{ c }[/math]作为其点区域的中点。在[math]\displaystyle{ P_{c} }[/math]的前提下迭代临界点0,在每一步检查轨道点的模量是否大于2。(即[math]\displaystyle{ x^2+y^2 }[/math]与[math]\displaystyle{ 2^2 }[/math]的关系)在这种情况下,我们知道[math]\displaystyle{ c }[/math]不属于曼德布洛特集。根据判断模量是否大于2所用的迭代次数为点区域着色。否则,我们将继续进行迭代直到遇到逃逸点,然后确定我们的参数“可能”在曼德布洛特集中,或者至少非常接近它,并将点区域涂成黑色。
在虚拟程序代码中,该算法如下所示。该算法不使用复数,对于那些没有复数类型的数据,使用两个实数手动对复数进行模拟。如果编程语言包含复杂数据类型的操作,则程序可以简化。
for each pixel (Px, Py) on the screen do
x0 = scaled x coordinate of pixel (scaled to lie in the Mandelbrot X scale (-2.5, 1))
y0 = scaled y coordinate of pixel (scaled to lie in the Mandelbrot Y scale (-1, 1))
x := 0.0
y := 0.0
iteration := 0
max_iteration := 1000
while (x×x + y×y ≤ 2×2 AND iteration < max_iteration) do
xtemp := x×x - y×y + x0
y := 2×x×y + y0
x := xtemp
iteration := iteration + 1
color := palette[iteration]
plot(Px, Py, color)
其中,虚拟程序代码中的[math]\displaystyle{ c }[/math], [math]\displaystyle{ z }[/math] 和[math]\displaystyle{ P_c }[/math]:
- [math]\displaystyle{ z = x + iy\ }[/math]
- [math]\displaystyle{ z^2 = x^2 +i2xy - y^2\ }[/math]
- [math]\displaystyle{ c = x_0 + i y_0\ }[/math]
因此,在计算“x”和“y”的虚拟程序代码中可以看出:
- [math]\displaystyle{ x = \mathop{\mathrm{Re}}(z^2+c) = x^2-y^2 + x_0 }[/math]
- [math]\displaystyle{ y = \mathop{\mathrm{Im}}(z^2+c) = 2xy + y_0 }[/math]
要获得曼德布洛特集的彩色图像,可以使用各种函数(线性、指数等)中的一个来为执行程序所需不同的迭代次数所对应的点区域分配颜色。
大众文化中的引用
曼德布洛特集合被许多人认为是最流行的分形[27][28],并在流行文化中被多次提及。
- 乔纳森·库尔顿 Jonathan Coulton的歌曲《曼德尔布洛特集》既是对分形本身的赞颂,也是对它的发现者 Benoît B. Mandelbrot的赞颂。[29]
- 皮尔斯·安东尼 Piers Anthony的《时尚》系列的第二本书《分形模式》描述了一个完美的三维模型世界。[30]
- 亚瑟·查理斯·克拉克 Piers Anthony的小说《来自大浅滩的幽灵》描绘了一个曼德布洛特集图形复刻版的人工湖。[31]
参阅
- 知乎 燃烧着船分形
- Buddhabrot 佛布罗特分形图案
- Collatz fractal 卡拉兹分形
- Fractint 命名由分形fractal和整数integer组合而来,是一种能够渲染和显示分形图案的计算机程序
- Gilbreath permutation 吉尔布雷斯排列
- Mandelbox 曼德尔箱
- Mandelbulb 曼德尔球
- Menger Sponge 门格海绵
- Newton fractal 牛顿分形
- Orbit portrait 轨道特征
- Orbit trap 轨道井
- Pickover stalk Pickover茎
参考文献
- ↑ 1.0 1.1 Adrien Douady and John H. Hubbard, Etude dynamique des polynômes complexes, Prépublications mathémathiques d'Orsay 2/4 (1984 / 1985)
- ↑ Robert Brooks and Peter Matelski, The dynamics of 2-generator subgroups of PSL(2,C), in Irwin Kra (1 May 1981). Irwin Kra. ed. Riemann Surfaces and Related Topics: Proceedings of the 1978 Stony Brook Conference. [Bernard Maskit]. Princeton University Press. ISBN 0-691-08267-7. http://www.math.harvard.edu/archive/118r_spring_05/docs/brooksmatelski.pdf.
- ↑ R.P. Taylor & J.C. Sprott (2008). "Biophilic Fractals and the Visual Journey of Organic Screen-savers" (PDF). Nonlinear Dynamics, Psychology, and Life Sciences, Vol. 12, No. 1. Society for Chaos Theory in Psychology & Life Sciences. Retrieved 1 January 2009.
- ↑ Benoit Mandelbrot, Fractal aspects of the iteration of [math]\displaystyle{ z\mapsto\lambda z(1-z) }[/math] for complex [math]\displaystyle{ \lambda, z }[/math], Annals of the New York Academy of Sciences 357, 249/259
- ↑ Peitgen, Heinz-Otto; Richter Peter (1986). The Beauty of Fractals. Heidelberg: Springer-Verlag. ISBN 0-387-15851-0.
- ↑ Frontiers of Chaos, Exhibition of the Goethe-Institut by H.O. Peitgen, P. Richter, H. Jürgens, M. Prüfer, D.Saupe. Since 1985 shown in over 40 countries.
- ↑ Gleick, James (1987). Chaos: Making a New Science. London: Cardinal. pp. 229.
- ↑ Dewdney, A. K. (1985). "Computer Recreations, August 1985; A computer microscope zooms in for a look at the most complex object in mathematics" (PDF). Scientific American.
- ↑ John Briggs (1992). Fractals: The Patterns of Chaos. p. 80.
- ↑ Pountain, Dick (September 1986). "Turbocharging Mandelbrot". Byte. Retrieved 11 November 2015.
- ↑ Lyubich, Mikhail (May–June 1999). "Six Lectures on Real and Complex Dynamics". Retrieved 2007-04-04.
{{cite journal}}
: Cite journal requires|journal=
(help) - ↑ Lyubich, Mikhail (November 1998). "Regular and stochastic dynamics in the real quadratic family" (PDF). Proceedings of the National Academy of Sciences of the United States of America. 95 (24): 14025–14027. Bibcode:1998PNAS...9514025L. doi:10.1073/pnas.95.24.14025. PMC 24319. PMID 9826646. Retrieved 2007-04-04.
- ↑ "Mandelbrot Set Explorer: Mathematical Glossary". Retrieved 2007-10-07.
- ↑ Kahn, Jeremy (8 August 2001). "The Mandelbrot Set is Connected: a Topological Proof" (PDF).
- ↑ The Mandelbrot set, theme and variations. Tan, Lei. Cambridge University Press, 2000. ISBN:978-0-521-77476-5. Section 2.1, "Yoccoz para-puzzles", p. 121
- ↑ Exploring the Mandelbrot set. The Orsay Notes by Adrien Douady and John H. Hubbard. page 12
- ↑ Wolf Jung, March 2002, Homeomorphisms on Edges of the Mandelbrot Set by Wolf Jung
- ↑ 18.0 18.1 Hubbard, J. H. (1993), "Local connectivity of Julia sets and bifurcation loci: three theorems of J.-C. Yoccoz" (PDF), Topological methods in modern mathematics (Stony Brook, NY, 1991), Houston, TX: Publish or Perish, pp. 467–511, MR 1215974. Hubbard cites as his source a 1989 unpublished manuscript of Yoccoz.
- ↑ Lei (1990). "Similarity between the Mandelbrot set and Julia Sets". Communications in Mathematical Physics. 134 (3): 587–617. Bibcode:1990CMaPh.134..587L. doi:10.1007/bf02098448.
- ↑ J. Milnor (1989). "Self-Similarity and Hairiness in the Mandelbrot Set". In M. C. Tangora. Computers in Geometry and Topology. New York: Taylor & Francis. pp. 211–257. https://books.google.com/books?id=wuVJAQAAIAAJ.)
- ↑ 21.0 21.1 Shishikura, Mitsuhiro (1998), "The Hausdorff dimension of the boundary of the Mandelbrot set and Julia sets", Annals of Mathematics, Second Series, 147 (2): 225–267, arXiv:math.DS/9201282, doi:10.2307/121009, JSTOR 121009, MR 1626737.
- ↑ Gary William Flake, The Computational Beauty of Nature, 1998. p. 125. ISBN:978-0-262-56127-3.
- ↑ Peitgen, Heinz-Otto; Saupe Dietmar (1988). The Science of Fractal Images. New York: Springer-Verlag. pp. 121, 196–197. ISBN 0-387-96608-0.
- ↑ Peitgen, Heinz-Otto; Richter Peter (1986). The Beauty of Fractals. Heidelberg: Springer-Verlag. pp. 85. ISBN 0-387-15851-0.
- ↑ Rudy Rucker's discussion of the CCM: CS.sjsu.edu
- ↑ http://archive.bridgesmathart.org/2010/bridges2010-247.pdf retrieved 19 August 2018
- ↑ Mandelbaum, Ryan F. (2018). "This Trippy Music Video Is Made of 3D Fractals." Retrieved 17 January 2019
- ↑ Moeller, Olga de. (2018)."what are Fratals?" Retrieved 17 January 2019.
- ↑ "Mandelbrot Set". JoCopeda. Retrieved 15 January 2015.
- ↑ Piers Anthony (1992). Fractal Mode. HarperCollins. ISBN 978-0-246-13902-3. https://books.google.com/books?id=XdUyAAAACAAJ.
- ↑ Arthur C. Clarke (29 September 2011). The Ghost From The Grand Banks. Orion. ISBN 978-0-575-12179-9. https://books.google.com/books?id=6ELsYigmXNoC.
进一步阅读
- John W. Milnor, Dynamics in One Complex Variable (Third Edition), Annals of Mathematics Studies 160, (Princeton University Press, 2006),ISBN 0-691-12488-4
- 约翰米尔诺,《一个复杂变量中的动力学》(第三版) ,数学纪事研究160,(普林斯顿大学出版社,2006) (首次出现于1990年,作为斯托尼布鲁克 IMS 预印本,可作为arXiV: math.DS / 9201272)
- Nigel Lesmoir-Gordon, The Colours of Infinity: The Beauty, The Power and the Sense of Fractals,ISBN 1-904555-05-5
- 奈杰尔·莱斯莫尔-戈登,《无限的颜色: 分形的美、力量和感觉》 ,包括一张 DVD,主角是亚瑟·查理斯·克拉克 Arthur C. Clarke 和大卫·吉尔摩 David Gilmour。
- Heinz-Otto Peitgen, Hartmut Jürgens, Dietmar Saupe, Chaos and Fractals: New Frontiers of Science (Springer, New York, 1992, 2004),ISBN 0-387-20229-3
- 海因茨-奥托·佩特根,哈特穆特·尤尔根斯,迪特马尔·索普,混沌与分形:科学的新前沿(施普林格,纽约,1992,2004),