Fourier series · complex plane

Fourier epicycles

A chain of circles, each turning at its own whole-number speed, draws a closed curve; each circle is one term of the curve’s complex Fourier series.

Fourier epicycles in this visualization

Walk once round a closed curve in the plane and read each point as a complex number z(t)z(t), with tt from 0 to 1 for one turn. Its complex Fourier series writes z(t)z(t) as a sum of terms cke2πiktc_k e^{2\pi i k t}, one for each whole number kk. Each term is a point going round a circle of radius ∣ck∣|c_k|, kk times per turn, starting at the angle arg⁡ck\arg c_k; negative kk turn the other way. Put the circles tip to tail and the last tip draws the sum. c0c_0 does not turn: it is where the chain starts.

z(t)=∑kck e2πikt,ck e2πikt: radius ∣ck∣, k turns per turn, starting at arg⁡ckz(t) = \sum_{k} c_k\, e^{2\pi i k t},\qquad c_k\, e^{2\pi i k t}:\ \text{radius } |c_k|,\ k \text{ turns per turn, starting at } \arg c_k

What the page computes

The page samples the curve at N=1024N = 1024 points, evenly along its length (or, for the heart, evenly in its formula’s tt), and computes the discrete coefficients c~k\tilde c_k. All NN terms together pass exactly through the NN samples. The error at the samples is ∑unused∣c~k∣2/∑∣c~k∣2\sqrt{\sum_\text{unused} |\tilde c_k|^2 / \sum |\tilde c_k|^2} (without c0c_0): by the discrete Parseval identity it is the root-mean-square distance between the drawn curve and the samples, relative to the shape’s size. It is not a comparison with your original stroke.

Build up adds circles level by level, M=1,2,4,8,…M = 1, 2, 4, 8, \ldots, one turn each, largest first (or lowest frequency first). Every finished level stays as a faint ghost: where the levels agree the light piles up, where they still disagree it spreads. One circle is a circle; k=1k = 1 and k=−1k = -1 together make an ellipse with semi-axes ∣c1∣+∣c−1∣|c_1| + |c_{-1}| and ∣∣c1∣−∣c−1∣∣\bigl||c_1| - |c_{-1}|\bigr|.

The ticks on the pen’s trace are 1/641/64 of a turn apart in time. Where they crowd, the pen went slowly over that stretch on average; between two ticks a fast circle can turn several times, so they do not show the speed at an instant.

Easy to misread

History

Ancient astronomers described the planets with circles riding on circles, the deferent and the epicycle. A chain of uniformly turning circles is exactly a sum of terms cke2πiktc_k e^{2\pi i k t}, which is why enough epicycles can follow any closed path [2]. Fourier introduced his series in 1807 for the flow of heat [1].

Related concepts

Complex plane · Unit circle · Fourier series (page planned) · Parametric equation (page planned)

References

  1. J. Fourier, Théorie analytique de la chaleur, Firmin Didot, Paris (1822); the 1807 memoir to the Institut de France.
  2. N. R. Hanson, “The mathematical power of epicyclical astronomy”, Isis 51 (1960) 150–158.

Further reading: Wikipedia: Fourier series; MacTutor History of Mathematics: Joseph Fourier.