В пуассоновском потоке событий (стационарном и нестационарном) число событий потока , попадающих на любой участок, распределено по закону Пуассона  

Таким образом, для исследуемой системы S с дискретными состояниями и непрерывным временем переходы из состояния в состояние происходят под действием пуассоновских потоков событий с определенной интенсивностью Я.  

Представим автомобиль как некоторую систему S с дискретными состояниями iSj,. 2. .... Sn, которая переходит из состояния S/ в состояние Sj(i - 1, 2,. .., n,j = I, 2,. .., и) под воздействием пуассоновских потоков событий (отказов) с интенсивностями Хд. Будем рассматривать следующие состояния автомобиля, в которых он может находиться в процессе эксплуатации и которые характеризуются целодневными простоями  

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

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

Между пуассоновскими потоками событий и дискретными марковскими процессами с непрерывным временем имеется тесная связь.  

Связь пуассоновских потоков событий с дискретными марковскими процессами с непрерывным временем  

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

Можно считать, что события, переводящие автомобиль из состояния в состояние, представляют собой потоки событий (например, потоки отказов). Если все потоки событий , переводящие систему (автомобиль) из состояния в состояние, пуассоновские (стационарные или нестационарные), то процесс, протекающий в системе, будет марковским, а плотности вероятности перехода Ху в непрерывной цепи Маркова представляют собой интенсивности потока событий, переводящего систему из состояния Si в состояние Sj. Например, Х03 - интенсивность потока отказов автомобиля, который переводит автомобиль из состояния исправен, работает в состояние находится в ТР.  

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

Пуассоновский стационарный (простейший) поток событий  

Пуассоновский стационарным (простейшим) поток событий  

Пуассоновский нестационарный поток событий  

Рассмотрим нестационарный пуассоновский поток с интенсивностью Mf), некоторый промежуток времени длиной г>0, начинающийся с момента t0 (и заканчивающийся, следовательно, в момент +г) и дискретную случайную величину Х р г) - число событий, наступающих в потоке за промежуток времени от ta до t0+r.  

