Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Моделирование по методу Монте-Карло





Моделирование по методу Монте-Карло определяется как процедура, в которой используются случайные числа, то есть случайные величины U (0, 1). Такая процедура предназначена для решения стохастических и детерминистических задач, в которых течение времени не играет особой роли. Следовательно, моделирование по методу Монте-Карло является скорее статическим, чем динамическим. Хотя иногда утверждают, что метод Монте-Карло применим к любомутипу моделирования, в котором используются случайные числа, лучше считать его использование более ограниченным. Название метода Монте-Карло появилось во время второй мировой войны, когда этот подход был применен к проблемам, связанным с разработкой атомной бомбы.

В качестве примера моделирования по методу Монте Карло приведём оценку интеграла

где g (x) является действительной функцией, которую нельзя интегрировать аналитически. На практике моделирование по методу Монте-Карло вряд ли будет использоваться для оценки одного единственного интеграла, так как для этого существуют более эффективные методы численного анализа. Вероятнее, что он будет применяться для решения задач на кратные интегралы с нерегулярной подынтегральной функцией.

Покажем, как детерминистическаязадача может быть решена с помощью моделирования методом Монте-Карло. Пусть Y будет случайной величиной (b - a) g (X), где X ‑ непрерывная случайная величина, распределенная равномерно в интервале [ а, b ] (плотность распределения вероятностей такой случайной величины обозначается как U [a, b ], U – от англ. uniform ‑ равномерный). Тогда математическое ожидание значения Y

где fX (x) = 1/(bа) ‑ плотность распределения вероятностей случайной величины X, равномерно распределённой в интервале [a, b ]. Таким образом, задача оценки интеграла упрощена до оценки ожидаемой величины E (Y). В частности, оценим E (Y) = I по выборочному среднему:

(n) =

где X 1, Х 2,..., Х n ‑ независимые и одинаково распределенные случайные величины U [ a, b ]. Выборочное среднее (n) можно рассматривать как площадь прямоугольника с длиной (b ‑ а) и высотой I /(b ‑ а), которое является непрерывным средним g (x) между [ а, b ]. Более того, можно показать, что Е [ (n) = 1], то есть (n) является несмещенной оценкой I, а дисперсия Var[ (n)] = Var (Y)/ n. Допустим, что дисперсия Var (Y) является конечной, из чего следует, что значение (n) будет сколь угодно близко к I для достаточно большого значения n с вероятностью, равной 1.

Для демонстрации описанной выше схемы в численном отношении предположим, что необходимо оценить интеграл

который, как показывают элементарные вычисления, может иметь значение 2. Результаты применения моделирования по методу Монте-Карло для оценки этого интеграла при разных значениях n приведены в табл. 4.2.

 

 

Таблица 4. 2. (n) при разных значениях n, полученные в результате применениямоделирования по методу Монте-Карло для оценки интеграла

n          
(n) 2, 213 1, 951 1, 948 1, 989 1, 993

В настоящее время моделирование по методу Монте-Карло широко применяется при решении определенных задач статистики, которые не поддаются аналитической обработке. Этот тип моделирования применялся для оценки критических значений или достоверности критерия проверки статистических гипотез. Примером такого применения является определение критических значений критерия Колмогорова-Смирнова для непараметрической проверки гипотезы нормальности.







Дата добавления: 2014-11-12; просмотров: 1214. Нарушение авторских прав; Мы поможем в написании вашей работы!




Композиция из абстрактных геометрических фигур Данная композиция состоит из линий, штриховки, абстрактных геометрических форм...


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


ТЕОРЕТИЧЕСКАЯ МЕХАНИКА Статика является частью теоретической механики, изучающей условия, при ко­торых тело находится под действием заданной системы сил...


Теория усилителей. Схема Основная масса современных аналоговых и аналого-цифровых электронных устройств выполняется на специализированных микросхемах...

Индекс гингивита (PMA) (Schour, Massler, 1948) Для оценки тяжести гингивита (а в последующем и ре­гистрации динамики процесса) используют папиллярно-маргинально-альвеолярный индекс (РМА)...

Методика исследования периферических лимфатических узлов. Исследование периферических лимфатических узлов производится с помощью осмотра и пальпации...

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

Признаки классификации безопасности Можно выделить следующие признаки классификации безопасности. 1. По признаку масштабности принято различать следующие относительно самостоятельные геополитические уровни и виды безопасности. 1.1. Международная безопасность (глобальная и...

Прием и регистрация больных Пути госпитализации больных в стационар могут быть различны. В цен­тральное приемное отделение больные могут быть доставлены: 1) машиной скорой медицинской помощи в случае возникновения остро­го или обострения хронического заболевания...

ПУНКЦИЯ И КАТЕТЕРИЗАЦИЯ ПОДКЛЮЧИЧНОЙ ВЕНЫ   Пункцию и катетеризацию подключичной вены обычно производит хирург или анестезиолог, иногда — специально обученный терапевт...

Studopedia.info - Студопедия - 2014-2025 год . (0.012 сек.) русская версия | украинская версия