Mehmet Emre Tiryaki
Back to Notes

Itô Calculus

Setting

L2(Ω)L^2(\Omega) is a Hilbert space with ⟨X,Y⟩=E[XY]\langle X, Y \rangle = E[XY]. A process f(t,ω)f(t, \omega) on [0,∞)×Ω[0, \infty) \times \Omega is progressively measurable if for every t≥0t \ge 0 the map (s,ω)↦f(s,ω)(s, \omega) \mapsto f(s, \omega) on [0,t]×Ω[0, t] \times \Omega is B([0,t])⊗Ft\mathcal{B}([0,t]) \otimes \mathcal{F}_t-measurable. Adapted processes with left- or right-continuous paths are progressively measurable. Lprog2([0,T]×Ω)L^2_{prog}([0,T] \times \Omega) is the space of progressively measurable ff with E∫0Tf2 dt<∞E\int_0^T f^2\, dt < \infty.

The Itô Integral

For a simple adapted process, constant Δ(tj)\Delta(t_j) on each [tj,tj+1)[t_j, t_{j+1}), and tk≤t<tk+1t_k \le t < t_{k+1},

I(t)=∫0tΔ(u) dW(u)=∑j=0k−1Δ(tj)[W(tj+1)−W(tj)]+Δ(tk)[W(t)−W(tk)]I(t) = \int_0^t \Delta(u)\, dW(u) = \sum_{j=0}^{k-1} \Delta(t_j)\big[W(t_{j+1}) - W(t_j)\big] + \Delta(t_k)\big[W(t) - W(t_k)\big]

For general Δ∈Lprog2\Delta \in L^2_{prog}, choose simple Δn\Delta_n with lim⁡n→∞E∫0T∣Δn(t)−Δ(t)∣2 dt=0\lim_{n \to \infty} E\int_0^T |\Delta_n(t) - \Delta(t)|^2\, dt = 0 (simple, then bounded continuous, then bounded, then general), and define ∫0TΔ dW\int_0^T \Delta\, dW as the L2L^2 limit of ∫0TΔn dW\int_0^T \Delta_n\, dW.

Properties

  • I(t)I(t) has continuous paths and is Ft\mathcal{F}_t-measurable.
  • It is linear in the integrand.
  • It is a martingale, E[I(t)∣Fs]=I(s)E[I(t)|\mathcal{F}_s] = I(s), so E[I(t)]=0E[I(t)] = 0.
  • Itô isometry E[I(t)2]=E∫0tΔ2(u) du=Var(I(t))E[I(t)^2] = E\int_0^t \Delta^2(u)\, du = \mathrm{Var}(I(t)), i.e. the integral is an isometry from Lprog2L^2_{prog} into L2(Ω)L^2(\Omega).
  • [I,I](t)=∫0tΔ2(u) du[I, I](t) = \int_0^t \Delta^2(u)\, du, i.e. dI dI=Δ2 dtdI\, dI = \Delta^2\, dt, and [I,J](t)=∫0tΔΓ du[I, J](t) = \int_0^t \Delta\Gamma\, du for J=∫Γ dWJ = \int \Gamma\, dW.
  • If Δ\Delta is deterministic, ∫0tΔ dW∼N(0,∫0tΔ2(s) ds)\int_0^t \Delta\, dW \sim \mathcal{N}\big(0, \int_0^t \Delta^2(s)\, ds\big).

Itô Processes

X(t)=X(0)+∫0tΔ(u) dW(u)+∫0tΘ(u) du,dX=Δ dW⏟diffusion+Θ dt⏟driftX(t) = X(0) + \int_0^t \Delta(u)\, dW(u) + \int_0^t \Theta(u)\, du, \qquad dX = \underbrace{\Delta\, dW}_{\text{diffusion}} + \underbrace{\Theta\, dt}_{\text{drift}}

