返回列表 回复 发帖
这个图形我两年前作过,你下载有分形图形学这本书没有?
xuefeiyang 发表于 2010-5-22 15:03
好象下载过,不是pdf格式。还没看。
利用119榕老师给的函数进行迭代发现迭代点的图像也很好看,有点像龙图。
changxde 发表于 2010-5-22 16:47
确实不错,这应该就是它的J集中的某一个图形了。如果各点让它颜色有变化就可能会更好看。
难道那个sqr是平方函数?我看到式子中有^2就想象那个sqr应该是开平方了。从那个软件的放大来看这个分形的结构是很丰富的。
看来是被那个sqr给误导了,定位、缩放、迭代次数一样扫出来的结果比较:其中281(1)是分形软件的结果

捕获280.JPG (34.84 KB)

捕获280.JPG

捕获281.JPG (109.5 KB)

捕获281.JPG

捕获281(1).gif (74.55 KB)

捕获281(1).gif

这两个分形就真的怪了,总是做不成功。只好再次请求协作:Magnet2Mandelbrot {
;
; Magnetic Mandelbrot set type 2. Use Switch Mode to select a
; magnetic Julia set.
;
init:
  z = p1
loop:
  z = sqr( (z^3 + 3 * (pixel-1) * z + (pixel-1) * (pixel-2))  /        \
           (3 * z^2 + 3 * (pixel-2) * z + (pixel-1) * (pixel-2) + 1) )
bailout:
  |z| < @bailout && |z - 1| > @lowerbailout
default:
  title = "Magnet 2 (Mandelbrot)"
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\magnet.html"
  center = (1.1, 0)
  magn = 0.77
  maxiter = 100
  periodicity = 0
  param p1
    caption = "Perturbation"
    default = (0, 0)
    hint = "Starting value for each point. You can use this to \
            'perturb' the fractal. Use (0, 0) for the classic set."
  endparam
  param bailout
    caption = "Bailout value"
    default = 100.0
    min = 1
$IFDEF VER40
    exponential = true
$ENDIF
    hint = "This parameter defines how soon an orbit bails out while \
            iterating. Larger values give smoother outlines; smaller values \
            generally produce more interesting shapes around the set."
  endparam
  param lowerbailout
    caption = "Convergent bailout value"
    default = 0.00005
    min = 0
$IFDEF VER40
    exponential = true
$ENDIF
    hint = "This parameter defines how soon a convergent orbit bails out while \
            iterating. Smaller values give more precise results but usually \
            require more iterations."
  endparam
switch:
  type = "Magnet2Julia"
  p1 = #pixel
  bailout = bailout
}

捕获.JPG (77.83 KB)

捕获.JPG

这是我做的结果,有一点不明确:式子中当Z=0时显然有意义,但初始点选原点就没有图象了。

捕获.JPG (48.57 KB)

捕获.JPG

谢谢feiyang老师,这个图形我已完成了,结果如:137#,我查了软件的函数列表:其中sqr确实是平方。我说的是138# :
z = sqr( (z^3 + 3 * (pixel-1) * z + (pixel-1) * (pixel-2))  /        \
           (3 * z^2 + 3 * (pixel-2) * z + (pixel-1) * (pixel-2) + 1) )
下午再做一次无果,得到不成样的结果:
捕获282.JPG

捕获283.JPG
再做一遍总算好了,因为项太多了可能某一项漏算了:

捕获284.JPG (23.59 KB)

捕获284.JPG

我怎么整不出整个图形外面的等势圈呢?
柳烟 发表于 2010-5-23 20:22
跟上次做简化牛顿是一样的。
返回列表