next up previous
Postscript version of these notes

STAT 380 Lecture 22

Properties of Brownian motion

align18

Notice the use of independent increments and of tex2html_wrap_inline180 .

Suppose t< s. Then tex2html_wrap_inline186 is a sum of two independent normal variables. Do following calculation:

tex2html_wrap_inline188 , and tex2html_wrap_inline190 independent. Z=X+Y.

Compute conditional distribution of X given Z:

align28

Now Z is tex2html_wrap_inline200 where tex2html_wrap_inline202 so

align37

for suitable choices of a and b. To find them compare coefficients of tex2html_wrap_inline208 , x and 1.

Coefficient of tex2html_wrap_inline208 :

displaymath216

so tex2html_wrap_inline218 .

Coefficient of x:

displaymath222

so that

displaymath224

Finally you should check that

displaymath226

to make sure the coefficients of 1 work out as well.

Conclusion: given Z=z the conditional distribution of X is tex2html_wrap_inline234 with a and b as above.

Application to Brownian motion:

The Reflection Principle

Tossing a fair coin:

HTHHHTHTHHTHHHTTHTH tex2html_wrap428
THTTTHTHTTHTTTHHTHT tex2html_wrap430

Both sequences have the same probability.

So: for random walk starting at stopping time:

Any sequence with k more heads than tails in next m tosses is matched to sequence with k more tails than heads. Both sequences have same prob.

Suppose tex2html_wrap_inline272 is a fair (p=1/2) random walk. Define

displaymath276

Compute tex2html_wrap_inline278 ? Trick: Compute

displaymath280

First: if tex2html_wrap_inline282 then

displaymath284

Second: if tex2html_wrap_inline286 then

displaymath288

Fix y < x. Consider a sequence of H's and T's which leads to say T=k and tex2html_wrap_inline294 .

Switch the results of tosses k+1 to n to get a sequence of H's and T's which has T=k and tex2html_wrap_inline302 . This proves

displaymath304

This is true for each k so

align90

Finally, sum over all y to get

align92

Make the substitution k=2x-y in the second sum to get

align96

Brownian motion version:

displaymath312

displaymath314

(called hitting time for level x). Then

displaymath318

Any path with tex2html_wrap_inline320 and X(t)= y< x is matched to an equally likely path with tex2html_wrap_inline320 and X(t)=2x-y> x.

So for y> x

displaymath330

while for y< x

displaymath334

Let tex2html_wrap_inline336 to get

align101

Adding these together gives

align103

Hence tex2html_wrap_inline338 has the distribution of |N(0,t)|.

On the other hand in view of

displaymath318

the density of tex2html_wrap_inline344 is

displaymath346

Use the chain rule to compute this. First

displaymath348

where tex2html_wrap_inline350 is the standard normal density

displaymath352

because P(N(0,1)> y) is 1 minus the standard normal cdf.

So

align113

This density is called the Inverse Gaussian density. tex2html_wrap_inline344 is called a first passage time

NOTE: the preceding is a density when viewed as a function of the variable t.

Martingales

A stochastic process M(t) indexed by either a discrete or continuous time parameter t is a martingale if:

displaymath364

whenever s< t.

Examples

Note: Brownian motion with drift is a process of the form

displaymath374

where B is standard Brownian motion, introduced earlier. X is a martingale if tex2html_wrap_inline380 . We call tex2html_wrap_inline382 the drift

Some evidence for some of the above:

Random walk: tex2html_wrap_inline400 iid with

displaymath402

and tex2html_wrap_inline404 with tex2html_wrap_inline406 . Then

align139

Things to notice:

tex2html_wrap_inline408 treated as constant given tex2html_wrap_inline410 .

Knowing tex2html_wrap_inline410 is equivalent to knowing tex2html_wrap_inline414 .

For j> k we have tex2html_wrap_inline418 independent of tex2html_wrap_inline414 so conditional expectation is unconditional expectation.

Since Standard Brownian Motion is limit of such random walks we get martingale property for standard Brownian motion.

Poisson Process: tex2html_wrap_inline422 . Fix t> s.

align144

Things to notice:

I used independent increments.

tex2html_wrap_inline426 is shorthand for the conditioning event.

Similar to random walk calculation.


next up previous



Richard Lockhart
Wednesday November 22 12:00:14 PST 2000