Research Article  Open Access
Wang Xingyuan, He Yijie, Sun Yuanyuan, "Accurate Computation of Periodic Regions' Centers in the General MSet with Integer Index Number", Discrete Dynamics in Nature and Society, vol. 2010, Article ID 653816, 12 pages, 2010. https://doi.org/10.1155/2010/653816
Accurate Computation of Periodic Regions' Centers in the General MSet with Integer Index Number
Abstract
This paper presents two methods for accurately computing the periodic regions' centers. One method fits for the general Msets with integer index number, the other fits for the general Msets with negative integer index number. Both methods improve the precision of computation by transforming the polynomial equations which determine the periodic regions' centers. We primarily discuss the general Msets with negative integer index, and analyze the relationship between the number of periodic regions' centers on the principal symmetric axis and in the principal symmetric interior. We can get the centers' coordinates with at least 48 significant digits after the decimal point in both real and imaginary parts by applying the Newton's method to the transformed polynomial equation which determine the periodic regions' centers. In this paper, we list some centers' coordinates of general Msets' periodic regions for the index numbers , all of which have highly numerical accuracy.
1. Introduction
According to the idea of complex dynamic system theory presented by Julia and Fatou, the famous mathematician Mandelbrot constructed and studied the Msets of complex mapping utilizing computer graphics technologies [1]. During the last 20 years, people have researched the embeddedlayer relationship and distribution of the bifurcation sequence and topological rule of periodic trajectories in the general Msets with and found there existed orderly structure within the Msets [2–16]. For example, Álvarez et al. studied the location and number of each periodic region in Msets [9]; Buchanan et al. studied the location of periodic region of the general Msets with [11]; Geum and Kim analyzed the quantitative relationship of each periodic region in the general Mandelbrot sets with positive integer index number, and calculated the coordinates of periodic regions’ centers [15]; The author studied the structure and distribution of the general Msets with integer index number [16].
Msets consist of differentperiod regions which constitute the fractal structures of the Msets. The analysis on the stability of the maps uncovered new and unexpected algebraic properties of the periodic regions. The centers of the periodic regions are determined by the transformed polynomials we worked on. The motivation for computing the periodic regions’ centers is provided by the need to consider the locations of the periodic regions and fractal structures of the Msets, which are useful to understand the inner infinite structures of the Msets.
On the basis of above research, we study the periodic region’s centers in the general Msets with integer index number, determine the relationship between each periodic region’s number in the general Msets with positive index number and negative index number, and then present a new method of calculating the coordinates of periodic regions’ centers in the general Msets with negative integer index number. Our research has good prospects in physics, information science, and other fields.
2. Periodic Region Theory of General MSets
Definition 2.1 (see [15]). Let for with , then the general Msets is defined to be the set
Definition 2.2 (see [15]). The sets defined by for are called the rays of symmetry. The set is called the principal ray of symmetry. As is shown in Figure 1(b).
(a)
(b) Primary symmetric axis and section
(c)
(d) Partial enlargement of little block in Figure 1(c)
(e)
(f) Partial enlargement of little block in Figure 1(e)
Definition 2.3. The set is called the principal symmetric sector, as is shown in Figure 1(b).
Theorem 2.4. In the parameter plane, M is symmetric about .
Proof. Let with . For all , there exists the following recursive relations: where . Using the mathematical induction, for all and , exists. So we obtain since with , we have whenever . This completes the proof.
Definition 2.5 (see [15]). An attracting periodic region is denoted by and is defined as a region of the set
Definition 2.6. When is a positive, if satisfying and , then is called the center of an attracting periodic region.
Definition 2.7. When is a negative, is the center of a stable periodic region only if belongs to this periodic region and satisfies .
3. Calculation of Periodic Regions’ Centers in General Msets with Integer Index Number
3.1. Calculation Method
The center of periodic region can be located by numerically solving the governing equation which is a polynomial of . The equation can be written as where is a recursive function defined as The in (3.2) has a degree of and thus will encounter a difficulty in obtaining accurate solutions as and increase. But the transformation reduces the degree by .
Let .
() When is a positive integer, (3.2) can be written as where and .
Definition 3.1 (see [15]). Let is a integer smaller than and satisfying , let
() When is a negative integer, (3.2) can be written as where and .
So we can solve the roots of (3.4) or (3.5) instead of (3.1). If is a positive integer, we can use (3.4) to calculate the center of periodic region. If is a negative integer, we can use (3.5) to calculate the center of periodic region.
3.2. Numerical Algorithm
Through solving the roots of or , we can solve the roots of . Let with be a root of or. The transformation with the symmetry of Msets yields where . If , then; If , then .
Among these values, we select the ones in the interior of the primary symmetric and the ones on the primary symmetric axis . Then by the rotation symmetry, we can get the center’s coordinate of each periodic region.
Using the Newton method, the following algorithm can locate the center of periodic region [15].()Set , and (number of precision digits).()Construct .()Construct SolveCenter(k) that does the following. () Using the function in to find the approximate roots of ; () Select the roots with nonnegative imaginary parts, and compute and .()Construct Newton that does the following. () Take the initial values of about three decimal digits of accuracy from the results of to precede the Newton sequence ; () Set the maximum iteration number ; () Set satisfying and .()Construct DefCenter(k) finding coordinate of the periodic region’s center: () Reset to a higher number NtDigits; () If , then ; If , then call Newton to get refined roots ; () Compute , , and .()Call SolveCenter(k), Compute residual error and .()If is not met, increase Digits and call SolveCenter ; If is met, then do the following. () Increase NtDigits; () Call DefCenter(k) to get refined roots , and calculate ; () If is not met, increase and call DefCenter(k); If is met, the check the convergence of the sequence , where : If is not convergent, increase NtDigits and call DefCenter ; If is convergent, then accept as the desired solution and terminate the entire procedure.
The above algorithm is appropriate for as a positive integer; if is a negative integer, we can take instead of . The algorithm is achieved by the . Tables 1 and 2 list typical coordinates of periodic region’ centers for and ; the accuracy is of 48 precision digits, but only the first 40 precision digits were printed for the tabulation. Table 3 shows the residual error defined by and the values of asymptotic error . The parameters are defined as ,, and in the experiment.



