返回列表 回复 发帖
pplane-4.gif
此分形我造了一遍,不对劲,大家想想办法.网址连接:
http://math.bu.edu/people/bob/papers.html
下了相应的文件,进行阅读,因为英文,所以阅读有些难度,算法不是太清楚.
111# 柳烟


算法应该没有什么特殊,只是Z的初始值是关键。
112# 榕坚
此类分形的J集问题不大,主要是M集C的定位.按照以前造复分形的经验,要得到标致的乌龟,必将Z值定位于二阶导数为0的z的方程的零值点,如经典M集,一些牛顿M集,均符合这一规律,否则得到的乌龟均走了样.我昨天干了一整天,均未找到这个点,怪了.可能有别样算法.
113# 柳烟
应该是Z的定位,请柳老师改一下。
找到定位点了:

捕获.JPG (24.66 KB)

捕获.JPG

定位点在那个地方?
116# 柳烟


还是一样的,定位于一阶导数的零点。
按榕坚老师提供的定位法,偷一下懒,将UF中已有的范例文件略改,迅速成就此分形,发张张图片。程序代码如下:
z^n+c/z^n的Mandelbrot {

init:
  z = #pixel^(1/(2*@power))
loop:
  z = z^@power + #pixel/(z^@power)
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF

  param power
    caption = "Power"
    default = (4,0)
    hint = "This parameter sets the exponent for the Mandelbrot formula. \
            Increasing the real part to 3, 4, and so on, will add discs to \
            the Mandelbrot figure. Non-integer real values and non-zero \
            imaginary values will create distorted Mandelbrot sets. Use (2, 0) \
            for the standard Mandelbrot set."
  endparam
  float param bailout
    caption = "Bailout value"
    default = 4.0
    min = 1.0
$IFDEF VER40
    exponential = true
$ENDIF
    hint = "This parameter defines how soon an orbit bails out while \
            iterating. Larger values give smoother outlines; values around 4 \
            give more interesting shapes around the set. Values less than 4 \
            will distort the fractal."
  endparam
switch:
  type = "Julia"
  seed = #pixel
  power = power
  bailout = bailout
}

大家可新建一个后缀名为ufm的文件,将此代码粘上去,关闭保存,再用UF打开,即可看到。用GSP造,迭代部分比较麻烦。
Fractal1.png
感谢老友的帮忙,终于解脱了。
再把边界给提取出来,最喜欢这样的边界了:

捕获.JPG (62.64 KB)

捕获.JPG

一个J_卷,其外部集三种不同算法于一身,不知如何实现?
03.jpg
返回列表