Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология
Образование Политология Производство Психология Стандартизация Технологии


Сущность метода статистических испытаний (метода Монте-Карло)



 

Датой рож­дения метода Монте-Карло принято считать 1949 г., когда появилась статья под названием «The Monte Carlo method». Создателями этого метода считают амери­канских математиков Дж. Неймана и С. Улама. В СССР первые статьи о методе Монте-Карло были опублико­ваны в 1955—1956гг.

Любопытно, что теоретическая основа метода была известна давно. Более того, некоторые задачи стати­стики рассчитывались иногда с помощью случайных вы­борок, т. е. фактически методом Монте-Карло. Однако до появления электронных вычислительных машин (ЭВМ) этот метод не мог найти сколько-нибудь широкого применения, ибо моделировать случайные величины вручную—очень трудоемкая работа. Таким образом, возникновение метода Монте-Карло как весьма универ­сального численного метода стало возможным только благодаря появлению ЭВМ. Само название «Монте-Карло» происходит от города Монте-Карло в княжестве Монако, знаменитого своим игорным домом.

Идея метода чрезвычайно проста и состоит она в следующем. Вместо того, чтобы описывать процесс с помощью аналитического аппарата (дифференциальных или алгебраических уравнений), производится «розыгрыш» случайного явления с помощью специально организованной процедуры, включающей в себя случайность и дающей случайный результат. В действительности конкретное осуществление случайного процесса складывается каждый раз по-иному; так же и в результате статистического моделирования мы получаем каждый раз новую, отличную от других реализацию исследуемого процесса. Что она может нам дать? Сама по себе ничего, так же как, скажем, один случай излечения больного с помощью какого-либо лекарства. Другое дело, если таких реализаций получено много. Это множество реализаций можно использовать, как некий искусственно полученный статистический материал, который может быть обработан обычными методами математической статистики. После такой обработки могут быть получены любые интересующие нас характеристики: вероятности событий, математические ожидания и дисперсии случайных величин и т. д. При моделировании случайных явлений методом Монте-Карло мы пользуемся самой случайностью как аппаратом исследования, заставляем ее «работать на нас».

Нередко такой прием оказывается проще, чем по­пытки построить аналитическую модель. Для сложных операций, в которых участвует большое число элемен­тов (машин, людей, организаций, подсобных средств), в которых случайные факторы сложно переплетены, где процесс — явно немарковский, метод статистиче­ского моделирования, как правило, оказывается проще аналитического (а нередко бывает и единственно воз­можным).

В сущности, методом Монте-Карло может быть ре­шена любая вероятностная задача, но оправданным он становится только тогда, когда процедура розыгрыша проще, а не сложнее аналитического расчета. Приведем пример, когда метод Монте-Карло возможен, но край­не неразумен. Пусть, например, по какой-то цели производится три независимых выстрела, из которых каж­дый попадает в цель с вероятностью 1/2. Требуется найти вероятность хотя бы одного попадания. Элементарный расчет дает нам вероятность хотя бы одного попадания 1 — (1/2)3 = 7/8. Ту же задачу можно решить и «розыгрышем», статистическим моделированием. Вместо «трех выстрелов» будем бросать «три монеты», считая, скажем, герб—за попадание, решку — за «промах». Опытсчитается«удачным», если хотя бы на одной из монетвыпадет герб. Произведем очень-очень много опытов, подсчитаем общее количество «удач» и разделим на число N произведенных опытов. Таким образом, мы получим частоту события, а она при большом числе опытов близка к вероятности. Применить такой прием мог бы разве человек, вовсе не знающий теории вероятностей, тем не менее, в принципе, он возможен.

Метод Монте-Карло - это численный метод решения математических задач при помощи моделирования случайных величин.

Пример 1. Предположим, что нам нужно вычислить площадь плоской фигуры S. Это может быть произвольная фигура с криволинейной границей, заданная графически или аналитически, связная или состоящая из нескольких кусков. Пусть это будет фигура изображенная на рис. 1, и предположим, что она вся расположена внутри единичного квадрата.

