Color normalization

In probability theory, reflected Brownian motion (or regulated Brownian motion, both with the acronym RBM) is a Wiener process in a space with reflecting boundaries.

RBMs have been shown to describe queueing models experiencing heavy traffic as first proposed by Kingman and proven by Iglehart and Whitt.

Definition

A d–dimensional reflected Brownian motion Z is a stochastic process on ${\mathbb{ℝ}}_{+}^d$ uniquely defined by

• a d–dimensional drift vector μ
• a d×d non-singular covariance matrix Σ and
• a d×d reflection matrix R.

where X(t) is an unconstrained Brownian motion and

$Z(t)=X(t)+RY(t)$

with Y(t) a d–dimensional vector where

• Y is continuous and non–decreasing with Y(0) = 0
• Y_j only increases at times for which Z_j = 0 for j = 1,2,...,d
• Z(t) ∈ S, t ≥ 0.

The reflection matrix describes boundary behaviour. In the interior of ${\mathbb{ℝ}}_{+}^d$ the process behaves like a Wiener process, on the boundary "roughly speaking, Z is pushed in direction R^j whenever the boundary surface $\{z\in {\mathbb{ℝ}}_{+}^d:z_j=0\}$ is hit, where R^j is the jth column of the matrix R."

Stability conditions

Stability conditions are known for RBMs in 1, 2, and 3 dimensions. "The problem of recurrence classification for SRBMs in four and higher dimensions remains open." In the special case where R is an M-matrix then necessary and sufficient conditions for stability are

1. R is a non-singular matrix and
2. R⁻¹μ < 0.

Stationary distribution

One dimension

The transient distribution of a one-dimensional Brownian motion starting at 0 restricted to positive values (a single reflecting barrier at 0) with drift μ and variance σ² is

$\mathbb{\mathbb{P}}(Z\le z)=\mathrm{\Phi }\left(\frac{z-\mu t}{\sigma t^{1/2}}\right)-e^{2\mu z/{\sigma }^2}\mathrm{\Phi }\left(\frac{-z-\mu t}{\sigma t^{1/2}}\right)$

for all t ≥ 0, (with Φ the cumulative distribution function of the normal distribution) which yields (for μ < 0) when taking t → ∞ an exponential distribution

$\mathbb{\mathbb{P}}(Z<z)=1-e^{2\mu z/{\sigma }^2}.$

Multiple dimensions

The stationary distribution of a reflected Brownian motion in multiple dimensions is tractable analytically when there is a product form stationary distribution, which occurs when the process is stable and

$2\mathrm{\Sigma }=RD+D{R}^{\prime }$

where D = diag(Σ). In this case the probability density function is

$p(z_1,z_2,\dots ,z_d)=\prod _{k=1}^d{\eta }_k{e}^{-{\eta }_k{z}_k}$

where η_k = 2μ_kγ_k/Σ_kk and γ = R⁻¹μ. Closed-form expressions for situations where the product form condition does not hold can be computed numerically as described below in the simulation section.

Hitting times

One dimension

Write T(y) for the first time a one dimensional RBM starting at 0 reaches the level y. Then

$\mathbb{\mathbb{P}}(T(y)>t)=\mathrm{\Phi }\left(\frac{y-\mu t}{\sigma t^{1/2}}\right)-e^{2\mu y/{\sigma }^2}\mathrm{\Phi }\left(\frac{-y-\mu t}{\sigma t^{1/2}}\right).$

Simulation

One dimension

In one dimension the simulated process is the absolute value of a Wiener process. The following MATLAB program creates a sample path.

%rbm.m
n=10^4; h=10^(-3); t=h.*(0:n); mu=-1;
X=zeros(1,n+1); M=X; B=X;
B(1)=3; X(1)=3;
for k=2:n+1
  Y=sqrt(h)*randn; U=rand(1);
  B(k)=B(k-1)+mu*h-Y;
  M=(Y + sqrt(Y^2-2*h*log(U)))/2;
  X(k)=max(M-Y,X(k-1)+h*mu-Y);
end
subplot(2,1,1)
plot(t,X,'k-');
subplot(2,1,2)
plot(t,X-B,'k-');

The error involved in discrete simulations has been quantified. An exact simulation method for time-dependent RBM is also proposed.

Multiple dimensions

QNET allows simulation of steady state RBMs.

Other boundary conditions

Feller described possible boundary condition for the process

• absorption or killed Brownian motion, a Dirichlet boundary condition
• instantaneous reflection, as described above a Neumann boundary condition
• elastic reflection, a Robin boundary condition
• delayed reflection (the time spent on the boundary is positive with probability one)
• partial reflection where the process is either immediately reflected or is absorbed
• sticky Brownian motion.

See also

• Skorokhod problem

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.