Solution to a stochastic differential equation
In probability theory and statistics , diffusion processes are a class of continuous-time Markov process with almost surely continuous sample paths. Diffusion process is stochastic in nature and hence is used to model many real-life stochastic systems. Brownian motion , reflected Brownian motion and Ornstein–Uhlenbeck processes are examples of diffusion processes. It is used heavily in statistical physics , statistical analysis , information theory , data science , neural networks , finance and marketing .
A sample path of a diffusion process models the trajectory of a particle embedded in a flowing fluid and subjected to random displacements due to collisions with other particles, which is called Brownian motion . The position of the particle is then random; its probability density function as a function of space and time is governed by a convection–diffusion equation .
Mathematical definition
A diffusion process is a Markov process with continuous sample paths for which the Kolmogorov forward equation is the Fokker–Planck equation .[ 1]
A diffusion process is defined by the following properties.
Let
a
i
j
(
x
,
t
)
{\displaystyle a^{ij}(x,t)}
be uniformly continuous coefficients and
b
i
(
x
,
t
)
{\displaystyle b^{i}(x,t)}
be bounded, Borel measurable drift terms. There is a unique family of probability measures
P
a
;
b
ξ
,
τ
{\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}
(for
τ
≥
0
{\displaystyle \tau \geq 0}
,
ξ
∈
R
d
{\displaystyle \xi \in \mathbb {R} ^{d}}
) on the canonical space
Ω
=
C
(
[
0
,
∞
)
,
R
d
)
{\displaystyle \Omega =C([0,\infty ),\mathbb {R} ^{d})}
, with its Borel
σ
{\displaystyle \sigma }
-algebra, such that:
1. (Initial Condition) The process starts at
ξ
{\displaystyle \xi }
at time
τ
{\displaystyle \tau }
:
P
a
;
b
ξ
,
τ
[
ψ
∈
Ω
:
ψ
(
t
)
=
ξ
for
0
≤
t
≤
τ
]
=
1.
{\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }[\psi \in \Omega :\psi (t)=\xi {\text{ for }}0\leq t\leq \tau ]=1.}
2. (Local Martingale Property) For every
f
∈
C
2
,
1
(
R
d
×
[
τ
,
∞
)
)
{\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times [\tau ,\infty ))}
, the process
M
t
[
f
]
=
f
(
ψ
(
t
)
,
t
)
−
f
(
ψ
(
τ
)
,
τ
)
−
∫
τ
t
(
L
a
;
b
+
∂
∂
s
)
f
(
ψ
(
s
)
,
s
)
d
s
{\displaystyle M_{t}^{[f]}=f(\psi (t),t)-f(\psi (\tau ),\tau )-\int _{\tau }^{t}{\bigl (}L_{a;b}+{\tfrac {\partial }{\partial s}}{\bigr )}f(\psi (s),s)\,ds}
is a local martingale under
P
a
;
b
ξ
,
τ
{\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}
for
t
≥
τ
{\displaystyle t\geq \tau }
, with
M
t
[
f
]
=
0
{\displaystyle M_{t}^{[f]}=0}
for
t
≤
τ
{\displaystyle t\leq \tau }
.
This family
P
a
;
b
ξ
,
τ
{\displaystyle \mathbb {P} _{a;b}^{\xi ,\tau }}
is called the
L
a
;
b
{\displaystyle {\mathcal {L}}_{a;b}}
-diffusion, where
L
a
;
b
=
L
a
;
b
+
∂
∂
t
{\displaystyle {\mathcal {L}}_{a;b}=L_{a;b}+{\frac {\partial }{\partial t}}}
is the time‐dependent infinitesimal generator.
Connection to Stochastic Differential Equations
If
(
X
t
)
t
≥
0
{\displaystyle (X_{t})_{t\geq 0}}
is an
L
a
;
b
{\displaystyle {\mathcal {L}}_{a;b}}
-diffusion, it satisfies the SDE
d
X
t
i
=
2
ν
∑
k
=
1
d
σ
k
i
(
X
t
)
d
B
t
k
+
b
i
(
X
t
)
d
t
{\displaystyle dX_{t}^{i}={\sqrt {2\nu }}\sum _{k=1}^{d}\sigma _{k}^{i}(X_{t})\,dB_{t}^{k}+b^{i}(X_{t})\,dt}
,
provided
a
i
j
(
x
)
=
∑
k
=
1
d
σ
k
i
(
x
)
σ
k
j
(
x
)
{\displaystyle a^{ij}(x)=\sum _{k=1}^{d}\sigma _{k}^{i}(x)\sigma _{k}^{j}(x)}
, and
σ
i
j
(
x
)
{\displaystyle \sigma ^{ij}(x)}
,
b
i
(
x
)
{\displaystyle b^{i}(x)}
are Lipschitz continuous. By Itô's lemma, for
f
∈
C
2
,
1
(
R
d
×
[
τ
,
∞
)
)
{\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times [\tau ,\infty ))}
we have
d
f
(
X
t
,
t
)
=
(
∂
f
∂
t
+
∑
i
=
1
d
b
i
∂
f
∂
x
i
+
ν
∑
i
,
j
=
1
d
a
i
j
∂
2
f
∂
x
i
∂
x
j
)
d
t
+
(martingale terms)
.
{\displaystyle df(X_{t},t)={\Bigl (}{\tfrac {\partial f}{\partial t}}+\sum _{i=1}^{d}b^{i}{\tfrac {\partial f}{\partial x_{i}}}+\nu \sum _{i,j=1}^{d}a^{ij}{\tfrac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}{\Bigr )}\,dt+{\text{(martingale terms)}}.}
Infinitesimal Generator
The infinitesimal generator
A
{\displaystyle {\mathcal {A}}}
of
X
t
{\displaystyle X_{t}}
is defined for
f
∈
C
2
,
1
(
R
d
×
R
+
)
{\displaystyle f\in C^{2,1}(\mathbb {R} ^{d}\times \mathbb {R} ^{+})}
by
A
f
(
x
,
t
)
=
∑
i
=
1
d
b
i
(
x
,
t
)
∂
f
∂
x
i
+
ν
∑
i
,
j
=
1
d
a
i
j
(
x
,
t
)
∂
2
f
∂
x
i
∂
x
j
+
∂
f
∂
t
.
{\displaystyle {\mathcal {A}}f(\mathbf {x} ,t)=\sum _{i=1}^{d}b^{i}(\mathbf {x} ,t)\,{\tfrac {\partial f}{\partial x_{i}}}+\nu \sum _{i,j=1}^{d}a^{ij}(\mathbf {x} ,t)\,{\tfrac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}+{\tfrac {\partial f}{\partial t}}.}
Transition Probability Density
For a diffusion process
(
X
t
,
P
a
;
b
ξ
,
τ
)
{\displaystyle (X_{t},\mathbb {P} _{a;b}^{\xi ,\tau })}
, the transition probability function is
H
a
;
b
(
τ
,
ξ
,
t
,
d
x
)
=
P
a
;
b
ξ
,
τ
[
ψ
:
ψ
(
t
)
∈
d
x
]
.
{\displaystyle H_{a;b}(\tau ,\xi ,t,\mathrm {d} x)=\mathbb {P} _{a;b}^{\xi ,\tau }[\psi :\psi (t)\in \mathrm {d} x].}
Under uniform ellipticity of
a
i
j
(
x
,
t
)
{\displaystyle a^{ij}(x,t)}
, this measure has a density
h
a
;
b
(
τ
,
ξ
,
t
,
x
)
{\displaystyle h_{a;b}(\tau ,\xi ,t,x)}
w.r.t. Lebesgue measure, satisfying
∂
h
∂
t
=
A
∗
h
,
{\displaystyle {\tfrac {\partial h}{\partial t}}={\mathcal {A}}^{*}h,}
where
A
∗
{\displaystyle {\mathcal {A}}^{*}}
is the adjoint of the infinitesimal generator.
See also
References
Concepts Applications Implementations
Audio–visual Text Decisional
People Architectures