Lévy construction of Brownian motion ($B_{(u+v)/2}=(B_u+B_v)/2+\sqrt{v-u}Z/2$)

The Lévy construction builds Brownian motion on $[0,1]$ by successive independent midpoint displacements on dyadic intervals. Start with $B_0=0$ and $B_1=Z_0$, where $Z_0$ has the standard normal distribution. On each dyadic interval $[u,v]$, choose a fresh independent standard normal random variable $Z$ and use the displayed midpoint rule. Conditional on the endpoints, this is the correct Brownian bridge midpoint law.

To verify the increment law, if $D=B_v-B_u\sim N(0,\ell )$ with $\ell =v-u$, the new increments are $D/2\pm \sqrt{\ell }Z/2$. They are Jointly Gaussian, each has variance $\ell /2$, and their covariance is zero, hence they are independent. Induction gives all dyadic Brownian increments. The change between consecutive linear interpolants has uniform norm $2^{-n/2-1}{\max }_{k<2^n}|Z_{n,k}|$. The Gaussian tail bound and union bound bound the chance that this maximum exceeds $n$ by $2^{n+1}{e}^{-n^2/2}$. The First Borel-Cantelli lemma and ${\sum }_{n}n{2}^{-n/2}<\infty$ prove uniform convergence almost surely. The continuous limit has the required Brownian increments at arbitrary times by dyadic approximation.

Let $H_{n,k}(s)=2^{n/2}(1_{(0,1/2)}-1_{(1/2,1)})(2^{n}s-k)$ for $n\ge 0$ and $0\le k<2^n$. Integrating these Haar wavelets gives tents of height $2^{-n/2-1}$ supported on the corresponding dyadic intervals. For independent standard normal random variables $Z_0,Z_{n,k}$, the displayed series, summed level by level, is exactly the succession of midpoint displacements in the Lévy construction of Brownian motion. Its almost sure convergence is uniform, and its limit is standard Brownian motion. This gives a continuous function from Gaussian coefficients even though the unintegrated unit-variance Haar wavelet series diverges in L2 space.