Жорж-Луи Бюффон (1707–1788) — французский биолог,
математик и автор огромной Естественной истории
Автор одной из самых известных задач
геометрической вероятности — об игле Бюффона
Будем бросать иглу на плоскость, разлинованную
параллельными прямыми.
Расстояние между линиями равно длине иглы.
Вопрос: с какой вероятностью игла пересечёт линию?
Сегодня мы ответим на этот вопрос.
Неожиданно, ответ оказывается связан с числом π
Дисклеймер: рассуждения в лекции будут достаточно неформальными
После каждого броска найдём долю игл, пересекших линию.
Во всех сериях она стабилизируется около одного и того же p.
Наша задача — узнать, чему равно p?
Начнём с другой задачи Бюффона. Будем бросать монетку диаметра 21.
С какой вероятностью она пересечёт линию?
Важно только положение центра: монетка пересекает линию, когда её центр
находится от неё на расстоянии меньше 41
Будем описывать положение монетки расстоянием от её
центра до верхней линии полосы, в которую он попал
Это число принимает случайные значения от 0 до 1
Монетка пересечёт линию, если число попадает в две
крайние четверти. Они занимают половину отрезка
Значит, вероятность равна 21
Чему равна вероятность, если диаметр d?
Тогда вероятность равна d
Разрежем иглу на три равные части. Сложим из
них треугольник и будем бросать на плоскость
Сколько (примерно) будет точек пересечения
треугольников и линий?
Покрасим стороны треугольников в зелёный,
жёлтый и оранжевый
Сначала мы найдём количество пересечений
линий со сторонами жёлтого цвета
Пересечений с зелёными и оранжевыми будет
столько же, а всех вместе в три раза больше
Как часто жёлтая сторона пересекает линию?
Сотрём все стороны, кроме жёлтой. У жёлтого отрезка совершенно
случайное (равновероятное) положение середины и направление
Жёлтая сторона имеет ту же статистику, что и отдельно брошенная
игла длины 31. Она пересекает линию в 3 раза реже, чем игла длины 1
Бросим треугольник N раз. Жёлтая сторона пересечёт линию ≈ 3p⋅N раз
Разрежем иглу на три равные части. Сложим из
них треугольник и будем бросать на плоскость
Сколько (примерно) будет пересечений?Сколько (примерно) будет пересечений?
За N бросков
• каждая сторона пересечёт линию ≈3p⋅N раз
• всего будет ≈p⋅N пересечений
Правильный треугольник со стороной 31
(как и игла длины 1) даёт ≈p⋅N пересечений
С какой вероятностью он пересекает линию?
Треугольник либо не пересекает линию, либо пересекает её дважды.
После N бросков ≈ 2p⋅N треугольников пересекает линию.
Вероятность пересечения равна 2p
А если вместо треугольника бросать квадратА если вместо треугольника бросать квадрат
периметра 1? пятиугольник? восьмиугольник?периметра 1? пятиугольник? восьмиугольник?
С какой вероятностью они пересекут линию?
Погнём иглу длины 1 на 8k равных частей и сложим из
них правильный многоугольник. Бросим его N раз:
• каждая сторона пересечёт линию ≈8p⋅Nkp⋅N раз,
• всего произойдёт ≈p⋅N пересечений,
• многоугольник пересекает линию 0 или 2 раза,
• ≈2p⋅N многоугольников пересекает линии.
Значит, вероятность пересечения равна 2p
Что произойдёт, если мы будем увеличивать
число углов «до бесконечности»?
Получится окружность длины 1
Она пересекает линию с вероятностью 2p
В начале лекции мы видели, что та же
вероятность равна d=π1. Следовательно
2p=π1⇒p=π2
Мы получили, что p=π2 ≈0.637
На самом деле по дороге мы доказали больше:
Если бросать на плоскость контур длины L,
среднее число пересечений равно π2L
Для выпуклого контура диаметра меньше 1 вероятность пересечения с линией равна πL
Мы получили ответ в задаче Бюффона очень
красивым, но довольно неожиданным путём.
Посмотрим на тот же вопрос другим способом —
напрямую, через описание нашего эксперимента
Будем бросать монетку радиуса 41 на плоскость,
разбитую на квадраты. С какой вероятностью она
пересечёт линии сетки? накроет центр?
Чтобы найти вероятность, надо понять, куда должен
попасть центр монетки внутри этого квадрата.
Вероятность равна площади этой области
В первом случае центр должен попасть в рамку
ширины 41 около границы, во втором в круг радиуса 41
Вероятность пересечь линии 43, накрыть центр 16π
Будем описывать положение иглы двумя числами.
Выберем нижний конец иглы. Пусть
• y — расстояние от конца до верхнего края полосы
• α — угол между иглой и вертикалью
Тогда все положения иглы — это точки прямоугольника0≤y≤1, −2π≤α≤2π
При каких y и α игла пересекает линию?
Игла пересекает верхнюю границу тогда и только тогда, когда ей
хватает вертикальной проекции, чтобы до неё дотянуться.
Расстояние до верхней границы y
Длина иглы равна 1, а её вертикальная проекция равна cosα
Значит, пересечение происходит тогда и только тогда, когда y≤cosα
Исходы эксперимента — это прямоугольник площади π
Благоприятные исходы — точки под y=cosα. Поэтому
p=πS⇒π2=πS⇒S=2
Можно ли найти S другим способом?
Точка P с единичной скоростью проходит четверть
единичной окружности. Параметр α — это время.
Когда P находится в положении α, скорость её
проекции x равна cosα. Значит, сверху изображён
график скорости x.
Площадь под графиком равна пути x, то есть 1
Про задачу Бюффона есть несколько хороших источников —
от популярных текстов до более обстоятельных книг