Itō's lemma

From Wikipedia Mirror

In mathematics, Itō's lemma is used in Itō stochastic calculus to find the differential of a function of a particular type of stochastic process. It is the stochastic calculus counterpart of the chain rule in ordinary calculus and is best memorized using the Taylor series expansion and retaining the second order term related to the stochastic component change. The lemma is widely employed in mathematical finance and its best known application is in the derivation of the Black-Scholes equation used to value options. Ito's Lemma is also referred to currently as the Itō–Doeblin Theorem in recognition of the recently discovered work of Wolfgang Doeblin.^[1]

Contents 1 Itō drift-diffusion processes 2 Continuous semimartingales 3 Poisson jump processes 4 Non-continuous semimartingales 5 Informal derivation 6 Examples 6.1 Geometric Brownian motion 6.2 The Doléans exponential 6.3 Black–Scholes formula 7 See also 8 References 9 Notes 10 External links

Itō drift-diffusion processes

In its simplest form, Itō's lemma states the following: for an Itō drift-diffusion process

<math> dX_t= \sigma_t\,dB_t + \mu_t\,dt</math>

where <math>dB_t</math> is a differential of Brownian motion, and for any twice continuously differentiable (having continuous first- and second-order partial derivatives) function ƒ(t, x) of two real variables t and x, letting

<math>

f'(t,x)=\frac{\partial f}{\partial x}(t,x),\quad f(t,x)=\frac{\partial^2f}{\partial x^2}(t,x),\quad \dot{f}(t,X_t)=\frac{\partial f}{\partial t}(t,x) </math>

then

<math>