Определение 6.2. Элементом вероятности появления события в нестационарном пуассоновском потоке называется вероятность >,(АО появления события за элементарный (достаточно малый) промежуток времени от t0 до t0+bt.  

Теорема 6.2. Для элемента вероятности появления события за элементарный промежуток времени от t0 до t0+Af в нестационарном пуассоновском потоке с интенсивностью A(t) имеет место приближенная формула  

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

Одной из основных стохастических характеристик нестационарного пуассоновского потока является дискретная случайная величина X(t т), представляющая собой случайное число событий, наступающих в потоке за промежуток [ t.+t.  

Другой основной стохастической характеристикой нестационарного пуассоновского потока является случайный интервал времени T(tB) между двумя соседними событиями, первое из которых наступило в момент t0.  

Доказательство Вероятность p (t At) того, что система S, находившаяся в момент времени t в состоянии sp за промежуток времени от t до t+Ы перейдет из него в состояние s (см. 4) равна элементу вероятности pfa t) появления события в пуассоновском потоке П.. на элементарном участке от t до +Д (см. Определение 5.11). Но (см. (4.3))  

Система, в которой протекает дискретный марковский процесс с непрерывным временем, перескакивает из одного состояния х в другое xj не самопроизвольно, а под воздействием определенного события, которое мы можем отнести к событиям некоторого пуассоновского потока П.. и считать, таким образом, что переход системы из состояния х в состояние х происходит под воздействием всего потока /L. Привлечение всего потока П.. дает нам возможность рассматривать интенсивность А() этого потока.  

Рассмотрим более подробно случай пуассоновского распределения спроса. Функция затрат будет иметь вид, аналогичный (5.6.18), с заменой интегрирования по х суммированием. Найдем плотность 1> (т) распределения времени дефицита. Распределение времени наступления k -го события пуассоновского потока подчинено закону Эрланга k -го порядка. Дефицит начинается при израсходовании всего запаса S и еще одной единицы, так что  

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

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

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

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

. (3.1)

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

, (3.2)

которая называется формулой Пуассона, а число - параметром Пуассона.

Задача 3.1. Вероятность брака при изготовлении некоторого изделия равна 0,008. Найти вероятность того, что при контроле среди 500 изделий будет не более двух бракованных.

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

Определение 3.1

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

Например , потоком событий будут вызовы, поступающие на АТС, сигналы при сеансе радиосвязи, сообщения, поступающие на сервер, и.т.д.

Определение 3.2

Поток событий называется пуассоновским (простейшим) если он обладает следующими свойствами:

1. Свойством стационарности , т.е. интенсивность потока - постоянная.

2. Свойством ординарности, т.е. появление двух или более событий за малый промежуток практически невозможно.

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

Если обозначить - вероятность появления событий пуассоновского потока c интенсивностью за время , то справедлива формула:

. (3.3)

Задача 3.2. Страховая компания обслуживает 10000 клиентов. Вероятность того, что в течение одного дня клиент обратится в компанию, равна 0,0003. Какова вероятность того, что в течение двух дней в нее обратятся 4 клиента?

Решение: Интенсивность потока клиентов в течение одного дня равна

Следовательно, .

Решение задач 3.1 и 3.2 в среде Mathcad показано на рис. 3.

Задача 3.3. Вероятность сбоя считывающего устройства турникета метрополитена в течение часа мала. Найти эту вероятность, если вероятность того, что за 8 часов будет хотя бы один сбой, равна 0,98, и если известно, что за час через турникет проходит в среднем 1000 человек?

Решение: По формулам (1.3) и (3.3) при вероятность того, что в течение 8 часов будет хотя бы один сбой, равна:

С помощью символьных команд, а затем определяется искомая вероятность .

Рассмотрим некоторую физическую систему S с дискретными состояниями которая переходит из состояния в состояние под влиянием каких-то случайных событий, например, вызовы на телефонной станции, выходы строя (отказы) элементов аппаратуры, выстрелы, направленные по цели и т. д.

Будем себе это представлять так, будто события, переводящие систему из состояния в состояние, представляют собой какие-то потоки событий (потоки вызовов, потоки отказов, потоки выстрелов и т. д.).

Пусть система S с графом состояний, показанным на рис. 4.27, в момент t находится в состоянии S; и может перейти из него в состояние под влиянием какого-то пуассоновского потока событий с интенсивностью как только появляется первое событие этого потока, система мгновенно переходит (перескакивает) из S в Как мы знаем, вероятность этого перехода за элементарный промежуток времени (элемент вероятности перехода) равна . Таким образом, плотность вероятности перехода в непрерывной цепи Маркова представляет собой не что иное, как интенсивность потока событий, переводящего систему по соответствующей стрелке.

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

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

Итак, рассматривается система S, в которой переходы из состояния в состояние происходят под действием пуассоновских потоков событий с определенными интенсивностями. Проставим эти интенсивности (плотности вероятностей переходов) на графе состояний системы у соответствующих стрелок.

Получим размеченный граф состояний (рис. 4.27); по которому, пользуясь правилом, сформулированным в § 3, можно сразу записать дифференциальные уравнения Колмогорова для вероятностей состояний.

Пример 1. Техническая система S состоит из двух узлов: I и II; каждый из них независимо от другого может отказывать (выходить из строя). Поток отказов первого узла - пуассоновский, с интенсивностью второго - также пуассоновский, с интенсивностью Каждый узел сразу после отказа начинает ремонтироваться (восстанавливаться). Поток восстановлений (окончаний ремонта ремонтируемого узла) для обоих узлов - пуассоновский с интенсивностью К.

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

Решение. Состояния системы:

Оба узла неправды,

Первый узел ремонтируется, второй исправен,

Первый узел исправен, второй ремонтируется,

Оба узла ремонтируются.

Размеченный граф состояний системы показан на рис. 4.28.

Интенсивности потоков событий на рис. 4.28 проставлены из следующих соображений. Если система S находится в состоянии то на нее действуют два потока событий: поток неисправностей узла I с интенсивностью X, переводящий ее в состояние и поток неисправностей узла II с интенсивностью переводящий ее в Пусть теперь система находится в состоянии (узел I ремонтируется, узел II - исправен). Из этого состояния система может, во-первых, вернуться в (это происходит под действием потока восстановлений с интенсивностью ); во-вторых, - перейти в состояние (когда ремонт узла I еще не закончен, а узел II тем временем вышел из строя); этот переход происходит под действием потока отказов узла II с интенсивностью Интенсивности потоков у остальных стрелок проставляются аналогично.

Обозначая вероятности состояний и пользуясь правилом, сформулированным в § 3, запишем уравнения Колмогорова для вероятностей состояний:

Начальные условия, при которых нужно решать эту систему: при

Заметим, что, пользуясь условием

можно было бы уменьшить число уравнений на одно. Действительно, любую из вероятностей можно выразить через остальные и подставить в уравнения (6.1), а уравнение, содержащее в левой части производную чтой вероятности - отбросить.

Заметим, кроме того, что уравнения (6.1) справедливы как для постоянных интенсивностей пуассоновских потоков X, так и для переменных:

Пример 2. Группа в составе пяти самолетов в строю «колонна» (рис. 4.29) совершает налет на территорию противника. Передний самолет (ведущий) является постановщиком помех; до тех пор, пока он не сбит, идущие за ним самолеты не могут быть обнаружены и атакованы средствами ПВО противника. Атакам подвергается только постановщик помех. Поток атак - пуассоновский, с интенсивностью X (атак/час). В результате атаки постановщик помех поражается с вероятностью р.

Если постановщик помех поражен (сбит), то следующие за ним самолеты обнаруживаются и подвергаются атакам ПВО; на каждый из них (до тех пор, пока он не поражен) направляется пуассоновский поток атак с интенсивностью X; каждой атакой самолет поражается с вероятностью р. Когда самолет поражен, атаки по нему прекращаются, но на другие самолеты не переносятся.

Написать уравнения Колмогорова для вероятностей состояний системы и указать начальные условия.

Решение. Будем нумеровать состояния системы соответственно числу сохранившихся самолетов в группе:

Все самолеты целы;

Постановщик помех сбит, остальные самолеты целы;

Постановщик помех и один бомбардировщик сбиты, остальные самолеты целы;

Постановщик помех и два бомбардировщика сбиты, остальные самолеты целы;

Постановщик помех и три бомбардировщика сбиты, один самолет цел;

Все самолеты сбиты.

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

Граф состояний системы показан на рис. 4 30. Чтобы разметить этот граф, определим интенсивности потоков событий, переводящих систему из состояния в состояние.

Из состояния систему переводит поток поражающих (или «успешных») атак, т. е. тех атак, которые приводят к поражению постановщика (разумеется, если он раньше не был поражен).

Интенсивность потока атак равна X, но не все они - поражающие: каждая из них оказывается поражающей только с вероятностью . Очевидно, интенсивность потока поражающих атак равна эта интенсивность и проставлена в качестве у первой слева стрелки на графе (рис. 4.30).

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

Здесь отброшены члены высшего порядка малости относительно Вычитая эту вероятность из единицы, получим вероятность перехода из за время (элемент вероятности перехода):

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

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

Система нескольких линейных дифференциальных уравнений такого типа только в редких случаях может быть проинтегрирована в квадратурах: обычно такую систему приходится решать численно - либо вручную, либо на аналоговой вычислительной машине (АВМ), либо, наконец, на ЭЦВМ. Все эти способы решения систем дифференциальных уравнений затруднений не доставляют; поэтому самое существенное - уметь записать систему уравнений и сформулировать для нее начальные условия, чем мы и ограничились здесь.

• Перевод

Введение

Одним из важнейших процессов, наблюдаемых в природе, является пуассоновский точечный процесс. Поэтому важно понять, как такие процессы можно моделировать. Методы моделирования различаются в зависимости от типа пуассоновского точечного процесса, т. е. пространства, в котором протекает процесс и однородности или неоднородности процесса. Мы не будем заинтересованы развитием пуассоновского точечного потока или с важными приложениями его в различных областях. Чтобы этот материал показался интересным, читателю настоятельно рекомендуется прочитать соответствующие разделы в Феллере (1965) и Синларе (1975) для основной теории и некоторые разделы в Триведи (1982) для приложений в ИТ.

На первом шаге мы определим пуассоновский процесс на = T
UNTIL Ложь (это бесконечный цикл; по желанию можно добавить критерий остановки)

Этот алгоритм просто реализовать, поскольку нет нужды генерировать пуассоновские случайные величины. Для других простых множеств A, существуют тривиальные обобщения теоремы 1.2. Например, когда A=x, где t может равняться бесконечности, 0 < T1 < T2 <… - равномерный пуассоновский процесс с интенсивностью λ и U1,U2,… - последовательность независимых одинаково распределённых равномерно на случайных величин, то (T1,U1),(T2,U2),… определяют пуассоновский процесс с интенсивностью λ на А.

Пример 1.1.

Равномерный пуассоновский процесс на единичной окружности
Если A - окружность с единичным радиусом, то разные свойства равномерного пуассоновского процесса можно использовать, чтобы получить несколько методов генерации (которые обобщаются на d-мерные сферы). Пусть λ - желаемая интенсивность.
Во-первых, мы просто могли бы сгенерировать случайную пуассоновскую величину N с параметром λπ, а затем вернуть последовательность N независимых одинаково распределённых равномерно на единично окружности векторов. Если мы применим метод порядковых статистик, предлагаемый теоремой 1.2, то пуассоновская случайная величина получается неявно. Например, перейдя в полярные координаты (R,φ) заметим, что для равномерного пуассоновского процесса R и φ независимы, и случайная величина R имеет плотность 2r, r меняется от 0 до 1, а φ равномерно распределена на . Таким образом, мы можем поступить следующим образом: Сгенерировать равномерный пуассоновский процесс 0 < φ1 < φ2 <… < φN с параметром интенсивности λ/(2π) на экспоненциальным методом и вернуть (φ1,R1),...,(φN,RN), где Ri - независимые одинаково распределённые случайные величины с плотностью 2r на , которые можно сгенерировать, взяв максимум из двух независимых равномерно распределённых на случайных величин. Особой причины применять эспоненциальный метод к углам нет. Таким же образом мы могли подобрать и радиусы. К сожалению, порядковые радиусы не формируют одномерный равномерный пуассоновский процесс на . Однако, тем не менее они образуют неоднородный пуассоновский процесс, и генерация таких процессов будет рассмотрена в следующем разделе.

Неоднородные пуассоновские процессы

Бывают такие ситуации, когда события происходят в «случайные моменты времени», но некоторые моменты более возможны, чем другие. Это случай прибытий в центры интенсивной терапии, предложений работ в компьютерных центрах и травмы игроков НХЛ. Для этих случаев очень хорошей моделью является модель неоднородного пуассоновского процесса, определённого здесь ради удобства на = T
UNTIL False

Пример 1.2. Однородный пуассоновский процесс

Для особого случая λ(t)=λ, Λ(t)=λt несложно видеть, что InvΛ(E+Λ(T))=T+E/λ, в результате чего мы снова получаем экспоненциальный метод.

Пример 1.3.

Для моделирования утреннего потока автомобилей перед часом пик, мы иногда можем взять λ(t)=t, тогда Λ(t)=t^2/2 и получим шаг

Если функцию интенсивности можно представить в виде суммы функций интенсивности, т. е. ,

0 < T i1 < T i2 <… T in - независимые реализации отдельных неоднородных пуассоновских процессов, то объединённая упорядоченная последовательность образует реализацию неоднородного пуассоновского процесса с функкцией интенсивности λ(t). Это относится к методу композиции, но разница теперь состоит в том, что нам нужны реализации всех компонентов процесса. Декомпозицию можно использовать, когда существует естественное разложение, продиктованное аналитической формой λ(t). Поскольку основная операция в слиянии процессов - взять минимальное значение из n процессов, для больших n преимущество может предоставить хранение моментов времени в куче из n элементов.

В итоге получим метод композиции:

Сгенерировать T,...,T для n пуассоновских процессов и хранить эти значения вместе с индексами соответствующих процессов в таблице
T = 0 (текущее время)
k = 0
REPEAT
Найти минимальный элемент T в таблице и удалить его
k = k + 1
T[k] = T
Сгенерировать T и вставить в таблицу
UNTIL False

Третий общий принцип - это принцип утоньшения (Льюис и Шедлер, 1979). Аналогично тому, что происходит в методе отклонения, предполагаем, что существует лёгкая доминирующая функция интенсивности λ(t) <= μ(t) для любого t.

Тогда идея состоит в том, чтобы сгенерировать однородный Пуассоновский процесс на части положительной полуплоскости между 0 и μ(t), затем рассмотреть однородный пуассоновский процесс под λ и, наконец, вернуть x-компоненты событий в этом процессе. Это требует следующей теоремы.

Теперь рассмотрим метод утоньшения Льюиса и Шедлера:

T = 0
k = 0
REPEAT
Сгенерировать Z, первое событие в неоднородном пуассоновском процессе с функцией интенсивности μ, который происходит после момента времени T. Присвоить T = Z
Сгенерировать равномерно распределённую на случайную величину U
IF U <= λ(Z)/μ(Z)
THEN k = k + 1, X[k] = T
UNTIL False

Утверждается, что последовательность X k так сгенерированная образует неоднородный пуассоновский процесс с функцией интенсивности λ. Заметим, что мы взяли неоднородный процесс 0 < Y1 < Y2 <… с функцией интенсивности μ и убрали некоторые точки. Насколько мы знаем, (Y i ,U i μ(Y i) - однородный пуассоновский процесс с единичной интенсивностью на кривой, если U i независимые одинаково распределённые равномерно на случайные величины в силу теоремы 1.3. Таким образом, подпоследовательность на кривой λ определяет однородный пуассоновский процесс с единичной интенсивностью на этой кривой (часть 3. теоремы 1.3). Наконец, взятие x-координат только этой подпоследовательности даёт нам неоднородный пуассоновский процесс с функцией интенсивности λ.
Неоднородный пуассоноский процесс с функцией интенсивности μ обычно моделируют методом инверсии.

Пример 1.4. Функция с циклической интенсивностью

Следующий пример также принадлежит Льюису и Шедлеру (1979). Рассмотрим функцию с циклической интенсивностью λ(t)= λ(1+cos(t)) с очевидным выбором доминирующей функции μ=2λ.

Тогда алгоритм моделирования примет вид:

T = 0
k = 0
REPEAT
Сгенерировать экспоненциальную случайную величину E c параметром 1
T = T + E/(2λ)
Сгенерировать равномерную на случайную величину U
IF U <= (1+cos(T))/2
THEN k = k + 1, X[k] = T
UNTIL False

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

Заключительное слово об эффективности алгоритма, когда моделируется неоднородный пуассоновский процесс на множестве . Среднее число событий, которое необходимо от доминирующего процесса, равно в то время, как среднее число возвращённых случайных величин равно
Отношение средних величин может быть рассмотрено как объективная мера эффективности, сравнимая в духе константы отклонения в стандартном методе отклонения. Заметим, что мы не можем использовать среднюю величину отношения, поскольку она, в общем случае, была бы равна бесконечности в силу положительной вероятности того, что ни одна величина не возвратится.

Информатика, кибернетика и программирование

Определение Пуассоновского потока. Пуассоновский поток это ординарный поток без последействия. Классической моделью трафика в информационных сетях является Пуассоновский простейший поток. Он характеризуется набором вероятностей Pk поступления k сообщений за временной интервал t: где k=01 число сообщений; λ интенсивность потока.

1. Определение Пуассоновского потока. Свойства.

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

Классической моделью трафика в информационных сетях является Пуассоновский (простейший) поток. Он характеризуется набором вероятностей P(k) поступления k сообщений за временной интервал t:

где k=0,1,… - число сообщений; λ - интенсивность потока.

Заметим, что интервал времени измерения количества сообщений t и интенсивность потока λ являются постоянными величинами.

Семейство Пуассоновских распределений P(k) в зависимости от λ изображено на рис.1. Большее значение λ соответствует более широкому и симметричному графику плотности вероятности.

Рис. 1. Пуассоновские распределения. Плотности вероятностей.

Математическое ожидание (среднее) и дисперсия Пуассоновского потока равны λ t .

Зная вероятность поступления данных за период, можно получить распределение интервала τ между соседними событиями:

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

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

При моделировании Пуассоновский поток можно получить мультиплексированием совокупности ON/OFF источников, которые называются Марковскими процессами (рис.2.).

Рис. 2. Получение Пуассоновского распределения

2. СМО с отказами (классическая система Эрланга)

Здесь мы рассмотрим одну из первых по времени, «классических» задач теории массового обслуживания; эта задача возникла из практических нужд телефонии и была решена в 1909 г. датским инженером-математиком А.К. Эрлангом. Задача ставится так: имеется n каналов (линий связи), на которые поступает поток заявок с интенсивностью λ. Поток обслуживаний каждого канала имеет интенсивность μ. Найти предельные вероятности состояний системы и показатели ее эффективности.

Система S (СМО) имеет следующие состояния (нумеруем их по числу заявок, находящихся в системе): S 0 , S 1 ,…, S n , где S k – состояние системы, когда в ней находится k заявок, т.е. занято k каналов.

Граф состояний СМО соответствует процессу гибели и размножения (рис. 3).

Рис. 3. Граф состояний СМО

Поток заявок последовательно переводит систему из любого левого состояния в соседнее правое с одной и той же интенсивностью λ. Интенсивность же потока обслуживаний, переводящих систему из любого правого состояния в соседнее левое, постоянно меняется в зависимости от состояния. Действительно, если СМО находится в состоянии S 2 (два канала заняты), то она может перейти в состояние S 1 (один канал занят), когда закончит обслуживание либо первый, либо второй канал, т.е. суммарная интенсивность их потоков обслуживаний будет 2μ . Аналогично суммарный поток обслуживаний, переводящий СМО из состояния S 3 (три канала заняты) в S 2 , будет иметь интенсивность 3μ , т.е. может освободиться любой из трех каналов, и т.д.

В формуле (1) для схемы гибели и размножения получим для предельной вероятности состояния:

(1)

где члены разложения - коэффициенты при p 0 в выражениях для предельных вероятностей p 1 , p 2 ,..., p n .

Заметим, что в формулу (1) интенсивности λ и μ входят не по отдельности, а только в виде отношения μ/λ. Обозначим: μ/λ = p , и будем называть величину ρ приведенной интенсивностью потока заявок или интенсивностью нагрузки канала. Она выражает среднее число заявок, приходящих за среднее время обслуживания одной заявки. Пользуясь этим обозначением, перепишем формулу (1) в виде:

(2)

При этом:

(3)

Формулы (2) и (3) для предельных вероятностей получили названия формул Эрланга в честь основателя теории массового обслуживания.

Вероятность отказа СМО есть предельная вероятность того, что все n каналов системы будут заняты, т.е.

Отсюда находим относительную пропускную способность – вероятность того, что заявка будет обслужена:

Абсолютную пропускную способность получим, умножая интенсивность потока заявок λ на Q:

(4)

Осталось только найти среднее число занятых каналов k. Эту величину можно было бы найти «впрямую», как математическое ожидание дискретной случайной величины с возможными значениями 0,1,..., n и вероятностями этих значений p 0 , p 1 , …, p n :

Подставляя сюда выражения (3) для p k и выполняя соответствующие преобразования, мы, в конце концов, получили бы формулу для k. Однако среднее число занятых каналов можно найти проще, если учесть, что абсолютная пропускная способность A системы есть не что иное, как интенсивность потока обслуженных системой заявок (в единицу времени). Так как каждый занятый канал обслуживает в среднем μ заявок (в единицу времени), то среднее число занятых каналов:

или, учитывая (4):

А также другие работы, которые могут Вас заинтересовать
43346.		Створення інформаційної бази даних служби продажу залізничних білетів	1.19 MB
	Курсова робота з дисципліни Організація баз даних та знань на тему: Створення інформаційної бази даних служби продажу залізничних білетів Курсова робота студента 3 курсу групи КН48 Нестеренка М. Проектування інформаційної бази даних Створення реляційної моделі бази даних Створення бази даних
43347.		Технологии аппаратной виртуализации	64.5 KB
	Аппаратная виртуализация - виртуализация с поддержкой специальной процессорной архитектуры. В отличие от программной виртуализации, с помощью данной техники возможно использование изолированных гостевых систем, управляемых гипервизором напрямую. Гостевая система не зависит от архитектуры хостовой платформы и реализации платформы виртуализации.
43349.		Публіцистика Уласа Самчука	108.5 KB
	Деяка молодь не знає своєї історії і не може відповісти на питання Хто були їх предки. А коли все одно то це значить що все одно для нас хто є ми самі Це значить що ми не нарід не якась спільна історична збірна сила а невиразна юрба сіра маса вічно принижена без всяких ідеалів чернь4. Замкнутість людини лише у сільському просторі неминуче призведе до відчуження її від міста.
43351.		РЕСУРСИ КОМЕРЦІЙНОГО БАНКУ, ЇХ ФОРМУВАННЯ ТА МЕНЕДЖМЕНТ	195 KB
	Ресурси комерційного банку їх склад та структура. Власний капітал комерційного банку та його формування. Залучений капітал комерційного банку та його характеристика.
43352.		Флористичні дані про медоносні рослини луків, які найбільш поширені на території України	210.5 KB
	Загальна характеристика Шишацького району та Полтавської області Розділ 2.6 Інші рослини Висновки Література ВСТУП Об"єктом нашого дослідження є медоносні рослини які поширені в Шишацько му районіПолтавської областіїх значення у житті людини та характеристика. В Шищацькому районіПолтавської області нараховується велика кількість медоносних рослин. Загальна характеристика Шишацького...
43354.		РЕАЛІЗАЦІЯ ЗМІШАНОЇ КРИПТОСИСТЕМИ ПОПЕРЕДНЬОГО ШИФРУВАННЯ	544 KB
	Цей пункт забезпечує разгортання системи конфіденційного звязку управління роботою системи конфіденційного звязку первісну генерацію та розповсюдження ключів управління персоналом. Для автентифікації відкритих ключів дозволяеться використовувати послуги регіонального центру сертифікації. Генерація сеансових ключів для ГОСТ 2814789 перешифрування в режими простої заміни алгоритма ГОСТ 2814789 послідовності що отримана відповідно до пункту 1.2 Розповсюдження ключів за допомогою асиметричних криптосистем } О сайте Контакты Реклама Написать нам