此分形在UF中的Pwc-convert.ufm中的Pokorny Julia Mod。
Mod2_Pokorny-bare {
; Original formula, Mod2_Pokorny_Rings, by Paul Carlson.
; Modified to separate coloring formula by Ken Childress.
; 15 June 1999
init:
complex c = @JuliaCoords
z = #pixel
loop:
z = 1.0 / (z * z + c) + c + 1.0
bailout:
|z| < @bailout
default:
title = "Pokorny Julia Mod"
maxiter = 5000
param JuliaCoords
caption = "Julia Set Coordinates"
default = (-2.25,0)
hint = "Julia Set Coordinates."
endparam
param bailout
caption = "Bailout Value"
default = 4.0
hint = "A bailout value of 4 is recommended to work with AtanMethods."
endparam
}作者: 柳烟 时间: 2010-10-27 23:28
上面这个分形,初看不起眼,但放大后,有意思。
迭代式简单: z = 1.0 / (z * z + c) + c + 1.0
C定位于(-2.25,0),阈值为4.
下面是我在面板上将c定位于(-2.8,0)后,在UF中放大后的图形。
位于mt.ufm
mt-attractor-julia {
; Mark Townsend, 29 Mar 1999
;
; Julia sets that bail out when the orbit reaches
; (hopefully) a periodic attractor. I iterate the
; critical point z = (0,0) before the main loop
; and test against the resulting value. In the
; hope of speeding the process up I include
; the option to give up the search when z
; grows too large.
;
init:
w = 0
u = #pixel
int iter = 0
if @no_bail
while iter < @max
w = w^2 + @c
iter = iter + 1
endwhile
else
while (iter <= @max) && (|u| < 10)
w = w^2 + @c
u = u^2 + @c
iter = iter + 1
endwhile
endif
z = #pixel
loop:
z = z^2 + @c
bailout:
|w - z| > @epsilon && |z| < @bailout
default:
title = "Attractor Julia"
maxiter = 1000
periodicity = 0
param c
caption = "Julia seed"
default = (-0.12,0.74)
endparam
param epsilon
caption = "Epsilon"
default = 1e-6
hint = "This is the bailout value for convergent orbits."
endparam
param bailout
caption = "Bailout value"
default = 100.0
hint = "This is the bailout value for divergent orbits."
endparam
param max
caption = "Search maxiter"
default = 100
hint = "The is the maxiter value for the attractor \
search."
endparam
param no_bail
caption = "Every pixel"
default = false
hint = "This determines whether to bailout of the attractor \
search when z grows too large. Enable this when the \
image is mainly inside."
endparam
switch:
type = "mt-attractor-m"
max = max
epsilon = epsilon
bailout = bailout
no_bail = no_bail
}
我老是作不出此图。UF效果是:
改一下代码,结果是一样的:
mt-attractor-julia {
; Mark Townsend, 29 Mar 1999
;
; Julia sets that bail out when the orbit reaches
; (hopefully) a periodic attractor. I iterate the
; critical point z = (0,0) before the main loop
; and test against the resulting value. In the
; hope of speeding the process up I include
; the option to give up the search when z
; grows too large.
;
init:
w = @c
z = #pixel
loop:
z = z^2 + @c
bailout:
|w - z| > @epsilon && |z| < @bailout
default:
title = "Attractor Julia"
maxiter = 1000
periodicity = 0
param c
caption = "Julia seed"
default = (-0.12,0.74)
endparam
param epsilon
caption = "Epsilon"
default = 1e-6
hint = "This is the bailout value for convergent orbits."
endparam
param bailout
caption = "Bailout value"
default = 100.0
hint = "This is the bailout value for divergent orbits."
endparam
param max
caption = "Search maxiter"
default = 100
hint = "The is the maxiter value for the attractor \
search."
endparam
param no_bail
caption = "Every pixel"
default = false
hint = "This determines whether to bailout of the attractor \
search when z grows too large. Enable this when the \
image is mainly inside."
endparam
switch:
type = "mt-attractor-m"
max = max
epsilon = epsilon
bailout = bailout
no_bail = no_bail
}