Stochastic Integration

Mathematics — Alex @ 5:04 pm

Integration is one of those mathematical operations that should conform to common-sense: no matter what method of integration you’re using, no matter how high-falutin’, the area of a square should be the square of the length of its sides. And that’s the way its been up to now: Riemann, Lebesgue-Stieltjes, and Lebesgue integration… they’re all cut from the same cloth, all succeedingly more careful generalizations of the intuitive notion of measure.

Yesterday we were briefly introduced to a form of stochastic integration, namely that of integration of L^2(\R) functions with respect to Brownian noise:

 \int_0^t f(s) dB_s = \sum_{i=1}^\infty \langle \chi_{[0,t]} f, \varphi_i \rangle \xi_i

where \varphi is an orthonormal basis of L^2 and \xi_i is a sequence of i.i.d. N(0,1) r.v.s, so the integral is a Gaussian r.v. for each t. This definition does not seem to have much to do with the concept of measurement. Is this even well-defined? By the Riemann-Lebesgue lemma the sum converges, but apparently any o.n.b. can be used, so does changing it affect the integral? Despite the obscurity of its origins, I like the simple form of this integral.

Today, I found an online book which does a great job of motivating stochastic integrals. To abbreviate the example he gives, consider the problem of sending a probe to the moon. We’ll describe the state of the system as X(t) = (\ell, v, \rho)^t \in \R^9, where \ell, v, \rho are the location, velocity, and momentum of the probe. In the ideal case, we know X’(t) = (v, d\rho/dt, F(\ell, v, \rho))^t where F is the force applied to the probe, so we can solve for

 X(t) = X(0) + \int_0^t X’(s) \: ds.

In the real world, there will be random effects we can’t account for deterministically, such as minor fluctuations in the gravitational field, collision with debris, etc., all of which we will call noise Z(t) (note that this model doesn’t account for noise which is partially caused by the state X(t), noise is a function of t– I don’t know if this is a limitation of the theory, or just this example). Then we’d like to say that how much the noise affects the state of the system varies depending on the state of the system, giving the new equation

 X(t) = X(0) + \int_0^t X’(s) \: ds + \int_0^t b(X(s)) \: dZ(s).

Since the noise is random, Z(s) = Z(s,\omega) is in fact a stochastic process. Each elementary event \omega determines a different sample path. Therefore if Z is such that the Lebesgue-Stieltjes integral is defined on each sample path (cadlag?), we can actually solve the above equation, and use it to calculate useful quantities. For instance, if T(\omega) represents the time when the probe reaches the moon (\infty if it misses), then we can calculate the expected distance from the target point s, given a particular noise model:

 \mathbb{E}[\|s - X(T)\|].

Pretty cool!

The problem is, apparently in general the sample paths of stochastic processes are not nice, so Lebesgue-Stieltjes integration cannot be done. Then you need to break out the power tools

Possibly relevant posts:

6 Comments »

  1. wow, alex, when reading your blog, your protriat shows up again in front of me, holding a library book stopping people in the hall way. i actually found some familiar concepts, like “gram schimidts” process, yes! ha, it’s not always easy to undertand what you are talking all about the math, but wish you every thing well in your grad study.

    Comment by ying hu — 11/17/2006 @ 2:19 pm
  2. Hi, Alex

    Really exciting formulation of stochastic integral. what book is it from ?

    thanks
    Cepreu

    Comment by Cepreu — 12/16/2006 @ 4:27 am
  3. Hi Alex,

    very nice and inspirational formulation of stochastic integral!
    What book/lectures is it from ?

    thanks,
    Ceprey

    Comment by Cepera — 12/18/2006 @ 7:44 am
  4. Cepera, I posted a link above to the book I saw that explanation in. It’s also available as a volume of the Encyclopedia of Mathematics and Its Applications, as “Stochastic Integration with Jumps”. Despite the very nice motivation of the concept of a stochastic integral, it’s not a good book for an introduction to the subject. I like Introduction to Stochastic Integration by Kuo better for that.

    Comment by Alex — 12/22/2006 @ 1:49 pm
  5. Alex,

    I meant first definition, where stochastic integral is represented by sum of function “f” projections onto basis multiplied by random Gaussian numbers.
    i didn’t find it in Bichteler’s book.

    thanks,
    C

    Comment by Cepreu — 12/30/2006 @ 6:34 am
  6. Ahh.. that was *given* to us in class, out of the blue, as our first working definition of a stochastic integral.

    Comment by Alex — 1/1/2007 @ 10:00 pm

RSS feed for comments on this post. TrackBack URI

Leave a comment

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 Unported License.
(c) 2008 ChapterZero | powered by WordPress with Barecity