If *X* is standard Brownian motion, what is the distribution of its absolute maximum |*X*|_{t}^{∗} = sup_{s ≤ t}|*X*_{s}| over a time interval [0, *t*]? Previously, I looked at how the reflection principle can be used to determine that the maximum *X*_{t}^{∗} = sup_{s ≤ t}*X*_{s} has the same distribution as |*X*_{t}|. This is not the same thing as the maximum of the absolute value though, which is a more difficult quantity to describe. As a first step, |*X*|_{t}^{∗} is clearly at least as large as *X*_{t}^{∗} from which it follows that it *stochastically dominates* |*X*_{t}|.

I would like to go further and precisely describe the distribution of |*X*|_{t}^{∗}. What is the probability that it exceeds a fixed positive level *a*? For this to occur, the suprema of both *X* and –*X* must exceed *a*. Denoting the minimum and maximum by

then |*X*|_{t}^{∗} is the maximum of *X*_{t}^{M} and –*X*_{t}^{m}. I have switched notation a little here, and am using *X*^{M} to denote what was previously written as *X*^{∗}. This is just to use similar notation for both the minimum and maximum. Using inclusion-exclusion, the probability that the absolute maximum is greater than a level *a* is,

As *X*_{t}^{M} has the same distribution as |*X*_{t}| and, by symmetry, so does –*X*^{m}, we obtain

This hasn’t really answered the question. All we have done is to re-express the probability in terms of *both* the minimum and maximum being beyond a level. For large values of *a* it does, however, give a good approximation. The probability of the Brownian motion reaching a large positive value *a* and then dropping to the large negative value –*a* will be vanishingly small, so the final term in the identity above can be neglected. This gives an asymptotic approximation as *a* tends to infinity,

(1) |

The last expression here is just using the fact that *X*_{t} is centered Gaussian with variance *t* and applying a standard approximation for the cumulative normal distribution function.

For small values of *a*, approximation (1) does not work well at all. We know that the left-hand-side should tend to 1, whereas 4ℙ(*X*_{t} > *a*) will tend to 2, and the final expression diverges. In fact, it can be shown that

(2) |

as *a* → 0. I gave a direct proof in this math.stackexchange answer. In this post, I will look at how we can compute joint distributions of the minimum, maximum and terminal value of Brownian motion, from which limits such as (2) will follow.

#### Distribution of the Minimum and Maximum

For *a* < 0 < *b*, we look at the event that standard Brownian motion *X* does not reach either of these levels before a fixed positive time *t*, and compute its probability. As explained in the post on the reflection principle, there are different ways of approaching such problems. While it is possible to apply the reflection principle here, giving an intuitive ‘pathwise’ approach to the solution, it can get a bit messy to go through the details. It is necessary to reflect through both levels, which can occur arbitrarily many times, and use an infinite inclusion-exclusion sequence to avoid over-counting any possibilities. Instead, I will take a more mathematically straightforward approach using the ‘square wave’ function *g*: ℝ → ℝ,

This is equal to 1 on the interval (*a*, *b*), and is antisymmetric about both *a* and *b*,

If *τ* is the first time at which *X* hits *a* or *b* then, by the strong Markov property *X*_{τ + s} – *X*_{τ} is standard Brownian motion independently of ℱ_{τ}. Consequently,conditioned on *τ* ≤ *t*, *X*_{t} – *X*_{τ} has a symmetric distribution and, by antisymmetry of *g* about the level *X*_{τ}, *g*(*X*_{t}) has zero expectation. On the other hand, if *τ* > *t* then *X*_{t} is in the interval (*a*, *b*), so *g*(*X*_{t}) = 1 giving,

As *τ* > *t* is equivalent to *a* < *X*_{t}^{m} and *X*_{t}^{M} < *b*, plugging in the definition of *g* above in the expectation on the left hand side expresses the joint distribution of *X*_{t}^{m}, *X*_{t}^{M} in terms of the normally distributed *X*_{t}.

Theorem 1IfXis standard Brownian motion,tis a fixed positive time anda< 0 <bare fixed levels then,

This idea can be taken further to compute the joint distribution of *X*_{t}^{m}, *X*_{t}^{M}, *X*_{t}. For any integrable function *f*: (*a*, *b*) → ℝ, there is a unique *g*: ℝ → ℝ equal to *f* on the interval (*a*, *b*) and antisymmetric about both *a* and *b*. Note that reflecting about both these levels in turn is the same as translating by 2(*b* - *a*), so *g* must be periodic with period 2(*b* - *a*). It can be explicitly constructed as

(3) |

where I am setting *u*_{n} = *n*(*b* - *a*) for *n* even and *u*_{n} = *n*(*b* - *a*) + *a* + *b* for *n* odd. Writing it in this way will be useful for our application, although it may not be immediately obvious that it satisfies the required conditions. To see that it is antisymmetric about *a*, it can be checked that reflecting *x* through this level exchanges terms *n* = 2*m* and *n* = 2*m* – 1 in the sum whilst flipping their sign, and reflecting through *b* exchanges *n* = 2*m* and *n* = 2*m* + 1 with opposite signs.

