Случайные события

Размещено на http://www.allbest.ru/

1. Поток случайных событий

Мы научились имитировать наступление случайных событий. То есть мы можем разыграть - какое из возможных событий наступит и в каком количестве. Чтобы это определить, надо знать статистические характеристики появления событий, например, такой величиной может быть вероятность появления события, или распределение вероятностей разных событий, если типов этих событий бесконечно много.

Но часто еще важно знать, когда конкретно наступит то или иное событие во времени.

Когда событий много и они следуют друг за другом, то они образуют поток. Заметим, что события при этом должны быть однородными, то есть похожими чем-то друг на друга. Например, появление водителей на АЗС, желающих заправить свой автомобиль. То есть, однородные события образуют некий ряд. При этом считается, что статистическая характеристика этого явления (интенсивность потока событий) задана. Интенсивность потока событий указывает, сколько в среднем происходит таких событий за единицу времени. Но когда именно произойдет каждое конкретное событие надо определить методами моделирования. Важно, что, когда мы сгенерируем, например, за 200 часов 1000 событий, их количество будет равно примерно величине средней интенсивности появления событий 1000/200 = 5 событий в час, что является статистической величиной, характеризующей этот поток в целом.

Интенсивность потока в некотором смысле является математическим ожиданием количества событий в единицу времени. Но реально может так оказаться, что в один час появится 4 события, в другой - 6, хотя в среднем получается 5 событий в час, поэтому одной величины для характеристики потока недостаточно. Второй величиной, характеризующей насколько велик разброс событий относительно математического ожидания, является, как и ранее, дисперсия. Собственно именно эта величина определяет случайность появления события, слабую предсказуемость момента его появления. Про эту величину мы расскажем в следующей лекции.

Итак.

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

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

Интенсивность потока л - это среднее число событий в единицу времени. Интенсивность потока можно рассчитать экспериментально по формуле:

л = N/T_н,

где N - число событий, произошедших за время наблюдения T_н.

Если интервал между событиями ф_j равен константе или определен какой-либо формулой в виде: t_j = f(t_j - 1), то поток называется детерминированным. Иначе поток называется случайным.

Случайные потоки бывают:

· ординарные: вероятность одновременного появления двух и более событий равна нулю;

· стационарные: частота появления событий л(t) = const(t);

· без последействия: вероятность появления случайного события не зависит от момента совершения предыдущих событий.

2 Пуассоновский поток

За эталон потока в моделировании принято брать пуассоновский поток.

Пуассоновский поток - это ординарный поток без последействия.

Как ранее было указано, вероятность того, что за интервал времени (t₀, t₀ + ф) произойдет m событий, определяется из закона Пуассона:

где a - параметр Пуассона.

Если л(t) = const(t), то это стационарный поток Пуассона (простейший). В этом случае a = л · t. Если л = var(t), то это нестационарный поток Пуассона.

Для простейшего потока вероятность появления m событий за время ф равна:

Вероятность непоявления (то есть ни одного, m = 0) события за время ф равна:

Очевидно, что чем больше время наблюдения, тем вероятность непоявления ни одного события меньше. Кроме того, чем более значение л, тем круче идет график, то есть быстрее убывает вероятность. Это соответствует тому, что если интенсивность появления событий велика, то вероятность непоявления события быстро уменьшается со временем наблюдения.

Вероятность появления хотя бы одного события (P_ХБ1С) вычисляется так:

так как P_ХБ1С + P₀ = 1 (либо появится хотя бы одно событие, либо не появится ни одного, - другого не дано).

Вероятность появления хотя бы одного события стремится со временем к единице, то есть при соответствующем длительном наблюдении события таковое обязательно рано или поздно произойдет. Чем дольше мы наблюдаем за событием (чем более t), тем больше вероятность того, что событие произойдет - график функции монотонно возрастает.

Чем больше интенсивность появления события (чем больше л), тем быстрее наступает это событие, и тем быстрее функция стремится к единице. На графике параметр л представлен крутизной линии (наклон касательной).

