返回列表 回复 发帖
“在同一逃逸区内的点离边界的远近是否能用什么加以区别。”是指着色吗?
一般来讲,每个逃逸区内的点按势值(等高线)具有相同的遍历性,所以我作的色因子常带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
}

Fractal1.jpg (56.39 KB)

Fractal1.jpg

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

未命名.JPG (56.61 KB)

未命名.JPG

真奇怪了,这个分形我破天荒的做了四次,每次得到的结果都不一样,第一次因为代码中s=s-1/s前有分号,而且是灰色的,我以为是注释(而且我试着把原软件中s=s-1/s删掉、保存。再打开结果没变。)就不算,后面算一次,算两次,可结果还是不行,晚上再来。就不信做不出来。
我又干了一遍,结果与前次效果一样.
我把阔值搞得无比大后,看起来还不错,效果:
未命名2.JPG
第五次了,真是气遍我了下面两幅是s=1-1/s算一遍与算两遍的结果:

捕获348.JPG (36.86 KB)

捕获348.JPG

捕获349.JPG (29.13 KB)

捕获349.JPG

原来这是M集的变异:

Fractal1.jpg (27.65 KB)

Fractal1.jpg

柳老师不要再试了,我重新跳过s=1-1/s做出来的形状与它一样,但是可能是软件有经过特殊的处理才会有那个结果。其实我第一次就做对了。只是结果与它不一样,外观是M集的外观。怪不得我删掉那两行没有影响图形。
Fractal1.jpg
这个分形也是累我干了好几次。
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
每次得到的结果是:
未命名.JPG
就是得不到UF中的帖图。
返回列表