If *τ* is the first time when *X* hits either *a* or *b* then, using exactly the same argument as above, antisymmetry about both *a* and *b* implies that the expectation of *g*(*X*_{t}) conditioned on *τ* ≤ *t* is zero, and *g*(*X*_{t}) = *f*(*X*_{t}) on *τ* > *t*,

Plugging in the expression above for *g* gives us the joint distribution of *X*_{t}^{m}, *X*_{t}^{M}, *X*_{t}.

Theorem 2IfXis standard Brownian motion,tis a fixed positive time,a< 0 <bare fixed levels andf: (a,b) → ℝis integrable then,

(4)

whereu_{n}=n(b-a)forneven andu_{n}=n(b-a) +a+bfornodd.

Here, the expectations on the right hand side are just the summand in (3) with *X*_{t} replacing *x*. By symmetry of the normal distribution, I also replaced (-1)^{n}*X*_{t} by *X*_{t} to simplify the expression.

The expectations on the right hand side of (4) can be rearranged. For any integrable function *h*, using the fact that *X*_{t} is normal with zero mean and variance *t*,

This is just the standard measure change for translating a normal density. Applying this to the expectations in the summation on the right hand side of (4) gives,

This gives the joint distribution of *X*_{t}^{m}, *X*_{t}^{M} conditioned on *X*_{t}.

Theorem 3IfXis standard Brownian motion,tis a fixed positive time,a< 0 <bare fixed levels then,

fora<X_{t}<b, whereu_{n}=n(b-a)forneven andu_{n}=n(b-a) +a+bfornodd.

In particular, taking *b* = –*a* gives the distribution of the absolute maximum conditioned on the terminal value.

Corollary 4IfXis standard Brownian motion,tis a fixed positive time anda> 0is a fixed level then,

for|X_{t}|<a, whereu_{n}= 2na.

Conditioning on *X*_{t} = 0 also gives the distribution of the absolute maximum of a Brownian bridge.

Corollary 5IfXis a standard Brownian bridge anda> 0is a fixed level then,

#### Transformed Sums

The sums given in theorems 1,2,3 and corollaries 4,5 all converge quickly in *n*. This is especially true for large absolute values of *a* and *b* where very few terms are needed to get an accurate approximation. On the other hand, for small values of *a* and *b*, a lot of terms may be needed before the convergence kicks in, so is not so useful in this case. I now look at how these can be transformed to give a sum where only a small number of terms will be needed for small absolute values of *a* and *b*. This is essentially a Jacobi identity, proved using Fourier series.

As the function *g* defined by (3) is periodic, it can be written as a Fourier series which consists of terms of the form exp(*i*π*nx*/(*b* - *a*)). Since expectations of exp(*i*π*nX*_{t}/(*b* - *a*)) have a simple closed form expression, this removes the summation from (4) replacing it, instead, with the Fourier expansion. Actually, antisymmetry says that *g* can be written as a sine series,

with coefficients,

(5) |

The expectation of the sine of a normal variable of mean *μ* and variance *ν* is

So, we can take expectations of the sine expansion of *g*(*X*_{t}) to obtain a description of the joint distribution of the maximum, minimum and terminal value of a Brownian motion in the same way as for theorem 2 above.

Theorem 6IfXis standard Brownian motion,tis a fixed positive time,a< 0 <bare fixed levels andf: (a,b) → ℝis integrable then,

(6)

wherec_{n}are given by (5).

This is a fast converging sequence and, when *b* – *a* is small very few terms will be required for an accurate approximation.

Computing the coefficients for *f* = 1,

gives an alternative expression for the distribution of the minimum and maximum from that in theorem 1.

Theorem 7IfXis standard Brownian motion,tis a fixed positive time anda< 0 <bare fixed levels then,

Applying this with *a* = –*b* gives an expression for the distribution of the absolute maximum.

Corollary 8IfXis standard Brownian motion,tis a fixed positive time anda> 0is a fixed level then,

Looking again at expression (6), substituting in the integral (5) for the sine components and commuting summation with the integration,

where

To compute the probability conditioned on *X*_{t}, we just divide by its probability density, giving a transformed version of theorem 3.

Theorem 9IfXis standard Brownian motion,tis a fixed positive time,a< 0 <bare fixed levels then,

(7)

fora<X_{t}<b.

Taking *b* = –*a* gives the distribution of the absolute maximum conditioned on the terminal value. Doing this, the product of sines on the right hand side of (7) is

We obtain the transformed version of corollary 4 for the distribution of the absolute maximum of *X* conditioned on its terminal value.

Corollary 10IfXis standard Brownian motion,tis a fixed positive time anda> 0is a fixed level then,

for|X_{t}|<a.

Finally, conditioning on *X*_{t} = 0 gives the transformed version of corollary 5 for the absolute maximum of a Brownian bridge.

Corollary 11IfXis a standard Brownian bridge anda> 0is a fixed level then,