柳烟

108# 柳烟
我看到你帖子中说，大家都忘了，以为可能你作出后，文件搞丢了，想不起来。有结误会。牛集p值闭上眼睛都能知呀，这是常识，要不，搞啥分形嘛。
我先按你的提示意思，理解成这样：
init:
  z =pixel
loop:
  z1=z
  z =(1-z^2)/(z-z^2*cos(z)) -1
bailout:
  |z-z1| >= @bailout
运行结果，不成。
然后，我想到，MJ集是|z|<= @bailout，既然按N集的p算，想了想，会不会是z|>= @bailout，所以成了前面我那帖子中的样子。

柳烟

此图来自网上：

下载 (28.35 KB)

2014-7-9 11:46

迭代公式：(1-z)/(8-4*z+2*z^2-z^3)+c.此图结构太美，但不明网上算法。我张帖在此，算是征解。两位老师思维活泛，看能否解决。

柳烟

上楼迭代公式的牛M集，要补等势圈。UF中自动生成等势圈，用GSP造，非补不可。分形很美，扫图一幅：

下载 (18.21 KB)

2014-7-9 11:46

UF中效果图：

下载 (10.74 KB)

2014-4-23 22:25

柳烟

125# 分形几何
从网上看，临界点象是一阶导数的零点。用maple计算结要是-.5459890945

柳烟

127# 分形几何
大体与原图吻合，你的思路应该没错，有何特别算法吗？c值可在UF中进行嗅探，找到与原图非常接近的c值应该没问题。我照普通的J集整，好象不行。

柳烟

130# 分形几何
请问胡老师，你的c定位值是多少，我造N集后，整出的图不对劲。尽得到些圆饼，啥原因？

柳烟

132# 分形几何
请问图片上的数据信息对造此分形有没有用？有用，有何作用？我用图片下的数据造。我曾用此复公式造M集，造出后，发现找小M集埋得太深，将图缩小后，发现与不要1000图大致相当，感觉到不要1000可有行，结果正如你指出的，肯定离原图差得很远。

柳烟

越整越不象样，用画板作，用迭代终点产生的em上色，结果逃逸点太多，图离正品差了十万八千里，只好放弃算了。 谢谢胡老师帮忙。

柳烟

141# 榕坚
将太靠近根的点与无穷的点给拖回来重新迭代，我理解成：根为z0=1,给定任意小的数, 比如0.00001,|z-z0|<0.00001(太靠近根的数)，给定任意大的数比如：10000，|z|>10000(趋向无穷的数)，并且按牛集的p，这三者的乘积作为现在的p进行点的迭代，不知我这理解正确与否？

我这理解，实践结果不成。这句技术上如何处理？

柳烟

z-(1000*(1-z)/(8-4*z+2*z^2-z^3)+C的N {
init:
  z =pixel
loop:
  z1=z
if |z|<1e20&&|z-1|>1e-12
  z = z-(1000*(1-z)/(8-4*z+2*z^2-z^3)+@c)/ (-1000/(8-4*z+2*z^2-z^3)-(1000*(1-z))*(-4+4*z-3*z^2)/(8-4*z+2*z^2-z^3)^2)
else
  z=z1
endif

bailout:
  |z-z1| >= @bailout
default:
  title = "z-(1000*(1-z)/(8-4*z+2*z^2-z^3)+cN"
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\julia.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param c
    caption = "c"
    default = (-.96164, .13619)
    hint = "Use this parameter to create many different Julia sets. A good \
            way to set this parameter is with the Switch, Eyedropper, or \
            Explore features."
  endparam

  param bailout
    caption = "Bailout value"
    default = 0.0001
    min = 0.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 = "Mandelbrot"
  power = power
  bailout = bailout
}
在画板中搞，又按胡老师与榕坚老师的提示，编写上面代码，不用et改用em作色，结果所得图皆如下图：

下载 (12.33 KB)

2014-4-26 20:33

问题出在那里呢？牛集不是照这个样子的吗？