Author Topic: can you identify this fractal ? (Read 770 times)

Adam Majewski · « **on:** June 04, 2022, 10:21:51 AM »

can you identify this fractal ? What function ?

It is fig 3 from paper : Near Parabolic Renormalization for Unicritical Holomorphic Maps by Arnaud Chéritat

It has degree 3. I do not know if it rational or polynomial.
It is parabolic.

TIA

Linkback: https://fractalforums.org/index.php?topic=4776.0

3DickUlus · « **Reply #1 on:** June 04, 2022, 06:06:39 PM »

"Quadratic Julia" from what I read in the paper, just skimming through, have not studied in depth.

Adam Majewski · « **Reply #2 on:** June 04, 2022, 06:18:56 PM »

Quote from: 3DickUlus on June 04, 2022, 06:06:39 PM

"Quadratic Julia" from what I read in the paper, just skimming through, have not studied in depth.

Thx for the answer, but : "Fig. 3 Another illustration of Theorem 6. This time, d = 3. "

Here d is a degree so IMHO it is not quadratic function, where degree d = 2.

What looks strange to me: holes between components. I have not see such Julia sets for polynomials. Am I right ?

marcm200 · « **Reply #3 on:** June 04, 2022, 06:46:39 PM »

Quote from: Adam Majewski on June 04, 2022, 06:18:56 PM

What looks strange to me: holes between components. I have not see such Julia sets for polynomials. Am I right ?

The basin of attraction of infinity (and I suspect those holes to be escaping points) for polynomial Julia sets is connected (cf Przytyck, Zdunik. On Hausdorff dimension of polynomial not totally disconnected Julia sets. Paragraph 2).

hgjf2 · « **Reply #4 on:** June 05, 2022, 05:57:05 AM »

This strange Julia set is from a Nova Mandelbrot set. Must check the Nova Mandelbrot sets by the fractional polynomial. Is obviously that isn't from a polynomial Mandelbrot set.

Adam Majewski · « **Reply #5 on:** June 05, 2022, 06:28:04 AM »

Quote from: hgjf2 on June 05, 2022, 05:57:05 AM

This strange Julia set is from a Nova Mandelbrot set. Must check the Nova Mandelbrot sets by the fractional polynomial. Is obviously that isn't from a polynomial Mandelbrot set.

from kf :
Nova
z := (2 z^3 / 3 - 2 z - 1) / (z + 1)^2 + c + 1
Nova is affine-conjugated from the usual form to move critical point from 1 to 0.

Newton Nova Mandelbrot (z^p - 1)
z := z - (z^p - 1) / (p z^(p-1)) + c
Halley Nova Mandelbrot (z^p - 1)
z := z - (2 (z^p - 1) (p z^(p-1))) / (2 (p z^(p-1))^2 - (z^p - 1) ((p (p-1)) z^(p-2))) + c
Schroder Nova Mandelbrot (z^p - 1)
z := z - ((z^p - 1) (p z^(p-1))) / ((p z^(p-1))^2 - (z^p - 1) ((p (p-1)) z^(p-2))) + c
Householder 3 Nova Mandelbrot (z^p - 1)
z := z - (6 (z^p - 1) (p z^(p-1))^2 - 3 (z^p - 1)^2 ((p (p-1)) z^(p-2))) / (6 (p z^(p-1))^3 - (z^p - 1) (p z^(p-1)) ((p (p-1)) z^(p-2))) + (z^p - 1)^2 (p (p-1) (p-2) z^(p-3)) + c
Householder Nova Mandelbrot (z^p - 1)
z := z - ((z^p - 1) (2(p z^(p-1))^2 + (z^p - 1) (p (p-1) z^(p-2))))/(2 (p z^(p - 1))^3) + c

from uf

"
Nova (Julia, Mandelbrot)

The Nova fractal is a modified Newton fractal. The Julia version can be used as a normal Newton fractal, but there are all kinds of other possibilities with intriguing spirals.

Both Mandelbrot and Julia versions of the Nova fractal are available as fractal formulas in Standard.ufm and as fractal formula plug-ins in Standard.ulb.

Use the Nova (Mandelbrot) formula to switch to interesting Nova (Julia) sets. The standard Newton fractals can be found in the empty circle to the right.

The formulas provide the following parameters:

Start Value
(Mandelbrot only)
For well-formed fractals, this should be set to (1, 0). Other values create distorted shapes that can be interesting, but they are usually not as well-formed.

Julia Seed
(Julia only)
This parameter specifies the point in the Mandelbrot version that corresponds to the current Julia version. It defines the shape and behavior of the fractal. Use the Switch feature to select good values.

Exponent
Increase the real value of the exponent to create more complex fractals. Non-integer real values and non-zero imaginary values create irregular fractals.

