Brownian bridge: Difference between revisions
→Relation to other stochastic processes: fix latex typo |
→Relation to other stochastic processes: The Brownian bridge is this context is defined for t\in[0,1] which is different from the Brownian motion for t\in[0,T] defined previously. |
||
(30 intermediate revisions by 20 users not shown) | |||
Line 1: | Line 1: | ||
{{Short description|A process in physics}} |
|||
[[Image:Brownian bridge.png|thumb|Brownian motion, pinned at both ends. This |
[[Image:Brownian bridge.png|thumb|Brownian motion, pinned at both ends. This represents a Brownian bridge.]] |
||
A '''Brownian bridge''' is a continuous-time [[ |
A '''Brownian bridge''' is a continuous-time [[gaussian process]] ''B''(''t'') whose [[probability distribution]] is the [[conditional probability distribution]] of a standard [[Wiener process]] ''W''(''t'') (a mathematical model of [[Brownian motion]]) subject to the condition (when standardized) that ''W''(''T'') = 0, so that the process is pinned to the same value at both ''t'' = 0 and ''t'' = ''T''. More precisely: |
||
:<math> B_t := (W_t\mid W_T=0),\;t \in [0,T] </math> |
:<math> B_t := (W_t\mid W_T=0),\;t \in [0,T] </math> |
||
The expected value of the bridge is zero, with variance <math>\textstyle\frac{t(T-t)}{T}</math>, implying that the most uncertainty is in the middle of the bridge, with zero uncertainty at the nodes. The [[covariance]] of ''B''(''s'') and ''B''(''t'') is ''s''(T − ''t'')/T if ''s'' |
The expected value of the bridge at any ''t'' in the interval [0,''T''] is zero, with variance <math>\textstyle\frac{t(T-t)}{T}</math>, implying that the most uncertainty is in the middle of the bridge, with zero uncertainty at the nodes. The [[covariance]] of ''B''(''s'') and ''B''(''t'') is <math>\min(s,t)-\frac{s\,t}{T}</math>, or ''s''(T − ''t'')/T if ''s'' < ''t''. |
||
The increments in a Brownian bridge are not independent. |
The increments in a Brownian bridge are not independent. |
||
== Relation to other stochastic processes == |
== Relation to other stochastic processes == |
||
If |
If <math display="inline">W(t)</math> is a standard Wiener process (i.e., for <math display="inline">t \geq 0</math>, <math display="inline">W(t)</math> is [[normal distribution|normally distributed]] with expected value <math display="inline">0</math> and variance <math display="inline">t</math>, and the [[Lévy process|increments are stationary and independent]]), then |
||
: <math> B(t) = W(t) - \frac{t}{T} W(T)\,</math> |
: <math> B(t) = W(t) - \frac{t}{T} W(T)\,</math> |
||
is a Brownian bridge for |
is a Brownian bridge for <math display="inline">t \in [0, T]</math>. It is independent of <math display="inline">W(T) </math><ref>Aspects of Brownian motion, Springer, 2008, R. Mansuy, M. Yor page 2</ref> |
||
Conversely, if |
Conversely, if <math display="inline">B(t)</math> is a Brownian bridge for <math display="inline">t \in [0, 1]</math> and <math display="inline">Z</math> is a standard [[normal distribution|normal]] random variable independent of <math display="inline">B</math>, then the process |
||
: <math>W(t) = B(t) + |
: <math>W(t) = B(t) + tZ\,</math> |
||
is a Wiener process for |
is a Wiener process for <math display="inline">t \in [0, 1]</math>. More generally, a Wiener process <math display="inline">W(t)</math> for <math display="inline">t \in [0, T]</math> can be decomposed into |
||
: <math>W(t) = B\left(\frac{t}{T}\right) + \frac{t}{\sqrt{T}} Z.</math> |
: <math>W(t) = \sqrt{T}B\left(\frac{t}{T}\right) + \frac{t}{\sqrt{T}} Z.</math> |
||
Another representation of the Brownian bridge based on the Brownian motion is, for |
Another representation of the Brownian bridge based on the Brownian motion is, for <math display="inline">t \in [0, T]</math> |
||
: <math> B(t) = \frac{ |
: <math> B(t) = \frac{T-t}{\sqrt T} W\left(\frac{t}{T-t}\right).</math> |
||
Conversely, for |
Conversely, for <math display="inline">t \in [0, \infty]</math> |
||
: <math> W(t) = \frac{T+t}{T} B\left(\frac{Tt}{T+t}\right).</math> |
: <math> W(t) = \frac{T+t}{T} B\left(\frac{Tt}{T+t}\right).</math> |
||
Line 32: | Line 33: | ||
The Brownian bridge may also be represented as a Fourier series with stochastic coefficients, as |
The Brownian bridge may also be represented as a Fourier series with stochastic coefficients, as |
||
: <math> B_t = \sum_{k=1}^\infty Z_k \frac{\sqrt{2} \sin(k \pi t / T)}{k \pi}</math> |
: <math> B_t = \sum_{k=1}^\infty Z_k \frac{\sqrt{2 T} \sin(k \pi t / T)}{k \pi}</math> |
||
where <math> Z_1, Z_2, \ldots </math> are [[independent identically distributed]] standard normal random variables (see the [[Karhunen–Loève theorem]]). |
where <math> Z_1, Z_2, \ldots </math> are [[independent identically distributed]] standard normal random variables (see the [[Karhunen–Loève theorem]]). |
||
A Brownian bridge is the result of [[Donsker's theorem]] in the area of [[empirical process]]es. It is also used in the [[Kolmogorov–Smirnov test]] in the area of [[statistical inference]]. |
A Brownian bridge is the result of [[Donsker's theorem]] in the area of [[empirical process]]es. It is also used in the [[Kolmogorov–Smirnov test]] in the area of [[statistical inference]]. |
||
Let <math>K=\sup_{t\in[0,1]}|B(t)|</math>, then the [[cumulative distribution function]] of <math display="inline">K</math> is given by<ref>{{Cite journal |vauthors=Marsaglia G, Tsang WW, Wang J |year=2003 |title=Evaluating Kolmogorov's Distribution |journal=Journal of Statistical Software |volume=8 |issue=18 |pages=1–4 |doi=10.18637/jss.v008.i18 |doi-access=free}}</ref><math display="block">\operatorname{Pr}(K\leq x)=1-2\sum_{k=1}^\infty (-1)^{k-1} e^{-2k^2 x^2}=\frac{\sqrt{2\pi}}{x}\sum_{k=1}^\infty e^{-(2k-1)^2\pi^2/(8x^2)}.</math> |
|||
== Intuitive remarks == |
== Intuitive remarks == |
||
A standard Wiener process satisfies ''W''(0) = 0 and is therefore "tied down" to the origin, but other points are not restricted. In a Brownian bridge process on the other hand, not only is ''B''(0) = 0 but we also require that ''B''(T) |
A standard Wiener process satisfies ''W''(0) = 0 and is therefore "tied down" to the origin, but other points are not restricted. In a Brownian bridge process on the other hand, not only is ''B''(0) = 0 but we also require that ''B''(''T'') = 0, that is the process is "tied down" at ''t'' = ''T'' as well. Just as a literal bridge is supported by pylons at both ends, a Brownian Bridge is required to satisfy conditions at both ends of the interval [0,''T'']. (In a slight generalization, one sometimes requires ''B''(''t''<sub>1</sub>) = ''a'' and ''B''(''t''<sub>2</sub>) = ''b'' where ''t''<sub>1</sub>, ''t''<sub>2</sub>, ''a'' and ''b'' are known constants.) |
||
Suppose we have generated a number of points ''W''(0), ''W''(1), ''W''(2), ''W''(3), etc. of a Wiener process path by computer simulation. It is now desired to fill in additional points in the interval [0,T], that is to interpolate between the already generated points ''W''(0) and ''W''(T). The solution is to use a Brownian bridge that is required to go through the values ''W''(0) and ''W''(T). |
Suppose we have generated a number of points ''W''(0), ''W''(1), ''W''(2), ''W''(3), etc. of a Wiener process path by computer simulation. It is now desired to fill in additional points in the interval [0,''T''], that is to interpolate between the already generated points ''W''(0) and ''W''(''T''). The solution is to use a Brownian bridge that is required to go through the values ''W''(0) and ''W''(''T''). |
||
== General case == |
== General case == |
||
For the general case when '' |
For the general case when ''W''(''t''<sub>1</sub>) = ''a'' and ''W''(''t''<sub>2</sub>) = ''b'', the distribution of ''B'' at time ''t'' ∈ (''t''<sub>1</sub>, ''t''<sub>2</sub>) is [[normal distribution|normal]], with [[expected value|mean]] |
||
:<math>a + \frac{t-t_1}{t_2-t_1}(b-a)</math> |
:<math>a + \frac{t-t_1}{t_2-t_1}(b-a)</math> |
||
and [[variance]] |
|||
:<math>\frac{(t_2-t)(t-t_1)}{t_2-t_1},</math> |
|||
and the [[covariance]] between ''B''(''s'') and ''B''(''t''), with ''s'' < ''t'' is |
and the [[covariance]] between ''B''(''s'') and ''B''(''t''), with ''s'' < ''t'' is |
||
Line 58: | Line 65: | ||
{{Stochastic processes}} |
{{Stochastic processes}} |
||
{{Authority control}} |
|||
[[Category:Wiener process]] |
[[Category:Wiener process]] |
Latest revision as of 02:35, 23 October 2024
A Brownian bridge is a continuous-time gaussian process B(t) whose probability distribution is the conditional probability distribution of a standard Wiener process W(t) (a mathematical model of Brownian motion) subject to the condition (when standardized) that W(T) = 0, so that the process is pinned to the same value at both t = 0 and t = T. More precisely:
The expected value of the bridge at any t in the interval [0,T] is zero, with variance , implying that the most uncertainty is in the middle of the bridge, with zero uncertainty at the nodes. The covariance of B(s) and B(t) is , or s(T − t)/T if s < t. The increments in a Brownian bridge are not independent.
Relation to other stochastic processes
[edit]If is a standard Wiener process (i.e., for , is normally distributed with expected value and variance , and the increments are stationary and independent), then
is a Brownian bridge for . It is independent of [1]
Conversely, if is a Brownian bridge for and is a standard normal random variable independent of , then the process
is a Wiener process for . More generally, a Wiener process for can be decomposed into
Another representation of the Brownian bridge based on the Brownian motion is, for
Conversely, for
The Brownian bridge may also be represented as a Fourier series with stochastic coefficients, as
where are independent identically distributed standard normal random variables (see the Karhunen–Loève theorem).
A Brownian bridge is the result of Donsker's theorem in the area of empirical processes. It is also used in the Kolmogorov–Smirnov test in the area of statistical inference.
Let , then the cumulative distribution function of is given by[2]
Intuitive remarks
[edit]A standard Wiener process satisfies W(0) = 0 and is therefore "tied down" to the origin, but other points are not restricted. In a Brownian bridge process on the other hand, not only is B(0) = 0 but we also require that B(T) = 0, that is the process is "tied down" at t = T as well. Just as a literal bridge is supported by pylons at both ends, a Brownian Bridge is required to satisfy conditions at both ends of the interval [0,T]. (In a slight generalization, one sometimes requires B(t1) = a and B(t2) = b where t1, t2, a and b are known constants.)
Suppose we have generated a number of points W(0), W(1), W(2), W(3), etc. of a Wiener process path by computer simulation. It is now desired to fill in additional points in the interval [0,T], that is to interpolate between the already generated points W(0) and W(T). The solution is to use a Brownian bridge that is required to go through the values W(0) and W(T).
General case
[edit]For the general case when W(t1) = a and W(t2) = b, the distribution of B at time t ∈ (t1, t2) is normal, with mean
and variance
and the covariance between B(s) and B(t), with s < t is
References
[edit]- ^ Aspects of Brownian motion, Springer, 2008, R. Mansuy, M. Yor page 2
- ^ Marsaglia G, Tsang WW, Wang J (2003). "Evaluating Kolmogorov's Distribution". Journal of Statistical Software. 8 (18): 1–4. doi:10.18637/jss.v008.i18.
- Glasserman, Paul (2004). Monte Carlo Methods in Financial Engineering. New York: Springer-Verlag. ISBN 0-387-00451-3.
- Revuz, Daniel; Yor, Marc (1999). Continuous Martingales and Brownian Motion (2nd ed.). New York: Springer-Verlag. ISBN 3-540-57622-3.