Escher_Julia 分形

changxde

下载 (107.8 KB)

2013-12-4 21:07

这是winfract中的一个分形图，名字叫escher_julia，我查了一下解释，就能得到下面的信息，
Escher_julia
Escher-like tiling of Julia sets from The Science of Fractal Images
   z(0) = pixel
  z(n+1) = z(n)^2 + (0, 0i)

The target set is a second, scaled, Julia set:

T = [ z: | (z * 15.0)^2 + c | < BAILOUT ]

Two parameters: real and imaginary parts of c
Iteration count and bailout size apply to both Julia sets.

c=0.32+0.043i


109#由给出了画板做法

http://www.inrm3d.cn/viewthread.php?tid=4591&page=14#pid43172
给出较干净的算法，调整了榕老师算法的小问题。

changxde

http://pan.baidu.com/s/1lDph7
我也忘了在那下的了。我把文件夹上传到网盘，看看能不能打开。

changxde

下载 (69.17 KB)

2013-12-5 22:12

c=-0.74+0.1i

changxde

我思考了一下，此分形的算法应该是：

第一步：z1=z^2
第二步：z= 15z1
第三步：z=z^2+c 迭代计算Julia集，
第四步：若z在Julia集内部，计算结束；若z在Julia集外部，z=z1回到第一步。

但还不知如何实现此算法。

changxde

UF试验

Escher {

init:
   z = #pixel
   ;z1=0
loop:
    z1 = z^2
    int k=0
    z = 15*z1

    while (k<@degree &&  |z| < @bailout)
        z = z^2 + @c
        k = k+1
    endwhile

    if |z| >= @bailout
       z = z1
    endif

bailout:
    |z| < @bailout

default:
  title = "Escher"
  int param degree
    caption="degree n"
    default=60
    min=1
  endparam
  complex param c
    default=(0.32,0.043)
  endparam
  float param bailout
    caption = "Bailout value"
    min=0
    default=4
  endparam
}

下载 (87.78 KB)

2014-7-9 11:06

选一个内部着色方法。

changxde

我也不能上传图片

changxde

31# 柳烟
上下对称，y加绝对值了。

changxde

这样做只是上面的几层，要体现极限的味道，须要增加层数。

changxde

简单修改carr199
Carr1997 {
init:
  z= 0
  c=pixel, c1=1.03125*c, c2=1.0625*c
  c3=1.125*c, c4=1.25*c , c5=1.5*c,  c6=2*c
  l1=round(real(p1)), l2=round(imag(p1))
  l3=round(real(p2)), l4=round(imag(p2))
  l5=round(real(p3)), l6=round(imag(p3))
  float bailout=16, int iter=0  ,int k = 256
loop:
  IF (iter==l1)
     c = c1 ,  k = 128
  ELSEIF (iter==l2)
     c = c2 ,  k = 64
  ELSEIF (iter==l3)
     c = c3 , k = 32
  ELSEIF (iter==l4)
     c = c4 , k = 16
  ELSEIF (iter==l5)
     c = c5 , k = 8
  ELSEIF (iter==l6)
     c = c6 , k = 4
  ENDIF
  z=z*z+c^k
  iter=iter+1
bailout:
  |real(z)| <= bailout
http://c.hiphotos.bdimg.com/album/s%3D1000%3Bq%3D90/sign=8102e445a38b87d65442af1f37381349/95eef01f3a292df5d2e0fa61be315c6034a873a6.jpg

最近网络不正常，老是不能上传图片。

changxde

http://c.hiphotos.bdimg.com/album/s%3D1000%3Bq%3D90/sign=e4902e72aeaf2eddd0f14de9bd203a98/79f0f736afc3793198ccd65de9c4b74543a91119.jpg

可以用右键目标另存观看