此篇文章是搬运自我自己的csdn博客。 2010年写的
其中因为图片失效了，幸好Evan_Gu在2015年转载过我的文章，图就用他的了
注意，我是原作者我，我有权转载
二维Mandelbrot集——Burning Ship
采用以下迭代公式
(x4-6x2y2+y4,4*|x|3*|y|-4*|y|3*|x|)
看不懂的，可以去学习深造了。。。\
Mandelbulb
这个3D的Mandelbrot集采用的是以下公式，这应该算是超复数的一种，人称“triplex”，三元复数
数学上的N次方
其中：
一般情况下，n取8。
简单说来，三元复数的平方的计算机表示为
xx = (xx+yy) * (1-zz/(xx+y*y));
yy = 2xy * (1-zz/(xx+y*y));
zz = -2zsqrt(xx+yy);
详情请参见en.wikipedia.org/wiki/Mandel…
以下内容摘自www.skytopia.com/project/fra…
{x,y,z}^n = r^n { sin(thetan) * cos(phin) ,sin(thetan) * sin(phin) , cos(thetan) }
...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(xx + yy + zz )
theta = atan2(sqrt(xx + yy) , z)
phi = atan2(y,x)
newx = r^n * sin(thetan) * cos(phin)
newy = r^n * sin(thetan) * sin(phin)
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
一个很奇妙的视频演示。“飞越曼德布罗盒子”
具体算法核心代码：
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;
详情参见“超复数分形”
www.bugman123.com/Hypercomple…
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}2 = {x2-y2-z2,2xy, -2xz}
上图是四元数Mandelbrot集，所用迭代公式： {x,y,z,w}2= {x2-y2-z2-w2, 2xy, 2xz, 2xw}
上图名叫Glynn Julia set
上面两幅图是4D Bicomplex Mandelbrot****集("Tetrabrot")
所用迭代公式：{x,y,z,w}2 = {x2-y2-z2+w2,2(xy-zw), 2(xz-yw), 2(xw+yz)}
** **
上面三幅图叫做NebulaBrot，是BuddhaBrot的三维形式
这些图均由20亿个点进行计算渲染而成
上图叫做3DChristmas Tree Mandelbrot Set，使用的迭代公式：
{x,y,z}n = rn{cos(θ)cos(φ),sin(θ)cos(φ), sin(φ)}
r=sqrt(x2+y2+z2), θ=n atan2(y,x), φ=natan2(z,x)
{x,y,z}n = rn{cos(θ)cos(φ),sin(θ)cos(φ), cos(θ)sin(φ)}
r=sqrt(x2+y2+z2), θ=n atan2(y,x), φ=n sin-1(z/r)
上图叫做 "Roundy" Mandelbrot Set 。迭代公式：
{x,y,z,w}2 = {x2-y2-z2-w2,2(xy+zw), 2(xz+yw), 2(xw+yz)}
所用迭代公式：{x,y,z}2 = {x2-y2-z2,y(2x-z),z(2x+y)}
**
**
上图所用迭代公式：{x,y,z}2 = {x2-y2-z2,2xy, 2(x-y)z}
