返回列表 回复 发帖
39# xiaongxp
经向老师这一说,懂了。谢。
38# xiaongxp


不过我再一次对比一下除了四个卷线以外接近卷线处的图形走向好象也不一样,那个螺纹线可能是某种逃逸条件造成的结果。
今天我弄榕坚推荐的软件中的这一范例,整出的图怪怪的,大家帮忙想想办法。原软件中的程序英文,由于我底子差,有些地方可猜出一二,有些地方不懂,如读天书。榕老师、梅老师英文棒。我把此摆在这公众出入多的地方:
原程序如下:
N_sin_M {
; newton's method applied to f(z) = sin(z) - c
; Q(z) = z - f(z)/f'(z) = z - ((sin(z)-c)/cos(z))

global:
   
init:
  z = 0
  complex delta_z = 0
  complex c = #pixel
  bool continue = true
  
loop:
  delta_z = ( sin(z) - c ) / cos(z)
  
  if |delta_z| < @delta
    continue = false
    if @zmode=="step"
        z = delta_z
    endif ; zmode
  else
    z = z - delta_z
  endif ; < @delta
  
bailout:
; -- this isn't the bailout but the continuing condition
;    false = bail; true = continue
  continue
  
default:
  title = "Newton-M, sin(z)-c" ; [v1.0]
  magn = 0.25
  periodicity = 0
  
  heading
    caption = "Bailout"
  endheading
  float param delta
    caption="Bailout delta"
    hint="Square of max step. Use small values."
    default=1e-8
    min=0
  endparam
  param zmode
    caption = "final z"
    enum = "z" "step"
    default = 1
  endparam
  
switch:
  type = "N_sin_J"
  cv = #pixel
  delta = delta
  zmode = zmode
  
}
制成的效果如下:
Fractal4.jpg
2010-6-1 20:36

图片看起来很夺人眼球的。我抵挡不住诱惑,一直在弄,均以失败告终。按牛顿迭代f(z) = sin(z) - c
; Q(z) = z - f(z)/f'(z) = z - ((sin(z)-c)/cos(z))的M集。谁先整出,请及时公开源文件。
32# 榕坚
我是错了,谢谢榕老师的指正。用了足足二个半小时才扫出这幅,效果也不怎么样,受不了 。迭代次数5000,放大12亿。还是用原来的那个文件,见8楼。
未命名.JPG
46# mjj_ljh


这个范例要求的迭代次数太高了,无法承受。不过它是钻到M集深层里去了。要想把那个纯色部分去掉迭代次数估计要10000次了。
45# 柳烟


这个应该与M集的做法一样啊,怎么会得不到结果呢?只是因为牛顿的逃逸条件比较特殊,你注意范例中的逃逸限制值及逃逸条件即可。
M集就是看不厌,我就是迷它。再发一个链接。
http://www.miqel.com/fractals_ma ... _fractal_guide.html
http://www.miqel.com/fractals_ma ... ndelbrot-magic.html
http://www.miqel.com/fractals_ma ... als_mandelbulb.html
以下是上面的一些图片,我们共同欣赏。
mardi-graz_mandelbrot_fractal.jpg
mandel_droplets_fractal.jpg
2009_orange_green.jpg
49# mjj_ljh


实在是太漂亮了,开眼。
返回列表