AMath 584, Autumn 2015
Sample Solutions for Assignment 4.
Reading: Lectures 9-11 in the text.
1. p. 68, Exercise 9.2.
(a) Since A is upper triangular, its eigenvalues are its diagonal entries, which
a
16
Period doubling route to chaos
We now study the "routes" or "scenarios" towards chaos.
We ask: How does the transition from periodic to strange attractor occur?
The question is analogous to the s
17
Intermittency (and quasiperiodicity)
In this lecture we discuss the other two generic routes to chaos, intermittency and quasiperiodicity. Almost all our remarks will be on intermittency; we close
AMath 584, Autumn 2015
Sample Solutions for Assignment 3.
Reading: Lectures 6-8 in the text.
1. p. 47, Exercise 6.1.
I 2P is unitary since (I 2P ) (I 2P ) = I 2P 2P + 4P P = I,
where the last equality
AMath 584, Autumn 2015
Sample Solutions for Assignment 2.
Reading: Lectures 3-5 in the text.
1. p. 24, Exercise 3.5.
The same is true for the Frobenius norm since
m
E
2
F
n
2
|Eij | =
=
i,j
m
|ui |2 |
HW #4 - Solutions of Selected Problems
5.1.9 (Page 141)
b) Follow the hint:
xx y y = x(y) y(x) = 0.
Notice that (x2 ) = 2xx, (y 2 ) = 2y y, integrating the above equation we obtain that
x2 y 2 = C,
wh
HW #2 - Solutions of The Rest of The Problems
Note: Whenever I say only one eigenvector, I mean (of course) up to scaling.
5.2.2 (Page 141)
a) The system can be written as follows:
x
y
=
1 1
1 1
x
y
T
HW 2 - Solutions of the Rest of the Problems
4.1.3 (Page 113)
The xed points are at sin2 = 0, i.e. at = 0, /2, , 3/2. The vector eld looks as shown
below:
4.1.4 (Page 113)
The xed points are at sin3 =
Solutions 1
2.2.5 There are no equilibria:
between 1/2 and 3/2.
1
2
f (c) 3 , so x(t) always moves to the right with speed
2
2.2.7 For x > 0, ex > 1 cos x, so f (x) = ex cos x > 0 has no equilibria i
HW #3 - Solutions of Selected Problems
3.6.2 (Page 86)
When h = 0 this system is the same as the one in Section 3.2. As h varies, the curves in
the following pictures move vertically.
r<0
r=0
r>0
V
V
HW #2 - Solutions of Selected Problems
3.1.2 (Page 79)
The graph of the function y = coshx is shown below on the right, with dotted lines
indicating the values of r. Its clear that x moves to the righ
15
Lyapunov exponents
Whereas fractals quantify the geometry of strange attractors, Lyaponov ex ponents quantify the sensitivity to initial conditions that is, in effect, their most salient feature. I
14
Fractals
We now proceed to quantify the "strangeness" of strange attractors. There are two processes of interest, each associated with a measurable quantity: sensitivity to initial conditions, quan
13
Experimental attractors
In this brief lecture we show examples of strange attractors found in experi ments.
13.1
Rayleigh-Bnard convection e
BPV, Figure VI.25 Two dynamical variables, (T ) and (T )
1
1.1
Pendulum
Free oscillator
To introduce dynamical systems, we begin with one of the simplest: a free oscillator. Specifically, we consider an unforced, undamped pendulum.
l
mg sin mg
The arc lengt
2
Stability of solutions to ODEs
How can we address the question of stability in general? We proceed from the example of the pendulum equation. We reduce this second order ODE, g + sin = 0, l to two f
3
Conservation of volume in phase space
We show (via the example of the pendulum) that frictionless systems conserve volumes (or areas) in phase space. Conversely, we shall see, dissipative systems co
Calculate
Pictorially
x1 x2 f = + = 0+0 x1 x2
x2
x1
Note that the area is conserved. Conservation of areas holds for all conserved systems. This is conventionally derived from Hamiltonian mechanics a
5
5.1
Forced oscillators and limit cycles
General remarks
How may we describe a forced oscillator? The linear equation + + 2 = 0 is in general inadequate. Why?
(3)
Linearity if (t) is a solution, then
6
6.1
Parametric oscillator
Mathieu equation
We now study a different kind of forced pendulum. Specifically, imagine subjecting the pivot of a simple frictionless pendulum to an alternating vertical m
7
Fourier transforms
Except in special, idealized cases (such as the linear pendulum), the precise oscillatory nature of an observed time series x(t) may not be identified from x(t) alone. We may ask
8
Poincar sections e
The dynamical systems we study are of the form d (t) = F ( , t) x x dt Systems of such equations describe a flow in phase space. The solution is often studied by considering the t
9
Fluid dynamics and Rayleigh-Bnard convection e
In these lectures we derive (mostly) the equations of viscous fluid dynamics. We then show how they may be generalized to the problem of RayleighBnard
10
Introduction to Strange Attractors
Thus far, we have studied only classical attractors such as fixed points and limit cycles. In this lecture we begin our study of strange attractors. We emphasize
11
Lorenz equations
In this lecture we derive the Lorenz equations, and study their behavior. The equations were first derived by writing a severe, low-order truncation of the equations of R-B convect
12
Hnon attractor e
The chaotic phenomena of the Lorenz equations may be exhibited by even simpler systems. We now consider a discrete-time, 2-D mapping of the plane into itself. The e points in R2 ar