Методы Монте-Карло

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

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

. (5.48)

Погрешность вычисления будет уменьшаться с ростом числа используемых узлов расчета функции по закону . Графическая иллюстрация метода представлена на рисунке 5.5

Рисунок 5.5. Статистический вариант метода прямоугольников

Формула (5.48) обобщается на случай кратных интегралов

(5.49)

Здесь – -мерный объем области интегрирования . Число узлов, в которых необходимо вычислять подынтегральную функцию, будет пропорционально .

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

,(5.50)

где находится в интервале . Тогда две случайные величины и можно рассматривать как координаты точек в единичном квадрате (рисунок 5.8). При равномерном распределении точек в квадрате за приближенное значение интеграла принимается отношение количества точек , попавших под кривую , к общему числу испытаний

.(5.51)

Этот алгоритм также обобщается на кратные интегралы.

Дата добавления: 2017-01-14