Если увеличивать л, то при наблюдении за событием в течение одного и того же времени ф, вероятность наступления события возрастает. Очевидно, что график исходит из 0, так как если время наблюдения бесконечно мало, то вероятность того, что событие произойдет за это время, ничтожна. И наоборот, если время наблюдения бесконечно велико, то событие обязательно произойдет хотя бы один раз, значит, график стремится к значению вероятности равной 1.

Изучая закон, можно определить, что: m_x = 1/л, у = 1/л, то есть для простейшего потока m_x = у. Равенство математического ожидания среднеквадратичному отклонению означает, что данный поток - поток без последействия. Дисперсия (точнее, среднеквадратичное отклонение) такого потока велика. Физически это означает, что время появления события (расстояние между событиями) плохо предсказуемо, случайно, находится в интервале

m_x - у < ф_j < m_x + у.

Хотя ясно, что в среднем оно примерно равно:

ф_j = m_x = T_н/N.

Событие может появиться в любой момент времени, но в пределах разброса этого момента ф_j относительно m_x на [-у; +у] (величину последействия).

По смыслу P равно r, поэтому, выражая ф из формулы, окончательно для определения интервалов между двумя случайными событиями имеем:

ф = -1/л · Ln(r),

где r - равномерно распределенное от 0 до 1 случайное число, которое берут из ГСЧ, ф - интервал между случайными событиями (случайная величина ф_j).

Рассмотрим поток изделий, приходящих на технологическую операцию. Изделия приходят случайным образом - в среднем восемь штук за сутки (интенсивность потока л = 8/24 [ед/час]). Необходимо промоделировать этот процесс в течение T_н = 100 часов. m = 1/л = 24/8 = 3, то есть в среднем одна деталь за три часа. Заметим, что у = 3.

3. Моделирование неординарных потоков событий

Если известно, что поток не является ординарным, то необходимо моделировать кроме момента возникновения события еще и число событий, которое могло появиться в этот момент. Например, вагоны на железнодорожную станцию прибывают в составе поезда в случайные моменты времени (ординарный поток поездов). Но при этом в составе поезда может быть разное (случайное) количество вагонов. В этом случае о потоке вагонов говорят как о потоке неординарных событий.

Допустим, что M_k = 10, у = 4 (то есть, в среднем в 68 случаях из 100 приходит от 6 до 14 вагонов в составе поезда) и их число распределено по нормальному закону. В место, отмеченное в предыдущем алгоритме, нужно вставить фрагмент,

Очень полезным в производстве является решение следующей задачи. Каково среднее время суточного простоя оборудования технологического узла, если узел обрабатывает каждое изделие случайное время, заданное интенсивностью потока случайных событий л₂? При этом экспериментально установлено, что привозят изделия на обработку тоже в случайные моменты времени, заданные потоком л₁ партиями по 8 штук, причем размер партии колеблется случайно по нормальному закону сm = 8, у = 2 (см. лекцию 25). До начала моделирования T = 0 на складе изделий не было. Необходимо промоделировать этот процесс в течение T_н = 100 часов.

Отмечая для обрабатывающего узла времена, когда он простаивал в ожидании очередной детали мы можем посчитать суммарное время простоев узла за все время наблюдения, а затем рассчитать среднее время простоя в течение суток. Для данной реализации это время вычисляется так:

T_пр.^ср. = 24 · (t₁^пр. + t₂^пр. + t₃^пр. + t₄^пр. + … + t_N^пр.)/T_н.

4. Моделирование нестационарных потоков событий

случайный пуассоновский неординарный простой

В ряде случаев интенсивность потока может меняться со временем л(t). Такой поток называется нестационарным. Например, среднее количество за час машин скорой помощи, покидающих станцию по вызовам населения большого города, в течение суток может быть различным. Известно, например, что наибольшее количество вызовов падает на интервалы с 23 до 01 часа ночи и с 05 до 07 утра, тогда как в остальные часы оно вдвое меньше.

Размещено на Allbest.ru