On the Boundedness and Symmetry Properties of the Fractal Sets Generated from Alternated Complex Map

A complex map can give rise to two kinds of fractal sets: the Julia sets and the parameters sets (or the connectivity loci) which represent different connectivity properties of the corresponding Julia sets. In the significative results of (Int. J. Bifurc. Chaos, 2009, 19:2123–2129) and (Nonlinear. Dyn. 2013, 73:1155–1163), the authors presented the two kinds of fractal sets of a class of alternated complex map and left some visually observations to be proved about the boundedness and symmetry properties of these fractal sets. In this paper, we improve the previous results by giving the strictly mathematical proofs of the two properties. Some simulations that verify the theoretical proofs are also included.


Introduction
Julia set [1], one of the most attractive fractal sets, is proposed by Gaston Julia when he studied the following complex quadratic polynomial Since it was proposed, Julia set has attracted significant interests in the topics of property analysis [2][3][4][5], superior iteration [6][7][8], time-delay [9], control [10,11] and so forth.Due to its good explanation for lots of complicated phenomena, Julia set has also been widely applied in the fields of physics [12], biology [13], image cryptography [14] and so on.
As mentioned in the research of [15,16], it may not be quite accurate to depict the evolution of the natural process mentioned above just by employing a unique system.In some cases, the alternated iterations can be used to have a better understanding of the natural process.In particular, for system (1), Danca et al. [17,18] proposed its alternated version as follows: where z n , c 1 = p + iq, c 2 = r + is ∈ C, C = (c 1 , c 2 ) T (T means transposition).
Based on the Fatou-Julia theorem [19] (see Theorem 2), the Julia set (see Definition 1) generated from P c 1 ,c 2 , which is denoted as J(P c 1 ,c 2 ), can be connected, disconnected and totally disconnected via choosing different parameters c 1 , c 2 [17].To test numerically which values of the parameters c 1 , c 2 give rise to connected, disconnected, or totally disconnected J(P c 1 ,c 2 ), the auxiliary complex quadratic polynomial Q c 1 ,c 2 was proposed in [17]: where w n = x n + iy n .Furthermore, they proved that the Julia sets generated from P c 1 ,c 2 and Q c 1 ,c 2 are the same with given c 1 , c 2 .Thus the connectivity of J(P c 1 ,c 2 ) relies on the critical orbits of system (3).In [18], the graphical results of the connectivity loci which represent the three connectivity states of the corresponding J(P c 1 ,c 2 ) were unveiled by analyzing the critical orbits of (3).
In concrete terms, the fractal sets which represent the connected, disconnected, totally disconnected Julia sets were respectively denoted as the Connectedness-Locus, Disconnectedness-Locus and Totally connectedness-Locus [17,18].Based on these graphical results, some observations about the boundedness and symmetry properties of these connectivity loci were summarized in Section 4 of [18].These unproved properties are just the original purpose and the starting point of this article.
Actually, there has been some research on the boundedness and symmetry properties for the fractal sets of complex maps which have only one critical point (see the classical M-J sets [2,20,21] and the generalized M-J sets [4,22,23]).For the systems like (3) which have more than one critical point, the results mainly focused on the properties of their Julia sets (see the Julia sets of cubic polynomials [3] and the superior Julia set [6]).To our knowledge, few research papers have addressed the problems about the boundedness and symmetry properties of their connectivity loci (or Mandelbrot set) since all the critical orbits should be considered.
Thus, as a supplementary research of [17,18], the purpose and achievement of this paper are no other than giving the mathematical proof about the two properties of the Connectedness-Locus and the Julia sets generated from (2) via analyzing all the critical orbits of system (3).
The remainder of this paper is outlined as follows.Section 2 recalls some definitions and theorems about the alternated complex map.In Section 3, the boundedness properties of the alternated Julia sets and the Connectedness-Locus are proposed and proved.In Section 4, the symmetry properties of the 3-D slice of Connectedness-Locus and the alternated Julia sets are analyzed in two theorems.Simulations are given to support the validity of the results.Finally, Section 5 concludes the paper.

