Math5846 Chapter10
Math5846 Chapter10
UNSW Sydney
OPEN LEARNING
Chapter 10
Brownian Motion
2 / 75
Outline:
10.1 Introduction
3 / 75
Outline - continued:
10.15 Martingales
4 / 75
10.1 Introduction
5 / 75
We will start with a very brief history of Brownian motion.
1827 Brown described the behaviour of random movements of particles
suspended in a medium such as a liquid or a gas.
6 / 75
Sample Paths of Brownian Motions
fractional Brownian motion − parameter: 0.5 fractional Brownian motion − parameter: 0.5
5 0
−2
−4
0 −6
−8
−10
−5 −12
−14
−16
−10 −18
100 200 300 400 500 600 700 800 900 1000 100 200 300 400 500 600 700 800 900 1000
7 / 75
Sample Paths of Brownian Motions
fractional Brownian motion − parameter: 0.5 fractional Brownian motion − parameter: 0.5
35 10
30
0
25
−10
20
15 −20
10 −30
5
−40
0
−50
−5
−10 −60
1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000
8 / 75
Now, we will discuss two mathematical definitions of Brownian motion (BM).
Definition
A continuous stochastic process {X(t), t ≥ 0} is said to be a Brownian
motion with drift coefficient µ if
(1) X(0) = 0
(2) {X(t), t ≥ 0} has stationary independent increments
(3) for every t ≥ 0, X(t) is normally distributed with mean µ t and variance
σ 2 t.
9 / 75
If we write X(t + s) = X(t + s) − X(t) + X(t), then
V ar(X(t + s)) = V ar(X(t + s) − X(t) + X(t))
= V ar(X(t + s) − X(t)) + V ar(X(t))
by independent increments
= V ar(X(s)) + V ar(X(t)) by stationary increments.
10 / 75
10.2 Standard Brownian Motion
11 / 75
We will write B(t) as a standard Brownian motion which means it has
zero drift (mean zero) and σ 2 = 1.
1 n x2 o
fB(t) (x) = √ exp − , t > 0, −∞ < x < ∞,
2π t 2t
12 / 75
Hence, we have
Z ∞
1 x2
P (B(t) ≥ a) = √ e− 2t dx
2 π t Za
y2 √ √
∞
1 x
= √ √
e− 2 t dy (by letting y = √ , so dx = t dy)
2 π t a/ t t
Z ∞
1 y 2
= √ √
e− 2 dy
2 π a/ t
!
a
= 1−Φ √ , a ∈ R, where Φ(x) ∼ N (0, 1).
t
and
a
P (B(t) ≤ a) = Φ √ , a ∈ R, where Φ(x) ∼ N (0, 1).
t
13 / 75
Note that X(t) = B(t) + µt is known as a Brownian motion with drift.
Any Brownian motion X(t) with mean µ t and variance σ 2 t can be converted
to a standard Brownian motion.
That is,
X(t) − µ t
σ
is a standard Brownian motion.
14 / 75
Now, consider a Brownian motion as a particle in a liquid medium.
Therefore, we have
Z ∞
ft (x|x0 ) ≥ 0 and ft (x|x0 ) dx = 1. (1)
−∞
15 / 75
In addition, we say that for small t, X(t + t0 ) is likely to be near X(t0 ) = x0 .
That is,
From the physical principles, Einstein (1905) showed that ft (x|x0 ) must
satisfy the following partial differential equation
∂f ∂ 2f
= D 2. (3)
∂t ∂x
The above is known as the diffusion equation, and D is called the diffusion
coefficient.
16 / 75
Actually, D = 2RT
Nh
, where R is the gas constant, T is temperature, N is
Avogadro’s number and h is a coefficient of friction.
In fact, it is the unique solution under the boundary conditions (1) and (2).
17 / 75
10.3 Another Approach to Defining Brownian
Motion
18 / 75
Another approach to Equation (3) is an approximation using a discrete
random walk.
Let pk (n) = P ( ni=1 ξi = k) be the probability that the random walk finds
P
itself k steps to the right of its starting point at time n.
19 / 75
Then by the Chapman-Kolmogorov equation and the Central Limit
Theorem , we have
1 1
pk (n + 1) = pk+1 (n) + pk−1 (n)
2 2
1
pk (n + 1) − pk (n) = [pk+1 (n) − 2pk (n) + pk−1 (n)]. (5)
2
Let the length between transitions be denoted by ∆ and the length of each
step η.
20 / 75
Now let ∆ and η shrink to zero preserving the relationship ∆ = η 2 and at the
same time let n and k increase to ∞ so kη → x while n ∆ → t.
21 / 75
Another limiting process for pk (n) requires the CLT. If we write
pk (n) = P (ξ1 + · · · + ξn = k), where
+1 with probability 1/2
ξi =
−1 with probability 1/2.
23 / 75
Before we continue, let us state some other invariance principles. Let X(t),
t ≥ 0 be a Brownian motion starting at x (i.e. X0 = x).
24 / 75
10.5 Other Properties of Brownian Motion
25 / 75
Without any proofs, some other properties of Brownian Motion are
sample paths of Brownian motion are nowhere differentiable.
26 / 75
(i)
Let Bt , i = 1, 2 be two independent Brownian motions and a1 , a2 ∈ R,
where both are non-zero.
Then,
(1) (2
B̃t (ω) ≡ (a21 + a22 )−1/2 (a1 Bt (ω) + a2 Bt (ω)), t ≥ 0,
is a Brownian motion.
The random variable σ B(t) + a t, t ≥ 0, is Gaussian random variable
N (a t, σ 2 t) for any constants σ and a.
This process is called a Brownian motion with linear drift or Brownian
motion with drift.
The constant σ is called the volatility and it contributes to the order of
the variance function σ(t) = σ 2 t.
27 / 75
10.6 Joint Distribution of Brownian Motion
28 / 75
Since we know the distribution of B(t), t ≥ 0, we can also determine the joint
distribution of
B(tk ) − B(tk−1 )
29 / 75
Notation
D
B(tk ) − B(tk−1 ) −
→ B(tk − tk−1 ).
D
Recall the definition of −
→ Convergence in Distribution .
30 / 75
It follows that the joint distribution is
31 / 75
It can be shown that the conditional distribution of B(s) given B(t) = a,
when s < t, is Normal with mean st a and variance st (t − s). That is,
s s
P B(s) B(t) = a ∼ N a, (t − s) .
t t
32 / 75
10.7 Hitting Times of Brownian Motion
33 / 75
Let τa be the first time at which the Brownian motion B(t) reaches a.
Suppose that a > 0.
We want to compute the distribution of τa , i.e., P (τa ≤ t), that the process
will hit level a before time t.
34 / 75
If τa ≤ t, the process hits a at some point in [0, t], and by symmetry, it is just
as likely to be above a or below a at time t:
1
P B(t) ≥ a τa ≤ t = .
2
35 / 75
Consider P (B(t) ≥ a) by conditioning on whether or not τa ≤ t (by the law
of total probability ).
Note that if τa > t, the process value cannot be greater than a without
having yet hit a (by continuity of sample paths), so that
P (B(t) ≥ a τa > t) = 0.
36 / 75
Distribution of τa
37 / 75
We obtain,
P (τa ≤ t) = 2 P (B(t) ≥ a)
Z ∞
2 x2
= √ e− 2t dx
2 π t Za
∞ 2 √ √
2 − y2 x
= √ e t dy (by letting y = √ , so dx = t dy)
2 π t a/√t t
Z ∞
2 y2
= √ √
e− 2 dy
2 π a/ t
!
a
= 2 1−Φ √ , a > 0, where Φ(x) ∼ N (0, 1).
t
38 / 75
For a < 0, the distribution of τa is the same as that of τ−a by symmetry. We
obtain Z ∞
2 2
− y2
P (τa ≤ t) = √ e dy, for a < 0.
2 π |a|/√t
39 / 75
10.8 Maximum Value of Brownian Motion
40 / 75
Another variable of interest is maximum value of a Brownian motion
attains in [0, t], max0≤s≤t B(s).
Note that B(0) = 0, so the maximum value of Brownian motion over [0, t]
must be non-negative.
For a > 0,
P max B(s) ≥ a = P ( τa ≤ t ).
0≤s≤t
41 / 75
Therefore, the cumulative distribution function of max0≤s≤t B(s) is
P max B(s) ≤ a = 1 − P max B(s) ≥ a
0≤s≤t 0≤s≤t
= 1 − P (τa ≤ t)
!
a
= 1−2 1−Φ √
t
a
= 2 Φ √ − 1, a > 0, Φ(x) ∼ N (0, 1).
t
42 / 75
10.9 Gambler’s Ruin Problem - Revisited
43 / 75
Let us now consider the probability that a Brownian motion hits a
before −b, where a, b > 0.
We will interpret a Brownian motion as being a limit of a symmetric random
walk.
Recall from the result of gambler’s ruin problem that the probability that
the symmetric random walk goes up to a steps before going down b steps
when each step is equally likely to be either up or down a distance ∆ x is (by
formula Pi = Ni with N = a+b
∆x
and i = ∆bx ) equal to
i b
Pi = = .
N a+b
Hence, letting ∆ x → 0, we see that
b
P (up a before down b) = .
a+b
44 / 75
10.10 Gaussian Processes
45 / 75
Definition
A random vector X1 , X2 , . . . , Xn is said to have a multivariate normal
distribution or a joint normal distribution, if every linear combination of
α 1 X1 + α 2 X2 + · · · + α n X n
46 / 75
Definition
A stochastic process {X(t), t ≥ 0} is called a Gaussian process if
X(t1 ), X(t2 ), . . . , X(tn ), with t1 < t2 < · · · < tn , has a multivariate normal
distribution.
47 / 75
10.11 Brownian motion is a Gaussian Process
48 / 75
A Brownian motion is a unique Gaussian process with zero mean (i.e.
E(B(t)) = 0) and covariance function, for 0 ≤ s ≤ t,
h i
Cov(B(s), B(t)) = E B(s) − E( B(s) ) B(t) − E( B(t) )
h i
= E B(s) − 0 B(t) − 0
h
= E B(s) B(t)
h i
= E B(s) (B(t) − B(s) + B(s))
h i
= E B 2 (s) + E B(s) (B(t) − B(s))
49 / 75
Similarly, if 0 ≤ t ≤ s, we see that Cov(B(s), B(t)) = σ 2 t.
50 / 75
Example
If B(t) is a standard Brownian motion, then X(t) = c B(t/c2 ) for a fixed
c > 0 is a version of standard Brownian motion.
51 / 75
Example
If B(t) is a standard Brownian motion, then X(t) = c B(t/c2 ) for a fixed
c > 0 is a version of standard Brownian motion.
Solution:
Assume that B(t) is a standard Brownian motion. Our goal is to show that
X(t) = c B(t/c2 ) for a fixed c > 0 is a version of standard Brownian motion.
51 / 75
Example
Solution - continued:
1 {X(t), t ≥ 0} is a Gaussian process because {B(t), t ≥ 0} is a standard
Brownian and hence Gaussian .
2 E(X(t)) = E c B(t/c2 ) = c E B(t/c2 ) = 0 since for any t ≥ 0,
B(t) ∼ N (0, t) and B(t/c2 ) ∼ N (0, t/c2 ).
52 / 75
10.12 Brownian Motion with Drift
53 / 75
Definition
A stochastic process {X(t), t ≥ 0} is said to be a Brownian Motion with
drift coefficient µ and variance parameter σ 2 if
1 X(0) = 0
2 {X(t), t ≥ 0} has stationary and independent increments
3 for every t ≥ 0, X(t) is normally distributed with mean µ t and variance
σ 2 t.
54 / 75
10.13 Geometric Brownian Motion
55 / 75
Definition
Let {Y (t), t ≥ 0} be a Brownian motion with drift coefficient µ and variance
parameter σ 2 . Then
X(t) = eY (t)
is called a geometric Brownian motion.
56 / 75
If Y (0) = 0, the X(0) = eY (0) = e0 = 1.
By setting u = 1 we have
1 2
E(X(t)) = E(eY (t) ) = eµ t+ 2 σ t
and
57 / 75
10.14 Application of Geometric Brownian
motion
58 / 75
Samuelson1 (1965) first introduced geometric Brownian motion to represent
the stock price changes.
1
Samuelson, Paul A. (1965) Rational Theory of Warrant Pricing, Industrial
Management Review, 6(2), 13-39.
59 / 75
Given the initial price of the stock price, S(0) = s0 , we can assume that the
price changes follow a geometric Brownian motion,
60 / 75
Example
Suppose that S(t) is a geometric Brownian motion with drift parameter
µ = 0.01 and volatility parameter σ = 0.2. If S(0) = 100, find
1 E(S(10))
2 P (S(10) > 100)
3 P (S(10) < 110)
61 / 75
Example
Solution:
From the question, we have S(t) = 100 e0.01 t+0.2 B(t) .
2
1 E(S(10)) = 100 exp{(0.01 + 0.22 ) 10} = 100 e0.3 .
Solution- continued
➌ Similarly,
−0.0234491
= Φ √
10
= 0.4970418
63 / 75
10.15 Martingales
64 / 75
Definition
A stochastic process {Y (t), t ≥ 0} is said to be a martingale process if, for
s<t
E Y (t) Y (u), 0 ≤ u ≤ s = Y (s). (9)
This definition tells us the best guess of the expected value of Y (t), at
time t, given the history (or information) up to time s is the current
value of Y (s) at time s.
65 / 75
There are other conditions to the definition of a martingale, but it is
beyond the level of this course.
Martingales are known as a type of betting strategy.
There exists a rich and growing theory on martingales. They are very
useful in studying stochastic processes, especially in mathematical
finance.
66 / 75
Example
Let {B(t), t ≥ 0} be a standard Brownian motion. Show that B(t) is a
martingale. That is, for s < t,
67 / 75
Example
Solution:
For s < t,
E(B(t) B(u), 0 ≤ u ≤ s) = E B(t) − B(s) + B(s) B(u), 0 ≤ u ≤ s
= E B(t) − B(s) B(u), 0 ≤ u ≤ s
+ E B(s) B(u), 0 ≤ u ≤ s
= E B(t) − B(s) by independent increments
+ B(s) expectation of a constant is a constant.
= B(s).
68 / 75
Supplementary Material
69 / 75
Supplementary Material - Chapman-Kolmogorov Equation
Write Pijn as the probability of the process going from state i to state j in n
steps. That is,
70 / 75
Supplementary Material - Chapman-Kolmogorov Equations
For any r ≤ n,
∞
X
n−r
Pijn = Pikr Pkj . (10)
k=0
71 / 75
Supplementary Material - Chapman-Kolmogorov Equations - Matrix Form
If we let P(n) denote the matrix of n-step transition probability Pijn , then
That is, the n-step transition matrix is the nth power of the one-step
transition matrix.
72 / 75
Supplementary Material - CLT
X̄n − µ D
√ − → Z,
σ/ n
where Z ∼ N (0, 1). It is common to write
X̄n − µ D
√ − → N (0, 1).
σ/ n
Memorise this! 73 / 75
Supplementary Material - Law of Total Probability
74 / 75
Supplementary Material - Convergence in Distribution
Definition
Let X1 , X2 , . . . be a sequence of random variables. We say that Xn
converges in distribution to X if
A common shorthand is
D
Xn −
→ X.