3. Моделирование случайной величины дискретного типа
А. Общий алгоритм моделирования.
Если случайная величина дискретная, то ее моделирование можно свести к моделированию независимых испытаний. В самом деле, пусть имеет место следующий ряд распределения:
| | | … | |
| | | … | |
Обозначим через событие, состоящее в том, что случайная величина
примет значение
, при этом
. Тогда нахождение значения, принятого случайной величиной
в результате испытания, сводится к определению того, какое из событий
появится. Так как события
несовместны и вероятность появления каждого из них не изменяется от испытания к испытанию, то для определения последовательности значений, принятых случайной величиной
можно использовать алгоритм моделирования последовательности независимых испытаний.
Б. Моделирование случайной величины с биномиальным распределением.
Случайная величина считается распределенной по биномиальному закону, если
где ;
— вероятность появления некоторого события
в каждом отдельно взятом испытании;
— вероятность появления события
в
независимых испытаниях
раз.
Введем случайную величину — число появлений события
в
-ом испытании,
Для этой величины имеет место:
,
. (1)
Тогда случайное число появлений события
в
испытаниях определяется по формуле
. (2)
Исходя из формул (1) и (2), значения случайной величины определяются следующим образом:
1) находят последовательность значений случайной величины
2) для каждого числа ,
проверяют, выполняется ли неравенство
если неравенство выполняется, то полагают
в противном случае считают
3) находят сумму значений случайных величин
которая совпадает со значением
Повторяя этот алгоритм, получим последовательность значений случайной величины с биномиальным законом распределения.
В. Моделирование случайной величины, распределенной по закону Пуассона.
Распределением Пуассона называется распределение вероятностей дискретной случайной величины, задаваемое формулой:
,
,
где — число событий простейшего потока, наступающих за некоторый промежуток времени. Распределение Пуассона применяется вместо биномиального распределения тогда, когда число
независимых испытаний велико (порядка нескольких сотен), а вероятность
появления события в каждом отдельно взятом испытании мала, при этом желательно, чтобы имело место
.
Алгоритм моделирования случайной величины , распределенной по закону Пуассона при заданном параметре
можно представить следующим образом:
1) выбираем таким образом, чтобы вероятность
была достаточно малой, например, меньше 0, 01;
2) получаем последовательность значений случайной величины
, равномерно распределенной на отрезке
;
3) для каждого числа ,
проверяем, выполняется ли неравенство
; если это неравенство выполняется, то полагают
, в противном случае считаем
;
4) вычисляем сумму которая совпадает со значением случайной величины
распределенной по закону Пуассона.
... 2. Таблица 2. Устранение недостатков анализа чувствительности и анализа сценариев при использовании для риск-анализа метода Монте-Карло Метод Недостаток Решение с помощью имитационного моделирования Анализ чувствительности Не учитывается наличие корреляции между различными составляющими проекта Корреляция моделируется различными методами и учитывается в модели Рассматривается влияние ...
... порыв ветра, интенсивность горения пороха, влажность и т.д. Из всех случайных величин рассмотрим только сопротивление воздуха. Учитывая то, что сопротивление воздуха направленно в противоположную сторону движения тела [4], то силу сопротивления можно представить в виде двух составляющих: вертикальной и горизонтальной, при этом, чем выше скорость, тем больше сопротивление и наоборот. Выделим ...
... как точки на временной оси. Для достижения основной цели моделирования достаточно наблюдать систему в моменты реализации основных событий. Рассмотрим пример одноканальной системы массового обслуживания. Целью имитационного моделирования подобной системы является определение оценок ее основных характеристик, таких, как среднее время пребывания заявки в очереди, средняя длина очереди и доля ...
... тех же условий анализа. Так, в точке эквивалентности при титровании 0,01 М AgNO3 глицином величина рAg оказалaсь: по потенциометрическим данным - 4,2, по результатам компьютерного моделирования - 3,5; соответственно при титровании 0,0010 М AgNO3 комплексоном III - 5,6 и 4,9. Невелики (менее 10%) и расхождения по объему титранта, соответствующему точке эквивалентности. Отмеченные расхождения можно ...
0 комментариев