Periodic points of complex quadratic mappings 869255 217149496 2008-06-04T20:36:53Z Mdd 113850 /* Further readings */ This article describes [[periodic point]]s of some [[Complex quadratic polynomial|complex quadratic map]]. This theory is applied in relation with the theories of [[Fatou set|Fatou]] and [[Julia set]]s. ==Definitions== Let :<math>f_c(z)=z^2+c\,</math> where <math>z</math> and <math>c</math> are [[Complex Number|complex-valued]]. (This <math>\ f</math> is the ''[[Complex quadratic polynomial|complex quadratic mapping]]'' mentioned in the title.) This article explores the ''[[periodic point]]s'' of this [[Map (mathematics)|map]]ping - that is, the points that form a periodic cycle when <math>\ f</math> is repeatedly applied to them. <math> \ f^{(k)} _c (z) </math> is the <math>\ k</math> -fold [[Function composition|compositions]] of <math>f _c\,</math> with itself = [[Iterated function|iteration of function]] <math>f _c\,</math> <math> \ f^{(k)} _c (z) = f_c(f^{(k-1)} _c (z))</math> then periodic points of [[Complex quadratic polynomial|complex quadratic mapping]] of [[Periodicity|period]] <math>\ p</math> are points <math> \ z</math> of [[phase space|dynamical plane]] such that : <math> \ z : f^{(p)} _c (z) = z</math> where <math>\ p</math> is the smallest positive integer. We can introduce new function: <math> \ F_p(z,f) = f^{(p)} _c (z) - z</math> so periodic points are zeros of function <math> \ F_p(z,f) </math> : <math> \ z : F_p(z,f) = 0</math> which is polynomial of degree <math> \ = 2^p</math> ==Stability of periodic points== The '''multiplier''' <math>m(f,z_0)=\lambda \,</math> of fixed point <math>z_0\,</math> is defined as <math> m(f,z_0)=\lambda = \begin{cases} f_c'(z_0), &\mbox{if }z_0\ne \infty \\ \frac{1}{f_c'(z_0)}, & \mbox{if }z_0 = \infty \end{cases} </math> where <math>f_c'(z_0)\,</math> is [[Complex quadratic polynomial| first derivative]] of <math> \ f_c</math> with respect to <math>z\,</math> at <math>z_0\,</math>. Multiplier is: *[[complex number]], *invariant under conjugation of any rational map at its fixed point<ref>Alan F. Beardon, Iteration of Rational Functions, Springer 1991, ISBN 0-387-95151-2, p. 41</ref> *used to check stability of periodic (also fixed) points. ==Period-1 points ([[Fixed point (mathematics)|fixed points]])== ===Finite fixed points=== {{Wikibooks|Fractals }} Let us begin by finding all [[finite]] points left unchanged by 1 application of <math>f</math>. These are the points that satisfy <math>\ f_c(z)=z</math>. That is, we wish to solve : <math>z^2+c=z\,</math> which can be rewritten : <math>\ z^2-z+c=0.</math> Since this is an ordinary quadratic equation in 1 unknown, we can apply [[Quadratic equation|the standard quadratic solution formula]]. Look in any standard mathematics textbook, and you will find that there are two solutions of <math>\ Ax^2+Bx+C=0</math> are given by : <math>x=\frac{-B\pm\sqrt{B^2-4AC}}{2A}</math> In our case, we have <math>A=1, B=-1, C=c</math>, so we will write : <math>\alpha_1 = \frac{1-\sqrt{1-4c}}{2}</math> and <math>\alpha_2 = \frac{1+\sqrt{1-4c}}{2}.</math> So for <math>c \in C \setminus [1/4,+\inf ]</math> we have two [[finite]] fixed points <math>\alpha_1 \,</math> and <math>\alpha_2\, </math>. Since : <math>\alpha_1 = \frac{1}{2}-m</math> and <math>\alpha_2 = \frac{1}{2}+ m</math> where <math>m = \frac{\sqrt{1-4c}}{2}</math> then <math>\alpha_1 + \alpha_2 = 1 \,</math>. It means that fixed points are symmetrical around <math>z = 1/2\,</math>. [[Image:fixed.gif|thumb|right|This image shows fixed points for c=i]] ====[[Complex dynamics]]==== Here different notation is commonly used: : <math>\alpha_c = \frac{1-\sqrt{1-4c}}{2}</math> and <math>\beta_c = \frac{1+\sqrt{1-4c}}{2}.</math> Using [[Viète's formulas]] one can show that: :<math> \alpha_c + \beta_c = -\frac{B}{A} = 1 </math> Since [[Complex_quadratic_polynomial#Derivative_with_respect_to_z|derivative with respect to z]] is : :<math>P_c'(z) = \frac{d}{dz}P_c(z) = 2z </math> then :<math>P_c'(\alpha_c) + P_c'(\beta_c)= 2 \alpha_c + 2 \beta_c = 2 (\alpha_c + \beta_c) = 2 \,</math> It implies that <math>P_c \,</math> can have at most one attractive fixed point. This points are distinguished by the facts that: * <math>\beta_c \,</math> is : **the landing point of [[external ray]] for angle=0 for <math>c \in M \setminus \left \{ \frac{1}{4} \right \}</math> **the most repelling fixed point, belongs to Julia set, ** the one on the right ( whenever fixed point are not symmetrical around the real axis), it is the extreme right point for connected Julia sets (except for cauliflower)<ref>[http://www.ibiblio.org/e-notes/MSet/Attractor.htm Periodic attractor by Evgeny Demidov]</ref>. * <math>\alpha_c \,</math> is: ** landing point of several rays ** is : ***attracting when c is in main cardioid of Mandelbrot set, then it is in interior of Filled-in Julia set, it means belongs to Fatou set ( strictly to basin of attraction of finite fixed point ) ***parabolic at the root point of the limb of Mandelbtot set ***repelling for other c values [[Image:Julia0bb.jpg|thumb| Fatou set for F(z)=z*z with marked fixed point]] ====Special cases==== An important case of the quadratic mapping is <math>c=0</math>. In this case, we get <math>\alpha_1 = 0</math> and <math>\alpha_2=1</math>. In this case, 0 is a superattractive [[Fixed point (mathematics)|fixed point]], and 1 belongs to the [[Julia set]]. ====Only one [[Fixed point (mathematics)|fixed point]]==== We might wonder what value <math>c</math> should have to cause <math>\alpha_1=\alpha_2</math>. The answer is that this will happen exactly when <math>1-4c=0</math>. This equation has 1 solution: <math>c=1/4</math> (in which case, <math>\alpha_1=\alpha_2=1/2</math>). This is interesting, since <math>c=1/4</math> is the largest positive, purely-real value for which a finite attractor exists. ===Infinite fixed point=== We can extend [[complex plane]] <math>\mathbb{C}</math> to [[Riemann sphere (extended complex plane)]] <math>\mathbb{\hat{C}}</math> by adding [[Point at infinity|infinity]] <math>\mathbb{\hat{C}} = \mathbb{C} \cup \{ \infty \}</math> and extend [[Complex quadratic polynomial| polynomial]] <math>f_c\,</math> such that <math>f_c(\infty)=\infty\,</math> Then [[Point at infinity|infinity]] is : *superattracting *fixed point <math>f_c(\infty)=\infty=f^{-1}_c(\infty)\,</math> of [[Complex quadratic polynomial| polynomial]] <math>f_c\,</math><ref>R L Devaney, L Keen (Editor}: Chaos and Fractals: The Mathematics Behind the Computer Graphics. Publisher: Amer Mathematical Society July 1989, ISBN-10: 0821801376 , ISBN-13: 9780821801376</ref>. ==Period-2 cycles== Suppose next that we wish to look at ''period-2 cycles''. That is, we want to find two points <math>\beta_1</math> and <math>\beta_2</math> such that <math>f_c(\beta_1) = \beta_2</math>, and <math>f_c(\beta_2) = \beta_1</math>. Let us start by writing <math>f_c(f_c(\beta_n)) = \beta_n</math>, and see where trying to solve this leads. : <math>f_c(f_c(z)) = (z^2+c)^2+c = z^4 + 2z^2c + c^2 + c.\,</math> Thus, the equation we wish to solve is actually <math>z^4 + 2cz^2 - z + c^2 + c = 0</math>. This equation is a polynomial of degree 4, and so has 4 (possibly non-distinct) solutions. ''However'', actually, we already know 2 of the solutions. They are <math>\alpha_1</math> and <math>\alpha_2</math>, computed above. It is simple to see why this is; if these points are left unchanged by 1 application of <math>f</math>, then clearly they will be unchanged by 2 applications (or more). Our 4th-order polynomial can therefore be factored in 2 ways : ===first method=== : <math>(z-\alpha_1)(z-\alpha_2)(z-\beta_1)(z-\beta_2) = 0.\,</math> This expands directly as <math>x^4 - Ax^3 + Bx^2 - Cx + D = 0</math> (note the alternating signs), where : <math>D = \alpha_1 \alpha_2 \beta_1 \beta_2\,</math> : <math>C = \alpha_1 \alpha_2 \beta_1 + \alpha_1 \alpha_2 \beta_2 + \alpha_1 \beta_1 \beta_2 + \alpha_2 \beta_1 \beta_2\,</math> : <math>B = \alpha_1 \alpha_2 + \alpha_1 \beta_1 + \alpha_1 \beta_2 + \alpha_2 \beta_1 + \alpha_2 \beta_2 + \beta_1 \beta_2\,</math> : <math>A = \alpha_1 + \alpha_2 + \beta_1 + \beta_2.\,</math> We already have 2 solutions, and only need the other 2. This is as difficult as solving a quadratic polynomial. In particular, note that : <math>\alpha_1 + \alpha_2 = \frac{1-\sqrt{1-4c}}{2} + \frac{1+\sqrt{1-4c}}{2} = \frac{1+1}{2} = 1</math> and : <math>\alpha_1 \alpha_2 = \frac{(1-\sqrt{1-4c})(1+\sqrt{1-4c})}{4} = \frac{1^2 - (\sqrt{1-4c})^2}{4}= \frac{1 - 1 + 4c}{4} = \frac{4c}{4} = c.</math> Adding these to the above, we get <math>D = c \beta_1 \beta_2</math> and <math>A = 1 + \beta_1 + \beta_2</math>. Matching these against the coefficients from expanding <math>f</math>, we get : <math>D = c \beta_1 \beta_2 = c^2 + c</math> and <math>A = 1 + \beta_1 + \beta_2 = 0.</math> From this, we easily get : <math>\beta_1 \beta_2 = c + 1</math> and <math>\beta_1 + \beta_2 = -1</math>. From here, we construct a quadratic equation with <math>A' = 1, B = 1, C = c+1</math> and apply the standard solution formula to get : <math>\beta_1 = \frac{-1 - \sqrt{-3 -4c}}{2}</math> and <math>\beta_2 = \frac{-1 + \sqrt{-3 -4c}}{2}.</math> Closer examination shows (the formulas are a tad messy) that : <math>f_c(\beta_1) = \beta_2</math> and <math>f_c(\beta_2) = \beta_1</math> meaning these two points are the two halves of a single period-2 cycle. ===Second method of factorization<ref>[http://www.ibiblio.org/e-notes/MSet/Attractor.htm Period 2 orbit by Evgeny Demidov]</ref>=== <math>(z^2+c)^2 + c -z = (z^2 + c - z)(z^2 + z + c +1 ) \,</math> The roots of the first factor are the two fixed points <math>z_{1,2}\,</math> . They are repelling outside the main cardioid. The second factor has two roots <math>z_{3,4} = -\frac{1}{2} \pm (-\frac{3}{4} - c)^\frac{1}{2} \,</math> These two roots form period-2 orbit. ====Special cases==== Again, let us look at <math>c=0</math>. Then : <math>\beta_1 = \frac{-1 - i\sqrt{3}}{2}</math> and <math>\beta_2 = \frac{-1 + i\sqrt{3}}{2}</math> both of which are complex numbers. By doing a little algebra, we find <math>| \beta_1 | = | \beta_2 | = 1</math>. Thus, both these points are "hiding" in the Julia set. Another special case is <math>c=-1</math>, which gives <math>\beta_1 = 0</math> and <math>\beta_2 = -1</math>. This gives the well-known superattractive cycle found in the largest period-2 lobe of the quadratic Mandelbrot set. ==Cycles for period>2== [[Abel–Ruffini theorem|There is no general solution]] in [[Radical (mathematics))|radicals]] to polynomial equations of degree five or higher, so it must be computed using [[Root-finding algorithm|numerical methods]]. ==References== {{Reflist}} == Further reading == *Alan F. Beardon, Iteration of Rational Functions, Springer 1991, ISBN 0-387-95151-2 *Michael F. Barnsley (Author), Stephen G. Demko (Editor), Chaotic Dynamics and Fractals (Notes and Reports in Mathematics in Science and Engineering Series) Academic Pr (April 1986), ISBN-10: 0120790602 *[http://www.math.sunysb.edu/cgi-bin/thesis.pl?thesis02-3 Wolf Jung : Homeomorphisms on Edges of the Mandelbrot Set. Ph.D. thesis of 2002] *[http://hdl.handle.net/10090/3895| The permutations of periodic points in quadratic polynominials by J Leahy] ==External links== *[http://cosinekitty.com/mandel_orbits_analysis.html ''Algebraic solution of Mandelbrot orbital boundaries'' by Donald D. Cross ] *[http://www.mrob.com/pub/muency/brownmethod.html ''Brown Method'' by Robert P. Munafo] *[http://arxiv.org/abs/hep-th/0501235 arXiv:hep-th/0501235v2] V.Dolotin, A.Morozov: ''Algebraic Geometry of Discrete Dynamics''. The case of one variable. [[Category:complex analysis]] [[Category:Fractals]] [[Category:Limit sets]] [[Category:Dynamical systems]]