Periodic points of complex quadratic mappings

Periodic points of complex quadratic mappings

This article describes periodic points of some complex quadratic map. This theory is applied in relation with the theories of Fatou and Julia sets.

Definitions

Let

:f_c(z)=z^2+c,

where z and c are complex-valued. (This f is the "complex quadratic mapping" mentioned in the title.) This article explores the "periodic points" of this mapping - that is, the points that form a periodic cycle when f is repeatedly applied to them.

f^{(k)} _c (z) is the k -fold compositions of f _c, with itself = iteration of function f _c,

f^{(k)} _c (z) = f_c(f^{(k-1)} _c (z))

then periodic points of complex quadratic mapping of period p are points z of dynamical plane such that :

z : f^{(p)} _c (z) = z

where p is the smallest positive integer.

We can introduce new function:

F_p(z,f) = f^{(p)} _c (z) - z

so periodic points are zeros of function F_p(z,f) :

z : F_p(z,f) = 0

which is polynomial of degree = 2^p

Stability of periodic points ( orbit )

The multiplier m(f,z_0)=lambda , of rational map f, at fixed point z_0, is defined as :

m(f,z_0)=lambda = egin{cases} f_c'(z_0), &mbox{if }z_0 e infty \ frac{1}{f_c'(z_0)}, & mbox{if }z_0 = infty end{cases}

where f_c'(z_0), is first derivative of f_c with respect to z, at z_0,.

Because multiplier is the same at all points of peiodic one can name it multiplier of periodic orbit.

Multiplier is:
*complex number,
*invariant under conjugation of any rational map at its fixed point [Alan F. Beardon, Iteration of Rational Functions, Springer 1991, ISBN 0-387-95151-2, p. 41]
*used to check stability of periodic (also fixed) points.


=Period-1 points (fixed points)=

Finite fixed points

Let us begin by finding all finite points left unchanged by 1 application of f. These are the points that satisfy f_c(z)=z. That is, we wish to solve

: z^2+c=z,

which can be rewritten

: z^2-z+c=0.

Since this is an ordinary quadratic equation in 1 unknown, we can apply the standard quadratic solution formula. Look in any standard mathematics textbook, and you will find that there are two solutions of Ax^2+Bx+C=0 are given by

