М_М_К_3

2.4. Вычисление значения числа 

Рассмотрим задачу вычисления значения числа . Эта задача интересна и с дидактической точки зрения, так как помогает совместить простые опыты с решением довольно сложной математической модели по определению числа .

1. Способ дождя. Идея определения числа  заключается в следующем: при радиусе круга r = 1, его площадь S равна . Рассмотрим четверть круга (рис.2.2) , тогда

. (2.35)

Если мы определим площадь четверти круга, то определим значение числа . Для определения площади круга используем метод Монте-Карло, который требует применения случайных чисел. Экспериментально случайные числа можно получить как при помощи рулетки, так и при помощи дождя. Для опыта приготовим кусок картона, нарисуем на нем квадрат и впишем в квадрат четверть круга.

Если такой чертеж некоторое время подержать под дождем, то на его поверхности останутся следы капель. Подсчитаем число следов внутри квадрата и внутри четверти круга. Очевидно, что их отношение будет приближенно равно отношению площадей этих фигур, так как попадание капель в различные места чертежа равновероятно, т.е. ,т.к. , то .

Таким образом, значение числа  можно найти, определив N_kr – число капель в кругу и N_kv– число капель в квадрате

. (2.36)

Фактически значение площади подынтегральной функции мы нашли, используя метод «зонтика» Неймана.

Генерация случайных чисел на отрезке. Дождь можно заменить таблицей случайных чисел, которая составляется с помощью компьютера по специальной программе (можно воспользоваться одной из многочисленных публикаций таких таблиц, например, в книгах[5,6,7,23,25]). Каждому следу капли поставим в соответствие два случайных числа, характеризующих его положение вдоль осей Ох и Оу (см. рис.2.2).

Случайные числа можно выбрать из таблицы в любом порядке, например, подряд. Пусть первое четырехзначное число в таблице 3265. Из него можно "приготовить" пару чисел, каждое из которых больше нуля и меньше единицы: х=0.32, у=0.65. Эти числа будем считать координатами капли, т. е. капля как будто попала в точку (0.32; 0.65).

Решение. Аналогично поступаем и со всеми выбранными случайными числами. Если окажется, что для точки (x_i; y_i) выполняется неравенство x_i^2 + y_i^2> 1, то, значит, она лежит вне круга. Если x_i^2 +y_i^2<=1, то точка лежит внутри круга.

Для подсчета значения  снова воспользуемся формулой (2.36). Ошибка вычислений по этому методу, как правило, пропорциональна , где С - некоторая постоянная, а N - число испытаний. В нашем случае N = N _kv. Из этой формулы видно: для того чтобы уменьшить ошибку в 10 раз, нужно увеличить N в 100 раз. Ясно, что широкое применение метода статистических испытаний (Монте-Карло) стало возможным только благодаря компьютерам. Программа 1 и 2 в приложении 3 реализуют описанный метод.

2. Способ "падающей иголки" или игла Бюффона. Идея метода также заключается в вычислении площади под некоторой кривой. В методе иглы Бюффона в качестве подынтегральной кривой используется тригонометрическая функция COS. Сущность метода заключается в том, чтобы экспериментально воспроизвести событие, вероятность которого выражается через число , и приближённо оценить эту вероятность.

Построение модели. Возьмем обыкновенную швейную иголку и лист бумаги. На листе проведем несколько параллельных, прямых линий так, чтобы расстояния между ними были равны а и не превышали длину иголки L.

Чертеж должен быть достаточно большим, чтобы случайно брошенная игла не упала за его пределами. Положение случайным образом брошенной на чертеж иглы (см. рис. 2.4) определяется расстоянием х от ее середины до ближайшей прямой, углом, который игла образует с перпендикуляром, опущенным из середины иглы на ближайшую прямую (см. рис. 2.3). Ясно, что 0<=x<= 1/2, -/2<= у <=/2.

На рис.2.5 изобразим графически функцию у=0.5*l*cos().

Всевозможные расположения иглы характеризуются точками с координатами (;y), расположенными на участке ABCD. Заштрихованный участок AED - это точки, которые соответствуют случаю пересечения иглы с прямыми линиями. Вероятность события - "игла пересекла прямую" - вычисляется по формуле

здесь , ,

т. е. . (2.37)

Рис. 2.3 Рис.2.4

Вероятность P(a) также можно приблизительно определить многократным бросанием иглы. Пусть иглу бросали на чертеж N раз и она упала k раз, пересекая одну из прямых, тогда при достаточно большом N имеем

. (2.38)

Приравнивая формулы (2.35) и (2.36) получаем формулу для расчета значения  (пи).

. (2.39)

При вычислении тригонометрических функций на языке Паскаль и Бейсике углы из градусной меры приходится переводить в радианную меру, где используется значение . Поэтому для построения нашей модели мы используем проекцию иглы (рис.2.4). Координаты проекции иглы на ось ОУ не влияют на наши расчеты. Координату Х1 зададим датчиком случайных чисел. Координата Х2 может изменяться в промежутке Х1<=X2<=X1+L, L - длина иглы. Условие пересечения иглы прямой Х1<=A<=X2, где а - координата х проекции прямой на ОХ. Для конкретности возьмем расстояние между прямыми – а = 2, длину иглы – L =1, количество линий на рис. 2.2. Описанная модель реализована в программах 3 и 4 (см. приложение 3.).

Задания на моделирование:

Составить программы определения числа  по рассмотренным методам.
Провести моделирование.
Оценить погрешность результатов и построить графики .

Содержание