返回列表 回复 发帖

UF分形图征解(phoenix)

Fractal2.png
2010-7-7 00:22

这是代码:
phoenix {
; Parameters 1 and 2 have been merged into one complex parameter
; because switching would not be possible otherwise.
init:
  complex z = pixel  
  complex y = 0
  complex oldz = 0
loop:
  oldz = z
  if @degree >= 2
    z = z^@degree + real(p1) * z^(@degree - 1) + imag(p1) * y
  elseif @degree <= -3
    z = z^(-@degree) + real(p1) * z^(-@degree - 2) + imag(p1) * y
  else
    z = z^2 + real(p1) + imag(p1) * y
  endif
  y = oldz
bailout:
  (@test == 0 && |z| <= @bailout) ||                                         \
  (@test == 1 && sqr(real(z)) <= @bailout) ||                                \
  (@test == 2 && sqr(imag(z)) <= @bailout) ||                                \
  (@test == 3 && (sqr(real(z)) <= @bailout && sqr(imag(z)) < @bailout)) ||   \
  (@test == 4 && (sqr(real(z)) <= @bailout || sqr(imag(z)) < @bailout)) ||   \
  (@test == 5 && (sqr(abs(real(z)) + abs(imag(z))) <= @bailout)) ||          \
  (@test == 6 && (sqr(real(z) + imag(z)) <= @bailout))
default:
  title = "Phoenix"
  maxiter = 149
  param p1
    caption = "Parameter"
    default = (0.56667, -0.5)
  endparam
  param degree
    caption = "Degree"
    default = 0
  endparam
  param test
    caption = "Bailout Test"
    default = 0
    enum = "mod" "real" "imag" "or" "and" "manh" "manr"
  endparam
  param bailout
    caption = "Bailout value"
    default = 4.0
    min = 1.0
  endparam
switch:
  type = "mandphoenix"
  degree = degree
  test = test
  bailout = bailout
}
此分形图,好象榕坚兄作过,最近不知榕坚兄那里去了。
榕坚已十多天没来了。
今天试了一下这个,外形出来了,还需要修改。
未命名.JPG
2010-7-10 13:28

未命名1.JPG
2010-7-10 13:28

UF-phoenix.gsp (26.19 KB)

changxde作得不错.想问两个问题:看了你的文件,按条件,迭代式应取代码中的z = z^2 + real(p1) + imag(p1) * y,其中y的初始值应为原点O,以后的y应为上一个z值,不知我这理解对不对?其次,你文件中的z'为自由点,z那里去了,若合并,合并到那里了,迭代时是如何迭代的?
我仍整来不对头.烦changxde老师指点一下,这分形程序就那么几句,却花了我不少时间,也浪费了不少时间.
这是一个二阶迭代关系,初值有两个,第一个合并到原点,对第二个着色。
这是一个二阶迭代关系,初值有两个,第一个合并到原点,对第二个着色。
changxde 发表于 2010-7-10 20:30
谢谢,懂了一些,但是仍糊涂.首先,请板友看看我楼上的理解是否有误,其次,z^2+P实+P虚.y中,z好象是你源文件中的z',文件中对这个进行了着色,是第二个初始点究竟如何处理?
其实这个分形与M集的J集是有关系的,把参数改一下:
snap1.JPG
返回列表