inRm3D: 画板论坛 - Powered by Discuz! Board

标题: UF分形图征解（phoenix） [打印本页]

作者: 柳烟 时间: 2010-7-7 00:22 标题: UF分形图征解（phoenix）

这是代码：
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
}

图片附件: Fractal2.png (2010-7-7 00:22, 50.16 KB) / 下载次数 1833
http://inrm3d.cn/attachment.php?aid=5915&k=a4ef11df6db878b7c8ccf96d09063019&t=1732395683&sid=9S3j3N

作者: 柳烟 时间: 2010-7-7 12:48

此分形图，好象榕坚兄作过，最近不知榕坚兄那里去了。

作者: inRm 时间: 2010-7-7 12:57

榕坚已十多天没来了。

作者: changxde 时间: 2010-7-10 13:28

今天试了一下这个，外形出来了，还需要修改。

图片附件: 未命名.JPG (2010-7-10 13:28, 14.95 KB) / 下载次数 1684
http://inrm3d.cn/attachment.php?aid=5989&k=4dcdf8eec86e27ec44e7e46a64c2fb31&t=1732395683&sid=9S3j3N

图片附件: 未命名1.JPG (2010-7-10 13:28, 19.29 KB) / 下载次数 1666
http://inrm3d.cn/attachment.php?aid=5990&k=a55222397af7a30a2ebac786b71a50f1&t=1732395683&sid=9S3j3N

附件: UF-phoenix.gsp (2010-7-10 13:28, 26.19 KB) / 下载次数 2046
http://inrm3d.cn/attachment.php?aid=5991&k=6ca9078ee44c67a78e9c9d52d03abf56&t=1732395683&sid=9S3j3N

作者: 柳烟 时间: 2010-7-10 17:33

changxde作得不错.想问两个问题:看了你的文件,按条件,迭代式应取代码中的z = z^2 + real(p1) + imag(p1) * y,其中y的初始值应为原点O,以后的y应为上一个z值,不知我这理解对不对?其次,你文件中的z'为自由点,z那里去了,若合并,合并到那里了,迭代时是如何迭代的?
我仍整来不对头.烦changxde老师指点一下,这分形程序就那么几句,却花了我不少时间,也浪费了不少时间.

作者: changxde 时间: 2010-7-10 20:30

这是一个二阶迭代关系，初值有两个，第一个合并到原点，对第二个着色。

作者: 柳烟 时间: 2010-7-10 21:28

这是一个二阶迭代关系，初值有两个，第一个合并到原点，对第二个着色。
changxde 发表于 2010-7-10 20:30

谢谢,懂了一些,但是仍糊涂.首先,请板友看看我楼上的理解是否有误,其次,z^2+P实+P虚.y中,z好象是你源文件中的z',文件中对这个进行了着色,是第二个初始点究竟如何处理?

作者: changxde 时间: 2010-7-15 18:35

作者: xuefeiyang 时间: 2010-7-15 19:00

http://www.inrm3d.cn/viewthread.php?tid=1448&page=1#pid13750
这个动画有意思。

作者: math 时间: 2010-7-15 20:48

其实这个分形与M集的J集是有关系的，把参数改一下：

图片附件: snap1.JPG (2010-7-15 20:48, 25.18 KB) / 下载次数 1534
http://inrm3d.cn/attachment.php?aid=6083&k=a0f76fbf75f366eb5b518f32af900d86&t=1732395683&sid=9S3j3N

欢迎光临 inRm3D: 画板论坛 (http://inrm3d.cn/)