Bailout
Specifies the magnitude of z at which the formula will stop iterating. Since z converges to a fixed value, smaller values will give more detailed images.

Relaxation
This can be used to influence the convergence of the fractal. Changing this parameter will twist and transform the fractal.
"

so it is rational function

Alef · « **Reply #6 on:** June 05, 2022, 04:10:33 PM »

Quote from: Adam Majewski on June 04, 2022, 06:18:56 PM

What looks strange to me: holes between components. I have not see such Julia sets for polynomials. Am I right ?

To have a holes formula must have some sort of division involved. Nova julia sets have holes but is different.

Alef · « **Reply #7 on:** June 05, 2022, 04:11:59 PM »

Quote from: Adam Majewski on June 04, 2022, 06:18:56 PM

What looks strange to me: holes between components. I have not see such Julia sets for polynomials. Am I right ?

To have a holes formula must have some sort of division involved. Nova julia sets have holes but is different. It is kind of curved. Here it is linear.

Adam Majewski · « **Reply #8 on:** June 10, 2022, 03:27:18 PM »

Arnaud Chéritat ( author) :

Quote

As far as I remember, this is the Blaschke fraction

z -> rho * z^2 * (z-3)/(1-3z)

where rho is a modulus one complex number chosen so that there is on the unit circle a parabolic period 3 cycle.

The initial picture was made in march 2000, the blue/yellow set was isolated 10 years later using filters in Gimp when I needed to make a slides presentation about parabolic points. Later I used it in the article you mention.
I had not noted the value of rho back then, but you can find it using algebra or numerical root finding.

marcm200 · « **Reply #9 on:** June 10, 2022, 05:28:23 PM »

First attempt, just browsing through some rho values. Here e^(2*i*pi*0.35). I haven't checked cycles or multipliers, and the holes between the three big components are currently filled (the circles, gray pixels, have NaN during orbit, white is bounded, black is escaping 10^25).

Adam Majewski · « **Reply #10 on:** June 11, 2022, 08:37:20 AM »

Code: [Select]

killa(all);
display2d:false;


period : 3;
t : 1.0/period ;
r: 1.0;
rho : r*exp(2*%pi*%i*t) ;
rho:float(rectform(rho));

rho: 0.8660254037844387*%i-0.4999999999999998

Code: [Select]

kill(all);
display2d:false;
ratprint : false; /* remove "rat :replaced " */


