Mandelbrot集的最新变化形态一览——MandelBox，Mandelbulb，Burning Ship，NebulaBrot

二维Mandelbrot集——Burning Ship

采用以下迭代公式

(x⁴-6*x²*y²+y⁴,4*|x|³*|y|-4*|y|³*|x|)

看不懂的，可以去学习深造了。。。

 

Mandelbulb

这个3D的Mandelbrot集采用的是以下公式，这应该算是超复数的一种，人称“triplex”，三元复数

 

数学上的N次方

其中：

一般情况下，n取8。

简单说来，三元复数的平方的计算机表示为

xx = (x*x+y*y) * (1-z*z/(x*x+y*y));

yy = 2*x*y * (1-z*z/(x*x+y*y));

zz = -2*z*sqrt(x*x+y*y);

详情请参见http://en.wikipedia.org/wiki/Mandelbulb

以下内容摘自http://www.skytopia.com/project/fractal/2mandelbulb.html

{x,y,z}^n = r^n { sin(theta*n) * cos(phi*n) ,sin(theta*n) * sin(phi*n) , cos(theta*n) }
...where:
r = sqrt(x^2 + y^2 + z^2)
theta = atan2( sqrt(x^2+y^2), z )
phi = atan2(y,x)

And the addition term in z -> z^n + c is similar tostandard complex addition, and is simply defined by:

{x,y,z}+{a,b,c} = {x+a, y+b, z+c}

The rest of the algorithm is similar to the 2D Mandelbrot!

Here is some pseudo code of the above:

r = sqrt(x*x + y*y + z*z )
theta = atan2(sqrt(x*x + y*y) , z)
phi = atan2(y,x)

newx = r^n * sin(theta*n) * cos(phi*n)
newy = r^n * sin(theta*n) * sin(phi*n)
newz = r^n * cos(theta*n)

...where n is the order of the 3D Mandelbulb. Use n=8to find the exact object in this article.

 

上图是16阶的MandelbrotBulb

 

放大之后的某部分

 

 

取n=8时候的图形

 

 

迭代次数取100时候的情形

“金色大峡谷”

 

“神秘洞穴”

“魔法球花甘蓝”

“Mandelbrot花园”

八阶的MandelBulb

圣诞节珊瑚球

冰激凌~~

贝壳~~

“冰封地狱”

 

 

 

 

 

 

 

 

 

 

 

3D的Mandelbrot集——Mandelbox

一个很奇妙的视频演示。“飞越曼德布罗盒子”

http://wimp.com/mandelboxflythrough/

具体算法核心代码：

for (each axis)

  if (v[axis]>1) v[axis] =2-v[axis];

  else if (v[axis]<-1) v[axis] =-2-v[axis];

if (v.magnitude() < 0.5) v *= 4;

else if (v.magnitude() < 1) v /= square(v.magnitude());

v = scale*v + c;

 

 

 

 

详情参见“超复数分形”

http://www.bugman123.com/Hypercomplex/index.html

 

n = 100; norm[x_] := x.x;

TriplexPow[{x_, y_, z_}, n_] := If[x == y == 0.0, 0.0, Module[{r =Sqrt[x^2 + y^2 + z^2], theta = n ArcTan[x, y], phi}, phi = n ArcSin[z/r];

r^n{Cos[theta]Cos[phi], Sin[theta]Cos[phi], -Sin[phi]}]];
Mandelbulb[c_] := Module[{p = {0, 0, 0}, i = 0}, While[i < 24 &&norm[p] < 4, p = TriplexPow[p, 8] + c; i++]; i];
image = Table[z = 1.1; While[z >= -0.1 && Mandelbulb[{x, y, z}] <24, z -= 2.2/n];

z, {y, -1.1, 1.1, 2.2/n}, {x, -1.1, 1.1, 2.2/n}];
ListDensityPlot[image, Mesh -> False, Frame -> False, PlotRange ->{-0.1, 1.1}]

 

 

上图是Lambdabulb

 

上图所用迭代公式：{x,y,z}² = {x²-y²-z²,2xy, -2xz}

 

 

上图是四元数Mandelbrot集，所用迭代公式：{x,y,z,w}²= {x²-y²-z²-w², 2xy, 2xz, 2xw}

 

 

上图名叫Glynn Julia set

 

 

 

上面两幅图是4D Bicomplex Mandelbrot集("Tetrabrot")

所用迭代公式：{x,y,z,w}² = {x²-y²-z²+w²,2(xy-zw), 2(xz-yw), 2(xw+yz)}

 

  

上面三幅图叫做NebulaBrot，是BuddhaBrot的三维形式

这些图均由20亿个点进行计算渲染而成

 

 

上图叫做3DChristmas Tree Mandelbrot Set，使用的迭代公式：

{x,y,z}ⁿ = rⁿ{cos(θ)cos(φ),sin(θ)cos(φ), sin(φ)}
r=sqrt(x²+y²+z²), θ=n atan2(y,x), φ=natan2(z,x)

 

 

 

{x,y,z}ⁿ = rⁿ{cos(θ)cos(φ),sin(θ)cos(φ), cos(θ)sin(φ)}
r=sqrt(x²+y²+z²), θ=n atan2(y,x), φ=n sin^-1(z/r)

 

 

 

 

 

上图叫做"Roundy" Mandelbrot Set 。迭代公式：

{x,y,z,w}² = {x²-y²-z²-w²,2(xy+zw), 2(xz+yw), 2(xw+yz)}

 

 

上图是Bristorbrot Set

所用迭代公式：{x,y,z}² = {x²-y²-z²,y(2x-z),z(2x+y)}

上图所用迭代公式：{x,y,z}² = {x²-y²-z²,2xy, 2(x-y)z}

作者DavidMakin