FractalPark
牛顿进阶

z³ − 1 牛顿分形

对 z³ − 1 使用牛顿法:平面分成三个收敛盆地,中间隔着一条永远细分下去的边界。

在探索器中打开
z³ − 1 的牛顿分形:三大盆地沿轮辐线相接,边界布满螺旋细节
FractalPark 渲染的 z³ − 1 牛顿分形:三个收敛盆地呈三分对称,由层层嵌套的螺旋边界相连。

概览

这个公式把熟悉的求根算法变成复平面的地图。它不像曼德布罗特集那样问轨道是否逃逸,而是问每个初值会被牛顿法带到 z³ − 1 的哪个根。三个单位立方根(1、e^{2πi/3}、e^{4πi/3})各有一个吸引盆:所有最后收敛到该根的初值。

三个盆地一眼能分开,共用边界却不简单。恰好落在边界上的点不会收敛到任何根,而会沿混沌、非重复的轨道走下去。这条边界就是牛顿映射——有理函数 N(z) = z − (z³ − 1)/(3z²)——的 Julia 集;放大后,新的细节还会继续冒出来。

这是收敛分形的经典例子:日常数值算法里,藏着一条不安分的边界。Arthur Cayley 在 1879 年提出牛顿盆地的一般分类问题,二次情形他解开了,三次情形却“似乎存在相当大的困难”。Fatou 与 Julia 的复动力系统理论让问题渐渐可说清,第一批计算机图像则到 20 世纪 80 年代初才出现。

数学原理

牛顿迭代

z(n+1) = z(n) - (z(n)^3 - 1) / (3 z(n)^2)

函数 f(z) 的牛顿迭代式为 N(z) = z − f(z) / f′(z)。这里 f(z) = z³ − 1,f′(z) = 3z²。每一步用多项式值除以导数,修正 z,把它推向某个根。对于 z³ − 1,迭代式可化简为:

N(z) = z − (z³ − 1) / (3z²) = (2z³ + 1) / (3z²)。

这是一个在黎曼球面上的二次有理映射。它的不动点恰好是三个单位立方根,且都是超吸引的——一旦轨道足够接近某个根,收敛就会极快。Julia 集就是这三个盆地相遇的边界,那里呈现混沌的动力学。

试一个具体初值:z₀ = 0.5。z₁ = (2·0.125 + 1) / (3·0.25) = (0.25 + 1) / 0.75 ≈ 1.667,z₂ ≈ 1.107,z₃ ≈ 1.001,轨道迅速收敛到根 1。再取 z₀ = 0.5i:z₁ ≈ −0.417 + 0.583i,z₂ ≈ −0.495 + 0.865i,轨道收敛到 e^{2πi/3} = −0.5 + 0.866i。再取 z₀ = −0.5 − 0.5i,轨道收敛到另一个根 e^{4πi/3} = −0.5 − 0.866i。三个盆地面积庞大、分界清晰,但在边界附近,初值轻微一动,最终的收敛方向就可能完全改变。

FractalPark 使用收敛判定阈值:当相邻两次迭代值的差小于给定容限时,即认为轨道已收敛。随后按收敛到哪个根以及所需步数染色。

历史

本页公式的核心迭代是牛顿法,由 Isaac Newton 于 1669 年发现,后由 Joseph Raphson 于 1690 年独立改进。对实函数而言它是实用的求根工具;对复多项式而言,它揭示了当时无人能预见的分形盆地。

对这些盆地的数学研究始于 1879 年 Arthur Cayley 的论文"The Newton–Fourier imaginary problem"(《American Journal of Mathematics》,第 2 卷,第 97 页)。对于复二次多项式,Cayley 证明吸引盆就是两个根连线的中垂线划分出的两个半平面。可当他转向三次多项式时,却发现该问题"似乎存在相当大的困难"。他未能将简洁的二次结论推广到三次情形,这一遗留问题后被称为 Cayley 问题。

1918 年前后,Pierre Fatou 与 Gaston Julia 发展出的复动力系统理论提供了描述盆地及其无限复杂边界所需的数学语言。但这些边界的第一批实际图像——z³ − 1 牛顿分形的图像——要等到 20 世纪 80 年代初计算机图形技术的出现。

Heinz-Otto Peitgen 与 Peter H. Richter 在 1986 年出版的《The Beauty of Fractals: Images of Complex Dynamical Systems》中收入专章"Newton's Method for Complex Polynomials: Cayley's Problem"(第 6 章,第 93–106 页)。它收录了其中最早一批公开发表的计算机渲染图,为广泛读者展示了隐藏在牛顿法中的分形结构。

视觉特征

三个大盆地各对应一个单位立方根,占据画面。它们沿中心放射的轮辐线相接,形成三分对称。盆地内部颜色平滑,收敛又快又稳。

任意两盆地的边界,尤其中心三岔处,挤满嵌套螺旋、三个盆地的微型副本和树突般细丝。沿螺旋往里看,三盆地结构会在更小尺度重现,只是被局部动力学拉扯变形。大块平静,边缘却一直不安。

再创作与示例

从 canonical Document 开始

打开本指南使用的同一份已审定公式状态,再在交互式探索器中调整视图、着色、变换或动画。

常见问题

主要色块分别表示什么?

每个有色区域都称为吸引盆,包含牛顿迭代最终收敛到 z³ − 1 同一个根的所有初值。三个单位立方根分别是 1、e^{2πi/3} = −0.5 + 0.866i 和 e^{4πi/3} = −0.5 − 0.866i。点的颜色记录它抵达哪个根,以及收敛有多快。

这是逃逸时间分形吗?

它主要是收敛分形。分类依据不是轨道是否逃向无穷(如曼德布罗特集那样),而是轨道趋近哪个根以及需要多少步。这里没有逃逸半径判断。渲染器检查相邻两次迭代值的差是否小于阈值——若已足够接近,便判定收敛并确认收敛到哪个根。

参考资料