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
}

收藏分享评分

回复引用

订阅 TOP

版主

Rank: 7 Rank: 7 Rank: 7

UID: 539
帖子: 87
精华: 0
积分: 338

2^#

柳烟发表于 2010-7-7 12:48 | 只看该作者

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

回复引用

TOP

版主

Rank: 7 Rank: 7 Rank: 7

UID: 3
帖子: 2358
精华: 4
积分: 2068

3^#

inRm发表于 2010-7-7 12:57 | 只看该作者

榕坚已十多天没来了。

回复引用

TOP

院士

Rank: 5 Rank: 5 Rank: 5

UID: 352
帖子: 1321
精华: 14
积分: 2599
来自: 河南永城

4^#

changxde发表于 2010-7-10 13:28 | 只看该作者

今天试了一下这个，外形出来了，还需要修改。
未命名.JPG

UF-phoenix.gsp (26.19 KB)

回复引用

TOP

版主

Rank: 7 Rank: 7 Rank: 7

UID: 539
帖子: 87
精华: 0
积分: 338

5^#

柳烟发表于 2010-7-10 17:33 | 只看该作者

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

回复引用

TOP