Построение имитационной модели методом Монте-Карло

 

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

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

Обозначим P{} теоретическое распределение, для которого мы хотим генерировать последовательность случайных чисел. Для любого отрезка [a,b] по определению P{[a,b]} равно вероятности того, что случайная величина, подчиняющаяся данному распределению, попадет в отрезок [a,b].

Пусть N количество чисел в последовательности, полученной с помощью датчика случайных чисел. Обозначим PN{} соответствующее эмпирическое распределение. ( По определению, , где N[a,b] количество чисел последовательности, попавших в отрезок [a,b].) Если последовательность случайных чисел распределена в соответствии с теоретическим распределением P{} , то для любого отрезка [a,b] при достаточно большом количестве N чисел в последовательности имеет место приблизительное равенство PN{[a,b]}= P{[a,b]}. (В пределе должно выполняться строгое равенство lim PN{[a,b]}= P{[a,b]}

N®¥

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

Пусть модель задана в виде Y=f (X1, X2,…. Xm). На основании известных теоретических распределений входных параметров X1, X2,…. Xm для каждого входного параметра Xk с помощью датчика случайных чисел строится последовательность чисел Xk1, Xk2,…. XkN, подчиняющаяся соответствующему теоретическому распределению. Затем с помощью последовательностей Xk1, Xk2,…. XkN , k=1,m ,строится последовательность чисел Y1,Y2, Y N для выходного параметра Y по формуле Yj=f(X1j, X2j,…. Xmj), где j=1,N.

Полученная последовательность Y1,Y2, Y N естественным образом задает эмпирическое распределение выходного параметра Y.

 


 








Дата добавления: 2019-04-03; просмотров: 233;


Поиск по сайту:

При помощи поиска вы сможете найти нужную вам информацию.

Поделитесь с друзьями:

Если вам перенёс пользу информационный материал, или помог в учебе – поделитесь этим сайтом с друзьями и знакомыми.
helpiks.org - Хелпикс.Орг - 2014-2024 год. Материал сайта представляется для ознакомительного и учебного использования. | Поддержка
Генерация страницы за: 0.003 сек.