Проблема с иглой Бюффона

В теории вероятностей проблема иглы Бюффона — это вопрос, впервые поставленный в 18 веке Жоржем -Луи Леклерком, графом де Бюффоном : [ 1 ]
- Предположим, у нас есть пол, сделанный из параллельных полос деревянных одинаковой ширины, и мы бросаем иголку на пол . Какова вероятность того , что игла окажется на линии между двумя полосками?
Игла Бюффона была первой проблемой геометрической вероятности , которую нужно было решить; [ 2 ] ее можно решить, используя интегральную геометрию . Решение для искомой вероятности p в случае, когда длина иглы l не превышает ширины t полосок, имеет вид
Это можно использовать для разработки метода Монте-Карло для аппроксимации числа π , хотя это не было первоначальной мотивацией вопроса де Бюффона. [ 3 ] Кажущееся необычным появление числа π в этом выражении происходит потому, что основная функция распределения вероятностей ориентации иглы является вращательно-симметричной.
Решение
[ редактировать ]Проблема в более математическом смысле такова: если игла длиной l упала на плоскость, разделенную параллельными линиями на расстоянии t единиц, какова вероятность того, что игла при приземлении окажется поперек этой линии?
Пусть x — расстояние от центра иглы до ближайшей параллельной линии, а θ — острый угол между иглой и одной из параллельных линий.
Равномерная функция плотности вероятности (PDF) x между 0 и т / 2 есть
Здесь x = 0 представляет собой иглу, центрированную прямо на линии, а x = t / 2 представляет собой иглу, идеально расположенную между двумя линиями. Единый PDF предполагает, что стрелка с одинаковой вероятностью упадет в любом месте этого диапазона, но не сможет выйти за его пределы.
Равномерная функция плотности вероятности θ между 0 и π / 2 это
Здесь θ = 0 представляет собой иглу, параллельную отмеченным линиям, а θ = π / 2 радиан представляет собой иглу, перпендикулярную отмеченным линиям. Любой угол в этом диапазоне считается равновероятным исходом.
Две случайные величины , x и θ , независимы, [ 4 ] поэтому совместная функция плотности вероятности представляет собой произведение
Стрелка пересекает линию, если
Теперь есть два случая.
Случай 1: Короткая игла ( l ≤ t )
[ редактировать ]Интегрирование функции плотности вероятности суставов дает вероятность того, что игла пересечет линию:
Случай 2: Длинная игла ( l > t )
[ редактировать ]Предположим, что l > t . В этом случае, интегрируя совместную функцию плотности вероятности, получаем:
где m ( θ ) — минимум между l / 2 sin θ и t / 2 .
Таким образом, выполняя описанное выше интегрирование, мы видим, что при l > t вероятность того, что игла пересечет хотя бы одну линию, равна
или
Во втором выражении первый член представляет вероятность того, что угол иглы будет таким, что она всегда пересечет хотя бы одну линию. Правильный член представляет собой вероятность того, что игла упадет под углом, где ее положение имеет значение, и пересечет линию.
В качестве альтернативы обратите внимание, что всякий раз, когда θ имеет такое значение, что l sin θ ≤ t , то есть в диапазоне 0 ≤ θ ≤ arcsin t / l , вероятность пересечения такая же, как и в коротком игольнице. Однако если l sin θ > t , то есть arcsin т / л < θ ≤ π / 2 вероятность постоянна и равна 1.
Использование элементарного исчисления
[ редактировать ]Следующее решение для случая «короткой иглы», хотя и эквивалентно приведенному выше, имеет более визуальный вид и позволяет избежать повторных интегралов.
Мы можем вычислить вероятность P как произведение двух вероятностей: P = P1 вероятность того · P2 , , где P1 — что центр иглы окажется достаточно близко к линии, чтобы игла могла пересечь ее, и P2 . — это вероятность того, что игла действительно пересечет линию, при условии, что центр находится в пределах досягаемости.
Глядя на иллюстрацию в приведенном выше разделе, становится очевидным, что игла может пересечь линию, если центр иглы находится в пределах л / 2 ед. с каждой стороны полосы. Добавление l / 2 + l / 2 с обеих сторон и разделив на всю ширину t , получим P 1 = l / t .
Теперь предположим, что центр находится в пределах досягаемости края полосы, и вычислим P 2 . Для упрощения расчета можно предположить, что .
Пусть x и θ такие, как на иллюстрации в этом разделе. Поместив центр иглы в точку x , игла пересечет вертикальную ось, если она попадает в диапазон 2 θ радиан из π радиан возможных ориентаций. Это представляет собой серую область слева от x на рисунке. Для фиксированного x мы можем выразить θ как функцию x : θ ( x ) = arccos( x ) . Теперь мы можем указать x в диапазоне от 0 до 1 и интегрировать:
Умножая оба результата, получаем P = P 1 · P 2 = l / t · 2 / π = 2 l / tπ как указано выше.
Есть еще более изящный и простой метод расчета «короткой иголки». Конец иглы, наиболее удаленный от любой из двух линий, граничащих с ее областью, должен располагаться на расстоянии по горизонтали (перпендикулярно граничащим линиям) l cos θ (где θ — угол между иглой и горизонталью) от этой области. линию, чтобы игла ее пересекла. Максимальное расстояние, на которое этот конец иглы может отойти от этой линии по горизонтали в ее области, равно t . Вероятность того, что самый дальний конец иглы находится не дальше, чем на расстоянии l cos θ от линии (и, таким образом, что игла пересекает линию) из общего расстояния t , на которое она может переместиться в своей области при 0 ≤ θ ≤ π / 2 определяется выражением
Без интегралов
[ редактировать ]Задачу о короткой игле можно также решить без какого-либо интегрирования таким образом, чтобы объяснить формулу для p тем геометрическим фактом, что круг диаметром t будет пересекать расстояние t полосок всегда (т. е. с вероятностью 1) ровно в двух точках. Это решение было предложено Жозефом-Эмилем Барбье в 1860 году. [ 5 ] Его еще называют « лапшой Бюффона ».
Оценка π
[ редактировать ]
2 л · п / ть = 2 × 9 × 17 / 9 × 11 ≈ 3,1 ≈ π .