GiveListOfCriticalPoints(fun):=
block(
  [d,s],
  /* derivative */
  d:diff(fun,z,1),
  /* critical points z: d=0 */
  s:solve(d=0,z),
  /* remove "z="  from list s */
  s:map('rhs,s),
  /* convert to form x+y*%i */
  s:map('rectform,s),
  s:map('float,s),
  return(s)
)$







period : 3; /* periodic */
t : 1/period ; /* parabolic */
r: 1.0; /* on the boundary */

/* multiplier of periodic point */
rho : r*exp(2*%pi*%i*t) ;
rho:float(rectform(rho));



f(z):= float(rectform((rho * z^2 * (z-3)/(1-3*z))))$



/* f(z) is used as a global function
   I do not know how to put it as a argument */

GiveOrbit(z0,OrbitLength):=
block(
 [z,Orbit],
 z:z0,
 Orbit:[z], 
 for i:1 thru OrbitLength step 1 do
        ( z:expand(f(z)),
         if cabs(z) > 2 then 
         	( print(z0), print("orbit of is escaping"),  
         	  return(Orbit)),
          Orbit:endcons(z ,Orbit)
          ),
 print(z0), print("orbit is not escaping"),  
 return(Orbit) 

)$


/* 
converts complex number z = x*y*%i 
to the list in a draw format:  
[x,y] 
*/
d(z):=[float(realpart(z)), float(imagpart(z))]$

ToPoint(z):= points([d(z)])$

/* give Draw List from one point*/
dl(z):=[d(z)]$

ToPoints(myList):= points(map(d,myList))$





compile(all);

/* ------------  */

iLength:1000;


ListOfCriticalPoints:GiveListOfCriticalPoints(f(z));
multiplicities;

Orbit1: GiveOrbit(ListOfCriticalPoints[1],iLength)$
Orbit2: GiveOrbit(ListOfCriticalPoints[2],iLength)$
Orbit3: GiveOrbit(ListOfCriticalPoints[3],iLength)$

Orbit1 : ToPoints(Orbit1)$
Orbit2 : ToPoints(Orbit2)$
Orbit3 : ToPoints(Orbit3)$


cr1: ToPoint(ListOfCriticalPoints[1])$
cr2: ToPoint(ListOfCriticalPoints[2])$
cr3: ToPoint(ListOfCriticalPoints[3])$

/*-----------------------------------------------------------------------*/
load(draw); /* ( interface to gnuplot ) by Mario Rodriguez Riotorto http://www.telefonica.net/web2/biomates */

draw2d(
    title = "All  critical orbits for f(z)= rho * z^2 * (z-3)/(1-3*z))",
    terminal  = png,
     user_preamble = "set size square; set key right top;", /* 360/26=13.85  ; 360/(2*26)=6,923 */
    file_name = concat("~/Dokumenty/julia_blaszke/period3/", string(iLength)),
    dimensions  = [600, 600],   /* Since Maxima 5.23, pic_width and pic_height are deprecated. */
    /*
    yrange = [-1.0,1.0],
    xrange = [-1.0,1.0],
    */
    xlabel     = "z.re ",
    ylabel     = "z.im",
    
    point_type    = filled_circle,
    points_joined = true,
    
    point_size    = 0.4,
    key="critical orbit 1",
    color             = red ,
    Orbit1,
    point_size    = 3.0,
    key= "critical point 1",
    points_joined = true,
    color           = red,
    cr1,
    
    
    
    key="critical orbit 2",
    point_size    = 0.4,
    points_joined = true,
    color             =blue ,
    Orbit2,
    point_size    = 1.0,
    key= "critical point 2",
    points_joined = false,
    color           = blue,
    cr2,
    
    key="critical orbit 3",
    point_size    = 0.4,
    points_joined = true,
    color             =green ,
    Orbit3,
    point_size    = 1.0,
    key= "critical point 3",
    points_joined = false,
    color           = green,
    cr3
    

 );

Critical points :

Code: [Select]

ListOfCriticalPoints:GiveListOfCriticalPoints(f(z))
(%o17) [1.193004794952656E-8*%i+1.000000003196646,
        0.999999996803353-1.193004772748196E-8*%i,0.0]

or maybe it is simply [1.0,0.0]

So there are 3 components:

critical point 0 is also fixed point
critical point 1 goes to parabolic period 3 orbit
other points escapes ( maybe infinity is also criticla point ? )

But now 3 points are not parabolic ( are inside components)

Adam Majewski · « **Reply #11 on:** June 11, 2022, 10:59:38 AM »

Second attempt
rho is a modulus one complex number chosen so that there is on the unit circle a parabolic period 3 cycle.

marcm200 · « **Reply #12 on:** June 11, 2022, 12:18:34 PM »

My 2nd attempt: Letting rho and z0 (one of the three periodic points) run independently over the unit circle in some small increments, looking for a strict period-3 cycle comprising z0 and checking the multiplier to be (near) parabolic, gave the following best brute-force hit (purely numerical, double precision, not rounding-controlled).

Code: [Select]

  rho [angle p/q]= 1092/3072                           (e^2*i*pi*rho)
  z0 [angle p/q]= 395/1024                             (e^2*i*pi*z0)
  distance p3= 5.9210162510539466761e-05               Euclidean distance last periodic point to first
  |mult|= 0.98813761                                   multiplier
  rho= -0.61523159058062670823,0.78834642762660633863  Cartesian 
  z0 = -0.7531867990436124094,0.65780669329707874837   Cartesion

orbit:
  -0.7531867990436124094,0.65780669329707874837
  -0.39147239529397764857,-0.92018985199946401465
  0.70360334925769740799,0.71059294037680265088
  -0.75314784888214114655,0.65785128845674756271

The image (4-square) shows the indentations, but the holes are still filled. And the p/q values seem rather arbitrary.

Adam Majewski · « **Reply #13 on:** June 11, 2022, 01:59:02 PM »

Adam Majewski · « **Reply #14 on:** June 12, 2022, 06:52:42 AM »

rho : 0.7883464276266063*%i - 0.6152315905806267;
f(z):= float(rectform((rho * z^2 * (z-3)/(1-3*z))))$

there are 3 basins
* basin of z=0 ( easy) with critical point z=0
* basin of infinity ( easy)
* basin of parabolic period 3 cycle ( hard), with critical point z=1

I found 2 period 3 cycles :
* cycle1; [-0.679019512966907,0.7341202224500998], [0.836586884826584,0.5478342670335168], [-0.5274543103865239,-0.849583398175056],
* cycle2; [[0.7032040964766633,0.7109880439912073], [-0.3910900888807991,-0.9203524011916365], [-0.7533753580242582,0.6575907313229275]

both cycles are on unit circle.

Critical point z=1 seems to land on cycle1.

I can't compute multiplier in Maxima ( long expression)

Code: [Select]



kill(all);
display2d:false;
ratprint : false; /* remove "rat :replaced " */



/* f(z) is used as a global function
   I do not know how to put it as a argument */

GiveOrbit(z0,OrbitLength):=
block(
 [z,Orbit],
 z:z0,
 Orbit:[z], 
 for i:1 thru OrbitLength step 1 do
        ( z:expand(f(z)),
         if cabs(z) > 2 then 
         	( print(z0), print("orbit of is escaping"),  
         	  return(Orbit)),
          Orbit:endcons(z ,Orbit)
          ),
 print(z0), print("orbit is not escaping"),  
 return(Orbit) 

)$



/* give Draw List from one point*/
/* 
converts complex number z = x*y*%i 
to the list in a draw format:  
[x,y] 
*/
d(z):=[float(realpart(z)), float(imagpart(z))]$
ToPoints(myList):= points(map(d,myList))$









rho : 0.7883464276266063*%i - 0.6152315905806267;
f(z):= float(rectform((rho * z^2 * (z-3)/(1-3*z))))$



iLength:10000;
Orbit: GiveOrbit(1.0,iLength)$



/* period 3 points */
e: f(f(f(z))) = z$

load (to_poly_solve);
s: to_poly_solve (e,z);

s: args(s); /*  https://stackoverflow.com/questions/12834709/create-a-union-into-a-list-in-maxima */
s:flatten(s);
s:map(rhs,s);

r:[];
for z in s do if (abs(abs(z) -1) < 0.001) then  r:cons(z,r);


cycle1:[];
z:r[1];
cycle1: cons(z, cycle1);
for i:1 thru 3 step 1 do (

	z:float(rectform(f(z))),
	cycle1:cons(z, cycle1)
	
);


cycle2:[];
z:r[2];
cycle2: cons(z, cycle2);
for i:1 thru 3 step 1 do (

	z:float(rectform(f(z))),
	cycle2:cons(z, cycle2)
	
);



Orbit : ToPoints(Orbit)$
r:ToPoints(r)$
cycle1 : ToPoints(cycle1)$
cycle2 : ToPoints(cycle2)$

load(draw); /* ( interface to gnuplot ) by Mario Rodriguez Riotorto http://www.telefonica.net/web2/biomates */

draw2d(
    title = "Two period 3 orbits for f(z)= rho * z^2 * (z-3)/(1-3*z))",
    terminal  = png,
     user_preamble = "set size square; set key right top;", /* 360/26=13.85  ; 360/(2*26)=6,923 */
    file_name = concat("~/Dokumenty/julia_blaszke/period3/maxima/roots/", concat("cycles_",string(iLength))),
    dimensions  = [600, 600],   /* Since Maxima 5.23, pic_width and pic_height are deprecated. */
    /*
    yrange = [-1.0,1.0],
    xrange = [-1.0,1.0],
    */
    xlabel     = "z.re ",
    ylabel     = "z.im",
    
    
    line_width        = 1,
    nticks = 50,
    color             = gray,
    transparent       = true,
    ellipse(0,0,1,1,0,360), /* unit circle */
    
    point_type    = filled_circle,
    
    
    
    point_size    = 0.9,
    key="critical orbit ",
    color             = black ,
    Orbit,
    
    point_size    = 1.2,
    key="critical point",
    color             = blue ,
    points([[1,0]]),
    
    point_type    = filled_circle,
    points_joined = true,
    point_size    = 0.9,
    key="cycle 1",
    color             = red ,
    cycle1,
    
    point_type    = filled_circle,
    points_joined = true,
    point_size    = 0.9,
    key="cycle 2",
    color             = blue ,
    cycle2
    

 );

Can you help me identify this (possibly) fractal? Started by username on Share a fractal	5 Replies 1073 Views	January 03, 2019, 12:02:36 PM by mrrudewords
when you're new to Mandelbulb3D and you are rendering your first animation Started by Nasko on Fractal Humor	10 Replies 3865 Views	April 29, 2020, 01:01:17 AM by Tas_mania
Article question: Julia sets you can trust Started by marcm200 on Fractal Mathematics And New Theories	13 Replies 1038 Views	August 04, 2019, 09:22:00 AM by marcm200
Can you run a Kalles Fraktaler job on a Linux Server? Started by Kartoffel on Noob's Corner	4 Replies 86 Views	April 21, 2022, 07:49:05 PM by Kartoffel
Have you seen this fractal before? Started by Chris_M_Thomasson on Fractal Mathematics And New Theories	7 Replies 285 Views	July 06, 2022, 01:25:05 AM by Chris_M_Thomasson

Author Topic: can you identify this fractal ? (Read 770 times)

Adam Majewski

can you identify this fractal ?

3DickUlus

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

marcm200

Re: can you identify this fractal ?

hgjf2

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

Alef

Re: can you identify this fractal ?

Alef

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

marcm200

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

marcm200

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

Adam Majewski

Re: can you identify this fractal ?

Similar Topics

Can you help me identify this (possibly) fractal?

when you're new to Mandelbulb3D and you are rendering your first animation

Article question: Julia sets you can trust

Can you run a Kalles Fraktaler job on a Linux Server?

Have you seen this fractal before?