The Azimuth Project
Itô formula

Idea
Details
Examples
- The Hopf-Bifurcation Example from TWF 308
References

Idea

The Itô formula is the chain rule for Itô stochastic calculus. The involvement of the white noise process leads to an inclusion of second order derivatives, instead of first order derivatives as in usual calculus.

Details

We will formulate the formula for a two dimensional process, the general formula as well as the one dimensional formula are simple to get from that.

Let $(W_{1}, W_{2})$ be a two dimensional Wiener process, and let an two dimensional Itô process be defined via

d X = a_{1} (t, x, y) d t + b_{1, 1} (t, x, y) d W_{1} + b_{1, 2} (t, x, y) d W_{2}

and

d Y = a_{2} (t, x, y) d t + b_{2, 1} (t, x, y) d W_{1} + b_{2, 2} (t, x, y) d W_{2}

A coordinate transformation is specified by two $C^{2}$ functions $U : = g_{1} (t, x, y)$ and $V : = g_{2} (t, x, y)$ , defining two new “coordinates” $U, V$ . Note that we need the functions to be twice continuous differentiable, since the transformation formula includes derivatives of the second order.

The Itô formula says that $U$ and $V$ are again the components of a two dimensional Itô process and satisfy the equations

d U = \frac{\partial g_{1}}{\partial t} dt + \frac{\partial g_{1}}{\partial x} d X + \frac{\partial g_{1}}{\partial y} d Y + \frac{1}{2} \frac{\partial^{2} g_{1}}{\partial x^{2}} d X d X + \frac{1}{2} \frac{\partial^{2} g_{1}}{\partial x \partial y} d X d Y + \frac{1}{2} \frac{\partial^{2} g_{1}}{\partial y^{2}} d Y d Y

and

d V = \frac{\partial g_{2}}{\partial t} dt + \frac{\partial g_{2}}{\partial x} d X + \frac{\partial g_{2}}{\partial y} d Y + \frac{1}{2} \frac{\partial^{2} g_{2}}{\partial x^{2}} d X d X + \frac{1}{2} \frac{\partial^{2} g_{2}}{\partial x \partial y} d X d Y + \frac{1}{2} \frac{\partial^{2} g_{2}}{\partial y^{2}} d Y d Y

Further, the Itô formula states that when one inserts the equations for $d X$ and $d Y$ in the above formula, products of “differentials” are all zero with the notable exception

d W_{i} d W_{i} = d t

So we see that, for example, a process without a drift (a is zero) could get transformed in a process with drift by a coordinate transformation.

Examples

The Hopf-Bifurcation Example from TWF 308

In this example we will take a look at the two dimensional system of SDE explained in This Weeks Finds 308.

Polar Coordinates

We start with the Euclidean plane $ℝ^{2}$ and Cartesian coordinates $x, y$ , we would like to do a transformation to polar coordinates $r, ϕ$ :

x = r \cos (ϕ)

y = r \sin (ϕ)

The inverse transformation is

r = (x^{2} + y^{2})^{\frac{1}{2}} = : g_{1} (x, y)

and

ϕ = \arctan (\frac{y}{x}) = : g_{2} (x, y)

for $x > 0$ .

Polar Coordinates, Two Dimensional Random Walk

A simple two dimensional random walk with independent but equally strong noise terms is given by

d X = λ d W_{1}

and

d Y = λ d W_{2}

We calculate the SDE in polar coordinates according to the Itô formula. For $r$ we get:

d r = \frac{\partial g_{1}}{\partial x} d X + \frac{\partial g_{1}}{\partial y} d Y + \frac{1}{2} \frac{\partial^{2} g_{1}}{\partial x^{2}} d X d X + \frac{1}{2} \frac{\partial^{2} g_{1}}{\partial y^{2}} d Y d Y

All the other terms in the Itô formula are zero. Inserting the formulae for the differentials of $g_{1}$ and the SDE we get:

d r = \frac{x}{(x^{2} + y^{2})^{\frac{1}{2}}} λ d W_{1} + \frac{y}{(x^{2} + y^{2})^{\frac{1}{2}}} λ d W_{2} + \frac{1}{2} (\frac{1}{(x^{2} + y^{2})^{\frac{1}{2}}}) λ^{2} d t

Replacing $x, y$ with polar coordinates on the right results in the SDE for $r$ :

d r = \frac{λ^{2}}{2 r} d t + λ \cos (ϕ) d W_{1} + λ \sin (ϕ) d W_{2}

We see that the radius $r$ picks up a drift term that pushes the process away from the origin, the nearer it comes the stronger. This accounts for the effect that the influence of noise tends to move the process away from any concentrically confined area.

Now we turn to the other coordinate $ϕ$ . We keep all nonzero terms in the Itô formula:

d ϕ = \frac{\partial g_{2}}{\partial x} d X + \frac{\partial g_{2}}{\partial y} d Y + \frac{1}{2} \frac{\partial^{2} g_{2}}{\partial x^{2}} d X d X + \frac{1}{2} \frac{\partial^{2} g_{2}}{\partial y^{2}} d Y d Y

As before we insert the formulae for the derivatives of $g_{2}$ and for $d X, d Y$ and get:

d ϕ = \frac{- y}{x^{2} + y^{2}} λ d W_{1} + \frac{x}{x^{2} + y^{2}} λ d W_{2}

The second derivatives cancel. Replacing $x, y$ with polar coordinates on the right results in the SDE for $ϕ$ :

d ϕ = \frac{- \sin (ϕ)}{r} λ d W_{1} + \frac{\cos (ϕ)}{r} λ d W_{2}

So the angle $ϕ$ does not pick up a drift term, but the effects of the noise decrease with the distance from the origin, as one would expect.

The Hopf Bifurcation Example System

There is no complication in applying our results for the “free” two dimensional random walk to the example of week 308, the system of equations in cartesian coordinates is:

d x = (- y + β x - x (x^{2} + y^{2})) d t + λ d W_{1}

and

d Y = (x + β y - y (x^{2} + y^{2})) d t + λ d W_{1}

Transforming to polar coordinates we get

d r = (β r - r^{3} + \frac{λ^{2}}{2 r}) d t + λ \cos (ϕ) d W_{1} + λ \sin (ϕ) d W_{2}

and

d ϕ = \frac{- \sin (ϕ)}{r} λ d W_{1} + \frac{\cos (ϕ)}{r} λ d W_{2}

Now we can calculate the Fokker-Planck operator for this system in polar coordinates, see over there.

References

Ito formula, Wikipedia

category: mathematical methods

Revised on November 7, 2011 14:30:00 by Andrew Stacey (129.241.15.200)

The Azimuth Project Itô formula

Contents