Introduction to Stochastic Differential Equations
Introduction to Stochastic Differential Equations
Equations
1 Introduction 1
2 Stochastic Integrals 2
2.1 Formalizing the ‘noise’ 2
2.2 Preparing the Terrain 3
2.3 Constructing the Itô Integral 5
2.4 Properties 9
2.5 Extensions 10
2.6 Stratonovich integral 10
A Measure Theory 14
A.1 Measurable Spaces 14
A.2 Measures 14
A.3 Measurable Functions 17
A.4 Integration 19
A.5 Product Spaces 23
A.6 The Monotone Class Lemma 24
A.7 Modes of Convergence 25
A.8 Conditional Expectation 26
B Probability Theory 28
B.1 Probability Spaces 28
B.2 Distribution and Density of a Random Variable 30
B.3 Expected Value and Variance 31
B.4 Stochastic Processes 34
B.5 Brownian Motion 35
B.6 Riemann-Stieltjes Integral 36
B.7 Modes of Convergence 37
References 39
2
1 Introduction
dX
= b(t, Xt ) + σ(t, Xt ) · ‘white noise’ (1)
dt
1
2 Stochastic Integrals
To study the equation (1), the first step is to give a formalization for the
‘noise’ that we are going to consider. In our case, we’ll consider the ‘white
noise’, which is a generalized stochastic process Wt with gaussian probability
distribution with mean zero, finite variance such that E[Ws Wt ] = 0 whenever
t ̸= s.
In other words, the ‘noise’ will be represented by some stochastic process
Wt and our goal is to satisfy the following properties:
• The collection {Wt } is stationary, i.e., the joint distribution of any collec-
tions of Wt does not depend on t.
However, this process Wt does not have reasonable properties (does not
have continuous paths, and it is not measurable for finite time). So our first
goal is to replace Wt with some convenient stochastic process.
Let 0 = t0 < t1 < . . . < tm = t and consider the discrete version
2
Notice that
Bt – Bs = Ws (t – s)
for all t ≥ s, i.e., the ‘white noise’ is the derivative, with respect to time, of the
Brownian motion.
0.5
–0.5
–1
0 50 100 150 200
Now we can start treating the equation. In order to do that, notice that in
the determinist case, in which there isn’t a ‘noise’, we have an equation of the
form
dXt
= b(t, Xt ), X0 = x0 (2)
dt
In this case, a solution is a function Xt such that
ZT
Xt – X0 = b(s, Xs ) ds (3)
0
3
where the left integral is a Lebesgue-Stieltjes’ integral, but the right integral is
what still needs to be formalized.
As said earlier, our task is to understand the following integral
ZT
f (t, ω) dBt (ω) (5)
S
where f : [0, ∞] × Ω −→ R.
The idea will be the following: we’ll first define the integral for the most
elementary functions and, after that, using approximations and convergence
criteria, we’ll define it for the other classes of functions.
Since it is natural to approximate a given function f (t, ω) by
X
f (tj∗ , ω)χ[tj ,tj+1 ) (t)
j
where tj∗ ∈ [tj , tj+1 ], we could define the integral (5) as follows
ZT
f (t, ω) dBt (ω) = lim f (tj∗ , ω)[Btj+1 – Btj ](ω)
n→∞
S
However, the choice of points tj∗ make a difference here. If we take tj∗ = tj ,
the left end point, then we obtain the Itô integral, denoted by
ZT
f (t, ω) dBt (ω)
S
If, on the other hand, we take tj∗ = (tj + tj+1 )/2, the mid point, then we
obtain the Stratonovich integral, denoted by
ZT
f (t, ω) ◦ dBt (ω)
S
4
Definition 2.1. Let Bt (ω) be a Brownian motion. We define Ft as the σ-algebra
generated by the random variables Bs , where s ≤ t. I.e., Ft is the smallest σ-
algebra containing all sets
The intuition behind this definition is that Ft can be thought as the ‘history
of Bs up to the time t’. A function h(ω) is Ft -measurable if h depends only on
the values from B0 up to Bt . I.e., it does not depend on the ‘future’.
Definition 2.2 (Nt -Adapted). Let {Nt }t≥0 be an increasing family of σ-algebras
of subsets of Ω. A process g(t, ω) : [0, ∞)×Ω −→ Rn is Nt -adapted if for each
t ≥ 0 the function ω −→ g(t, ω) is Nt -measurable.
Example 2.1. The process h1 (t, ω) = Bt/2 (ω) is Ft -measurable and, hence,
Ft -adapted. Notice that this process only depends on previous information.
However, the process h2 (t, ω) = B2t (ω) is not Ft -measurable and, hence,
it is not Ft -adapted. That’s because the process h2 depends on future informa-
tion.
The class of functions for which the Itô integral will be first defined is the
following.
f (t, ω) : [0, ∞) × Ω −→ R
such that
5
Given functions f ∈ V, in order to define the Ito integral
ZT
I[f ](ω) = f (t, ω) dBt (ω)
S
we are going to define I[ϕ] for a simple class of functions ϕ and then show
how each f ∈ V can be approximated by ϕ.
To further our development of the Ito integral, we need the following result.
1. Let g(ω) ∈ V be bounded and continuous for each ω. Then there exists
elementary functions ϕn ∈ V such that
Z T
E (g – ϕn )2 dt −→ 0
S
as n → ∞.
6
3. Let f ∈ V. Then there exists a sequence of functions (hn ) ⊆ V such that
hn is bounded for each n and
Z T
E (f – hn )2 dt −→ 0
S
as n → ∞.
Please note that we started from bounded and continuous functions, gen-
eralized into bounded functions and then generalized further into any function
in our class V.
These steps can be summarized in the following result, which guarantees
that the elementary functions are dense in V(S, T).
Theorem 2.2. If g ∈ V(S, T), then there exists a sequence of elementary func-
tions ϕn ∈ V(S, T) such that ϕn → g in L2 (P)
where the limit is in L2 (P) and (ϕn ) is a sequence of elementary functions such
that
T
Z
E (f (t, ω) – ϕn (t, ω))2 dt −→ 0
S
as n → ∞.
7
as n → ∞, then
ZT ZT
fn (t, ω) dBt (ω) −→ f (t, ω) dB)t(ω)
S S
in L2 (P) as n → ∞.
Example 2.2. We wish to compute the integral
ZT
B(t) dB(t)
0
Step 1. Let
n–1
X
ϕn (t) = Bn (tin )[Bn (ti+1
n
) – Bn (tin )]
i=0
be a sequence of elementary functions. Then,
Z T
X Z tj+1 X Z tj+1
2 2
E (ϕn – Bs ) ds = E (Bj – Bs ) ds = (s – tj )2 ds
0 j tj j tj
X1 2
= tj+1 – tj −→ 0 when ∆tj → 0
2
j
n – tn ) → 0 by the conti-
I.e., ϕn (t) converges to B(t) almost surely as maxi (ti+1 i
nuity of B(t).
Step 2. Now notice that
Bn (ti )[Bn (ti+1 ) – Bn (ti )] = Bn (ti+1 )Bn (ti ) – B2n (ti ) + B2n (ti+1 ) – B2n (ti+1 )
1
= [B2n (ti+1 ) – B2n (ti ) – (Bn (ti+1 ) – Bn (ti ))2 ]
2
Step 3. With these two results, we can write the original integral as
ZT n–1 n–1
1 X 2 1X
Xn (t) dB(t) = Bn (ti+1 ) – B2n (ti ) – [Bn (ti+1 ) – Bn (ti )]2
0 2 2
i=0 i=0
Step 4. Since the first sum is a telescopic sum and the second one converges
in probability to T by the quadratic variation of Brownian motion, we have
ZT
1 1 1 1 1
Xn (t) dB(t) = B2 (t) – B2 (0) – T = B2 (t) – T
0 2 2 2 2 2
8
Hence, the integral converges in L2 (P) to
ZT ZT
1 1
B(t) dB(t) = lim Xn (t) dB(t) = B2 (t) – T
0
n→∞
0 2 2
2.4 Properties
Before heading to more theoretical and important results, let us notice some
natural facts of the Itô integral.
With the Doob’s martingale inequality, it can be proved that the Itô integral
can be chosen to depend continuously on t. A proof of this fact is on the third
chapter of [Oks13].
Now an important result is that the Itô integral is a martingale.
9
2.5 Extensions
Using the concept of martingales, we can generalize the Itô integral for a
larger class of functions than V.
Considering the Definition 2.3, we can relax the condition 2. into
The idea here is that Bt must remain a martingale with respect to the history
of fs .
Another way of extending the Itô integral definition is by weaking the con-
dition 3. into
R T
3.’ P f (s, ω) 2 ds < ∞ = 1
S
Too understand why this works, let WH (S, T) denote the class of processes
satisfying conditions 1, 2’ and 3’ above. Then, in the 1-dimensional Brownian
motion, for all t there exists step functions fn ∈ WH such that
Zt
|fn – f |2 ds −→ 0
0
in probability.
Rt
Now, since 0 fn (s, ω) dBs (ω) converges in probability to a random variable
and the limit depends only on f , we define
Zt Zt
f (s, ω) dBs (ω) = lim fn (s, ω) dBs (ω)
n→∞
0 0
10
The Stratonovich integral is another possibility and, in general, leads to
different results (except when the integrating function has a derivative). In
some cases, the Stratonovich definition may be adequate.
In other cases, Itô’s feature of ‘not looking in the future’ (as [Oks13] puts
it), justifies its use in biology and finance, for example. Moreover, Itô integral is
a martingale, and Stratonovich’s is not. This gives an important computational
advantage to Itô’s definition.
And, to quote [S+ 04],
11
3 The Ito Formula and the Martingale Representation
Theorem
12
4 Stochastic Differential Equations
13
A Measure Theory
1. ;, S ∈ S.
2. If A ∈ S, then Ac ∈ S.
S∞
3. If A1 , A2 , . . . ∈ S, then i=1 Ai ∈ S.
For any set S and any collection A of subsets of S there is at least one σ-
algebra containing A: the family of all subsets of S. Taking the intersection of
all the σ-algebras containing A, we obtain the smallest σ-algebra containing
A, which is called the σ-algebra generated by A.
Particularlly, the smallest σ-algebra containing all of the open sets of R is
called the Boreal algebra for R and is denoted by B. Any set in B is called a
Borel set.
A.2 Measures
Now, how can we assign a size (or a probability) to all the sets in S?
1. µ(;) = 0;
14
The triple (S, S, µ) is called a measure space.
1. If (An )∞
n=1 is an increasing sequence in S (i.e. if An ⊆ An+1 for all n),
then
∞
[
µ An = lim µ(An )
n→∞
n=1
2. If (Bn )∞
n=1 is a decreasing sequence in S (i.e. if Bn+1 ⊆ Bn for all n) and
if µ(Bm ) < ∞ for some m, then
∞
\
µ Bn = lim µ(Bn )
n→∞
n=1
1. ;, S ∈ A;
2. A ∈ A implies Ac ∈ A;
15
Sn
3. A1 , A2 , . . . , An ∈ A implies i=1 Ai ∈ A.
1. µ(;) = 0;
Notice that last item is different from the definition on a σ-algebra. Here,
we require that the union is contained in the algebra.
Although it is easier to define measures on algebras, it is more convenient
to work with σ-algebras.
Theorem A.2 (Caratheodory Extension Theorem). Let S be a set, A an algebra
of its subsets, and µ a measure on A. Let S be the smallest σ-algebra contain-
ing A. Then there exists a measure µ∗ on S such that µ(A) = µ∗ (A), for all
A ∈ A.
To rule out the possibility of more than one extension of µ to S, we’ll need
the following definition
Definition A.5. Let S be a set, A an algebra of its subsets, and µ a measure on
A. If there is a countable sequence of sets (Ai )∞
i=1
∈ A with µ(Ai ) < ∞ for all
S∞
i, and S = i=1 Ai , then µ is called σ-finite.
16
The completion of S is the family S′ constructed by starting with any set
B ∈ S, and then adding and subtracting from it sets in C. That is
S′ = {B′ ⊆ S : B′ = (B ∪ C1 ) \ C2 , B ∈ S, C1 , C2 ∈ C}
{s ∈ S : f (s) ≤ a} ∈ S, ∀a ∈ R
17
where (Ai )ni=1 is a sequence of subsets of S, and (ai )ni=1 is a sequence of real
numbers. Functions like this are called simple functions.
If the sets (Ai ) form a partition of S and if all of the values ai are distinct,
then the previous equation is the standard representation of the function. If
this is the case, then ϕ is measurable iff. each Ai ∈ S.
As we’ll see in the next two theorems, the set of all measurable functions
consists of those that are pointwise limits of measurable simple functions.
First, any function f that is pointwise limit of a sequence (fn ) of measurable
function is itself measurable.
Theorem A.4 (Pointwise convergence preserves measurability). Let (S, S) be
a measurable space, and let (fn ) be a sequence of S-measurable functions con-
verging pointwise to f , i.e.,
lim f (s) = f (s), ∀s ∈ S
n→∞ n
With these two theorems, we see that a function is measurable iff. it is the
pointwise limit of a sequence of measurable simple functions. The standard
way to prove that a function is measurable is to find such an approximation.
Some important properties follow. Suppose that f and g are S-measurable
functions on S and c ∈ R. Then,
18
• The functions f + g, f · g, |f | and c · f are S-measurable.
Definition A.8. Let (S, S) and (T, T) be measurable spaces. Then the function
f : S −→ T is measurable if the inverse image of every measurable set is
measurable, i.e.,
{s ∈ S : f (s) ∈ A} ∈ S, ∀A ∈ T
Hence, if (S, S), (T, T), and (U, U) are measurable spaces, and f : S −→ T
and g : T −→ U are measurable functions, then h : S −→ U defined by h(s) =
g ◦ f (s) is a measurable function.
A.4 Integration
Let (S, S, µ) be a fixed measure space and M(S, S) be the space of measur-
able, extended real-valued functions on S.
For a nonnegative, measurable and simple function, we have the following
definition.
19
To extended from simple functions to all M+ (S, S), we define
Definition A.11. For f ∈ M+ (S, S), the integral of f with respect to µ is
Z Z
f (s) µ(ds) = f (s)χA (s)µ(ds)
A S
R R
We also denote the integrals above as f dµ and A f dµ.
and Z Z
c · f dµ = c f dµ
20
3. If f ∈ M+ (S, S) and A, B ∈ S with A ⊆ B, then
Z Z
f dµ ≤ f dµ
A B
µ(A) = 0.
21
In probability, if (S, S, µ) is a probability space and f is a measurable func-
tion, then f is a random variable. And if f is integrable, the integral of f with
respect to µ is called the expected value of f and is denoted by
Z
E[f ] = f dµ
Now let L(S, S, µ) denote the set of all S-measurable, real-valued functions
on S that are integrable with respect to µ. How can we determine whether a
function is integrable?
3. If f isRS-measurable,
R g is µ-integrable, and |f | ≤ |g|, then f is µ-integrable
and |f | dµ ≤ |g| dµ.
µ(A) = 0 =⇒ λ(A) = 0 ∀A ∈ S
22
Theorem A.11 (Radon-Nikodym Theorem). Let λ and µ be σ-finite positive
measures on (S, S) with λ << µ. Then there exists an integrable function h
such that Z
λ(A) = h(s) µ(ds), all A ∈ S
A
And the function h is called the Radon-Nikodym derivative of λ with respect
to µ.
Lemma A.12. Let λ1 and λ2 be finite measures on (S, S). Then there exists a
triple of measures γ, α1 , α2 such that
λi = γ + α1
where i = 1, 2, and α1 ⊥ α2 .
We denote by C the set of all measurable rectangles and E the set of all finite
unions of measurable rectangles. Notice that E is an algebra and that every set
in E can be written as the finite union of disjoint measurable rectangles.
Theorem A.13. Let (X, X), (Y, Y), C, and E be as specified above. Let µ : C −→
R+ have the following properties:
• µ(;) = 0;
S∞
• If (Ci ) = ((Ai × Bi ))∞ is a sequence of disjoint sets in C and i=1 Ci is
S∞ i=1 P∞
in C, then µ i=1 Ci = i=1 µ(Ci ).
23
This theorem can be naturally extended to any space that is the product of
a finite number of measurable spaces.
Definition A.16. Let (X, X) and (Y, Y) be measurable spaces, and let (Z, F) be
the product space. And let E ⊆ Z and x ∈ X. Then the x-section of E is the set
Ex = y ∈ Y : (x, y) ∈ E. And the y-section of E is the set Ey = x ∈ X : (x, y) ∈ E.
Let f : Z −→ R and let x ∈ X. The x-section of f is the function fx : Y −→ R
defined by fx (y) = f (x, y). And the y-section of f is the function fy : X −→ R
defined by fy (x) = f (x, y).
Theorem A.14. Let (X, X) and (Y, Y) be measurable spaces, and let (Z, F) be
the product space.
The main question that this section aims to answer is whether a propery P
holds at a given σ-algebra.
24
3. The σ-algebra generated by A contains the monotone class generated by
A.
Using this lemma, our original task in this section is solved. To establish
that P holds for all sets in a product σ-algebra, it suffices to show that
3. almost uniformly if, for every ϵ > 0, there exists a set E ∈ B of measure
µ(E) ≤ ϵ such that fn converges uniformly to f on the complement of E.
25
5. in measure if, for every ϵ > 0, the measures
converge to zero as n → ∞.
1. Aη ∈ F, ∀η ∈ H.
2. η∈H Aη = H.
S
3. Aη ∩ Aη′ = 0, ∀η ̸= η′ .
If the index set H is countable, then we call the measurable partition countable.
26
To formalize it better, extending the concept to σ-algebras not generated
by a countable partition, we’ll introduce the following.
Definition A.21 (Conditional Expectation). Let A ⊂ F be a σ-algebra and let
f : Ω −→ R be an integrable function. Then the conditional expectation of f
relative to A is an A-measurable function E[f |A] : Ω −→ R such that
Z Z
E[f |A](ω)µ (dω) = f (ω)µ (dω), ∀C ∈ A
C C
It can be proved that this function exists and is unique in the sense that if
g also satisties the condition above, then E[f |A] = g.
Notice that to compute conditional probability P(B|Aη ) we may take f to
be the indicator function χB . Then
Z
E[χB |A](ω̂)µ(Aη ) = χb dµ = µ(B ∩ Aη )
Aη
27
B Probability Theory
A probability space starts with a set Ω called the sample space, which is,
intuitively, a list of all possible outcomes of an experiment. Each ω ∈ Ω is a
sample point, and each subset A ⊆ Ω is an event.
To filter ‘well-behaved’ subsets, where it will be possible to measure a prob-
ability, the following definition is necessary.
1. ; ∈ F.
2. If A ∈ F, then Ac ∈ F.
S∞
3. If A1 , A2 , . . . ∈ F, then i=1 Ai ∈ F.
The pair (Ω, F) is said to be a measurable space, and any subset B ⊆ Ω that
also belongs to F is called a measurable set.
1. The family {;, Ω} is called the trivial σ-algebra, and is the smallest one
possible.
2. The power set P(Ω) is called the discrete σ-algebra, and is the largest
one possible.
But the most important σ-algebra for probability theory is the Borel σ-
algebra, which will be denoted by B. Taking Ω = R, the Borel σ-algebra is
generated by the intersection of all σ-algebras containing the real line inter-
vals.
Notice that the Borel σ-algebra contains all open sets, closed sets, and all
their countable operations with union ∪, intersection ∩, and their complements
c . This is the smallest σ-algebra containing all open subsets.
28
More generally, if U is the collection of all open subsets of a topological
space Ω, then B is called the Borel σ-algebra on Ω. The elements B ∈ B are
called Borel sets.
Given a measurable space, it’s possible to assign each outcome to a proba-
bility.
P : F −→ [0, 1]
satisfying
for all open sets U ∈ Rn (i.e. for all Borel sets U ⊆ Rn ). In other words, the
inverse image of U is in the σ-algebra.
This definition means that if we know which event U in the F has occurred,
then we know which value of X has occurred.
29
Consider the measurable space (Ω, P(Ω)), and let X be a random variable
with values xi , i = 1, 2, . . . , k. The sets
Ai = {ω ∈ Ω : X(ω) = xi } ⊆ Ω
Definition B.7 (Joint Distribution). If X and Y are random variables, then their
joint distribution function is
where x, y ∈ R.
The distribution function gives the probability that the random variable X
is on the interval (–∞, x]. This also allows to look for an interval [a, b] using
FX (b) – FX (a).
While the distribution of X returns the probability of an event {ω : X(ω) ∈
B}.
30
Definition B.8 (Density function). A random variable X has a probability den-
sity function (PDF) f : Rn −→ R if f is a measurable function and
Zx
FX (x) = P[X ≤ x] = f (y) dy
–∞
Another useful tool for any given random variable is to know its mean value
and how much it varies. This motivates the following.
Notice that the expectation is linear. I.e., if X and Y are integrable and a
and b are constants, then
1
P[|X| ≥ λ] ≤ p E[|X|p ]
λ
for all λ ≥ 0.
Definition B.10 (Variance). Let µ := E[X], i.e., the expected value of a random
variable X. Then the variance of X is given by
Z
Var[X] = E[∥X – µ∥2 ] = |X – µ|2 dP(ω)
Ω
31
Definition B.11 (Covariance). Let X and Y be integral random variables. Let
µX = E[X] and µY = E[Y]. If XY is integrable, then the covariance of X and Y
is:
Cov[X, Y] = E[(X – µX )(Y – µY )] = E[XY] – µX µY
Notice that Var[X] = Cov[X, X].
The variance of X can be computed more simply using the following theo-
rem.
2. Var[aX + b] = a2 Var[X].
f (x, y)
f (x|y) =
fY (y)
32
The intuition behind this definition is to build an estimate of the random
variable X given the information available in Y.
Some important and very useful properties of conditional expectation are
listed below.
2. E[E[X|H]] = E[X].
3. E[X|H] = X if X is H-measurable.
6. If G ⊆ H is a σ-algebra, then
E[X|G] = E[E[X|H]|G]
ϕ(E[X|H]) ≤ E[ϕ(X)|H]
The next lemma helps us check if some sequence of events occurs infinitely
often.
P∞
Lemma B.4 (Borel-Cantelli). If n=1 P(An ) < ∞, then
∞ [∞
\
P(An i.o.) = P Am = 0
n=1 m=n
The most important distribution in this text is the normal (or gaussian)
distribution.
33
Definition B.15 (Normal Distribution). If random variable X has mean µ, vari-
ance σ2 and a density function of the form
–(x – µ)2
1
f (x) = p exp
2πσ 2σ2
F = {F0 , F1 , . . . , Ft , . . . , FT }
34
1. Mt is Mt -measurable for all t.
In 1828 the Scottish botanist Robert Brown described the irregular motion
of pollen grains suspended in fluid, that is now known as the Brownian Motion.
To describe this mathematically, we use a stochastic process Bt (ω), which is the
position at time t of the particle ω.
This can also be understood as the model for the cumulative effect of ‘noise’.
And is also called Wiener process, after N. Wiener, who formalized mathemat-
ically the Brownian motion.
35
Theorem B.6 (Quadratic variation). If t ≥ s,
E[(Bt – Bs )2 ] = t – s
36
B.7 Modes of Convergence
I.e., Z Z
lim f (Xn ) dPn = f (X) dP
n→∞
Ωn Ω
37
Definition B.26 (Lp -Convergence). A sequence of random variables {Xn } con-
verges in Lp to X if {Xn } ⊆ Lp , p ∈ [1, ∞), and
38
References
[Dur19] Rick Durrett. Probability: theory and examples, volume 49. Cam-
bridge university press, 2019.
[S+ 04] Steven E Shreve et al. Stochastic calculus for finance II: Continuous-
time models, volume 11. Springer, 2004. 11
39