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

Метод Монте-Карло является весьма популярным методом в математическом моделировании физических задач и основывается на проведении серии статистических испытаний. Метод Монте-Карло является самостоятельным методом и не относится к числу методов раздела математики, называемого “Численные методы”.

Проиллюстрируем метод на примере следующей задачи. Каким образом можно измерить площадь пруда, имея только кучу камней? Предположим, что пруд расположен на поле известной площади S.

Sdf

Рис.6

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

S пруда = S*m /n.

Рассмотрим данный метод применительно к вычислению определенного интеграла.

y S S=H(b-a)

H

Y=f(x)

a b x

Рис.7

Площадь прямоугольника, ограниченного линиями y=H, y=0, x=a, x=b, равна S. Площадь под кривой y=f(x) на интервале (a,b) равна F.

Будем генерировать n пар случайных чисел x_i, y_i (I=1, 2,…n), таких, что:

A £x_i£ b, 0 £y_i£ H. (19)

Доля m пар точек (x_i,y_i), которые удовлетворяют условию

y_i£ f(x_i), (20)

представляют собой отношение интеграла F к площади прямоугольника S. Отсюда оценка интеграла методом Монте-Карло определяется выражением

F_n=A*m/n» , (21) где A=H(b-a), m – доля точек, удовлетворяющих условию (20); n – общее количество сгенерированных пар точек (x_i,y_i), удовлетворяющих условию (15).

Другая разновидность метода Монте-Карло основывается на теореме, гласящей, что оценка определенного интеграла определяется следующей формулой:

» < f > *(b-a), (22)

где < f > - среднее значение функции f(x) на интервале (a,b).

Для вычисления < f > будем генерировать случайные числа x_i, удовлетворяющие условию a £x_i £b, и вычислять значения функции f(x_i) в этих точках. Тогда оценка одномерного интеграла будет иметь вид

» F_n=< f > (b-a)= (b-a)*1/n (23)

где x_i – случайные числа; n- число испытаний.

Данный метод называется методом выборочного среднего.

Рассмотрим вопрос точности вычисления интегралов методом Монте-Карло. Понятно, что чем больше случайных точек мы сгенерируем, тем точнее будет оценка интеграла, причем можно показать, что точное значение интеграла F лежит в некотором интервале с центром в точке F_n: FÎ(F_n-s_n, F_n+s_n).

Здесь s_n=s/n,

где s - это дисперсия; n –число испытаний.

Дисперсия s определяется по формулам:

s²=< f ²> - < f >², (24)

< f > = , (24а)

< f >² = . (24б)

Дата добавления: 2015-05-26