复动力系统 · z → z² + c

朱利亚集

固定c,反复做z → z² + c:轨道保持有界的起点组成填充朱利亚集,它的边界是朱利亚集。

朱利亚集与本页画面

对固定的复数cc,填充朱利亚集KcK_c是让下面这个数列保持有界的起点z0z_0组成的集合;一旦∣zn∣>max⁡(2,∣c∣)|z_n| > \max(2, |c|),数列就会跑到无穷远。朱利亚集JcJ_c是KcK_c的边界:它的每一点附近,既有逃逸的起点,也有不逃逸的起点。

zn+1=zn2+c,Kc={ z0∈C:∣zn∣≤max⁡(2,∣c∣)  ∀n≥0 },Jc=∂Kcz_{n+1} = z_n^2 + c,\qquad K_c = \{\, z_0 \in \mathbb{C} : |z_n| \le \max(2, |c|) \ \ \forall n \ge 0 \,\},\qquad J_c = \partial K_c

两个画面上是什么

大图是z平面:每一点是一个起点z0z_0,颜色表示它的轨道多快逃逸。平滑的颜色只是显示效果。黑色表示“迭代上限NN步内没有逃逸”——这是数值判定。尘埃那几个视角和c=ic = i的树枝面积为0,那里几乎每个像素都会逃逸,所以按距离估计离KcK_c不到半个像素的点也涂深,好让人看见它:这同样只是显示加工。

面板里的小屏是c平面:曼德博集合,红色十字标着本页的cc。拖动十字或点一下,就换一个cc。大图里黄色的点是0的轨道(或者你点的那个起点的轨道)。

为什么是这个样子

红键逐级画出满足∣zn∣≤R|z_n| \le R的起点,R=max⁡(2,∣c∣)R = \max(2, |c|),n=1,2,3,…n = 1, 2, 3, \ldots:每一层都在上一层里面,一层层收紧到KcK_c。每一层是上一层在z↦z2+cz \mapsto z^2 + c下的原像;圆盘的原像连成一片,当且仅当圆盘包含临界值cc。所以只要0的轨道还在半径RR以内,各层都连成一片;从它跑出去的那一步起,区域分成2块、4块、8块……

这就是Fatou和Julia关于z2+cz^2 + c的定理 [1, 2]:0的轨道有界时KcK_c连通;逃逸时KcK_c是完全不连通的“尘埃”(康托集)。曼德博集合正是让KcK_c连通的那些cc。

吸引循环一定吸引0的轨道,所以z2+cz^2 + c至多有一个吸引循环:KcK_c黑色的内部都流向它。在尖点c=1/4c = 1/4和c=−3/4c = -3/4,循环是抛物的;在主心形边界上黄金分割角的那一点,不动点周围是西格尔圆盘,轨道在里面永远转圈 [3]。这几种情形都没有吸引循环;默认设置下,页面上的周期读数是“未判定”(NN很大时,数值检测可能把缓慢的趋近误认成循环)。

容易误读的地方

来历

1918–1920年,Gaston Julia和Pierre Fatou在还没有计算机图像的年代建立了有理函数迭代的理论 [1, 2]。1942年,Carl Ludwig Siegel证明了有些中性不动点周围有一圈轨道永远转圈的圆盘 [3]。1980年以后的计算机图像和Douady、Hubbard的工作,让这些集合广为人知 [4]。

相关概念

曼德博集合 · 牛顿分形(页面筹备中)· 复数

参考文献

  1. G. Julia, “Mémoire sur l’itération des fonctions rationnelles”, Journal de Mathématiques Pures et Appliquées, série 8, 1 (1918) 47–245.
  2. P. Fatou, “Sur les équations fonctionnelles”, Bulletin de la Société Mathématique de France 47 (1919) 161–271; 48 (1920) 33–94, 208–314.
  3. C. L. Siegel, “Iteration of analytic functions”, Annals of Mathematics 43 (1942) 607–612. doi:10.2307/1968952
  4. J. Milnor, Dynamics in One Complex Variable, 3rd ed., Annals of Mathematics Studies 160, Princeton University Press (2006).

延伸阅读:维基百科:朱利亚集合;MacTutor数学史:Gaston Julia(英文);MacTutor数学史:Pierre Fatou(英文)。