\begin{align} df(t,X_t)&=\dot{f}(t,X_t)\,dt+f'(t,X_t)\,dX_t+\frac{1}{2}f(t,X_t)\sigma^2_t\,dt,\\ &=\dot{f}(t,X_t)\,dt+f'(t,X_t)(\mu_t\,dt + \sigma_t\,dB_t)+\frac{1}{2}f(t,X_t)\sigma^2_t\,dt,\quad\mbox{rearranging terms}\\ &=\left(\dot{f}(t,X_t)+\mu_tf'(t,X_t)+\frac{\sigma_t^2}{2}f(t,X_t)\right)dt+f'(t,X_t)\sigma_t\,dB_t \end{align} </math>

In multi-dimension,

<math>

df\left(t,X_{t}\right)=\dot{f}_{t}\left(t,X_{t}\right)dt+\nabla_{X_{t}}^{T}f\cdot dX_{t}+\frac{1}{2}dX^{T}_{t}\cdot\nabla_{X_{t}}^{2}f\cdot dX_{t} </math>

where <math>X_{t}=\left(X_{t,1},X_{t,2},\cdots,X_{t,n}\right)^{T}</math> is a vector of Itō process, <math>\dot{f}_{t}\left(t,X\right)</math> is the partial differential w.r.t. t, <math>\nabla_{X}^{T}f</math> is the gradient of ƒ w.r.t. X, and <math>\nabla_{X}^{2}f</math> is the Hessian matrix of ƒ w.r.t. X.

Continuous semimartingales

More generally, the above formula also holds for any continuous d-dimensional semimartingale X=(X¹,X²,…,X^d), and twice continuously differentiable and real valued function f on R^d. And some people prefer to presenting the formula in another form with cross variation shown explicitly as follows, f(X) is a semimartingale satisfying

<math>df(X_t) = \sum_{i=1}^d f_{,i}(X_t)\,dX^i_t + \frac{1}{2}\sum_{i,j=1}^df_{,ij}(X_t)\,d[X^i,X^j]_t.</math>

In this expression, the term f_,i represents the partial derivative of f(x) with respect to xⁱ, and [Xⁱ,X^j ] is the quadratic covariation process of Xⁱ and X^j.

Poisson jump processes

We may also define functions on discontinuous stochastic processes.

Let h be the jump intensity. The Poisson process model for jumps is that the probability of one jump in the interval <math>[t,t+{\Delta}t]</math> is <math>h {\Delta}t</math> plus higher order terms. h could be a constant, a deterministic function of time or a stochastic process. The survival probability <math>p_s(t)</math> is the probability that no jump has occurred in the interval <math>[0,t]</math>. The change in the survival probability is

<math>d p_s(t) = -p_s(t) h(t) \, dt.</math>

So

<math>p_s(t) = \exp \left(-\int_0^t h(u) \, du \right).</math>

Let <math>S(t)</math> be a discontinuous stochastic process. Write <math>S(t^-)</math> for the value of S as we approach t from the left. Write <math>d_j S(t)</math> for the non-infinitesimal change in <math>S(t)</math> as a result of a jump. Then

<math>d_j S(t)=\lim_{{\Delta}t \to 0}(S(t+{\Delta}t)-S(t-))</math>

Let z be the magnitude of the jump and let <math>\eta(S(t^-),z)</math> be the distribution of z. The expected magnitude of the jump is

<math>E[d_j S(t)]=h(S(t^-)) dt \int_z z \eta(S(t^-),z) dz \, .</math>

Define <math>d J_S(t)</math>, a compensated process and martingale, as

<math>d J_S(t)=d_j S(t)-E[d_j S(t)]=S(t)-S(t^-)-(h(S(t^-)) \int_z z \eta(S(t^-),z) dz) dt \,.</math>

Then

<math>d_j S(t)=E[d_j S(t)]+d J_S(t)=h(S(t^-)) \int_z z \eta(S(t^-),z) dz + d J_S(t) \, .</math>

Consider a function <math>g(S(t),t)</math> of jump process <math>d S(t)</math>. If <math>S(t)</math> jumps by <math>{\Delta}s</math> then <math>g(t)</math> jumps by <math>{\Delta}g</math>. <math>{\Delta}g</math> is drawn from distribution <math>\eta_g()</math> which may depend on <math>g(t^-)</math>, dg and <math>S(t^-)</math>. Itō's Lemma for <math>g(S(t),t)</math> is then

<math> d g(t) = g(t)-g(t^-) \,,</math>

<math>=h(t) dt \int_{\Delta}g {\Delta}g \eta_g(.) d{\Delta}g + d J_g(t) \,,</math>

<math>=\frac{\partial g}{\partial t}+\mu \frac{\partial g}{\partial S}+\frac{1}{2} \sigma^2 \frac{\partial g}{\partial S^2}+h(t)\int_{{\Delta}g} ({\Delta}g \eta_g(.) d{\Delta}g)dt+\frac{\partial g}{\partial S} \sigma d W(t)+ d J_g(t) \,.</math>

Itō's lemma for a process which is the sum of a drift-diffusion process and a jump process is just the sum of the Itō's lemma for the individual parts.

Non-continuous semimartingales

Itō's lemma can also be applied to general d-dimensional semimartingales, which need not be continuous. In general, a semimartingale is a càdlàg process, and an additional term needs to be added to the formula to ensure that the jumps of the process are correctly given by Itō's lemma. For any cadlag process Y_t, the left limit in t is denoted by Y_t-, which is a left-continuous process. The jumps are written as ΔY_t = Y_t - Y_t-. Then, Itō's lemma states that if X = (X¹,X²,…,X^d) is a d-dimensional semimartingale and f is a twice continuously differentiable real valued function on R^d then f(X) is a semimartingale, and

<math>

\begin{align} f(X_t)= & f(X_0)+\sum_{i=1}^d\int_0^t f_{,i}(X_{s-})\,dX^i_s + \frac{1}{2}\sum_{i,j=1}^d \int_0^t f_{,ij}(X_{s-})\,d[X^i,X^j]_s\\ &+\sum_{s\le t}\left(\Delta f(X_s)-\sum_{i=1}^df_{,i}(X_{s-})\Delta X^i_s-\frac{1}{2}\sum_{i,j=1}^d f_{,ij}(X_{s-})\Delta X^i_s\Delta X^j_s\right). \end{align} </math>

This differs from the formula for continuous semimartingales by the additional term summing over the jumps of X, which ensures that the jump of the right hand side at time t is Δf(X_t).

Informal derivation

A formal proof of the lemma requires us to take the limit of a sequence of random variables, which is not done here. Instead, we can derive Itō's lemma by expanding a Taylor series and applying the rules of stochastic calculus.

Assume the Itō process is in the form of

<math> dx= a\,dt + b\,dB.\!</math>

Expanding f(x, t) in a Taylor series in x and t we have

<math> df = \frac{\partial f}{\partial x}\,dx + \frac{\partial f}{\partial t}\,dt + \frac{1}{2}\frac{\partial^2 f}{\partial x^2}\,dx^2 + \cdots </math>

and substituting a dt + b dB for dx gives

<math> df = \frac{\partial f}{\partial x}(a\,dt + b\,dB) + \frac{\partial f}{\partial t}\,dt + \frac{1}{2}\frac{\partial^2 f}{\partial x^2}(a^2\,dt^2 + 2ab\,dt\,dB + b^2\,dB^2) + \cdots. </math>

In the limit as dt tends to 0, the dt² and dt dB terms disappear but the dB² term tends to dt. The latter can be shown if we prove that

<math> dB^2 \rightarrow E(dB^2), </math> since <math> E(dB^2) = dt. \, </math>

The proof of this statistical property is however beyond the scope of this article.

Deleting the dt² and dt dB terms, substituting dt for dB², and collecting the dt and dB terms, we obtain

<math> df = \left(a\frac{\partial f}{\partial x} + \frac{\partial f}{\partial t} + \frac{1}{2}b^2\frac{\partial^2 f}{\partial x^2}\right)dt + b\frac{\partial f}{\partial x}\,dB </math>

as required.

The formal proof is beyond the scope of this article.

Examples

Geometric Brownian motion

A process S is said to follow a geometric Brownian motion with volatility σ and drift μ if it satisfies the stochastic differential equation dS = S(σdB + μdt), for a Brownian motion B. Applying Itō's lemma with f(S) = log(S) gives

<math>

\begin{align} d\log(S) &= f^\prime(S)\,dS + \frac{1}{2}f^{\prime\prime}(S)S^2\sigma^2\,dt \\ &= \frac{1}{S}\left( \sigma S\,dB + \mu S\,dt\right) - \frac{1}{2}\sigma^2\,dt \\ &= \sigma\,dB +(\mu-\sigma^2/2)\,dt. \end{align} </math> It follows that log(S_t) = log(S₀) + σB_t + (μ - σ²/2)t, and exponentiating gives the expression for S,

<math>S_t=S_0\exp\left(\sigma B_t+(\mu-\sigma^2/2)t\right).</math>

The Doléans exponential

The Doléans exponential (or stochastic exponential) of a continuous semimartingale X can be defined as the solution to the SDE dY = YdX with initial condition Y₀ = 1. It is sometimes denoted by Template:Unicode(X). Applying Itō's lemma with f(Y)=log(Y) gives

<math>

\begin{align} d\log(Y) &= \frac{1}{Y}\,dY -\frac{1}{2Y^2}\,d[Y] \\ &= dX - \frac{1}{2}\,d[X]. \end{align} </math> Exponentiating gives the solution

<math>Y_t = \exp\left(X_t-X_0-[X]_t/2\right).</math>

Black–Scholes formula

Itō's lemma can be used to derive the Black–Scholes formula for an option. Suppose a stock price follows a Geometric Brownian motion given by the stochastic differential equation dS = S(σdB + μdt). Then, if the value of an option at time t is f(t,S_t), Itō's lemma gives

<math>df(t,S_t) = \left(\frac{\partial f}{\partial t} + \frac{1}{2}\left(S_t\sigma\right)^2\frac{\partial^2 f}{\partial S^2}\right)\,dt +\frac{\partial f}{\partial S}\,dS_t.</math>

The term (∂f/∂S) dS represents the change in value in time dt of the trading strategy consisting of holding an amount ∂f/∂S of the stock. If this trading strategy is followed, and any cash held is assumed to grow at the risk free rate r, then the total value V of this portfolio satisfies the SDE

<math> dV_t = r\left(V_t-\frac{\partial f}{\partial S}S_t\right)\,dt + \frac{\partial f}{\partial S}\,dS_t.</math>

This strategy replicates the option if V = f(t,S). Combining these equations gives the celebrated Black-Scholes equation

<math>\frac{\partial f}{\partial t} + \frac{1}{2}\sigma^2S^2\frac{\partial^2 f}{\partial S^2} + rS\frac{\partial f}{\partial S}-rf = 0.</math>

See also

• Wiener process
• Itō calculus
• Feynman–Kac formula

References

• Kiyoshi Itō (1951). On stochastic differential equations. Memoirs, American Mathematical Society 4, 1–51.
• Hagen Kleinert (2004). Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets, 4th edition, World Scientific (Singapore); Paperback ISBN 981-238-107-4. Also available online: PDF-files. This textbook also derives generalizations of Itō's lemma for non-Wiener (non-Gaussian) processes.
• Bernt Øksendal (2000). Stochastic Differential Equations. An Introduction with Applications, 5th edition, corrected 2nd printing. Springer. ISBN 3-540-63720-6. Sections 4.1 and 4.2.
• Domingo Tavella (2002). Quantitative Methods in Derivatives Pricing: An Introduction to Computational Finance, John Wiley and Sons. ISBN 9780471274797. Pages 36–39.

Notes

External links

• Derivation, Prof. Thayer Watkins
• Discussion, quantnotes.com
• Informal proof, optiontutorde:Lemma von Itō

fr:Lemme d'Itô it:Lemma di Itô he:הלמה של איטו ja:伊藤の補題 pt:Lema de Itō sv:Itōs lemma zh-yue:伊藤引理