Theme

Sensitive to everything — page 2

Systems where a difference too small to draw becomes the whole difference, and the reason that is a property of the rule rather than of the measurement.
A quadrilateral the Hénon map carries strictly inside itself. Hénon's trapping quadrilateral with corners (-1.33, 0.42), (1.32, 0.133), (1.245, -0.14), (-1.06, -0.5), the image of its boundary under the Hénon map lying strictly inside it, and an orbit on the attractor. Dynamics

A region the orbit cannot leave

That the Hénon map and the Lorenz flow have attractors takes two lines of arithmetic each — a region mapped strictly inside itself, a quantity that falls outside an ellipsoid. That the attractors are strange took a computer and twenty years for Lorenz, and for Hénon's classical parameters it has never been done. The difference between the two questions is the difference between a region and what lives in it.

A strand of the logistic map's diagram between r = 3.8284 and 3.8572, upside down. The bifurcation diagram of the logistic map cropped to r in [3.8284, 3.8572] and x in [0.44, 0.565], with the value axis inverted. Dynamics

Every window is the whole diagram again

Zoom into one strand of the period-three window of the logistic map and the whole bifurcation diagram is there again — a single value, a fork, a cascade of doublings, chaos, bands and windows — upside down and ninety-seven times narrower. The doublings inside it shrink by Feigenbaum's 4.669, exactly as the originals do, because near the top of its hump the map's third iterate is itself a one-humped map, turned over.

What a rule collects when every value comes from the same distribution. uniform on [0, 1]: best rule 0.995, one threshold 0.635 of E[max] at n = 200; exponential: best rule 0.905, one threshold 0.678 of E[max] at n = 200; Pareto, tail exponent 2: best rule 0.802, one threshold 0.714 of E[max] at n = 200; Pareto, tail exponent 1.2: best rule 0.802, one threshold 0.682 of E[max] at n = 200. Probability

When every value comes from the same hat

A rule that sees values one at a time and must keep or discard each on the spot can guarantee half of what a prophet collects, and no more, when the values come from different distributions. When they all come from the same one, the guarantee rises to 0.745 — and a single fixed threshold, set so that each value crosses it with chance 1/n, already secures 1 − 1/e. For bounded values the best rule collects nearly everything; only a heavy tail, where one enormous value carries the prize, keeps the gap open.

The expected rank a threshold rule achieves with the values shown, against what is known. n=1: 1.0000 (c=4.000); n=2: 1.2500 (c=1.000); n=3: 1.4009 (c=1.124); n=5: 1.5868 (c=1.257); n=10: 1.8141 (c=1.416); n=20: 1.9950 (c=1.555); n=50: 2.1557 (c=1.701); n=100: 2.2284 (c=1.781); n=200: 2.2725 (c=1.838); n=400: 2.2983 (c=1.877); n=800: 2.3130 (c=1.902). Probability

The rank that remembers every value

Values arrive one at a time, each must be kept or discarded on the spot, and the aim is to keep one whose rank among all of them is low on average. Told only who is leading, the best rule gets 3.87. Shown the values, a rule gets below 2.33 — and how much lower the best possible rule goes is not known, because the rank of what is kept depends on every value seen, and the best rule may need to remember all of them.

Traffic with random dawdling: jams from nowhere at density 0.18. Nagel–Schreckenberg traffic, 29 cars on 160 cells, top speed 5, dawdling probability 0.25; up to 11 cars stopped at once. Dynamics

A jam that comes from nowhere

Give cars on a ring road a top speed of five cells a step, let each slow to the gap ahead, and add one more rule: now and then, at random, a driver eases off by one. That is the whole of the Nagel–Schreckenberg model, and it produces what the exactly solvable rule 184 could not — jams that form in free traffic with no obstacle, drift backwards against the flow, and cost the road more than a third of its capacity. Set the top speed to one and remove the chance, and it is rule 184 again, cell for cell.

The medial axis of a wobbly rectangle, whole and trimmed. A rectangle with sides wobbling by ±0.02: its medial axis (654 Voronoi vertices, Hausdorff distance 0.499 from the clean axis) and its λ-medial axis at λ = 0.2 (63 vertices, distance 0.078). Geometry

A skeleton that ignores small teeth

The medial axis of a shape grows a branch to every bump on its outline, however small, so a scanned shape has a skeleton made mostly of noise. Keep only the points of the axis whose nearest boundary points are spread across a circle of radius at least λ, and the branches to bumps smaller than about 2λ disappear while the real skeleton stays. That λ-medial axis moves only a little when the outline wobbles by much less than λ — and the price is λ itself, a size below which a feature is declared not to exist, which can even cut a shape's skeleton in two.

A random sequence that shrinks, grows or wanders, by one parameter. Runs of tₖ₊₁ = tₖ ± β·tₖ₋₁ at β = 0.6, 0.70258, 0.8 on a logarithmic scale: shrinking, wandering and growing. Dynamics

A threshold no average can see

Damp the random Fibonacci sequence — add or subtract β times the term before last, by a coin — and somewhere between β = 0.6 and β = 0.8 almost every run stops shrinking and starts growing. The switch is at β ≈ 0.70258. The average of the squares, which can be computed exactly, grows for every β there is, so the one number that decides the sequence's fate is invisible to the one number anybody can compute.

Two chaotic maps, coupled weakly and strongly. The gap between two coupled copies of the logistic map at full chaos: at ε = 0.15 it stays large; at ε = 0.3 it falls to nothing. Dynamics

Two chaotic orbits made to agree

Chaos means two copies of a system started a hair apart soon have nothing in common. Couple them — let each take a fraction ε of its next step from the other — and above a definite strength they agree to every digit for ever, while each stays exactly as chaotic as before. The strength needed is read off the same number that measured the chaos: the Lyapunov exponent.

Newton's method for z² − 1: two half-planes and the line between. The complex plane split down the imaginary axis into two shaded halves, the roots ±1 marked, and five orbits of Newton's method each converging to the root on its own side. Dynamics

Chaos on the line between two roots

Newton's method for a quadratic is the one case solved completely: every starting point goes to the nearer root, and the plane splits into two half-planes. Cayley did it in 1879. What his solution also contains is a line on which the method never settles — and a change of variable shows that on that line Newton's method is exactly the doubling map, the simplest chaos there is. Run on the real line for x² + 1, which has no real root to find, it becomes a source of Cauchy-distributed randomness.

All themes