advanced mathematical perspectives 1 · 2019-03-26 · galton board the galton board or quincunx...

17
Advanced Mathematical Perspectives 1 Lecture 13: Random Walks Matthew Roughan <[email protected]> www.maths.adelaide.edu.au/matthew.roughan/notes/AMP1/ School of Mathematical Sciences, University of Adelaide

Upload: others

Post on 15-Jul-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Advanced Mathematical Perspectives 1Lecture 13: Random Walks

Matthew Roughan<[email protected]>

www.maths.adelaide.edu.au/matthew.roughan/notes/AMP1/

School of Mathematical Sciences,University of Adelaide

Page 2: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Probability is the very guide of life.Leonard Mlodinow, The Drunkard’s Walk: HowRandomness Rules Our Lives

It takes only one drink to get me drunk. The trouble is, Ican’t remember if it’s the thirteenth or the fourteenth.

George Burns

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 2 / 16

Page 3: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Section 1

The Drunkard’s Walk

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 3 / 16

Page 4: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

A Simple Random Walk

i-2 i-1 i i+1 i+2

1/2 1/2

Imagine moving around the integers using the following rule:I starting at 0I at each time step, toss a fair coinI if heads, move left by 1I if tails, move right by 1

This is a simple random walk

Questions:I If we did this lots of times, what is the probability distribution across

the set of points?I If we had barriers, how long would it take to hit them?I How often does the walk cross over itself (particularly in 2D)?I ...

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 4 / 16

Page 5: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Random Walk Generalisations

We can generalise this in various ways

allow a biased coin

allow larger jumps

jump around on a 2D lattice of points

But for the moment let’s keep it simple

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 5 / 16

Page 6: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Random Walk Mathematics

Take a series of random variables {Xi} for i = 1, 2, . . . defined by

Xi =

{1, with probability 1/2,−1, with probability 1/2.

Now we could describe the state of our random walk at time n as arandom variable Sn, defined by S0 = 0 and

Sn =n∑

i=1

Xi

This is a very common type of random process, and often analysed.

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 6 / 16

Page 7: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Bernoulli distribution

A Bernoulli random variable is X such that

Xi =

{0, with probability 1− p,1, with probability p.

Think of it as flipping a biased coin, with probability p of heads (or a“success”)

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 7 / 16

Page 8: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Binomial distribution

The Binomial random variable Y ∼ B(n, p) is what we get when we sum nBernoulli random variables. It has probability distribution:

Prob (Y = k) =

(n

k

)pk(1− p)n−k .

It has mean E [Y ] = np and variance Var (Y ) = np(1− p).

The random-walk random variable is almost the same as the Binomial,except the components take values -1 and 1 instead of 0 and 1.

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 8 / 16

Page 9: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Random Walk Distribution

The random-walk random variable is almost the same as the Binomial,except the components take values -1 and 1 instead of 0 and 1.

The main difference isI at time 1, we can only be in state -1 or 1I at time 2, we can only be in state -2, 0, or 2

so at odd times = odd state, and even times = even state

So, combining insight from Binomial, and the above we get adistribution following Pascal’s triangle, i.e.,

(nk

)k -4 -3 -2 -1 0 1 2 3 4

Prob (Sn = 0) 12Prob (Sn = 1) 1 1

22Prob (Sn = 2) 1 2 123Prob (Sn = 3) 1 3 3 124Prob (Sn = 4) 1 4 6 4 1

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 9 / 16

Page 10: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Tricks for understanding Random Walk Distribution

Expectation is the average or mean

Expectation is a linear operator which means

E [A + B] = E [A] + E [B]

And so

E [Sn] = E

[n∑

i=1

Xi

]=

n∑i=1

E [Xi ] = 0

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 10 / 16

Page 11: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Tricks for understanding Random Walk Distribution

Variance is a measure of the spread of variability

Variance of independent random variables add

Var (A + B) = Var (A) + Var (B)

And so

Var (Sn) = Var

(n∑

i=1

Xi

)=

n∑i=1

Var (Xi ) =n∑

i=1

1 = n

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 11 / 16

Page 12: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Tricks for understanding Random Walk Distribution

The Central Limit Theorem says that for any sum like this1

√n

[Snn− µ

]→ N(0, σ2),

where N(0, σ2) denotes the normal or Gaussian distribution with mean 0and variance σ2 (where this is the variance of the Xi ).

Here, Xi has µ = 0 and σ2 = 1, so

Sn√n→ N(0, 1),

1There are some conditions, and we need to define the notion of limit forprobabilities more carefully.

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 12 / 16

Page 13: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Galton Board

The Galton board or Quincunx demonstratesthe Central Limit Theorem

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 13 / 16

Page 14: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Tricks for understanding Random Walk Distribution

These three together (or just the last) give us a pretty clear picture thatthe distribution evolves towards a Gaussian (normal) distribution with wellunderstood parameters.

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 14 / 16

Page 15: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Example

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 15 / 16

Page 16: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Activities

Simulate and play with random walks

Start thinking more seriously about your project

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 16 / 16

Page 17: Advanced Mathematical Perspectives 1 · 2019-03-26 · Galton Board The Galton board or Quincunx demonstrates the Central Limit Theorem Matthew Roughan (School of Mathematical Sciences,

Further reading I

Paul C. Bressloff, Stochastic processes in cell biology, ch. Diffusion in Cells:Random Walks and Brownian Motion, Springer, 2014,http://www.springer.com/gp/book/9783319084879.

Sheldon Ross, Introduction to probability models, Academic Press, 2010.

Matthew Roughan (School of Mathematical Sciences, University of Adelaide)AMP1 17 / 16