Buffon's needle problem
By Yatir Halevi on Tuesday, July 30, 2002
- 06:51 pm:
I read once that the probability that a pin dropped from a certain height from
a plain that consists of parallel lines will not fall on any line has it in
p.
Does any body know anything about this?
Yatir
By Andre Rzym on Wednesday, July 31, 2002
- 08:30 pm:
Yatir,
Lets take the case when the length of the needle=separation of
the lines, d. Hopefully you can then solve the general
case.
Suppose the lines are at y=0, y=+/-d, y=+/2d etc.
Consider a square of side d centred on the origin. If we can
figure out the probability of the pin touching a line
given that the centre of the pin falls in the square, then
can you see that this is the result for the entire plane?
What does it mean for the pin to fall 'randomly'? The centre of
the pin can lie anywhere in the square with equal probability and
at any angle to the x axis with equal probability. Therefore the
probability of the pin being in
[x,x+dx; y,y+dy; q+dq] is dxdydq/2d2p.
As a check,
|
|
ó õ
|
d/2
x=-d/2
|
|
ó õ
|
d/2
y=-d/2
|
|
ó õ
|
2p
q = 0
|
dxdydq/2d2p = 1 |
|
which is a good start.
Now define 1(x,y,q) to be a function that is 1 if the pin
touches a line and 0 otherwise. Then the probability of the pin touching is
|
|
ó õ
|
d/2
x=-d/2
|
|
ó õ
|
d/2
y=-d/2
|
|
ó õ
|
2p
q = 0
|
1(x,y,q )dxdydq/2d2p |
|
Now 1 is not actually dependent on y, so
|
|
ó õ
|
d/2
x=-d/2
|
|
ó õ
|
d/2
y=-d/2
|
|
ó õ
|
2p
q = 0
|
1(x,y,q )dxdydq/2d2p = |
ó õ
|
d/2
x=-d/2
|
|
ó õ
|
2p
q = 0
|
1(x,q) dxdq/2dp |
|
Symmetry for positive and negative x reduces this to
|
|
ó õ
|
d/2
x=0
|
|
ó õ
|
2p
q = 0
|
1(x,q)dxdq/dp = |
ó õ
|
2p
q = 0
|
|
ó õ
|
d/2
x=0
|
1(x,q)dxdq/dp |
|
Now the pin touches the line (i.e. 1=1) if |d/2.cos(q)| > x,
so 1 is a function of
By Yatir Halevi on Wednesday, July 31,
2002 - 09:16 pm:
I don't know much probability, so can you elaborate on how
exactly you reached the integral. I guess it has to do with the
fact you added an infinite number of small probabilities.
(definition of an integral)...
Yatir
By Andre Rzym on Thursday, August 01, 2002
- 05:12 pm:
Yatir, you are correct, it is a question of sums converging to
an integral.
First, a few bits of basic probability. If there are exactly 4 possible
outcomes (a,b,c,d), with probability Pa, ..., Pd then we know
Pa+Pb+Pc+Pd=1
If we have a function f which has values depending on the outcomes, then
the expected value of f is
f(a).Pa+f(b).Pb+f(c).Pc+f(d).Pd
Now suppose our trial yields not 4 discrete outcomes but x, a continuous
variable for example in the range [0,1].
We could denote the prob. that x is in [0,0.5], [0.5,1.0] by p(0) and
p(0.5). Similarly if we divided the range into 4 parts the probabilities
would be q(0/4), q(1/4), q(2/4), q(3/4). Now our ps will sum up to 1,
just as our qs will. However, [ignoring pathogenic distributions] the more
we divide the interval, the smaller the probabilities will become (think,
for example, of x being uniformly distributed in [0,1]. Then each
p=1/2, each q=1/4).
Instead, therefore, we denote the probability of being in an interval as being
a function multiplied by the width of the interval. Therefore the analogue to
p above would be a function P, such that the probability of being in
[0,0.5], [0.5,1] is P(0)Dx, P(0.5).Dx respectively, where
Dx=0.5. Ditto the analogue of q is Q, the probabilities being
Q(0/4).Dx, Q(1/4).Dx, Q(2/4).Dx, Q(3/4).Dx where
Dx=0.5.
As you might guess, in the limit as Dx® 0 we get the probability of
x being in the range [x,x+dx]=Pr(x)dx. We call Pr(x) the
probability density.
Since the sum of probabilities add up to 1, we must have
òPr(x)dx=1
Similarly the expected value of some function, f(x), is just the sum of the
product of f(x) and Pr(x)dx, i.e.
òf(x)Pr(x)dx
Andre
By Yatir Halevi on Friday, August 09, 2002
- 01:30 pm:
why is the probability:
1/2d2p?
Yatir
By Andre Rzym on Friday, August 09, 2002
- 07:41 pm:
I take it you are referring to the probability density 1/2pd2 in paragraph
4.
Consider, first, the x coordinate of the pin centre. Since we assume that
the distribution is uniform, the probability of it being in {x,x+D}
is proportional to D, i.e.
pr[pin centre in {x,x+D}]=kD
Hence for a small interval, dx,
pr[pin centre in {x,x+dx}]=k. dx
Since we are assuming that the centre of the pin is in [-d/2,d/2] we require
1=ò-d/2d/2k. dx
Þ k=1/dÞ probability density =1/d
The same argument applied to the y coordinate gives
pr[pin centre in {y,y+dy}]=1/d.dy However since we assume that there is
no correlation between the x and y coord,
pr[pin centre in {x,x+dx;y,y+dy}]=pr[pin centre in {x,x+dx} .pr[pin centre in {y,y+dy}]=[1/d.dx][1/d.dy]=dxdy/d2
Now the angle of the pin to the x axis is also assumed to be uniformly
distributed, so
pr[pin centre in {q,q+D}]=k.D
where q is in radians. Since the range of angles is [0,2p] we must
have
1=ò02pk.dqÞ k=1/2p
Using the independence argument again,
pr[pin centre in {x,x+dx;y,y+dy;q,q+dq}]=pr[pin centre in {x,x+dx;y,y+dx}].pr[pin centre in {q,q+dq}] = [dxdy/d2].[dq/2p]
Andre
By Yatir Halevi on Saturday, August 10,
2002 - 10:32 pm:
Thanks a lot Andre, I get it now!
Yatir