MAT 335 Homework #6

* * * This assignment does not have to be handed in * * *
Solutions will be available Wednesday, April 24.



Last updated: April 16, pm.

Some useful formula



  • (1*)
    Use the applet "Logistic Movie" to determine the period-doubling points b1 (period 6), b2 (period 12) and b3 (period 24) in the period 3 window near a = 3.83 (see Figure 11.41). Use these values to estimate the Feigenbaum constant delta (see pages 611 and 636) for this period-doubling scenario. (You will have to estimate the bi to 4 decimal places.)
  • (2*)
    A dynamical system depends on a parameter a. Initially, you observe a steady state (i.e., a period 1 orbit). As a increases you observe a period 2 oscillation appearing at a = a1 = 7. Then at a = a2 = 10 you observe that the period 2 orbits splits into a period 4 orbit. As a continues to increase a series of period-doublings occurs. Assuming Universality, at what a value would you expect to observe the onset of chaos? ('Assuming Universality' means assuming that the system will go through a series of period-doubling bifurcations as the parameter a changes, and that the distance (in a) between bifurcations is given by the Feigenbaum constant delta.)
  • (3) (hand in part (d) )
    Answer the following questions about the final state diagram of the logistic equation f = f a (see for example, Figure 11.5 in the text).
  • (4*)
    Write down the relative ordering of the following numbers using Charkovsky's ordering of the integers: 56, 31, 128, 160. That is, which one of these numbers is the 'smallest', the next 'smallest', ... , the 'largest' according to Charkovsky's ordering?
  • (5)
    Prove that if f(x) (here f is an arbitrary continuous function) has an orbit of prime period 2 m, then f(x) has an orbit of period 2kfor all k = 0,...,m. Prove this without using Charkovsky's Theorem (just sketch the graph of f k(x)) (i.e., use graphical analysis).
  • (6)
  • (7*)
    Consider the two dimensional discrete dynamical system defined by the function F(x,y) = r2*A(x,y), where r2 = (x2 + y2) and A is the matrix that is rotation by pi/4 counterclockwise.
  • (8*)
    Sketch the orbits of 0, 3, i, -i, (2+i), (2-i) under iteration by the function g(z) = (1/2)z + 1. What do you conclude about the dynamics of this function?
  • (9)
    Sketch the orbits of the points 1/2, 1/2*exp(i*pi), 2*exp(i*pi), and 2*exp(i*pi/4) under the quadratic function q0(z) = z2. (Here, exp(x) is the exponential function ex.)
  • (10*)
    Show that for each n = 1,2,..., the points z(n,m) = exp(itm) , m = 0,1,... 2n-1, where tm = m*(2*pi)/(2n-1), are periodic points of q0(z) of period n. From this deduce that periodic points of q0(z) "fill out" J0. Note that all these periodic points are repelling (i.e., unstable) periodic points of q0(z).
  • (11)
    (z is a complex number)
  • (12*)
    Verify that there are two fixed points of qc(z) = z2 + c for c = -(1/2) + i(1/2) and that one fixed point is attracting (stable) and the other is repelling (unstable).
  • (13*)
    Let c=0 and let W = W1 U W2 be the IFS for q0(z) (as described in Section 13.7).



    The following questions are for practice:


  • Explain why W^s(z_2) is a subset of P_c, for c = -1/2 + 1/2*i. Here, z_2 is the stable fixed point of q_c(z) for c = -1/2 + 1/2*i, and W^s(z_2) is the stable set of z_2.


  • Show that the prisoner set of the Cantor tent function, T(x) = 3x if x less than 1/2, 3-3x if x greater than 1/2, is the Cantor set. (See pages 826-829.)


  • Explain why Q_c^(-n) >> P_c as n >> infinity. In other words, show that if z is in P_c, then z is in Q_c^(-n) for all n, and conversely, if z is in P_c then z is in Q_c^(-n) for all n. Here Q_c^(-n) is as defined in the text on page 796.


  • Show that P_c and E_c are invariant sets of q_c(z).


  • Consider the quadratic function q(z) = z^2. Recall that the circle of radius 1 is invariant under q(z). Show that periodic points of q(z) are dense on the circle, i.e., let t_1 < t_2 be any two angles, then there is a t, t_1 < t < t_2, such that the complex number z = e^(it) is a periodic point of q(z). (Hint: z_n = q^n(z) = e^(i2^nt). Find a t, t_1 < t < t_2, such that 2^nt = t + 2k*pi for some integers n and k.)


  • Let (c) be the complex conjugate of c. Is J_(c) = J_c?


  • Why is the Mandelbrot set M symmetrical with respect to the (real) x-axis? That is, why is (M) = M? (Here, like in the question above, (z) denotes the complex conjugate of z.) Hint: Use an appropriate definition of M.


  • Consider Figure 14.27 in the text. Explain the relationships between the two objects in the figure (the final state diagram of the real quadratic function x^2 +c and the Mandelbrot set).


  • Explain the shapes of the main bodies of the "Mandelbrot" sets of the functions z^3 +c and z^4 + c (cf. page 8 of the handout). Also explain the symmetries of the sets.


  • Suppose q_c has a stable period 3 orbit. Make a sketch of what the Julia set J_c looks like in general. Explain the shape of J_c (or P_c) in terms of the stable set ( = basin of attraction) of the perioodic orbit. (See April 9 notes and Figure 14.18).


  • What to the Julia sets look like for those values of c where the periodic orbit is parabolic (indifferent, neutral)? Where are these c-values located on the Mandelbrot set?


  • Derive the formula for the threshold radius r(c) of q_c.


  • Derive the formulae describing the main cardioid region of the Mandelbrot set and the period 2 bulbs (including an explanation of why we are considering stable periodic orbits of q_c).


  • In what sense does the Mandelbrot set contain copies of Julia sets?


  • Explain how one could use the chaos game to draw Julia sets. Is this really a chaos game?


  • Explain why the Julia set is a fixed point of an "IFS". Begin by defining the "IFS".


  • Sketch the encirclements Q_o^(-n) for the Julia set q_o(z) = z^2.


  • Suppose P_c has zero area. Is c necessarily in the Mandelbrot set M or is c necessarily not in M, or could c be either in M or not in M?


  • Calculate the Jacobian matrix of the transformation W(w) = sqrt(w-c) (i.e., one of the inverse transformations of q_c(z) = z^2 + c = w) and determine the region in the complex plane where W contracts area.
    Correction: It is easier to calculate the determinant of the Jacobian matrix of q_c(z) and determine the regions where q_c expands area. Then the regions where W(w) contracts areas would be those (since W(w) is an inverse of q_c). To do this calculation write q_c(z) as a two-dimensional function F(x,y) where z = x + iy.



  • Back to the MAT 335 Home Page