Solutions To Oksendal
Solutions To Oksendal
Solutions To Oksendal
Solutions Manual
Christopher Kennedy
May 3, 2021
Contents
1 Introduction 2
3 Itô Integrals 9
1
Chapter 1
Introduction
This is a solutions manual for Stochastic Differential Equations by Bernt Øksendal. This is a
working document last updated May 3, 2021. Progress to date:
2
Chapter 2
Hence,
Z Z ∞
X
E(f (X)) = f (x) dPX = S∞
f (x) dPX = f (ak )P(X = ak )
R k=1 {ak } k=1
3
CHAPTER 2. SOME MATHEMATICAL PRELIMINARIES 4
5. Since the measures are σ-finite, f (x, y) = xy is PX ⊗ PY measurable and E(|XY |) < ∞,
apply Fubini-Tonelli and compute
Z
E(XY ) = xy dPXY (x, y)
R2
Z
= xy dPX (x) ⊗ dPY (y)
R2
Z Z
= y x dPX (x) dPY (y)
R R
Z
= E(X) y dPY (y)
R
= E(X)E(Y ).
P∞
6. (Borel-Cantelli) Let {Ak }∞
k=1 ⊆ F and suppose k=1 P(Ak ) < ∞. Then
P(∩∞ ∞
m=1 ∪k=m Ak ) ≤ lim sup P(Ak ) = 0
m→∞ k≥m
by dominated convergence.
7. Let Ω = ni=1 Gi .
F
(a) Note ∅ ∈ G and G is closed under unions by construction. It is also closed under
complements as Ω \ Gi = ∪j6=i Gj ∈ G.
(b) Write a new sequence defined by Fi = Gi \ ∪j≤i Fj and {Fi } will satisfy (a).
(c) Note that {X −1 (x ∈ R)} ⊆ F is disjoint. So, by (a) and (b), F is finite if and only if
all but finitely many X −1 (x ∈ R) are empty.
Since Wt is the sum of two normal distributions, it is also normal and Wt ∼ N (0, t).
13. Compute
ρ2
Z Z ρ Z
1 − |~x|2 2 2π 2
− r2t
2t ρ2
P0 (Bt ∈ Dρ ) = e 2t d ~x = re dr = e−u du = 1 − e− 2t .
|x|<ρ 2πt 2πt 0 0
14. Compute
Z Z
Ex χK (Bt ) dt = P(Bt ∈ K) dt
[0,∞] [0,∞]
Z Z
1 − 2t
|~
x−~y |2
n
= n/2
e d ~x dt
[0,∞] K (2πt)
Z
1 |~
x−~y |2
− 2t
≤
n/2
e
µ(K) dt
[0,∞] (2πt)
∞
=0
15. Note that U U T = I, whence | det U | = 1 and the probability measures are identical by
change of variables. It follows that both are Brownian.
P0 (Wt ∈ U ) = P0 (Bc2 t ∈ cU )
Z
= p(c2 t, 0, y) dy
ZcU
1
= p(t, 0, y/c) dy
cU c
Z
1
= p(t, 0, y 0 )(cdy 0 )
U c
= P0 (Bt ∈ U ),
(2)
So hB, Bit (w) = t.
(b) Note that the Brownian motion has positive quadratic variation t on [0, t]. So
(2)
(1) hB, Bit (w)
hB, Bit (w) ≥ lim = ∞.
k∆Bk k→0+ k∆Bk k
Chapter 3
Itô Integrals
1. Compute
dnte
Z t t
−1
X jt
s dBs = lim (B (j+1)t − B jt )
0 n→∞
j=0
n n n
dnte
t
−1
dnte t X t
= lim B dnte − lim B jt + lim (B0 − B dnte )
n→∞ n n n→∞ n n n→∞ n n
j=0
Z t
= tBt − Bs ds.
0
2. Compute
dnte
Z t t
−1
X
Bs2 dBs = lim B 2jt (B (j+1)t − B jt )
0 n→∞ n n n
j=0
dnte
t
−1
X 1 1 3 1
3 2 3
= lim B (j+1)t − B j − B jt (B (j+1)t − B j ) − (B (j+1)t − B j )
n→∞
j=0
3 n 3 n n n n 3 n n
dnte
t
−1
1 X t
= Bt3 − lim B jt + O(t2 /n)
3 n→∞
j=0
n n
Z t
1 3
= Bt − Bs ds.
3 0
(X)
3. Let {Nt } be some filtration and let {Ht } be the filtration of process Xt .
(a) Compute
E(Xt | Hs(X) ) = E E(Xt | Ns ) | Hs(X) = E(Hs | Hs(X) ) = Hs .
9
CHAPTER 3. ITÔ INTEGRALS 10
(b) Compute
(X)
E(Xt ) = E(E(Xt | H0 )) = E(X0 ).
4. Compute
E(Bt + 4t | Fs ) = Bs + 4t 6= Bs + 4s
E(Bt2 | Fs ) = E((Bt − Bs )2 + 2Bs (Bt − Bs ) + Bs2 | Fs ) = Bs2 + t − s 6= Bs2
Z t Z s Z t Z s
2 2 2
E t Bt − 2 uBu du | Fs = t Bs − 2 uBu du − 2 uBs du = s Bs − 2 uBu du
0 0 s 0
(1) (2) (1) (2)
E(Bt Bt | Fs ) = E(Bt | Fs )E(Bt | Fs ) = Bs(1) Bs(2) ,
R that {0 ≤ u1 · · · ≤ un } is Borel
(a) Note measurable and χ0≤u1 ···≤un is Ft -adapted. Finally
T
E 0 f (t1 , . . . tn , ω) dt1 . . . dtn ≤ T n < ∞.
2
(c) Deduce that d(Bt3 − 3tBt ) = 3(Bt2 − t) dBt and so Yt := Bt3 − 3tBt is a martingale.
8. There exists continuous martingale Mt iff there exists Y ∈ L1 such that Mt = E(Y | Ft ).
(a) Verify that E(|E(Y | Ft )|) ≤ E(E(|Y | | Ft ) = E(|Y |) < ∞ and
E(Mt | Fs ) = E(E(Y | Ft ) | Fs ) = E(Y | Fs ) = Ms .
(b) If Mt is a continuous martingale such that supt>0 E(|X|p ) < ∞ for p ∈ (1, ∞), then
∃M such that kMt − M kL1 → 0 as t → ∞. So let Y = M and
Z Z
lim |Ms − E(M | Fs )| dP = lim |E(Ms − M | Fs )| dP
s→∞ Ω s→∞ Ω
s s
Z
≤ lim E(|Ms − M | | Fs ) dP
s→∞ Ω
s
Z
= lim |Ms − M | dP
s→∞ Ωs
= 0.
9. Compute
bntc
Z T t
−1
X 1
Bt ◦ dBt = lim (B jt + B (j+1)t )(B (j+1)t − B jt )
0 n→∞
j=0
2 n n n n
bntc bntc
t
−1 t
−1
X X 1
= lim B jt (B (j+1)t − B jt ) + lim (B (j+1)t − B jt )2
n→∞
j=0
n n n n→∞
j=0
2 n n
1 t t
= Bt2 − +
2 2 2
1 2
= Bt .
2
10. If f (t, ω) varies smoothly in t, then the Itô and Stratonovich integrals coincide. Compute
Z T Z T
1
f (t, ω) ◦ dBt = f (t, ω) dBt + hf (t, ω), Bt i(2)
0 0 2
and
E(hf (t, ω), Bt i(2) )2 ≤ E(hBt , Bt i(2) E(hf (t, ω), f (t, ω)i(2)
T
≤ T lim + sup (K|∆tk |1+ε )
k∆tk k→0 |∆tk | |∆tk |
= KT 2 lim k∆tk kε
k∆tk k→0+
= 0.
CHAPTER 3. ITÔ INTEGRALS 12
(N )
11. Define white noise Wt = max{−N, min{Wt , N }}. Since Wt and Ws are independent
(N ) (N )
and identically distributed, it follows that Wt and Ws are as well. If Wt is continuous,
(N )
then since |Wt | ≤ N and by bounded convergence
(N ) 2 (N )
lim 2E(Wt ) = lim E(|Wt − Ws(N ) |2 ) = 0.
t→s t→s
a.s.
But then Wt = 0, which is a contradiction.
α2
dXt = (γ + )Xt dt + αXt dBt .
2
Since (t2 + cos(Xt )) ◦ dBt = − sin(X
2
t)
(t2 + cos(Xt )) dt + (t2 + cos(Xt ) dBt ,
sin(Xt )
dXt = (cos(Xt ) − t2 ) dt + (t2 + cos(Xt )) dBt .
2
α2
(ii) Since αXt dBt = αXt ◦ dBt − 2
Xt dt,
α2
dXt = (r − )Xt dt + αXt ◦ dBt .
2
Since Xt2 dBt = Xt2 ◦ dBt − Xt3 dt,
= 0.
CHAPTER 3. ITÔ INTEGRALS 13
14. Show that h(ω) is Ft measurable if and only if it is the pointwise limit of a sum-product of
bounded continuous functions g(Btj ).
(a) Assume that h is bounded since {hn (ω) := h(ω)1{|h(ω)|<n} } converges pointwise to h.
(b) Let Hn be the σ-algebra generated by B(tj ) for tj = 2jn ≤ t. Then Ft = σ (∪n Hn )
and so by Corollary (C.9), h = E[h|Fn ] = lim E[h|Hn ].
n→∞
(c) By Doob-Dynkin, E[h|Hn ](ω) = g Bt1 , . . . B(tb2n tc ) . Since C(Rk ) is dense in
L1 (Rk ) and by Stone-Weierstrass P (Rk ) is dense in C(Rk ), a limiting sequence must
exist.
RT RT
15. Suppose C + S f (t, ω) dBt (ω) = D + S g(t, ω) dBt (ω). Then we have that
Z T Z T
C − D = E[C − D] = E g(t, ω) dBt (ω) − f (t, ω) dBt (ω) = 0 =⇒ C = D,
S S
1. Compute
3. Let Xt and Yt be Itô processes. Then, letting f (t, x, y) = xy and by Itô’s formula
14
CHAPTER 4. THE ITÔ FORMULA 15
< ∞.
(a) Calculate
1
dXt = cect+αBt dt + αect+αBt dBt + α2 ect+αBt d[B, B]t
2
α2
= Xt (c + ) dt + α dBt .
2
(b) Calculate
n n
!
X (j) 1X
dYt = Yt c dt + αj dBt + αi αj d[B (i) , B (j) ]t
j=1
2 i,j=1
n n
!
1X 2 X (j)
= Yt (c + α ) dt + αj dBt .
2 j=1 i j=1
and so
Z t Z t Z t
1
f (Bt ) − f (B0 ) = d(f (Bs )) = h∇f (Bs ), dBs i + ∆f (Bs ) ds.
0 0 2 0
(b) Assume that g is of class C 1 everywhere, as well as C 2 and uniformly bounded outside
of finitely many points with |g 00 (z)| ≤ M for z ∈/ {z1 , . . . zk }. Then the set of functions
{f } of class C 2 uniformly bounded with |f 00 (z)| ≤ M are C k -dense. So we can extract
a sequence {fk } such that fk ⇒ g, fk0 ⇒ g 0 as well as fk00 → g 00 and |fk00 | ≤ M on
R \ {z1 , . . . zk }. So
Z t
1 t 00
Z
0 0 00
lim (fk − g)(Bt ) + (fk − g)(0) +
(fk − g ) dBs + (fk − g ) ds
k→∞ 0 2 0
1 t 00
Z
0 0
≤ lim |(fk − g)(Bt )| + |(fk − g)(0)| + tkfk − g k∞ + |f − g 00 | ds
k→∞ 2 0 k
= 0,
9. Clearly
Z t Z t∧τn
∂gn ∂g
v (s, Xs )χs≤τn dBs = v (s, Xs ) dBs
0 ∂x 0 ∂x
and the result follows by Itô’s lemma where dXt = u dt + v dBt . Since E(|Xt |) < ∞, it
follows that lim P(τn > t) = lim P(Xt < n) = 1 and so the identity holds almost surely.
n→∞ n→∞
10. (Tanaka) In this problem, Tanaka’s formula for Brownian motion is derived.
Then
u(t, ω) = E(u(t, ω) | Ft ) = lim− E(u(t, ω) | Fs ) = 0.
s→t
13. Let dXt = u(t, ω) dt + dBt where u(t, ω) ∈ V([0, T ]). Then Yt = Xt Mt is a martingale,
where
Z t
1 t 2
Z
Mt = exp − u(r, ω) dBr − u (r, ω) dr
0 2 0
√ qR t
since E(|Mt |) < ∞ (see question 4b), E(|Xt |) ≤ t 0
u2 (r, ω) dr + 1 < ∞ and
14. In this problem, the martingale representation of stochastic processes is explicitly shown.
Deduce that
Z T
Bt −t/2
e BT
=e T /2
1+ e dBt =⇒ f (t, ω) = eBt +(T −t)/2 .
0
Deduce that
Z T
−T /2
sin(BT ) = e et/2 cos(Bt ) dBt =⇒ f (t, ω) = e−(T −t)/2 cos(Bt ).
0
1. Compute
B
(a) dXt = d(eBt ) = eBt dBt + 12 t d[B, B]t = 12 Xt dt + Xt dBt
Bt 1 Bt 1 1
(b) dXt = d 1+t = 1+t dBt − (1+t) 2 dt = 1+t dBt − 1+t Xt dt
(1) (2)
2. Let Xt = a cos(Bt ) and Xt = b sin(Bt ). Then
is
Z t
(2) (2) (1) (2)
Xt = X0 + X0 Bt + (s + Bs(1) ) dBs(2) .
0
(c) Using integrating factors, solve dXt = −Xt dt + e−t dBt for
Z t
t
e Xt − X0 = dBs
0
Xt = e−t (X0 + Bt ).
6. Suppose Yt is given by
dYt = r dt + αYt dBt .
Using integrating factors, solve for
α2
−αBt −αBt
d(e Yt ) = e Yt r− dt
2
and
2
Z t
α2
−αBt + α2 t
e Yt − Y0 = re−αBs + 2
s
ds.
0
Deduce that
2
Z t
α2
αBt − α2 t
Yt = e Y0 + r eα(Bt −Bs )− 2
(t−s)
ds.
0
9. Let dXt = ln(1 + Xt2 ) dt + χ{Xt >0} Xt dBt . It suffices to check that
2
|b(t, x)| + |σ(t, x)| = ln(1 + x2 ) + χ{x>0} |x| ≤ (|x| + 1) + |x| ≤ 2(|x| + 1),
e
2 2
E(|X0 | ) = α < ∞, and
|b(t, x) − b(t, y)| + |σ(t, x) − σ(t, y)| ≤ | ln(x2 ) − ln(y 2 )| + |x − y| ≤ 3|x − y|.
Hence, by Theorem 5.2.1, there is a unique strong solution to the stochastic differential
equation.
10. Calculate
Z t Z t 2
E(Xt2 ) =E Z+ b(s, Xs ) ds + σ(s, Xs ) dBs
0 0
Z t 2 Z t 2 !
≤ 3 E(Z 2 ) + E b(s, Xs ) ds +E σ(s, Xs ) dBs
0 0
Z t Z t
2 2 2
≤ 3 E(Z ) + T E b(s, Xs ) ds + E σ(s, Xs ) ds
0 0
Z t
2 2 2
≤ 3E(Z ) + 6C T + E(|Xs | ) ds (T + 1)
0
Z t
2 2 2
= (3E(Z ) + 6C T (T + 1)) + 6C (T + 1) E(|Xs |2 ) ds.
0
and apply Gronwall to derive the result.
CHAPTER 5. STOCHASTIC DIFFERENTIAL EQUATIONS 24
(a) Rewrite
dyt 0 1 yt 0 0 yt
= dt + dBt .
dy˙t −1 0 y˙t −ε 0 y˙t
Rt Rt
(b) Check that, if y(t) = y(0) + y 0 (0)t + 0
(r − t)y(r) dr + 0
ε(r − t)y(r) dBr , then
Z t Z t Z t
0 0 0
y (t) = y (0) − y(r) dr − εy(r) dBr = y (0) − y(r)(1 + εWr ) dr
0 0 0
13. Let x00t + a0 x0t + w2 xt = (T0 − α0 x0t )ηWt where Wt is 1-dimensional white noise. Then
dxt 0 1 xt 0 0 xt 0
= 2 dt + dBt + dBt
dẋt −w −a0 ẋt 0 −α0 η ẋt T0 η
α0 ηe−λ(t−s) 0
A(t−s) sin(ξ(t − s)) xs −α0 ηys gt−s
e KXs = − =
ξ 0 ξ cos(ξ(t − s)) − λ sin(ξ(t − s)) x˙s −α0 ηys ht−s
and
T0 ηe−λ(t−s)
A(t−s) sin(ξ(t − s)) ηT0 gt−s
e M= = .
ξ ξ cos(ξ(t − s)) − λ sin(ξ(t − s)) ηT0 ht−s
It follows that
Z t
xt = η (T0 − α0 ys )gt−s dBs
0
and
Z t
yt = η (T0 − α0 ys )ht−s dBs .
0
(1) (2)
14. Letting Zt = F (Bt ), where Bt = Bt + iBt , calculate
(1) (2)
dZt = Fx (Bt ) dBt + Fy (Bt ) dBt
1
+ Fxx (Bt ) d[B (1) , B (1) ]t + Fxy (Bt ) d[B (1) , B (2) ]t + Fyy (Bt ) d[B (2) , B (2) ]t
2
(1) (2) 1
= (ux + ivx ) dBt + (uy + ivy ) dBt + (uxx + ivxx + uyy + ivyy ) dt
2
1
= hF 0 (Bt ), dBt i + (vxy − iuxy + uyy + ivyy ) dt
2
1
= hF 0 (Bt ), dBt i + (−uyy − ivyy + uyy + ivyy ) dt
2
= hF 0 (Bt ), dBt i.
CHAPTER 5. STOCHASTIC DIFFERENTIAL EQUATIONS 26
and calculate
2
dZt = −βZt dBt + re(rk−β )t dt
Z t
−βBt −1 (rk−β 2 )s+βBs
=⇒ Zt = e x +r e ds .
0
Conclude that
2 )t 2
e(rk−β e(rk−β )t+βBt
Xt = = Rt .
Zt x−1 + r 0 e(rk−β 2 )s+βBs ds
which implies
s Z t
Yt = Y02 +2 e−2αBt +α2 s ds
0
and
s
2
Z t
αBt − α2 t
Xt = e x2 +2 e−2αBt +α2 s ds.
0
and
Z t 1
1−γ
2 2
αBt − α2 t 1−γ −(1−γ)Bs +(1−γ) α2 s
Xt = e x + (1 − γ) e ds .
0
Rt Rt
17. Let v ≥ 0 satisfy v(t) ≤ C + A 0
v(s) ds and consider quantity w(t) = 0
v(s) ds. Then
Z t
0
w (t) = v(t) ≤ C + A v(s) ds = C + Aw(t).
0
and
Z t
−At C
w(t)e ≤ Ce−As ds = (1 − e−At )
0 A
C At
=⇒ w(t) ≤ (e − 1).
A
Deduce that
28
Chapter 7
29
Chapter 8
30
Chapter 9
31
Chapter 10
32
Chapter 11
33
Chapter 12
34