3.3. Numerical Results
Now, we study the relationship between the numbers of roots on and in the interior of . According to the Rotation symmetry, we consider that is and only is on and in the interior of .
Let , we select roots having only nonnegative imaginary part and suppose is the number of all such roots. Let denotes the number of roots for in the complex plane , denotes the number of centers lying in the interior of , denotes the number of centers lying on , and denotes the total number of periodic regions’ centers lying in the complex plane . Then, the following relation holds [15]:
Next we discuss the calculation of periodic regions’ centers and the relation among ,, and when is a positive integer and a negative integer, respectively.
() If is a positive integer,
The general Msets with positive integer index number are similar to flowers combined with major petals, and stable region is embedded in unstable region. If , the complex mapping degenerates into linear mapping, the general Msets is only a circular which is a trivial structure without complex boundary and selfsimilar structure. Once deviating from , the complex structure appear immediately [3].
Theorem 3.2. If is a positive integer, then the relation (or ) yields the following [15].()If , then and for any ;()If , then and for all odd ;()If , then and for all even .()If , then and for all even ;()If , then and for all even ;()If , then , respectively, for all even .
We can get the following from Theorem 3.2. Statement () implies that a twoperiodic region’s center lies on ; Statements () and () show that only a twoperiodic region’s center lies on when is odd. Statements () and () show that every periodic region’s center with lies on when is even. Statements () and () show that for any even integer one threeperiodic and two fourperiodic periodic region’s centers lie on , respectively.
() If is a negative integer
The general Msets with negative integer index number have planetary configuration consisting of a central planet with major satellite structures, and the unstable region is embedded in the stable region, as is shown in Figure 1. The numbers and locations of the periodic region’s centers in the general Msets with negative integer index number can be calculated using the same method as being a positive integer.
Theorem 3.3. If is a negative integer, then the relation (or )) yields the following:()If and , then ;()If and , then , namely with negative integeris equal to with ;()If or , then for all ;()If , then and for being even; , for being odd;()If , then for all .
We can get the following from Theorem 3.2. Statement () implies that one twoperiodic and one threeperiodic region’s center lies on ; Statement () state that two fourperiodic region’s center lies on when is even, as is shown in Figure 1(d); When is odd, only one fourperiodic region’s center lies on , as is shown in Figure 1(f). Table 4 shows the relationship among ,, and when is a negative integer, which testifies Theorem 3.3 and (3.7). The numbers in Table 4 are . The corresponding and satisfy Theorem 3.2.

