柳烟

29楼，将阀值搞得很巨大，圆柄萎缩至几近无了。J集还真是无法。

2015-1-28 22:48

lnszdzg

270# 柳烟


能分享一下UF代码不？J集我就是弄不出来啊！

柳烟

272# lnszdzg
[1-z2]/[z-z2cosZ]+CJ{

init:
  c =@e
  z=#pixel
loop:
  z1=z
  z2 = (1-z1^2)/(z1-z1^2*cos(z1))+c
  z=(1-z2^2)/(z2-z2^2*cos(z2))+c
bailout:
|z-z1| >= 0.00001
  default:
  title = "J"
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\julia.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param e
    caption = "M seed"
    default = (-1, 0)
    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 = 120.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

}
bailout后面，这样的逃逸条件也行：|z|>10000&&|z-z1| >= 0.00001.J集里的圆柄不知如何让其消失，可能有高手知道，但不想教，只好在黑暗中再探索一阵子看能否碰巧了。

柳烟

画板玩分形有几大毛病：1.迭代次数太受限；2.成图慢，特别是迭代次数大时的放大图，纯粹是与耐力比拼；3.搞一个分形化时多，调色太费周章……用程序造作分形，似乎优点多多。看到各位老师用编程法玩分形玩得风吼，也有点心动，但有点畏难情绪……

2015-1-29 08:39

lnszdzg

邀请大家到Mathcad吧做客。
http://tieba.baidu.com/f?kw=mathcad&fr=ala0

lnszdzg

273# 柳烟

用z2和z迭代了两次？我说我怎么弄不出来呢！

lnszdzg

柳老师，次分形对于阀值和次数十分敏感，稍有不同便大相径庭

2015-1-29 21:44

这是根据您的程序做的，这是UF代码及其参数，您看看
Mandelbrot5 {
init:
  z = #pixel
  c = @start
  z2 = 10
  z1 = 20
loop:
  z1 = z
  z2= (1 - z1^@power)/(z1 - z1^@power*cos(z1))+ c
  z =(1 - z2^@power)/(z2 - z2^@power*cos(z2))+ c
bailout:
   ;|z| < @bailout1
  |z-z1| > @bailout && |z| < @bailout1
default:
  title = "Mandelbrot5"
  center = (0, 0)
param start
    caption = "Starting point"
    default = (-1.05445654699195010595,0)
  endparam
   param power
    caption = "Power"
    default = (2,0)
  endparam
float param bailout
    caption = "Bailout value"
    default = 1E-6
  endparam
param bailout1
    caption = "Bailout1"
    default = 1200000.0
  endparam
}

lnszdzg

2015-1-29 21:45

c=-1.2543/0.0

lnszdzg

2015-1-29 21:46

-1.2443/0.00

lnszdzg

-1.2443/0.00037

2015-1-29 21:46