Stochastic calculus

Overview

When you have a function X, interpreted as a table

t X(t)
0 1
1 2
2 5

you can think of the "discrete derivative" ΔX as the extended table

t ΔX(t) X(t)
0 - 1
1 1 2
2 3 5

For the understanding of stochastic calculus, it is better to think that the usual derivative is not a "real" thing: it emerges when you take increasingly smaller steps Δt for the t column and you consider the ratio ΔX(t)/Δt.
If we stick to real life, with real measurements, we only have this kind of discrete tables, not functions. The problem of recovering X from ΔX, or even from dX/dt, is called integration. And it is written as

dXdt=g(t),

or

dX=g(t)dt,

where g(t) is the known data. With this setup,

X=X0+g(t)dtX0+ΔX=X0+g(t)Δt.

So far, so good. But, what if ΔX(t) does not correspond to a fixed list of values but, instead, to a random list?
For example, suppose that the contents of the list are generated by multiplying a certain quantity H times a random number from a Gaussian. This is usually represented as H(t)ΔW, where W stands for Wiener process. We enter into the realm of stochastic calculus. Indeed, the reconstruction of a function from this noisy derivative could be done in several different ways, for example, Itô integral,

Xt=X0+0tH(s)dWsX0+iH(ti)ΔWi.

or Stratonovich integral

Xt=X0+0tH(s)dWsX0+iH(ti+ti+12)ΔWi.

The resulting object is a stochastic process, which can be interpreted as a curve-valued random variable, instead of a single curve. The canonical example is Brownian motion as a limit of random walks.

Coming back to classical analysis, it may happen that we want to recover a function X from a table in which the column doesn't solely depend on t, but also on X(t). That is, we would have a kind of 2D table with t free in the first column, X(t) free in the first row, and the values of ΔX in the interior of the table:

tX(t) 0 1 2 3
0 0 1 2 3
1 1 2 3 4
2 2 3 4 5
3 3 4 5 6

This problem is called a first-order ODE. It can be expressed as dXdt=F(t,X(t)) or

(*)dX=F(t,X(t))dt,

and the idea is the same: recover X from the known data F.

Now, like above, the content of the table can also be randomized. For example, it could be a Gaussian, multiplied by t, or by X(t), or in general H(t,Xt)ΔWt. This is known as stochastic ordinary differential equation, denoted as

(**)dXt=H(t,Xt)dWt.

The solution is again a stochastic process, but now is more difficult to find it out. In the same way that equation () must be seen as encoding "any function/table X(t) such that its increments correspond to F(t,X(t))dt", equation () means "any stochastic process such that its increments are H(t,Xt)dWt, with dWt the increments of a Brownian motion".

Finally, we have more general stochastic ordinary differential equation, in which the increments have a deterministic part and a random part:

dXt=F(t,Xt)dt+H(t,Xt)dWt.

Main definitions

%%{init: {'themeVariables': { 'fontSize': '13px', 'nodeSpacing': 30, 'rankSpacing': 25 }}}%%
graph TD
    A["Stochastic Process"]
    subgraph Discrete-Time
        C["Markov Chain"]
        D["Random Walk"]
    end
    subgraph Continuous-Time
        G["Semimartingale"]
        E["Brownian Motion"]
    end
    A -->|Markov property| B["Markov Process"]
    B --> C
    C --> D
    B --> G
    G --> E
    %% Random Walk approximation to Brownian
    D -.->|scaling limit| E
    %% Itô Integral as operator
    G -.->|Itô integral operator| H["New Processes"]
    %% Outputs: Martingales, SDEs, Diffusions
    H --> F["Martingale"]
    H --> I["Itô Diffusion"]
    %% Style
    classDef operator fill:#f9f,stroke:#333,stroke-width:2px;
    class G,H,J operator

Solutions

See solutions of an SDE.

Symmetries

See overview on symmetries of SDEs.