4. Calculation of Periodic Regions’ Centers in General Msets with Negative Integer Index Number
Let , then where and . The coordinates of periodic regions’ centers can be obtained by solving the equation with satisfying .
Tables 5 and 6 list the expressions of centers’ coordinates from twoperiodic to fiveperiodic region.


For the Mset constructed from the complex mapping (shown in Figure 1(a)), oneperiodic region’s center is for . Twoperiodic region’s center reaches the oneperiodic region’s center by one iteration, so the origin is the center of twoperiodic region, that is, in Table 5 deduces . Similarly, threeperiodic region’s centers can be located by the roots of .
For the Mset constructed from the complex mapping (shown in Figure 1(c)), oneperiodic region is a huge area whose boundary is defined by . Twoperiodic stable region divides the area outside of oneperiodic region into three parts. Oneperiodic region’s center is . Twoperiodic region’s center is the origin. Similarly, threeperiodic region’s centers can be located by the roots of .
The discussion above indicates that periodic region’s centers can be located by the roots of . The results of centers’ coordinates are almost the same as the results obtained by the first method described in Section 2 when is a negative integer.
5. Conclusions
() In this paper, we proposed two methods for calculating the periodic regions’ centers of the general Msets. The first method fits for calculating the periodic regions’ centers in the general Msets with integer index number, which is to transform the polynomial equation that governs the periodic regions’ centers, obtain high precision of the coordinates by the simple method, and analyze the relationship between the number of each periodic regions’ centers on the principal symmetric axis and in the principal symmetric interior, then comparatively analyze the relation of periodic regions’ number in general Msets with opposite integer index number. The second methods as discussed in Section 3 suits for calculating the periodic regions’ centers in the general Msets with negative integer index number, which also transforms the polynomial equation. The results of centers’ coordinates obtained by the second method are almost the same as that of the first method when is a negative integer which is described in Section 2.
() The investigation of the periodic regions’ centers of the Msets can help us explore the distribution of periodic regions of the Msets, which can further help us to study the fractal structures of the Msets. The centers of the periodic regions are located as the roots of certain polynomials, which are shown to coincide with solutions of the Douady and Hubbard formula [17]. In addition, the methods we proposed are helpful for solving polynomial equations, especially of high degree.
() This research is some inspiration for the people studying on the difficult problems in their professional and interdisciplinary fields.
As a classical example of physics, Brownian movement is the most simple and typical random movement. The Langevin equation can depict the rule of a charged particle under the circularly successive influence of the impulse functions. However, it is difficult to visually depict the trajectory and dynamics of these systems with many random variables. If we construct the complex general Msets using the rules of Langevin equation, the fractal structure characteristics of the general Msets can reveal the changing rule of the particle velocity visually [18]. This study makes it possible to depict complex Brownian movement more accurately.
In addition, the theories of Msets have potential applications on image processing. We have known that general Msets are illustrated dictionary of the corresponding general Julia sets [16], which means a single point on the Msets can represent the huge amount image data of the general Julia sets with manifold shapes and complicated structures. This research provides the technology support for determining rapidly the coordinates of the points on the Msets by the Julia images. On the basis of the above research results, future work includes establishing dictionary of fractal compression and studying the corresponding coding algorithm to improve the transmission and memory of the information, which could provide the new thoery for fractal compression technology.
Acknowledgment
This research is supported by the National Natural Science Foundation of China (nos. 60573172, 60973152), the Superior University Doctor Subject Special Scientific Research Foundation of China (no. 20070141014), and the Natural Science Foundation of Liaoning province (no. 20082165).
References
 B. B. Mandelbrot, The Fractal Geometry of Nature, Freeman, San Francisco, Calif, USA, 1982. View at: Zentralblatt MATH  MathSciNet
 H. O. Peitgen and D. Saupe, The Science of Fractal Images, Springer, New York, NY, USA, 1988. View at: Zentralblatt MATH  MathSciNet
 X. Y. Wang, The Fractal Mechanism of Generalized MJ Set, Dalian University of Technology Press, Dalian, China, 2002. View at: Zentralblatt MATH
 A. K. Dewdney, “A computer microscope zooms in fro a look at the most complex object in mathematics,” Scientific American, vol. 255, no. 7, pp. 16–32, 1985. View at: Google Scholar
 C. A. Pickover and E. Khorasani, “Computer graphics generated from the iteration of algebraic transformations in the complex plane,” Computers & Graphics, vol. 9, no. 1, pp. 147–151, 1985. View at: Publisher Site  Google Scholar
 A. K. Dewdney, “Beauty and profundity, the Mandelbrot set and a flock of its cousins called Julia,” Scientific American, vol. 257, no. 5, pp. 118–122, 1987. View at: Publisher Site  Google Scholar
 M. Romera, G. Pastor, and F. Montoya, “Misiurewicz points in onedimensional quadratic maps,” Physica A, vol. 232, no. 12, pp. 517–535, 1996. View at: Publisher Site  Google Scholar
 G. Pastor, M. Romera, and F. Montoya, “On the calculation of Misiurewicz patterns in onedimensional quadratic maps,” Physica A, vol. 232, no. 12, pp. 536–553, 1996. View at: Publisher Site  Google Scholar
 G. Álvarez, M. Romera, G. Pastor, and F. Montoya, “Determination of Mandelbrot set's hyperbolic component centres,” Chaos, Solitons & Fractals, vol. 9, no. 12, pp. 1997–2005, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. Y. Wang, X. Liu, W. Zhu, and S. Gu, “Analysis of cplane fractal images from $z\leftarrow {z}^{a}+c$ for $a<0$,” Fractals, vol. 8, no. 3, pp. 307–314, 2000. View at: Google Scholar  MathSciNet
 W. Buchanan, J. Gomatam, and B. Steves, “Generalized Mandelbrot sets for meromorphic complex and quaternionic maps,” International Journal of Bifurcation and Chaos, vol. 12, no. 8, pp. 1755–1777, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. C. Sasmor, “Fractals for functions with rational exponent,” Computers & Graphics, vol. 28, no. 4, pp. 601–615, 2004. View at: Publisher Site  Google Scholar
 M. Romera, G. Pastor, G. Álvarez, and F. Montoya, “External arguments of Douady cauliflowers in the Mandelbrot set,” Computers & Graphics, vol. 28, no. 3, pp. 437–449, 2004. View at: Publisher Site  Google Scholar
 G. Pastor, M. Romera, G. Álvarez, and F. Montoya, “Chaotic bands in the Mandelbrot set,” Computers & Graphics, vol. 28, no. 5, pp. 779–784, 2004. View at: Publisher Site  Google Scholar
 Y. H. Geum and Y. I. Kim, “Accurate computation of component centers in the degree$n$ bifurcation set,” Computers & Mathematics with Applications, vol. 48, no. 12, pp. 163–175, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 X. Y. Wang and P. J. Chang, “Research on fractal structure of generalized MJ sets utilized Lyapunov exponents and periodic scanning techniques,” Applied Mathematics & Computation, vol. 175, no. 2, pp. 1007–1025, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Douady and J. H. Hubbard, “On the dynamics of polynomiallike mappings,” Annales Scientifiques de l'École Normale Supérieure, vol. 18, no. 4, pp. 287–343, 1985. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 X. Y. Wang, W. Liu, and X. Yu, “Research on brownian movement based on generalized MandelbrotJulia sets from a class complex mapping system,” Modern Physics Letters B, vol. 21, no. 20, pp. 1321–1341, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2010 Wang Xingyuan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.