牛顿迭代
函数 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 使用收敛判定阈值:当相邻两次迭代值的差小于给定容限时,即认为轨道已收敛。随后按收敛到哪个根以及所需步数染色。

