Math 4/5027 - Nonlinear Dynamics and Chaos
Solution Set #2 - Spring 1999

Problems = 5 points each (T2.7 is 10 points). Undergraduate total = 25. Graduate total = 40.

Problem T2.2. If we view the map in components as (f(x,y),g(x,y)), then we must solve

eqnarray7

simultaneously for the fixed points. The first fixed point (0,0) is evident; but don't forget the second fixed point (-0.6,-0.6).

*Problem T2.5. To find the period two points, we must solve the equations tex2html_wrap_inline144 or

displaymath100

Solving the second equation for y and substituting into the first equation, we have (eventually)

displaymath101

Notice that we have factored out the known period-1 points, leaving a quadratic for the period-2 points. Solving the quadratic, the period-2 points are

displaymath102

The condition for existence of period-2 points is tex2html_wrap_inline148 or tex2html_wrap_inline150 .

*Problem T2.7. The real goal of this problem is to determine, not just verify, the given intervals in a.

(a) The fixed points of the Henon map are the solutions of the system

eqnarray24

Setting y=x in the first equation, the first component satisfies tex2html_wrap_inline156 which has roots

displaymath103

The problem considers the case b=0.4, which gives the fixed points tex2html_wrap_inline160 where

displaymath104

Notice that the fixed points exist only for tex2html_wrap_inline162 . Also note that tex2html_wrap_inline164 for all tex2html_wrap_inline162 , while tex2html_wrap_inline168 for tex2html_wrap_inline170 and tex2html_wrap_inline172 for a >0. Furthermore, tex2html_wrap_inline176 as tex2html_wrap_inline178 . The diagram, below and left, shows how the fixed points vary with the parameter a.

Now we need to investigate the stability of these fixed points. The Jacobian is given by

displaymath105

A quick calculation shows that the eigenvalues satisfy tex2html_wrap_inline182 which has roots

displaymath106

Some sleuthing is now required. Let's consider tex2html_wrap_inline184 first. Note that tex2html_wrap_inline186 for all values of tex2html_wrap_inline162 . Therefore

displaymath107

For the behavior of tex2html_wrap_inline190 , notice that when a=-0.09, we have tex2html_wrap_inline194 and tex2html_wrap_inline196 . As a increases, tex2html_wrap_inline184 decreases, and tex2html_wrap_inline190 approaches 0. Therefore, for tex2html_wrap_inline162 , tex2html_wrap_inline206 and tex2html_wrap_inline208 , making the fixed point tex2html_wrap_inline210 a saddle for all tex2html_wrap_inline162 .

Now consider tex2html_wrap_inline214 . It helps to determine for what values of a and tex2html_wrap_inline214 we have tex2html_wrap_inline220 . A little calculation shows that when a=-0.09, tex2html_wrap_inline224 , tex2html_wrap_inline226 and tex2html_wrap_inline196 . Similarly, when a=0.27, tex2html_wrap_inline232 , tex2html_wrap_inline234 and tex2html_wrap_inline236 . As a increases, tex2html_wrap_inline214 increases, and tex2html_wrap_inline242 decreases from 1 to 0, while tex2html_wrap_inline190 decreases from -0.4 to tex2html_wrap_inline248 (see figure above and right). We see that the fixed point tex2html_wrap_inline250 changes from a sink to a source as tex2html_wrap_inline190 passes through -1, which occurs at a=0.27.

(b) As noted above, when a=0.27, tex2html_wrap_inline232 and the eigenvalues of the Jacobian are tex2html_wrap_inline260 and 0.4.

(c) We found the period-2 points of the Henon map in T2.5 above. They are the roots of the quadratic

displaymath108

With b=0.4, the period two points are tex2html_wrap_inline250 and tex2html_wrap_inline210 (recall that x=y for all of the periodic points), where

displaymath109

Clearly, these period-2 points exist only for a>0.27. To determine the stability of these points, we need to compute tex2html_wrap_inline274 , the Jacobian of tex2html_wrap_inline276 , which is given by tex2html_wrap_inline278 . Using the expression for the Jacobian given above, we see that

displaymath110

The eigenvalues satisfy the quadratic equation

displaymath111

which has roots

displaymath112

The discriminant in this expression is negative when tex2html_wrap_inline280 or when 0.36<a<0.76. Thus for values of a in this interval, the eigenvalues are complex. Furthermore, for values of a in this interval, the eigenvalues have a constant absolute value of 0.4. So for 0.36<a<0.76, the period-2 pair is a sink. But there is a larger interval in which the period-2 pair is a sink.

tex2html_wrap382

(d) As described above, when a=0.85, we have tex2html_wrap_inline236 .

Problem 2.1a. To find the eigenvalues, we must solve the equation

displaymath113

The roots are tex2html_wrap_inline318 and tex2html_wrap_inline320 , which says that (0,0) is a source (because both eigenvalues are greater than 1 in absolute value).

Problem 2.1c. To find the eigenvalues, we must solve the equation

displaymath114

The roots are tex2html_wrap_inline324 and tex2html_wrap_inline326 , which says that (0,0) is a sink (because both eigenvalues are less than 1 in absolute value).

Problem 2.3. The fixed points satisfy

displaymath115

Solving this system, we find two fixed points, (0,0) and (3,9). The Jacobian calculation gives us

displaymath116

The eigenvalues corresponding to the point (0,0) are tex2html_wrap_inline336 , meaning that (0,0) is a saddle. The eigenvalues corresponding to the point (3,9) are tex2html_wrap_inline342 , which says (3,9) is a source (disagreeing with the book).

Problem 2.8a. Following the procedure outlined in Theorem 2.24 of the book, we must compute the eigenvalues and eigenvectors of tex2html_wrap_inline346 . We find that

displaymath117

The eigenvalues are tex2html_wrap_inline348 and tex2html_wrap_inline350 , with corresponding eigenvectors (1,1) and (1,-1). Thus the image of the unit disk is an ellipse with a semi-major axis in the (1,1) direction of length tex2html_wrap_inline358 and a semi-minor axis in the (1,-1) direction of length tex2html_wrap_inline362 . The area of the image ellipse is tex2html_wrap_inline364 .

Problem 2.8b. We must compute the eigenvalues and eigenvectors of

displaymath118

The eigenvalues are tex2html_wrap_inline366 , with corresponding eigenvectors tex2html_wrap_inline368 and tex2html_wrap_inline370 . This says that the image of the unit disk is an ellipse with a semi-major axis in the (1,-2) direction of length tex2html_wrap_inline374 and a semi-minor axis in the (2,1) direction of length tex2html_wrap_inline378 . The area of the image ellipse is tex2html_wrap_inline380 .



Mon Feb 22 05:11:04 MST 1999