In probability theory, Maxwell's theorem (known also as Herschel-Maxwell's theorem and Herschel-Maxwell's derivation) states that if the probability distribution of a random vector in
is unchanged by rotations, and if the components are independent, then the components are identically distributed and normally distributed.
If the probability distribution of a vector-valued random variable X = ( X1, ..., Xn )T is the same as the distribution of GX for every n×n orthogonal matrix G and the components are independent, then the components X1, ..., Xn are normally distributed with expected value 0 and all have the same variance.
This theorem is one of many characterizations of the normal distribution.
The only rotationally invariant probability distributions on Rn that have independent components are multivariate normal distributions with expected value 0 and variance σ2In, (where In = the n×n identity matrix), for some positive number σ2.
James Clerk Maxwell proved the theorem in Proposition IV of his 1860 paper.
[1] Ten years earlier, John Herschel also proved the theorem.
[2] The logical and historical details of the theorem may be found in.
[3] We only need to prove the theorem for the 2-dimensional case, since we can then generalize it to n-dimensions by applying the theorem sequentially to each pair of coordinates.
Since rotating by 90 degrees preserves the joint distribution, both
has the same probability measure.
is a Dirac delta distribution at zero, then it's a gaussian distribution, just degenerate.
By Lebesgue's decomposition theorem, we decompose it to a sum of regular measure and an atomic measure:
contains an atomic part, then there exists some
, and so the double ray
has nonzero probability.
Now by rotational symmetry of
, any rotation of the double ray also has the same nonzero probability, and since any two rotations are disjoint, their union has infinite probability, contradiction.
(As far as we can find, there is no literature about the case where
is singularly continuous, so we will let that case go.)
have probability density function
, and the problem reduces to solving the functional equation
ρ ( x ) ρ ( y ) = ρ ( x cos θ + y sin θ ) ρ ( x sin θ − y cos θ )