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

Арзамасский государственный педагогический институт

имени А.П.Гайдара

Кафедра математического анализа

Зубанов М. А., студент

3 курса очного отделения

 физико-математического

факультета

КУРСОВАЯ РАБОТА Метод Монте-Карло и его применение

Научный руководитель:

канд. тех. наук, доцент

Потехин В.А.

Арзамас-2002 г.

Содержание

 

Введение……………………………………………………………..3

Глава 1. Некоторые сведения теории вероятностей ………….5

§1. Математическое ожидание, дисперсия……………………..5

§2. Точность оценки, доверительная вероятность. Доверительный

интервал……………………………………………………….6

§3. Нормальное распределение…………………………………..6

Глава 2. Метод Монте-Карло……………………………………...8

§1. Общая схема метода Монте-Карло……………………….….8

§2. Оценка погрешности метода Монте-Карло…………………8

Глава 3. Вычисление интегралов методом Монте-Карло…….12

§1. Алгоритмы метода Монте-Карло для решения

интегральных уравнений второго рода………………….…12

§2. Способ усреднения подынтегральной функции………….…13

§3. Способ существенной выборки, использующий

«вспомогательную плотность распределения»…………… .16

§4. Способ, основанный на истолковании интеграла как

площади……………………………………………………. ..19

§5. Способ «выделения главной части»……………………… ...21

§6. Программа вычисления определенного интеграла методом

Монте-Карло…………………………………………………..23

§7. Вычисление кратных интегралов методом Монте-Карло.…25

Заключение…………………………………………………………..28

Приложение……………………………………………………….. ..29

Литература…………………………………………………………...30

Введение.

Метод Монте-Карло можно определить как метод моделирования случайных величин с целью вычисления характеристик их распределений.

Возникновение идеи использования случайных явлений в области приближённых вычислений принято относить к 1878 году, когда появилась работа Холла об определении числа p с помощью случайных бросаний иглы на разграфлённую параллельными линиями бумагу. Существо дела заключается в том, чтобы экспериментально воспроизвести событие, вероятность которого выражается через число p, и приближённо оценить эту вероятность. Отечественные работы по методу Монте-Карло появились в 1955-1956 годах. С того времени накопилась обширная библиография по методу Монте-Карло. Даже беглый просмотр названий работ позволяет сделать вывод о применимости метода Монте-Карло для решения прикладных задач из большого числа областей науки и техники.

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

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

Глава 1. Некоторые сведения теории вероятностей

 

§1. Математическое ожидание, дисперсия.

Дискретной называют случайную величину, которая принимает отдельные, изолированные возможные значения с определёнными вероятностями. Число возможных значений дискретной случайной величины может быть конечным или бесконечным.

Математическим ожиданием дискретной случайной величины называют сумму произведений всех её возможных значений на их вероятность.

,

где Х – случайная величина,  - значения, вероятности которых соответственно равны .

Математическое ожидание приближённо равно (тем точнее, чем больше число испытаний) среднему арифметическому наблюдаемых значений случайной величины.

Дисперсией (рассеянием) случайной величины называют математическое ожидание квадрата отклонения случайной величины от её математического ожидания: .

Средним квадратичным отклонением случайной величины Х называют квадратный корень из дисперсии: .

§2. Точность оценки, доверительная вероятность. Доверительный интервал.

Точечной называют оценку, которая определяется одним числом.

Интервальной называют оценку, которая определяется двумя числами – концами интервала. Интервальные оценки позволяют установить точность и надёжность оценок.

Пусть, найденная по данным выборки, статистическая характеристика  служит оценкой неизвестного параметра . Ясно, что  тем точнее определяет параметр , чем меньше абсолютная величина разности . Другими словами, если d>0 и , то , чем меньше d, тем оценка точнее. Положительное число d характеризует точность оценки.

Надёжностью (доверительной вероятностью) оценки  по  называют вероятность g, с которой осуществляется неравенство .

Доверительным называют интервал , который покрывает неизвестный параметр с заданной надёжностью g.

§3. Нормальное распределение.

Нормальным называют распределение вероятностей непрерывной

случайной величины, которое описывается дифференциальной функцией

.

а - математическое ожидание, s - среднее квадратичное отклонение нормального распределения.

 

Глава 2. Метод Монте-Карло

§1. Общая схема метода Монте-Карло.

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

Практически же поступают так: производят n испытаний, в результате которых получают n возможных значений Х; вычисляют их среднее арифметическое  и принимают x в качестве оценки (приближённого значения) a^* искомого числа a:

.

Поскольку метод Монте-Карло требует проведения большого числа испытаний, его часто называют методом статистических испытаний. Теория этого метода указывает, как наиболее целесообразно выбрать случайную величину Х, как найти её возможные значения. В частности, разрабатываются способы уменьшения дисперсии используемых случайных величин, в результате чего уменьшается ошибка, допускаемая при замене искомого математического ожидания а его оценкой а^*.

§2. Оценка погрешности метода Монте-Карло.

 Пусть для получения оценки a^* математического ожидания а случайной величины Х было произведено n независимых испытаний (разыграно n возможных значений Х) и по ним была найдена выборочная средняя , которая принята в качестве искомой оценки: . Ясно, что если повторить опыт, то будут получены другие возможные значения Х, следовательно, другая средняя, а значит, и другая оценка a^*. Уже отсюда следует, что получить точную оценку математического ожидания невозможно. Естественно возникает вопрос о величине допускаемой ошибки. Ограничимся отысканием лишь верхней границы d допускаемой ошибки с заданной вероятностью (надёжностью) g: .

Интересующая нас верхняя грань ошибки d есть не что иное, как «точность оценки» математического ожидания по выборочной средней при помощи доверительных интервалов. Рассмотрим следующие три случая.