В первом, более простом случае, приведенном выше, формулу, полученную для вероятности P, можно переписать в виде
Таким образом, если мы проведем эксперимент для оценки P , у нас также будет оценка для π .
Предположим, мы уронили n игл и обнаружили, что h из этих игл пересекают линии, поэтому P аппроксимируется дробью час / п . Это приводит к формуле:
В 1901 году итальянский математик Марио Лаццарини провел эксперимент Бюффона с иглой. Подбросив иголку 3408 раз, он получил известное приближение 355/113 с точностью до шести для π знаков после запятой. [ 6 ] «Эксперимент» Лаццарини является примером предвзятости подтверждения , поскольку он был создан для того, чтобы воспроизвести уже известную аппроксимацию 355/113 , Milü (фактически, не существует лучшего рационального приближения с менее чем пятью цифрами в числителе и знаменателе, см. также ) , дающее более точный «предсказание» π чем можно было бы ожидать по количеству испытаний, следующее: [ 7 ]
Лаццарини выбрал иглы длиной 5/6 планок . ширины деревянных В этом случае вероятность того, что иглы пересекут линии, равна 5 / 3 π . Таким образом, если бы нужно было опустить n игл и получить x пересечений, можно было бы оценить π как
Так что если Лаццарини стремился к результату 355 / 113 , ему нужны были n и x такие, что
или эквивалентно,
Для этого нужно выбрать n кратным 213, так как тогда 113 n / 213 — целое число; Затем роняют n иголок и надеются ровно на x = 113 п / 213 успехов. Если вы уроните 213 иголок и получите 113 успехов, то сможете с триумфом сообщить об оценке π с точностью до шести знаков после запятой. Если нет, то можно просто провести еще 213 попыток и надеяться на 226 успехов; если нет, просто повторите при необходимости. Лаццарини выполнил 3408 = 213 × 16 испытаний, поэтому вполне вероятно, что именно эту стратегию он использовал для получения своей «оценки».
Приведенное выше описание стратегии можно даже считать милосердием по отношению к Лаццарини. Статистический анализ промежуточных результатов, о которых он сообщил для меньшего числа бросков, приводит к очень низкой вероятности достижения столь близкого соответствия ожидаемому значению на протяжении всего эксперимента. Это делает вполне возможным, что сам «эксперимент» никогда не проводился физически, а был основан на числах, выдуманных в воображении, чтобы соответствовать статистическим ожиданиям, но, как оказалось, слишком хорошо. [ 7 ]
Голландский научный журналист Ханс ван Маанен, однако, утверждает, что статью Лаццарини никогда не предполагалось воспринимать слишком серьезно, поскольку для читателей журнала (ориентированного на школьных учителей) было совершенно очевидно, что аппарат, который, по словам Лаццарини, построил, не может возможно, работает так, как описано. [ 8 ]
Расширение Лапласа (короткая игольница)
[ редактировать ]Теперь рассмотрим случай, когда плоскость содержит два набора параллельных линий, ортогональных друг другу, образующих стандартную перпендикулярную сетку. Наша цель — найти вероятность того, что игла пересечет хотя бы одну линию сетки. Пусть a и b — стороны прямоугольника, в котором находится середина иглы длиной l . Поскольку это короткая игольница, l < a , l < b . Пусть ( x , y ) обозначают координаты средней точки иглы, а φ обозначает угол, образованный иглой и осью x . Подобно примерам, описанным выше, мы рассматриваем x , y , φ как независимые однородные случайные величины в диапазонах 0 ≤ x ≤ a , 0 ≤ y ≤ b , − π / 2 ≤ φ ≤ π / 2 .
Чтобы решить такую задачу, мы сначала вычисляем вероятность того, что игла не пересечет ни одной линии, а затем принимаем ее дополнение. Мы вычисляем эту первую вероятность, определяя объем области, в которой игла не пересекает ни одной линии, а затем делим ее на объем всех возможностей V . Легко видеть, что V = πab .
Теперь пусть V * — объём возможностей, при котором игла не пересекает ни одной линии. Разработан СП Успенским , [ 9 ]
где F ( φ ) — область, в которой игла не пересекает ни одну линию под углом φ . Чтобы определить F ( φ ) , давайте сначала рассмотрим случай горизонтальных краев ограничивающего прямоугольника. Общая длина стороны равна a , а средняя точка не должна находиться в пределах l / 2 cos φ любой конечной точки ребра. Таким образом, общая допустимая длина без пересечений равна a − 2( л / 2 потому что φ ) или просто а - л потому что φ . Эквивалентно, для вертикальных ребер длиной b имеем b ± l sin φ . ± учитывает случаи, когда φ положителен или отрицателен. Беря положительный случай и затем добавляя знаки абсолютных значений в окончательный ответ для общности, мы получаем
Теперь мы можем вычислить следующий интеграл:
Таким образом, вероятность того, что игла не пересечет ни одну прямую, равна
И, наконец, если мы хотим вычислить вероятность P того, что игла пересечет хотя бы одну линию, нам нужно вычесть приведенный выше результат из 1, чтобы вычислить его дополнение, что даст
- .
Сравнение оценок π
[ редактировать ]Как упоминалось выше, эксперимент Бюффона с иглой можно использовать для оценки π . Этот факт справедлив и для расширения Лапласа, поскольку π также присутствует в этом ответе. Тогда естественным образом возникает следующий вопрос, который обсуждался Э. Ф. Шустером в 1974 г. [ 10 ] Какой эксперимент Бюффона или Лапласа лучше позволяет оценить значение π ? Поскольку в расширении Лапласа имеется два набора параллельных прямых, мы сравниваем N капель при наличии сетки (Лаплас) и 2 N капель в исходном эксперименте Бюффона.
Пусть A — это событие, когда игла пересекает горизонтальную линию (параллельную оси x ).
и пусть B — это событие, когда игла пересекает вертикальную линию (параллельную оси y ).
Для простоты дальнейшей алгебраической формулировки пусть a = b = t = 2 l такое, что исходный результат в задаче Бюффона равен P ( A ) = P ( B ) = 1 / π . Далее, пусть N = 100 капель.
Теперь давайте рассмотрим P ( AB ) для результата Лапласа, то есть вероятности пересечения иглой как горизонтальной, так и вертикальной линии. Мы знаем, что
Из приведенного выше раздела P ( A ′ B ′) , или вероятность того, что игла не пересечет ни одной линии, равна
Мы можем найти P ( A ′ B ) и P ( AB ′), используя следующий метод:
Решив для P ( A ′ B ) и P ( AB ′ ) и подставив это в исходное определение P ( AB ) несколькими строками выше, мы получим
Хотя это и не обязательно для решения проблемы, теперь можно увидеть, что P ( A ′ B ) = P ( AB ′) = 3 / 4 π . Имея приведенные выше значения, мы теперь можем определить, какая из этих оценок является лучшей оценкой для π . Для варианта Лапласа пусть p̂ будет оценкой вероятности того, что существует пересечение линий такое, что
- .
Нас интересует дисперсия такой оценки, чтобы понять ее полезность или эффективность. Чтобы вычислить дисперсию p̂ , мы сначала вычисляем Var( x n + y n ) , где
Решение для каждой части индивидуально,
Из предыдущего раздела мы знаем, что
уступчивость
Таким образом,
Возвращаясь к исходной задаче этого раздела, дисперсия оценки p̂ равна
Теперь давайте посчитаем количество капель M , необходимое для достижения той же дисперсии, что и 100 капель по перпендикулярным линиям. Если M < 200 , то можно заключить, что установка только с параллельными линиями более эффективна, чем схема с перпендикулярными линиями. И наоборот, если M равно или больше 200, то эксперимент Бюффона соответственно одинаково или менее эффективен. Пусть q будет оценкой исходного эксперимента Бюффона. Затем,
и
Решение для M ,
Таким образом, для того чтобы иметь ту же уверенность, что и 100 капель в случае Лапласа, требуется 222 капли только с параллельными линиями. На самом деле это неудивительно, учитывая наблюдение, что Cov( x n , y n ) < 0 . Поскольку x n и y n являются отрицательно коррелированными случайными величинами, они уменьшают общую дисперсию оценки, которая является средним значением двух из них. Этот метод уменьшения дисперсии известен как метод антитетических вариаций .
См. также
[ редактировать ]Ссылки
[ редактировать ]- ^ История акад. Рой. принадлежащий. Наука (1733), 43–45; Естественная, всеобщая и частная история Приложение 4 (1777), с. 46.
- ^ Сенета, Евгений ; Паршалл, Карен Хангер; Йонгманс, Франсуа (2001). «Развитие геометрической вероятности в девятнадцатом веке: Дж. Дж. Сильвестр, М. В. Крофтон, Ж.-Э. Барбье и Ж. Бертран» . Архив истории точных наук . 55 (6): 501–524. дои : 10.1007/s004070100038 . ISSN 0003-9519 . JSTOR 41134124 . S2CID 124429237 .
- ^ Берендс, Эрхард. «Бюффон: Он палками бросал или нет?» (PDF) . Проверено 14 марта 2015 г.
- ^ Формулировка задачи здесь позволяет избежать необходимости работать с регулярными условными плотностями вероятности .
- ^ Айгнер, Мартин; Циглер, Гюнтер М. (2013). Доказательства из КНИГИ (2-е изд.). Springer Science & Business Media. стр. 189–192.
- ^ Лаццарини, М. (1901). «Применение теории вероятностей к экспериментальному исследованию π » числа приближения . Периодический журнал по математике для среднего образования (на итальянском языке). 4 : 140–143.
- ^ Jump up to: а б Ли Бэджер, Лаццарини «Удачное приближение числа π » , Mathematics Magazine 67, 1994, 83–91.
- ↑ Ханс ван Маанен, 'Палка Лазарини' (Палка Лазарини) , «Скипетр» 31.3, 2018.
- ^ Ю. В. Успенский, «Введение в математическую вероятность» , 1937, 255.
- ^ Э. Ф. Шустер, Эксперимент Бюффона с иглой , The American Mathematical Monthly, 1974, 29-29.
Библиография
[ редактировать ]- Бэджер, Ли (апрель 1994 г.). Лазарини «Удачное приближение числа π ». Журнал «Математика» . 67 (2). Математическая ассоциация Америки: 83–91. дои : 10.2307/2690682 . JSTOR 2690682 .
- Рамали, Дж. Ф. (октябрь 1969 г.). «Задача Бюффона о лапше». Американский математический ежемесячник . 76 (8). Математическая ассоциация Америки: 916–918. дои : 10.2307/2317945 . JSTOR 2317945 .
- Матай, AM (1999). Введение в геометрическую вероятность . Ньюарк: Гордон и Брич. п. 5. ISBN 978-90-5699-681-9 .
- Делл, Закари; Франклин, Скотт В. (сентябрь 2009 г.). «Проблема иглы Бюффона-Лапласа в трех измерениях». Журнал статистической механики: теория и эксперимент . 2009 (9): 010. Бибкод : 2009JSMTE..09..010D . дои : 10.1088/1742-5468/2009/09/P09010 . S2CID 32470555 .
- Шредер, Л. (1974). «Проблема иглы Бюффона: захватывающее применение многих математических концепций». Учитель математики , 67 (2), 183–6.
- Успенский, Джеймс Виктор. «Введение в математическую вероятность». (1937).
Внешние ссылки
[ редактировать ]
- Проблема с иглой Бюффона при разрезании узла
- Математические сюрпризы: лапша Буффона в самом разгаре
- MSTE: Игла Бюффона
- Java-апплет «Игла Буффона»
- Оценка визуализации PI (Flash)
- Игла Бюффона: веселье и основы (презентация) на слайдшере
- Анимации для моделирования иглы Бюффона, выполненные Ихуэй Се с использованием R. пакета анимации
- 3D-физическая анимация Джеффри Вентреллы
- Падилья, Тони. «Пи Пи и игла Бюффона» . Числофил . Брэйди Харан . Архивировано из оригинала 17 мая 2013 г. Проверено 9 апреля 2013 г.