xiaongxp

“在同一逃逸区内的点离边界的远近是否能用什么加以区别。”是指着色吗？
一般来讲，每个逃逸区内的点按势值（等高线）具有相同的遍历性，所以我作的色因子常带0.5(et±log(abs(em))) ，以0.5et区别不同逃逸区，以±0.5log(abs(em))) 来反映同一逃逸区的明暗过渡效果。这种做法对经典M集和经典J集很有效。但对其他复分形，要协调好两者的系数才能达到好的视觉效果，这是一个需要耐心、充满期待、成后无比欣慰的过程。我们都非常享受这一过程，不痴迷复分形的人是不能理解个中乐趣、享受到这种快乐的。

榕坚

20# xiaongxp


我之前也是采取类似做法，但没办法耐心扫描就放弃了。因为它要设置较慢的速度进行扫描才会使那些线相对连续。等过一阵子再来研究一下UF的做法。谢谢向老师。

榕坚

这个奇怪的UF中的范例，我做了两遍都与结果相去甚远：
Critper3 {
; Generic Mandelbrot set
init:
  complex s = #pixel
  ;s=1-1/s
  ;s=1-1/s
  complex a3 = 0.3333333333333333333
  complex pp = (1/(s - s*s) - s)*3/s
  complex a0 = sqrt(pp)
  complex a2 = a0*((s - 1)/pp - a3)
  z = -2*a2
loop:
  z = (a3*z + a2)*sqr(z) + a0
bailout:
  |z| <= @bailout
default:
  title = "Critper3"
  center = (0, 0)
  param start
    caption = "Starting point"
    default = (1,0)
    hint = "Perturbation. Use (1,0) for the standard Mandelbrot set."
  endparam
  param power
    caption = "Power"
    default = (1,0)
    hint = "This defines the power of the Mandelbrot set. Use (1,0) \
            for the standard Mandelbrot set."
  endparam
  param bailout
    caption = "Bailout value"
    default = 128.0
    min = 1.0
    hint = "Defines how soon an orbit bails out, i.e. doesn't belong \
            to the Mandelbrot set anymore."
  endparam
switch:
  type = "Jritper3"
  seed = #pixel
  power = power
  bailout = bailout
}

柳烟

我弄出来与榕兄的帖图相近了,我的图与原图差了180度,.但与原图有点差别,仔细检查了计算,应该说没有问题.我仍不明原因何在.我作此图时,先是对c着色,再合并z,结果作出的图乱糟糟,我又采取先合并z,再对c进行着色,结果得下图:

榕坚

真奇怪了，这个分形我破天荒的做了四次，每次得到的结果都不一样，第一次因为代码中s=s-1/s前有分号，而且是灰色的，我以为是注释（而且我试着把原软件中s=s-1/s删掉、保存。再打开结果没变。）就不算，后面算一次，算两次，可结果还是不行，晚上再来。就不信做不出来。

柳烟

我又干了一遍,结果与前次效果一样.
我把阔值搞得无比大后,看起来还不错,效果:

下载 (56.89 KB)

2010-6-17 20:31

榕坚

第五次了，真是气遍我了下面两幅是s=1-1/s算一遍与算两遍的结果：

榕坚

原来这是M集的变异：

榕坚

柳老师不要再试了，我重新跳过s=1-1/s做出来的形状与它一样，但是可能是软件有经过特殊的处理才会有那个结果。其实我第一次就做对了。只是结果与它不一样，外观是M集的外观。怪不得我删掉那两行没有影响图形。

柳烟

下载 (71.44 KB)

2010-6-17 22:48

这个分形也是累我干了好几次。
complex st = #pixel
  if @insout
    st = 1/st
  endif
  complex sqst = sqrt(sqr(st) - 0.33333333333333333333)
  complex s = sqst + st
  complex t = sqst - st
  complex a3 = 0.33333333333333333333
  complex a2 = -0.5*(s + t)
  complex a1 = -0.33333333333333333333
  complex q = ((a3*s + a2)*s + a1)*s
  complex r = ((a3*t + a2)*t + a1)*t
  complex pn2 = 2/(q - r)
  complex u = 1 - q*pn2
  complex p = sqrt(pn2)
  a2 = a2*p
  a1 = a1*pn2
  complex a0 = u*p
  z = p
loop:
  z = ((a3*z + a2)*z + a1)*z + a0
每次得到的结果是：

下载 (35.23 KB)

2010-6-17 22:48

就是得不到UF中的帖图。