Welcome to Fractal Forums

Fractal Art => Images Showcase (Rate My Fractal) => Topic started by: Pauldelbrot on February 21, 2014, 05:18:07 PM




Title: Siegels and Seahorses Blended
Post by: Pauldelbrot on February 21, 2014, 05:18:07 PM
Siegels and Seahorses Blended

(http://nocache-nocookies.digitalgott.com/gallery/15/511_21_02_14_5_18_06.jpeg)

http://www.fractalforums.com/index.php?action=gallery;sa=view;id=15629

A Triskelion Julia set with parameters very close to "Siegels and Seahorses Forever". But all three attractors have been destabilized, so the Julia set becomes spacefilling. Colored using curvature average or something similar, it actually still gives the three quadratic structures the same colors as before; but now they are disconnected quadratic Julia shapes entwined infinitely about every single point.


Title: Re: Siegels and Seahorses Blended
Post by: knighty on February 22, 2014, 05:11:29 PM
That's awesome. Any info about the method for obtaining the formulas?  There are some docs about mating julia sets but I couldn't find anything about (ahem!) mating 3 or more julias (:evil1:  :tongue1:)


Title: Re: Siegels and Seahorses Blended
Post by: Pauldelbrot on February 23, 2014, 02:36:10 AM
That's awesome. Any info about the method for obtaining the formulas?  There are some docs about mating julia sets but I couldn't find anything about (ahem!) mating 3 or more julias (:evil1:  :tongue1:)

Thanks.

I started with the form f(z) = \frac{1}{z^2} + \frac{z^3 - 1}{g(z)}. For most choices of g that will have fixed points at the cube roots of unity, since the left hand term fixes the cube roots and the right hand term is zero there.

The expression for the derivative of the above is complicated, but when z is a cube root of unity it simplifies significantly (because all instances of the right hand term's numerator go to zero) and is just f'(z) = -\frac{2}{z^3} + \frac{3z^2}{g(z)}.

That suggested making g the sum of three polynomials, each with zeros at two of the three cube roots of 1. I defined h(z) = z^2 + z + 1, which is zero at the two primitive cube roots of unity and 3 at 1, and g(z) = \frac{h(z)}{a} + w^2\frac{h(w^2z)}{b} + w\frac{h(wz)}{c}. Note how the three terms differ by being rotated by 0, 1/3, or 2/3 of a circle about the origin (via multiplication by w where w is a primitive cube root of 1) and scaled by a different parameter a, b, or c. Consider the effect of the three terms on g(1). The second and third terms are both zero there, because they contain a factor of h(w) and h(w^2), respectively, and h is zero at both points. The first term is 3/a there, and so f'(1) = -\frac{2}{z^3} + \frac{3}{g(1)} = a - 2. So by controlling a one controls the derivative at the fixed point 1, and so the stability of that fixed point, and indeed can essentially directly input that fixed point's dynamics and "dial-a-Julia".

The rotations of the other two terms of g produce similar results for w and w2, and b and c can be used to dial Julias at those two sites, separately and independently from one another and from the one at 1.

The attachment shows a prototype implemented in UF's formula language being configured. You can essentially directly read off the stability parameters (magnitude and internal angle) from the display, and see how they correspond to components of the fractal. For example, seed a expressed in polar form (relative to the point 2) as magnitude 1, angle 0.61803398 (where the angle is scaled to go from 0 to 1, not 0 to 2pi or 0 to 360) corresponds to a Golden Mean Siegel disk Julia, and indeed that's what has its Siegel disk centered on the point 1, with the disk's preimages clustered all around it. (You might recognize the fractal as "Siegels and Seahorses Forever".)

I've made further refinements to the formula. Polar can be toggled on and off (rectangular is much nicer for importing points picked from Mandelbrot views, while polar is preferable for direct-dialing a particular Julia) and there are some additional options. For example one can add a term t(z^3 - 1) to g without perturbing the derivatives at the cube roots of 1, while adding another parameter t that varies the fractal, and one can change \frac{h(z)}{a} + w^2\frac{h(w^2z)}{b} + w\frac{h(wz)}{c} to z(\frac{h(z)}{a} + w\frac{h(w^2z)}{b} + w^2\frac{h(wz)}{c}) without perturbing the derivatives while altering the dynamics in various ways. Playing around with these I've gotten as many as five attracting basins (or "near-misses").


Title: Re: Siegels and Seahorses Blended
Post by: knighty on February 24, 2014, 03:31:03 PM
Thank you very much :)
Does it works for any parameter set (a,b,c) ?
Is it possible to generalize it by using f(z)=1/z^(n-1) + (z^n - 1) / g_n(z) where g_n(z) fixes the n-th roots of unity?


Title: Re: Siegels and Seahorses Blended
Post by: Pauldelbrot on February 25, 2014, 08:15:08 AM
Thank you very much :)
Does it works for any parameter set (a,b,c) ?
Is it possible to generalize it by using f(z)=1/z^(n-1) + (z^n - 1) / g_n(z) where g_n(z) fixes the n-th roots of unity?

In theory. The method of constructing g should also generalize (with h being a product of z - ɑ terms as ɑ varies over all but one of the nth roots of unity).


Title: Re: Siegels and Seahorses Blended
Post by: knighty on February 27, 2014, 07:52:48 PM
Thank you. I think this fractal type deserves a thread of it's own.
A -quite abstract math- question: Conjugating f(z) with a mobius transform this way: M-1(f(M(z)); should give a similar but deformed fractal. Is there a set of Mobius transforms (that would depend on (a,b,c) paramters) that keeps the fractal invariant? if true, what is the structure of that set if any? I gess it's a -kleinian- group but I'm not sure at all.


Title: Re: Siegels and Seahorses Blended
Post by: Pauldelbrot on February 27, 2014, 08:32:37 PM
Thank you. I think this fractal type deserves a thread of it's own.
A -quite abstract math- question: Conjugating f(z) with a mobius transform this way: M-1(f(M(z)); should give a similar but deformed fractal. Is there a set of Mobius transforms (that would depend on (a,b,c) paramters) that keeps the fractal invariant? if true, what is the structure of that set if any? I gess it's a -kleinian- group but I'm not sure at all.

I'm really not sure. Sorry...


Title: Re: Siegels and Seahorses Blended
Post by: cKleinhuis on February 27, 2014, 08:40:26 PM
Thank you. I think this fractal type deserves a thread of it's own.
A -quite abstract math- question: Conjugating f(z) with a mobius transform this way: M-1(f(M(z)); should give a similar but deformed fractal. Is there a set of Mobius transforms (that would depend on (a,b,c) paramters) that keeps the fractal invariant? if true, what is the structure of that set if any? I gess it's a -kleinian- group but I'm not sure at all.
lol, me neither ;)


Title: Re: Siegels and Seahorses Blended
Post by: knighty on March 01, 2014, 03:46:01 PM
:D
The question was about the existance of "simple" symmetries in this kind of fractals.
After "googling" a little I found these:
http://www.math.dartmouth.edu/~doyle/docs/icos/icos/node1.html
http://www.math.univ-toulouse.fr/~buff/Symmetries/Patterns.html
Well... it seems the set of symmetric rational mappings is very small.