5

In Ito's lemma for a Brownian motion $B_t$ the term in $dB_t^2$ is replaced with $dt$ without any averaging. It seems that higher moments are an order $dt$ smaller and that the term $dB_t^2$ is dominated by its expectation and it becomes deterministic and equal to $dt$.

I would like to demonstrate this by doing a numerical simulation. Now since $dB_t = g_t \sqrt{dt}$ where $g_t \simeq N(0,1)$ I could simply draw a value of $g_t$ and calculate $dB_t^2$. But such a term is random for all values of $dt$ and so does not demonstrate what I want to show.

Alternatively I can split the period $dt$ up into $N$ periods of length $dt/N$. Will this work ?

Dom
  • 256
  • 3
  • 9
  • 2
    Well you don't seem to have a rigorous definition of brownnian motion and of "quadratic variation", for example. $dB^2_t = dt$ a just a symbol, what is actually meant is that $<B_t> = t$. Look up quadratic variation to get more insights, and a correct version of Ito will make use of $< B_t> $, not $dB^2_t$ – Ant Aug 19 '16 at 10:22
  • I am not sure I understand your notation. Is $<B_t>$ an average of $B_t$. In which case it should be zero. – Dom Aug 19 '16 at 10:53
  • @Dom Quote: "Look up quadratic variation". – Did Aug 19 '16 at 11:16
  • I would search for Levy's characterisation of Brownian Motion for example in Philip E. Protter: Stochastic Integration and Differential Equations in order to understand this term better; for numeric interpretation you can use https://www.statistik.tu-dortmund.de/fileadmin/user_upload/Lehrstuehle/Ingenieur/Mueller/Lehre/Seminar/SS2012/GriesingP2.pdf page 14,15 Maybe it helps... – user361972 Aug 19 '16 at 10:58

2 Answers2

5

You seem to be confused by some terminology; I'll try to explain briefly.

Ito's lemma says that (in the context of a one-dimensional brownian motion; it is actually more general, see below):

$$f(B_t) = f(B_0) + \int_0^t f'(B_s) dB_s + \frac 12 \int_0^t f''(B_s) ds$$ for $f \in C^2(\mathbb R, \mathbb R)$. As you can see, the term $dB_s^2$ is nowhere to be seen.

It is useful now to define the quadratic variation of certain stochastic processes:

Let $X_t$ be a square integrable local martingale. Then the quadratic variation $\langle X\rangle_t$ is the unique, right-continuos, increasing adapted process starting at $0$ such that $X^2_t - \langle X\rangle_t$ is a local martingale. It can be shown that $\langle X\rangle_t$ arises as the quadratic variation of the trajectory of $X_\cdot$ along certain partitions, hence the name. One can additionally show that $\langle B \rangle_t = t$. Finally we also define $\langle M, N \rangle $ in a similar way; such that $M_tN_t - M_0 N_0 - \langle M, N \rangle_t $ is a local martingale. Note that $\langle X, X \rangle_t = \langle X \rangle _t$

Why is this important? Because the of the version of Ito's lemma for continuos semimartingales: Let $X_t$ be a continuos $d$-dimensional semimartingale, then

$$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_{i,j}(X_{s})\,d\langle X^i,X^j\rangle_s $$

for $f \in C^2(\mathbb R^d, \mathbb R)$

If $d=1$ and $X_t = B_t$ a brownian motion, then using the fact that $\langle B\rangle_t = t$ one recovers the previous formula.

Another use of quadratic variation is the Ito isometry:

$$E\left[\left(\int_0^t Y_s dX_s\right)^2\right] = E\left[\int_0^t Y_s^2 d\langle X\rangle _s \right]$$. Note that for $X = B$ we get

$$E\left[\left(\int_0^t Y_s dB_s\right)^2\right] = E\left[\int_0^t Y_s^2 ds \right]$$

So as you can see $dB^2_s$ is not present in the above discussion.


Why did you hear then that $dB_s^2 = ds$? This is a simplification and stems from different facts (namely, that $E[B^2_s] = s$ and $\langle B\rangle _s = s$). Let $X_t$ be a Ito diffusion process:

$$dX_t = b(t) dt + \sigma(t) dB_t$$ Using $dB^2_s = ds$ and $dB_s ds = 0$ you get

$$dX^2_t = \sigma^2(s) ds$$

And then you can "say" that the variance of $X^2_t$ is $\int_0^t \sigma^2(s) ds$, which is correct (as you can see properly from Ito's formula). So $dB_s^2 = ds$ can be useful but shouldn't be confused with the rigorous definition of quadratic variation.

Ant
  • 21,098
  • OK. Suppose you wanted to use a Monte-Carlo simulation to verify that $dB_s^2=dt$. What would your code do ? Can you sketch a simple algorithm ? – Dom Aug 19 '16 at 13:36
  • @Dom Well you could work on the fact that $$\lim_{n \to \infty} \sum_{k \ge 0, \frac{k+1}{2^n} \le t} \left| B_{\frac{k+1}{2^n}} - B_{\frac{k}{2^n}} \right|^2 = \langle B \rangle _t = t$$ almost surely and in $L^2$. So you could simulate a trajectory of the brownian motion up until time $t$, divide the interval into pieces of length $2^{-n}$ for $n$ very big, calculate the squared difference at the endpoints, sum everything. You should get $t$ ;-) – Ant Aug 19 '16 at 13:41
  • OK. I have done that. And I get $t$. My issue is that it seems to work for any $t$ provided $n$ is large. It does not require $t \rightarrow 0$. – Dom Aug 19 '16 at 13:53
  • @Dom Did you read my answer? $dB$ does not exist and you should forget about it, think about $\langle B \rangle _t$ instead; and the fact that $\langle B \rangle _t = t$ holds for every $t$ and does not require small $t$ – Ant Aug 19 '16 at 14:17
  • @Dom What I meant was that the equation $dB^2_t = t$ does not exists, does not make sense, you should forget it and think about $\langle B \rangle_t$ instead. If you insist, you could look it another way and say that $$dB^2_t = (B_{t + dt} - B_{t})^2$$; taking expectation on both sides you have $$E[dB^2_t] = dt$$. Then you could say that "since $dt$ is small, $dB^2_t$ is almost constant and its expectation is itself,to find $$dB^2_t = dt$$. But 1) It does not make sense and 2) The only numerical verification to do is to show that the variance of $B_t $ is $t$, but that's trivial – Ant Aug 19 '16 at 17:13
0

Expanding a bit on Ant's last point so that your intuition might be satisfied:

We know that $E[dB_t^2] = dt$, but you can also easily show that $\text{var}[dB_t^2] = 2dt^2$ by noting that $\frac{B_t^2}{t}$~$\chi^2(1)$. Any RV with zero variance is a constant, so we have $dB_t^2$~$dt$, in exactly the same way that you have $dt^2$~$0$ (i.e. in the limit as $dt \to 0$).