Выберем внутри квадрата N случайных точек. Обозначим через F число точек, попавших при этом внутрь S. Геометрически очевидно, что площадь S приближенно равна отношению F/N. Чем больше N, тем больше точность этой оценки.

Две особенности метода Монте-Карло.

Первая особенность метода - простая структура вычислительного алгоритма.

Вторая осо­бенность метода - погрешность вычислений, как правило, пропорциональна D/N2, где D - неко­торая постоянная, N - число испытаний. Отсюда видно, что для того, чтобы уменьшить по­грешность в 10 раз (иначе говоря, чтобы получить в ответе еще один верный десятичный знак), нужно увеличить N (т. е. объем работы) в 100 раз. Ясно, что добиться высокой точности таким путем невозможно. Поэтому обычно говорят, что метод Монте-Карло особенно эффективен при решении тех задач, в которых результат ну­жен с небольшой точностью (5-10%). Способ применения метода Монте-Карло по идее доволь­но прост. Чтобы получить искусственную случайную выборку из совокупности величин, опи­сываемой некоторой функцией распределения вероятностей, следует:

1. Построить график или таблицу интегральной функции распределения на основе ряда чи­сел, отражающего исследуемый процесс (а не на основе ряда случайных чисел), причем значе­ния случайной переменной процесса откладываются по оси абсцисс (х), а значения вероятности (от 0 до 1) - по оси ординат (у).

2. С помощью генератора случайных чисел выбрать случайное десятичное число в преде­лах от 0 до 1 (с требуемым числом разрядов).

3. Провести горизонтальную прямую от точки на оси ординат соответствующей выбран­ному случайному числу, до пересечения с кривой распределения вероятностей.

4. Опустить из этой точки пересечения перпендикуляр на ось абсцисс.

5.Записать полученное значение х. Далее оно принимается как выборочное значение.

6. Повторить шаги 2-5 для всех требуемых случайных переменных, следуя тому порядку, в котором они были записаны.

В задачах исследования операций метод Монте-Карло применяется в трех основных направлениях:

1) при моделировании сложных, комплексных операций, где присутствует много взаимодействующих случайных факторов;

2) при проверке применимости более простых, аналитических методов и выяснении условий их применимости;

3) в целях выработки поправок к аналитическим формулам типа «эмпирических формул» в технике.

Сущность метода Монте-Карло состоит в следующем: требуется найти значение а некоторой изучаемой величины. С этой целью выбирают такую случайную величину Х, математическое ожидание которой равно а: М(Х)=а.

Практически же поступают так: вычисляют (разыгрывают) n возможных значений xi случайной величины Х, находят их среднее арифметическое

и принимают в качестве оценки (приближенного значения) а* искомого числа а: a~a*= . Таким образом, для применения метода Монте-Карло необходимо уметь разыгрывать случайную величину.

Пусть требуется разыграть дискретную случайную величину Х, т.е. вычислить последовательность ее возможных значений хi (i=1, 2, …), зная закон распределения Х. Введем обозначения: R- непрерывная случайная величина, распределенная равномерно в интервале (0, 1); ri (j=1, 2, …) – случайные числа (возможные значения R).

Правило: Для того чтобы разыграть дискретную случайную величину Х, заданную законом распределения

Х х1 х2 … хn

P p1 p2 … pn

надо:

1. разбить интервал (0, 1) оси or на n частичных интервалов:

Δ 1=(0; р1), Δ 2 =1; р1+р2), …, Δ n =12+…+рn-1; 1).

2. выбрать случайное число rj.

Если rj попало в частичный интервал Δ i, то разыгрываемая величина приняла возможное значение хi.


Поделиться:



Популярное:

Последнее изменение этой страницы: 2016-03-15; Просмотров: 2144; Нарушение авторского права страницы


lektsia.com 2007 - 2024 год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! (0.015 с.)
Главная | Случайная страница | Обратная связь