The Fractal Sets Generated from Alternated Complex Map
In this section, some necessary definitions and theorems are recalled and given.Definition 1. [17] The filled Julia set of the system (2), denoted as K(P c 1 ,c 2 ), is the set of all the values of initial conditions z 0 such that lim where P n (•) represents the n-th iteration of the initial point z 0 .The Julia set of P c 1 ,c 2 which is the boundary of the filled Julia set is denoted as Theorem 1. [17,18] J(P c 1 ,c 2 ) and J(Q c 1 ,c 2 ) are the same for given c 1 and c 2 parameter values.
Theorem 2. [19] The connectivity properties of the Julia set for a complex polynomial of degree d > 2 can be identified based on the following cases: • The Julia set is connected if and only if all the critical orbits are bounded.
• The Julia set is totally disconnected, a cantor set, if (but not only if) all the critical orbits are unbounded.
• For a polynomial with at least one critical orbit unbounded, the Julia set is totally disconnected if and only if all the bounded critical orbits are aperiodic.
Then from Theorem 1 and Theorem 2, the connectivity of the J(P c 1 ,c 2 ) or K(P c 1 ,c 2 ) (J and K have the same connectivity properties [18]) relies on the boundedness of the three critical points of The following definition is given.Definition 2. [18] The Connectedness-Locus of the alternated complex map (2), denoted as M e (P c 1 ,c 2 ), is composed by all the initial constant c 1 , c 2 such that the orbits of 0 and √ −c 1 are bounded (since ± √ −c 1 behave the same boundedness properties).
In [17,18], M e (P c 1 ,c 2 ) is named as "Connectedness-Locus" or "Connected zone".The reason we denote it as M e (P c 1 ,c 2 ) in this paper is that it has the same efficacy with the classical Mandelbrot set [20] which gives rise to the connectedness property.Thus, in this paper, we could also call M e (P c 1 ,c 2 ) the Mandelbrot-efficacy set.
As shown in Definition 2, M e (P c 1 ,c 2 ) is formed by the initial constants c 1 , c 2 , which illustrates that it belongs to a four-dimensional space (p, q, r, s).By fixing q 0 = 0, the 3-D slice of M e (P c 1 ,c 2 ) is shown in Figure 1a (also see the red part of Figure 1 in [18]).The internal slice of M e (P c 1 ,c 2 ) q 0 =0 with p 0 = 0 is revealed in Figure 1b.The partial enlarged views of Figure 1b are successively shown in Figure 1c,d which fully display the beautiful and complex fractal structure of M e (P c 1 ,c 2 ).

Boundedness of the Fractal Sets Generated from Alternated Complex Map
In this section, by analysing the orbits of the three critical points of system (3), the boundedness properties of J(P c 1 ,c 2 ) and M e (P c 1 ,c 2 ) are given in the next two theorems.
where R is the same as Theorem 3.
Proof.Derive the max{|c 1 |, |c 2 |, R} in Theorem 3 into three cases as follows: ( 1 ): For case 1 , one knows that |c 1 |>|c 2 | and |c 1 |>T, the orbit of the critical value 0 satisfies that Then, |Q(0)| satisfies the condition of Theorem 3, it is clear that the orbit of 0 is diverging.
• Theorem 4 becomes: M(P c ) ⊂ {c |c| ≤ 2}.J(P c ) and M(P c ) mean the Julia set and Mandelbrot set of P c .The conclusions in this remark are consistent with the classical results [20].

Symmetry of the Fractal Sets Generated from Alternated Complex Map
Symmetry [2,4,23] is an important property of fractal sets.This section focuses on the symmetry properties of the Julia set J(P c 1 ,c 2 ) and the 3-D slice of the Mandelbrot-efficacy set M e (P c 1 ,c 2 ) q 0 =0 .The following two theorems are given and proved.
Theorem 5.The Mandelbrot-efficacy set M e (P c 1 ,c 2 ) q 0 =0 is symmetric around the pr-plane.
Proof.For the critical points 0, ± √ −c 1 , we give them a unitary expression w 0 .Then by choosing any two parameters values (c 1 , c 2 ) T = (p 0 , 0, r 0 , s 0 ) T , (c 1 , c 2 ) T = (p 0 , 0, r 0 , −s 0 ) T , the orbits of all the three critical points w 0 preserve complex conjugacy when n > 1.That is: Thus, the critical orbits w n behave the same boundedness property by choosing any two parameters which are symmetric around the pr-plane.Then, based on Definition 2, the two parts of M e (P c 1 ,c 2 ) q 0 =0 divided by the pr-plane are symmetric.
Proof.When c 1 , c 2 are given, one gets a corresponding Julia set.For any two initial points w 0 and −w 0 , they map to the same w 1 = (w 2 0 + c 1 ) 2 + c 2 .Thus for any two points which are centrally symmetric around the origin, their orbits behave the same boundedness properties.From Definition 1, one knows that J(P c 1 ,c 2 ) is centrally symmetric around the origin.
With the help of the graphical method proposed in [18], we supplement the Disconnectedness-Locus (gray zone) and the Totally disconnectedness-Locus (white zone) to M e (P c 1 ,c 2 ) q 0 =0,s 0 =−0.75 .Then three different points are picked out and their corresponding filled Julia sets are shown in Figure 4.
The simulations of both Figure 4 and Figure 2 verify the correctness of Theorem 6.

Conclusions
Fractals is one of the hottest topics in the fields of nonlinear dynamics.The properties of fractal sets are necessary to be investigated and favorable for the practical applications.This work focus on the boundedness and symmetry properties of the fractal sets generated from a class of alternated complex map.This is an important supplementary theoretical research for the existing related research and could represent a guide for the future research on the properties analysis about the complex maps with higher order or more complicated alternate rules.

Thus, |Q 2
( √ −c 1 )| satisfies the condition of Theorem 3, one can know that the orbit of √ −c 1 is diverging.