在二维平面中,复数可表示为z=x+yi,x和y是实数,i2=-1。自平方分形的变换函数是z2=λz(1-z),λ是一个常量复数。
由z2=λz(1-z)得z=f-1(z2)=1/2(1+sqrt(1-4z2/λ),
所以令二次判别式1-(4z2)/λ=discr,设z=re(z)+im(z),re(z),im(z)分别为实部和虚部。
则:
re(z)=1/2*(1+sqrt((abs(discr)+re(discr))/2),
im(z)=1/2*sqrt((abs(discr)+re(discr))/2。
当λ=1.5+0.8i时的图形如下:
程序如下:
若修改一下自平方函数,令:
z0=c,c为常量复数,
zk=z2k-1+z0;
也即f(z)=z2+c。则可以得到Mandelbrot集。
如图: