Question 1 · The doubling map and chaotic dynamics
Define the doubling map $T: [0,1) \to [0,1)$ by
$T(x) = 2x \bmod 1 = \begin{cases} 2x & \text{if } 0 \le x < \tfrac{1}{2} \\ 2x - 1 & \text{if } \tfrac{1}{2} \le x < 1 \end{cases}$
The orbit of a starting point $x_0$ is the sequence $x_0, x_1, x_2, \ldots$ where $x_{n+1} = T(x_n)$. Throughout this question you may express numbers in binary using the notation $0.b_1 b_2 b_3 \ldots_2$ where each $b_i \in \{0,1\}$.
(a) [3 marks]
Compute the first five terms of the orbit starting at $x_0 = \tfrac{1}{3}$.
Reveal worked solution
Step 1. Apply $T$ repeatedly.
Observation: the orbit is a period-2 cycle $\tfrac{1}{3} \to \tfrac{2}{3} \to \tfrac{1}{3} \to \ldots$ (A1 · A1 · A1)
(b) [4 marks]
Let $x_0 = 0.b_1 b_2 b_3 \ldots_2$ be a binary expansion. Prove that $x_1 = T(x_0) = 0.b_2 b_3 b_4 \ldots_2$.
Reveal worked solution
Setup. By definition, $x_0 = \sum_{i=1}^{\infty} b_i \, 2^{-i}$.
Case split.
- If $b_1 = 0$: $x_0 < \tfrac12$, so $x_1 = 2x_0 = 0.b_2 b_3 \ldots_2$. (M1)
- If $b_1 = 1$: $x_0 \ge \tfrac12$, so $x_1 = 2x_0 - 1 = b_1 + 0.b_2 b_3 \ldots_2 - 1 = 0.b_2 b_3 \ldots_2$. (M1)
In both cases $x_1 = 0.b_2 b_3 b_4 \ldots_2$ — the doubling map is a binary shift. (R1 · AG1)
(c) [3 marks]
Deduce that every rational $x_0 \in (0,1)$ has an eventually periodic orbit.
Reveal worked solution
Argument. A rational has a binary expansion that is either terminating or eventually periodic. (R1)
By part (b), $T$ shifts this expansion left by one bit each step, so after finitely many applications the tail is purely periodic; the orbit then repeats every $p$ steps where $p$ is the binary period. (R1 · A1)
(d) [5 marks]
Two starting points $x_0 = 0.0101\overline{01}_2$ and $y_0 = 0.0101\overline{10}_2$ differ from the 5th binary digit onwards. Prove that after $n$ applications of $T$, $|x_n - y_n| \le 2^{-(4-n)}$ for $n \le 4$, and that the orbits diverge to be at least $\tfrac{1}{4}$ apart by step 4.
Reveal worked solution
Setup. Because the doubling map is a binary shift (part b), the difference in the first surviving bit is what matters after $n$ shifts.
Initial difference: agreeing on bits 1–4, differing at bit 5. So $|x_0 - y_0| \le 2^{-4}$. (M1)
After $n$ applications of $T$, the leading disagreeing bit is at position $5 - n$. In general, $|x_n - y_n|$ can double each step (approximately). (M1)
- $n=1$: bits 4 disagree ⇒ $|x_1 - y_1| \le 2^{-3} = \tfrac{1}{8}$
- $n=2$: bits 3 disagree ⇒ $|x_2 - y_2| \le 2^{-2} = \tfrac{1}{4}$
- $n=3$: bits 2 disagree ⇒ $|x_3 - y_3| \le \tfrac{1}{2}$
- $n=4$: bits 1 disagree; the numbers are macroscopically apart. (A1 · A1)
Sensitivity to initial conditions: tiny changes in the starting point ($\le 2^{-4}$) grow to differences of $\ge \tfrac{1}{4}$ within 4 steps — a signature of chaos. (R1)
(e) [5 marks]
Prove by mathematical induction that for any $x_0 \in [0,1)$ whose binary expansion starts with $k$ zeros ($b_1 = b_2 = \ldots = b_k = 0$), $x_k = 2^k x_0$.
Reveal worked solution
Base case ($k = 0$): $x_0 = 2^0 x_0 = x_0$. Trivially true. (M1)
Inductive step. Assume $x_j = 2^j x_0$ for some $0 \le j < k$. Because bits $b_1, \ldots, b_k$ are all $0$, $x_j$ has its bit-$1$ position equal to $b_{j+1} = 0$, so $x_j < \tfrac12$. (M1)
Conclusion. By induction, $x_k = 2^k x_0$ for all $0 \le k \le$ (first non-zero bit position). (R1 · AG1 · A1)
(f) [4 marks]
The Lyapunov exponent of the doubling map at $x_0$ is defined as $\lambda(x_0) = \lim_{n \to \infty} \tfrac{1}{n} \sum_{i=0}^{n-1} \ln |T'(x_i)|$ wherever $T'$ exists. Show that $\lambda(x_0) = \ln 2$ for every $x_0$ where the derivative is defined.
Reveal worked solution
Compute $T'$. On each of $[0, \tfrac12)$ and $[\tfrac12, 1)$, $T(x)$ is linear with slope $2$. So $T'(x) = 2$ wherever it exists. (M1 · A1)
Apply the definition.
Interpretation. A positive Lyapunov exponent ($\ln 2 \approx 0.693$) confirms the chaos observed in (d): nearby orbits diverge exponentially fast, doubling their separation each step. (R1 · AG1)
(g) [4 marks]
Explain, using parts (b) and (f), why despite being deterministic the doubling map is unpredictable in practice.
Reveal worked solution
Argument. By part (b), each application of $T$ discards one bit of the initial condition. After $n$ steps, the current state depends only on bits $b_{n+1}, b_{n+2}, \ldots$ of $x_0$. (R1)
In practice we know $x_0$ only to finite precision — say $B$ bits. After $B$ iterations, all known information about $x_0$ has been "shifted out"; the future orbit is determined by bits we never measured. (R1)
Equivalently, from part (f), errors grow by a factor of $2$ per step ($e^{\lambda} = 2$). Any starting-value uncertainty of $2^{-B}$ becomes $O(1)$ after $B$ steps. (R1 · A1)