with Δ,Θ\Delta, \Theta adapted. Its quadratic variation is [X,X](t)=∫0tΔ2(u) du[X, X](t) = \int_0^t \Delta^2(u)\, du, i.e. dX dX=Δ2 dtdX\, dX = \Delta^2\, dt, since dW dW=dtdW\, dW = dt and dW dt=dt dt=0dW\, dt = dt\, dt = 0.

Itô's Lemma

If f(t,x)f(t, x) has continuous ftf_t, fxf_x, fxxf_{xx}, then

f(T,X(T))=f(0,X(0))+∫0Tft dt+∫0Tfx dX+12∫0Tfxx d[X,X](t)f(T, X(T)) = f(0, X(0)) + \int_0^T f_t\, dt + \int_0^T f_x\, dX + \frac{1}{2}\int_0^T f_{xx}\, d[X, X](t)

or in differential form

df=ft dt+fx dX+12fxx dX dX=(ft+Θfx+12Δ2fxx) dt+Δfx dWdf = f_t\, dt + f_x\, dX + \frac{1}{2} f_{xx}\, dX\, dX = \Big(f_t + \Theta f_x + \frac{1}{2}\Delta^2 f_{xx}\Big)\, dt + \Delta f_x\, dW

For X=WX = W this is df=ft dt+fx dW+12fxx dtdf = f_t\, dt + f_x\, dW + \frac{1}{2} f_{xx}\, dt.

Multidimensional df=ft dt+∇f⋅dX+12 dXTHf dXdf = f_t\, dt + \nabla f \cdot dX + \frac{1}{2}\, dX^T H_f\, dX, with HfH_f the Hessian in xx.

Product rule d(XY)=X dY+Y dX+dX dYd(XY) = X\, dY + Y\, dX + dX\, dY.

Examples

∫0TW dW\int_0^T W\, dW Approximate WW by the simple process Δn(t)=W(jT/n)\Delta_n(t) = W(jT/n) on [jT/n,(j+1)T/n)[jT/n, (j+1)T/n). Writing Wj=W(jT/n)W_j = W(jT/n),

∑j=0n−1Wj(Wj+1−Wj)=12Wn2−12∑j=0n−1(Wj+1−Wj)2  ⟶  12W(T)2−12T\sum_{j=0}^{n-1} W_j (W_{j+1} - W_j) = \frac{1}{2}W_n^2 - \frac{1}{2}\sum_{j=0}^{n-1} (W_{j+1} - W_j)^2 \;\longrightarrow\; \frac{1}{2}W(T)^2 - \frac{1}{2}T

since the sum of squared increments converges to [W,W](T)=T[W, W](T) = T.

Others

∫0Tt dW=TW(T)−∫0TW(t) dt,∫0TW2 dW=13W(T)3−∫0TW(t) dt\int_0^T t\, dW = TW(T) - \int_0^T W(t)\, dt, \qquad \int_0^T W^2\, dW = \frac{1}{3}W(T)^3 - \int_0^T W(t)\, dt

from the product rule and from Itô's lemma with f(x)=x3/3f(x) = x^3/3.

Geometric Brownian motion For dS=μS dt+σS dWdS = \mu S\, dt + \sigma S\, dW and f(x)=ln⁡xf(x) = \ln x,

d(log⁡S)=dSS−12S2dS dS=μ dt+σ dW−12σ2 dt=(μ−12σ2)dt+σ dWd(\log S) = \frac{dS}{S} - \frac{1}{2S^2} dS\, dS = \mu\, dt + \sigma\, dW - \frac{1}{2}\sigma^2\, dt = \Big(\mu - \frac{1}{2}\sigma^2\Big) dt + \sigma\, dW

Ornstein-Uhlenbeck For dY=−aY dt+σ dWdY = -aY\, dt + \sigma\, dW, let Z=YeatZ = Ye^{at}. Then dZ=σeat dWdZ = \sigma e^{at}\, dW, so

Y(t)=e−atY(0)+σ∫0te−a(t−s) dW(s)Y(t) = e^{-at}Y(0) + \sigma\int_0^t e^{-a(t-s)}\, dW(s)