OCT21
FRI2022

Stochastic Differential Equations & Geometric Brownian Motion

Tracking how ΔXtN(0,Δt)\Delta X_t \sim \mathcal{N}(0,\Delta t) becomes dSt=μStdt+σStdWtdS_t = \mu S_t dt + \sigma S_t dW_t.
stochastic differential equationsgeometric brownian motionjump diffusionsquare root diffusion

Attention Conservation Notice: Read the full warning

This is the 2nd post in a series of five (first post here) where we build up to the Black-Scholes-Merton model for pricing European put and call options.

  1. Brownian Motion as a Symmetric Random Walk Limit
  2. Stochastic Differential Equations & Geometric Brownian Motion (you are here)
  3. Ito's Lemma
  4. The Black-Scholes Partial Differential Equation
  5. (Coming soon) - tying the pieces together for the pricing formula

This track closely follows Stephen Blyth's An Introduction to Quantitative Finance with a few Python experiments sprinkled in. Expect:

  • not-very-practical advice if you simply want to price an option,
  • a heavy dose of math if you're allergic,
  • and perhaps gaps if you're already deep into stochastic calculus.

Preamble

Last time we built Brownian motion by taking the limit of a symmetric random walk as the step size shrinks to zero. With that in hand, we can start to sketch the idea of a stochastic differential equation (SDE).

Difference equations

For Δt>0\Delta t > 0 we can write

St+ΔtSt=μΔt+σ(Wt+ΔtWt)=μΔt+σΔtW,WN(0,1).S_{t + \Delta t} - S_t = \mu \Delta t + \sigma \bigl(W_{t + \Delta t} - W_t\bigr) = \mu \Delta t + \sigma \sqrt{\Delta t}\, W,\qquad W \sim \mathcal{N}(0,1).

The first term μ\mu is a deterministic drift; the second scales a Brownian increment. Holding one term fixed while varying the other gives intuition for their roles:

00.510.60.811.21.400.5100.511.5200.510.811.21.400.5100.511.52
mu = 0, sigma = 0.2How drift and volatility shape stochastic pathsTimeTimeTimeTimeValueValueValueValuemu = 0, sigma = 0.2mu = 0, sigma = 0.6mu = 0.1, sigma = 0.2mu = 0.1, sigma = 0.6

μ\mu sets the slope-add or subtract a little every timestep and the path tilts accordingly. σ\sigma dials how violently Brownian motion kicks the process around; larger σ\sigma makes the path spikier unless a strong drift overpowers it.

Stochastic differential equations

Letting Δt0\Delta t \to 0 gives the differential form

dSt=μdt+σdWt,dS_t = \mu\, dt + \sigma\, dW_t,

or more generally

dSt=μ(St,t)dt+σ(St,t)dWt.dS_t = \mu(S_t, t)\, dt + \sigma(S_t, t)\, dW_t.

The discrete intuition (ΔtW\sqrt{\Delta t} W) still helps, but dWtdW_t is subtle-it represents an infinitesimal Brownian increment.

Why SDEs show up everywhere

Stochastic differential equations take an ordinary differential equation and inject randomness. Use a normal draw, you get diffusions; use something else and you can model jumps, regime switches, and more.

  • Short rates. The Cox-Ingersoll-Ross model adds a square-root diffusion term to produce mean-reverting short rates.

    01234524681012
    Mean pathCIR square-root diffusion (rates mean-revert)YearsRate (bps)Kappa=1.1, Theta=6%, Sigma=15%, Start=4%
  • Jump diffusion. Swap the Gaussian draw for a Poisson jump and you reach the Merton jump-diffusion model, capturing sudden moves from earnings releases or breaking news.

    01230.60.811.21.41.601230.60.811.21.41.61.82
    JumpsJump diffusion: contrasting low vs high jump intensityYearsYearsPrice (normalized)Price (normalized)Low jump intensity (lambda = 0.2)High jump intensity (lambda = 1.2)

In both cases, parameter choices (e.g., the Poisson rate λ\lambda) change the texture of the path dramatically.

Choosing μ(St,t)\mu(S_t, t) and σ(St,t)\sigma(S_t, t) for equity prices

A simple, finance-flavored specification is:

  • Drift is proportional to price: μ(St,t)=μSt\mu(S_t, t) = \mu S_t. Without randomness this solves to St=S0eμtS_t = S_0 e^{\mu t}.
  • Volatility scales with price: for small Δt\Delta t, assume Var(ΔSt/St)σ2Δt\operatorname{Var}(\Delta S_t / S_t) \approx \sigma^2 \Delta t, so Var(ΔSt)σ2St2Δt\operatorname{Var}(\Delta S_t) \approx \sigma^2 S_t^2 \Delta t.

That gives the familiar geometric Brownian motion (GBM):

dSt=μStdt+σStdWt.dS_t = \mu S_t\, dt + \sigma S_t\, dW_t.
00.511.522.530.511.522.5
80% bandMedian pathGeometric Brownian motion sample pathsYearsPrice (normalized)

Why GBM is beloved (and where it falls short)

GBM has the right vibe for equities-multiplicative growth, log-normal distributions, paths that "look" like prices. But it's not quite a duck:

  • Volatility is fixed. Real markets exhibit volatility smiles and clustering. Allowing σ\sigma to depend on (St,t)(S_t, t) leads to local volatility models; letting σ\sigma follow its own SDE yields stochastic volatility models.
  • No jumps. Earnings surprises, takeovers, and crises produce discontinuities. Jump-diffusion blends or Levy processes push past GBM's smoothness.

Next time

We'll need Ito's Lemma to manipulate functions of diffusions like GBM. That's next on the path toward Black-Scholes.


Extra notes

  • Most surfaces look cooler in 3D. Here's a Brownian surface for eye candy:

    Brownian surface (toy rendering)
  • Why stochastic volatility matters. Black-Scholes assumes constant volatility; letting volatility evolve stochastically repairs that simplification and better fits observed option prices.