: x=frac{-Bpmsqrt{B^2-4AC{2A}

In our case, we have A=1, B=-1, C=c, so we will write

: alpha_1 = frac{1-sqrt{1-4c{2} and alpha_2 = frac{1+sqrt{1-4c{2}.So for c in C setminus [1/4,+inf ] we have two finite fixed points alpha_1 , and alpha_2, .

Since : alpha_1 = frac{1}{2}-m and alpha_2 = frac{1}{2}+ m where m = frac{sqrt{1-4c{2}

then alpha_1 + alpha_2 = 1 ,.

It means that fixed points are symmetrical around z = 1/2,.

Complex dynamics

Here different notation is commonly used:

: alpha_c = frac{1-sqrt{1-4c{2} and eta_c = frac{1+sqrt{1-4c{2}.

Using Viète's formulas one can show that:

: alpha_c + eta_c = -frac{B}{A} = 1

Since derivative with respect to z is :

:P_c'(z) = frac{d}{dz}P_c(z) = 2z

then

:P_c'(alpha_c) + P_c'(eta_c)= 2 alpha_c + 2 eta_c = 2 (alpha_c + eta_c) = 2 ,

It implies that P_c , can have at most one attractive fixed point.

This points are distinguished by the facts that:
* eta_c , is :
**the landing point of external ray for angle=0 for c in M setminus left { frac{1}{4} ight }
**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) [ [http://www.ibiblio.org/e-notes/MSet/Attractor.htm Periodic attractor by Evgeny Demidov] ] .
* alpha_c , 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

pecial cases

An important case of the quadratic mapping is c=0. In this case, we get alpha_1 = 0 and alpha_2=1. In this case, 0 is a superattractive fixed point, and 1 belongs to the Julia set.


=Only one fixed point=

We might wonder what value c should have to cause alpha_1=alpha_2. The answer is that this will happen exactly when 1-4c=0. This equation has 1 solution: c=1/4 (in which case, alpha_1=alpha_2=1/2). This is interesting, since c=1/4 is the largest positive, purely-real value for which a finite attractor exists.

Infinite fixed point

We can extend complex plane mathbb{C} to Riemann sphere (extended complex plane) mathbb{hat{C by adding infinity

mathbb{hat{C = mathbb{C} cup { infty }

and extend polynomial f_c, such that f_c(infty)=infty,

Then infinity is :
*superattracting
*fixed point f_c(infty)=infty=f^{-1}_c(infty),of polynomial f_c, [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] .

Period-2 cycles

Suppose next that we wish to look at "period-2 cycles". That is, we want to find two points eta_1 and eta_2 such that f_c(eta_1) = eta_2, and f_c(eta_2) = eta_1.

Let us start by writing f_c(f_c(eta_n)) = eta_n, and see where trying to solve this leads.

: f_c(f_c(z)) = (z^2+c)^2+c = z^4 + 2z^2c + c^2 + c.,

Thus, the equation we wish to solve is actually z^4 + 2cz^2 - z + c^2 + c = 0.

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 alpha_1 and alpha_2, computed above. It is simple to see why this is; if these points are left unchanged by 1 application of f, then clearly they will be unchanged by 2 applications (or more).

Our 4th-order polynomial can therefore be factored in 2 ways :

first method

: (z-alpha_1)(z-alpha_2)(z-eta_1)(z-eta_2) = 0.,

This expands directly as x^4 - Ax^3 + Bx^2 - Cx + D = 0 (note the alternating signs), where

: D = alpha_1 alpha_2 eta_1 eta_2,

: C = alpha_1 alpha_2 eta_1 + alpha_1 alpha_2 eta_2 + alpha_1 eta_1 eta_2 + alpha_2 eta_1 eta_2,

: B = alpha_1 alpha_2 + alpha_1 eta_1 + alpha_1 eta_2 + alpha_2 eta_1 + alpha_2 eta_2 + eta_1 eta_2,

: A = alpha_1 + alpha_2 + eta_1 + eta_2.,

We already have 2 solutions, and only need the other 2. This is as difficult as solving a quadratic polynomial. In particular, note that

: alpha_1 + alpha_2 = frac{1-sqrt{1-4c{2} + frac{1+sqrt{1-4c{2} = frac{1+1}{2} = 1

and

: 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.

Adding these to the above, we get D = c eta_1 eta_2 and A = 1 + eta_1 + eta_2. Matching these against the coefficients from expanding f, we get

: D = c eta_1 eta_2 = c^2 + c and A = 1 + eta_1 + eta_2 = 0.

From this, we easily get :eta_1 eta_2 = c + 1 and eta_1 + eta_2 = -1.

From here, we construct a quadratic equation with A' = 1, B = 1, C = c+1 and apply the standard solution formula to get

: eta_1 = frac{-1 - sqrt{-3 -4c{2} and eta_2 = frac{-1 + sqrt{-3 -4c{2}.

Closer examination shows (the formulas are a tad messy) that :

f_c(eta_1) = eta_2 and f_c(eta_2) = eta_1

meaning these two points are the two halves of a single period-2 cycle.

econd method of factorization

(z^2+c)^2 + c -z = (z^2 + c - z)(z^2 + z + c +1 ) ,

The roots of the first factor are the two fixed points z_{1,2}, . They are repelling outside the main cardioid.

The second factor has two roots

z_{3,4} = -frac{1}{2} pm (-frac{3}{4} - c)^frac{1}{2} ,

These two roots form period-2 orbit. [ [http://www.ibiblio.org/e-notes/MSet/Attractor.htm Period 2 orbit by Evgeny Demidov] ]

pecial cases

Again, let us look at c=0. Then

: eta_1 = frac{-1 - isqrt{3{2} and eta_2 = frac{-1 + isqrt{3{2}

both of which are complex numbers. By doing a little algebra, we find | eta_1 | = | eta_2 | = 1. Thus, both these points are "hiding" in the Julia set.

Another special case is c=-1, which gives eta_1 = 0 and eta_2 = -1. This gives the well-known superattractive cycle found in the largest period-2 lobe of the quadratic Mandelbrot set.

Cycles for period>2

There is no general solution in radicals to polynomial equations of degree five or higher, so it must be computed using numerical methods.

References

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.
* [http://arxiv.org/abs/0802.2565 Gvozden Rukavina : Quadratic recurrence equations - exact explicit solution of period four fixed points functions in bifurcation diagram]


Wikimedia Foundation. 2010.

Игры ⚽ Нужно решить контрольную?

Look at other dictionaries:

  • Complex quadratic polynomial — A complex quadratic polynomial is a quadratic polynomial whose coefficients are complex numbers. Contents 1 Forms 2 Conjugation 2.1 Between forms 2.2 With doubling map …   Wikipedia

  • Quadratic equation — This article is about quadratic equations and solutions. For more general information about quadratic functions, see Quadratic function. For more information about quadratic polynomials, see Quadratic polynomial. In mathematics, a quadratic… …   Wikipedia

  • Quadratic function — A quadratic function, in mathematics, is a polynomial function of the form f(x)=ax^2+bx+c ,!, where a e 0 ,!. The graph of a quadratic function is a parabola whose major axis is parallel to the y axis.The expression ax^2+bx+c in the definition of …   Wikipedia

  • Quadratic — In mathematics, the term quadratic describes something that pertains to squares, to the operation of squaring, to terms of the second degree, or equations or formulas that involve such terms. Quadratus is Latin for square . Mathematics Algebra… …   Wikipedia

  • Periodic point — In mathematics, in the study of iterated functions and dynamical systems, a periodic point of a function is a point which returns to itself after a certain number of function iterations or a certain amount of time. Iterated functions Given an… …   Wikipedia

  • List of complex analysis topics — Complex analysis, traditionally known as the theory of functions of a complex variable, is the branch of mathematics that investigates functions of complex numbers. It is useful in many branches of mathematics, including number theory and applied …   Wikipedia

  • Orbit portrait — In mathematics, an orbit portrait is a combinatorial tool used in complex dynamics for understanding the behavior of one complex dimensional quadratic maps. In simple words one can say that it is : a list of external angles for which rays… …   Wikipedia

  • List of mathematics articles (P) — NOTOC P P = NP problem P adic analysis P adic number P adic order P compact group P group P² irreducible P Laplacian P matrix P rep P value P vector P y method Pacific Journal of Mathematics Package merge algorithm Packed storage matrix Packing… …   Wikipedia

  • Orbit (dynamics) — In mathematics, in the study of dynamical systems, an orbit is a collection of points related by the evolution function of the dynamical system. The orbit is a subset of the phase space and the set of all orbits is a partition of the phase space …   Wikipedia

  • External ray — In complex analysis, particularly in complex dynamics and geometric function theory, external rays are associated to a compact, full, connected subset K, of the complex plane as the images of radial rays under the Riemann map of the complement of …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”