Игла Бюффона
Про меня
  • Закончил Матфак ВШЭ
  • Работаю в 179 школе
  • Преподаю математику
    школьникам, студентам и
    взрослым, организую
    математические события
  • Веду телеграм-канал
    Кроссворд Тьюринга
    @turings_crossword
  • Телеграм @d1_d57

Жорж-Луи Бюффон (1707–1788) — французский биолог,
математик и автор огромной Естественной истории

Автор одной из самых известных задач
геометрической вероятности — об игле Бюффона

Будем бросать иглу на плоскость, разлинованную
параллельными прямыми.
Расстояние между линиями равно длине иглы.
Вопрос: с какой вероятностью игла пересечёт линию?

Сегодня мы ответим на этот вопрос.
Неожиданно, ответ оказывается связан с числом π\boldsymbol{\pi}
Дисклеймер: рассуждения в лекции будут достаточно неформальными

После каждого броска найдём долю игл, пересекших линию.
Во всех сериях она стабилизируется около одного и того же p\boldsymbol{p}.
Наша задача — узнать, чему равно p\boldsymbol{p}?

Начнём с другой задачи Бюффона. Будем бросать монетку диаметра 12\boldsymbol{\tfrac12}.
С какой вероятностью она пересечёт линию?
Важно только положение центра: монетка пересекает линию, когда её центр
находится от неё на расстоянии меньше 14\boldsymbol{\tfrac14}

Будем описывать положение монетки расстоянием от её
центра до верхней линии полосы, в которую он попал

Это число принимает случайные значения от 0\boldsymbol{0} до 1\boldsymbol{1}
Монетка пересечёт линию, если число попадает в две
крайние четверти. Они занимают половину отрезка

Значит, вероятность равна 12\boldsymbol{\tfrac12}

Чему равна вероятность, если диаметр d\boldsymbol{d}?

Тогда вероятность равна d\boldsymbol{d}

Доказательство I
Магическое

Разрежем иглу на три равные части. Сложим из
них треугольник и будем бросать на плоскость
Сколько (примерно) будет точек пересечения
треугольников и линий?

Покрасим стороны треугольников в зелёный,
жёлтый и оранжевый

Сначала мы найдём количество пересечений
линий со сторонами жёлтого цвета
Пересечений с зелёными и оранжевыми будет
столько же, а всех вместе в три раза больше

Как часто жёлтая сторона пересекает линию?
Сотрём все стороны, кроме жёлтой. У жёлтого отрезка совершенно
случайное (равновероятное) положение середины и направление

Жёлтая сторона имеет ту же статистику, что и отдельно брошенная
игла длины 13\boldsymbol{\tfrac13}. Она пересекает линию в 3\boldsymbol{3} раза реже, чем игла длины 1\boldsymbol{1}
Бросим треугольник N\boldsymbol{N} раз. Жёлтая сторона пересечёт линию ≈ p3N\boldsymbol{\tfrac{p}{3}\cdot N} раз

Разрежем иглу на три равные части. Сложим из
них треугольник и будем бросать на плоскость
Сколько (примерно) будет пересечений?Сколько (примерно) будет пересечений?

За N\boldsymbol{N} бросков
• каждая сторона пересечёт линию ≈p3N\boldsymbol{\tfrac{p}{3}\cdot N} раз
• всего будет ≈pN\boldsymbol{p\cdot N} пересечений

Правильный треугольник со стороной 13\boldsymbol{\tfrac13}
(как и игла длины 1\boldsymbol{1}) даёт ≈pN\boldsymbol{p\cdot N} пересечений
С какой вероятностью он пересекает линию?

Треугольник либо не пересекает линию, либо пересекает её дважды.
После N\boldsymbol{N} бросков ≈ p2N\boldsymbol{\tfrac{p}{2}\cdot N} треугольников пересекает линию.
Вероятность пересечения равна p2\boldsymbol{\tfrac{p}{2}}

А если вместо треугольника бросать квадратА если вместо треугольника бросать квадрат
периметра 1\boldsymbol{1}? пятиугольник? восьмиугольник?периметра 1\boldsymbol{1}? пятиугольник? восьмиугольник?
С какой вероятностью они пересекут линию?

