changxde

看了xuefeiyang老师的提示，把开方换成平方，图形出来了。

下载 (26.33 KB)

2010-5-22 17:30

柳烟

根号的那个没作成功，飞扬老师推出的平方的那个M集，作出来了，将Z初始位置定在原点。

下载 (41.64 KB)

2010-5-22 22:34

下载 (38.24 KB)

2010-5-22 23:16

榕坚

难道那个sqr是平方函数？我看到式子中有＾2就想象那个sqr应该是开平方了。从那个软件的放大来看这个分形的结构是很丰富的。

柳烟

问题是SQR确实是平方根的记号，程序者不可能犯这个低级的错误。如果是平方函数，为何没有用到那两个r1，r2那两个返回条件而造出图来了呢？而不用平方根，反而弄出榕兄所发图片中间的部分了呢？迷糊。我搞了几个钟点，始终不能出现外面的等势圈。

柳烟

看来xuefeiyang老师的提示，把开方换成平方，图形出来了。
4558
changxde 发表于 2010-5-22 17:30

changxde老师此图的白色部分有点逸出，有问题。

xuefeiyang

仅仅是算法的不同，一个用了逃逸时间算法，一个没有用逃逸时间算法。不用逃逸时间算法就会出现空白。

榕坚

看来是被那个sqr给误导了，定位、缩放、迭代次数一样扫出来的结果比较：其中281（1）是分形软件的结果

榕坚

这两个分形就真的怪了，总是做不成功。只好再次请求协作: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
}

榕坚

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

xuefeiyang

你肯定是哪里计算有了问题。我试过z=0没有问题。图形如下：

下载 (37.98 KB)

2010-5-23 14:10

另外这个分形可以看作是一个特殊的N集，如果把阀值设值设为很小时，你会看到在这里看不到的一些现象！