柳烟

好象N集光滑作色，不是完全抹平等势圈，如果完全抹平后，每块缺乏立体感，成平面状，边界看起来也不太美。

下载 (59.36 KB)

2012-7-27 20:33

这是UF中的效果，各块看起来光滑，仔细一瞧，仍有等势圈，只不过等势圈十分摸糊，晃眼一看，很是光滑。不知能否象圆陷阱那样，使每块成为 伪3D.

柳烟

4# changxde
在ahm.ucl中，里面只有一个特效文件。打开任一牛集，再点选outside选项卡，找到ahm.ucl中的文件Extrapolated Bailout，打开后，再在面板中选中四个开关项的newton，即得此图形。代码可单独将此开关项抽出来。这个效果我干了几次，整不出这效果，可能我的造法有问题，代码有些地方费解。

柳烟

我将那个特效文件，进行大幅度缩水，仅保留newton开关项部分，并动了手脚后，代码及图形效果如下：
VaryingBailout {

final:
complex x = #pixel
If Imag(#pixel)==0
  x=#pixel + 1e-30i
Endif

int iter = 0
repeat

   xold=x
   x=x -( x^3-1)/(3*x^2)

  iter=iter+1
until (|x-xold|<=@SmallBailout )

  float bailTest=|x-xold|
  Power=2
If iter!=1000
    #index=cabs(sqrt(log(Power)*(#maxiter-iter) + log(log(bailTest))))
Else
   #index=sqrt(log(log(bailTest)))
Endif
default:
title="Extrapolated Bailout"

float param SmallBailout
  caption="Small Bailout"
  hint="Should be positive. Smaller values mean more time, but more accurate \
  pictures (assuming you increase the iterations)."

  default=0.0001
endparam

}

下载 (196.95 KB)

2012-7-28 07:48

柳烟

我今天试着模糊等势圈，用常老师法，我不得法，结果等势圈越整越浓，用|z-z0|的终象除任意小正数ε(阈值)与et搭配，我是分二个阵营，一是||dM||>ε，二是||dM||<ε分别调色。如何搭配才得当？
我试着翻译8楼光滑N集的代码，结果#index=cabs(sqrt(log(Power)*(#maxiter-iter) + log(log(bailTest))))不好驾驭，弄出无意义。又until (|x-xold|<=@SmallBailout )是啥意思？与后文如何串联？丈二和尚摸不着头脑。再精简代码如下，我在UF中验证了，得到的图形完全一样。个别句子再动手脚。
VaryingBailout {

final:
complex x = #pixel


int iter = 0
repeat

   xold=x
   x=x -( x^3-1)/(3*x^2)

  iter=iter+1
until (|x-xold|<=@SmallBailout )

  float bailTest=|x-xold|
  Power=2
If iter-#maxiter<0
    #index=cabs(sqrt(log(Power)*(#maxiter-iter) + log(abs(log(bailTest)))))
Else
   #index=sqrt(log(abs(log(bailTest))) )
Endif
default:
title="Extrapolated Bailout"

float param SmallBailout
  caption="Small Bailout"
  hint="Should be positive. Smaller values mean more time, but more accurate \
  pictures (assuming you increase the iterations)."

  default=0.0001
endparam

}
短短的代码，怎么那么麻烦，怪得出奇。

柳烟

13# changxde
漂亮！看了你这帖，我仅保留代码IF至endif中的第一个#index,将判断的其余部分删去，仍不影响效果。上午，我仍失败，看了这帖，很受启发，我再试试看。

柳烟

15# 榕坚
如果将红色部分去除，在UF中可发现等势圈大显。

柳烟

请教一个问题x-xold|>=@SmallBailout 实际上就是用GSP造牛集时的p，而代码中的：|x-xold| 〈=@SmallBailout ，其实就相当于1-p，要造出牛集，缩放点z-f(z)/f'(z),用的是p,而非1-p，否则造不出链条。决定是否终止循环的应该是p,而非1-p,我没有考虑红色部分，造出的效果如下：

下载 (91.11 KB)

2012-7-31 17:53

请问红色部分究竟在画板中，如何关联？

柳烟

按UF的光滑作色法,扫一M集:

下载 (75.89 KB)

2012-8-2 23:25

按#index=sqrt(ln(2)*(n-et) + log(abs(|zn|)))),这里|zn|为模的迭代终值.我弄出后,扫出的M集的肚皮上,发现有部分白色溢出,采取了一些补洞措施,其中:迭代次数为n=100时,我用的是n+2-et,且在根号内加上了一个正整数.不知有无更好办法.这是在前面说的那个光滑代码中,找到的并进行了一些改造,以达填洞的目的.大家看看有无更好办法.