Погнём иглу длины 1\boldsymbol{1} на 8\boldsymbol{8}k\boldsymbol{k} равных частей и сложим из
них правильный многоугольник. Бросим его N\boldsymbol{N} раз:
• каждая сторона пересечёт линию ≈p8N\boldsymbol{\tfrac{p}{8}\cdot N}pkN\boldsymbol{\tfrac{p}{k}\cdot N} раз,
• всего произойдёт ≈pN\boldsymbol{p\cdot N} пересечений,
• многоугольник пересекает линию 0\boldsymbol{0} или 2\boldsymbol{2} раза,
• ≈p2N\boldsymbol{\tfrac{p}{2}\cdot N} многоугольников пересекает линии.

Значит, вероятность пересечения равна p2\boldsymbol{\tfrac{p}{2}}

Что произойдёт, если мы будем увеличивать
число углов «до бесконечности»?

Получится окружность длины 1\boldsymbol{1}
Она пересекает линию с вероятностью p2\boldsymbol{\tfrac{p}{2}}

В начале лекции мы видели, что та же
вероятность равна d=1π\boldsymbol{d=\tfrac1\pi}. Следовательно

p2=1πp=2π\boldsymbol{\tfrac{p}{2}=\tfrac1\pi\Rightarrow p=\tfrac2\pi}

Мы получили, что p=2π\boldsymbol{p=\tfrac{2}{\pi}} 0.637\boldsymbol{\approx 0.637}

На самом деле по дороге мы доказали больше:
Если бросать на плоскость контур длины L\boldsymbol{L},
среднее число пересечений равно 2Lπ\boldsymbol{\tfrac{2L}{\pi}}

Для выпуклого контура диаметра меньше 1\boldsymbol{1} вероятность пересечения с линией равна Lπ\boldsymbol{\tfrac{L}{\pi}}

Мы получили ответ в задаче Бюффона очень
красивым, но довольно неожиданным путём.
Посмотрим на тот же вопрос другим способом
напрямую, через описание нашего эксперимента

Доказательство II
Геометрическое

Будем бросать монетку радиуса 14\boldsymbol{\tfrac14} на плоскость,
разбитую на квадраты. С какой вероятностью она
пересечёт линии сетки? накроет центр?

Чтобы найти вероятность, надо понять, куда должен
попасть центр монетки внутри этого квадрата.
Вероятность равна площади этой области

В первом случае центр должен попасть в рамку
ширины 14\boldsymbol{\tfrac14} около границы, во втором в круг радиуса 14\boldsymbol{\tfrac14}
Вероятность пересечь линии 34\boldsymbol{\tfrac34}, накрыть центр π16\boldsymbol{\tfrac{\pi}{16}}

Будем описывать положение иглы двумя числами.
Выберем нижний конец иглы. Пусть
y\boldsymbol{y} — расстояние от конца до верхнего края полосы
α\boldsymbol{\alpha} — угол между иглой и вертикалью
Тогда все положения иглы — это точки прямоугольника0y1, π2απ2\boldsymbol{0\leq y\leq 1,\ -\tfrac\pi2\leq\alpha\leq\tfrac\pi2}

При каких y\boldsymbol{y} и α\boldsymbol{\alpha} игла пересекает линию?
Игла пересекает верхнюю границу тогда и только тогда, когда ей
хватает вертикальной проекции, чтобы до неё дотянуться.

Расстояние до верхней границы y\boldsymbol{y}
Длина иглы равна 1\boldsymbol{1}, а её вертикальная проекция равна cosα\boldsymbol{\cos\alpha}
Значит, пересечение происходит тогда и только тогда, когда ycosαy\,\boldsymbol{\leq\cos\alpha}

Исходы эксперимента — это прямоугольник площади π\boldsymbol{\pi}
Благоприятные исходы — точки под y=cosα\boldsymbol{y=\cos\alpha}. Поэтому

p=Sπ2π=SπS=2\boldsymbol{p=\tfrac{S}{\pi}\Rightarrow\tfrac2\pi=\tfrac{S}{\pi}\Rightarrow S=2}

Можно ли найти S\boldsymbol{S} другим способом?

Точка P\boldsymbol{P} с единичной скоростью проходит четверть
единичной окружности. Параметр α\boldsymbol{\alpha} — это время.

Когда P\boldsymbol{P} находится в положении α\boldsymbol{\alpha}, скорость её
проекции x\boldsymbol{x} равна cosα\boldsymbol{\cos\alpha}. Значит, сверху изображён
график скорости x\boldsymbol{x}.

Площадь под графиком равна пути x\boldsymbol{x}, то есть 1\boldsymbol{1}

Что почитать?

Про задачу Бюффона есть несколько хороших источников —
от популярных текстов до более обстоятельных книг

Спасибо за внимание!
← → листать · F — во весь экран · B — чёрный экран · O — обзор