Поток событий называется потоком с ограниченным последействием, если случайные величины T1,T2,...,Tn..., представляющие собой интервалы времени между соответственно 1-м и 2-м, 2-м и 3-м и т.д., n-м и (п+ 1)-м событиями и т.д., независимы.
Поток событий называется стационарным, если вероятность наступления того или иного числа событий за какой-либо промежуток времени зависит только от длины этого промежутка и не зависит от момента его начала. Стационарность потока означает, что его вероятностные характеристики не зависят от времени, т.е. не изменяются с течением времени.
Стационарный поток с ограниченным последействием называется потоком Пальма.
У потока Пальма случайные величины T1,T2,...,Tn... имеют один и тот же закон распределения (имеют одну и ту же функцию распределения). Пуассоновский процесс является потоком Пальма, поскольку он стационарен, случайные величины T1,T2,...,Tn...распределены по показательному закону (докажите это) и независимы в силу независимости приращений.
Важными специальными случаями потока Пальма являются потоки Эрланга. Потоком Эрланга k-го порядка с параметром 0 называется поток, получающийся из пуассоновского с интенсивностью 0 сохранением в нем
каждого k-то события.
Например, поток Эрланга 1-го порядка совпадает с исходным пуассоновским потоком и, следовательно, пуассоновский поток является потоком Эрланга 1-го порядка.
Для случайной величины Tk - промежутку времени между двумя соседними
событиями в потоке Эрланга k-го порядка, порожденном Пуассоновским процессом с интенсивностью 0 можно найти плотность:
p |
t |
t k 1 |
e t ,t 0, |
|
|||
T |
|
k 1 ! |
|
k |
|
||
|
|
|
|
а также функцию распределения
|
|
k 1 |
x |
m |
FT k |
1 e t |
|
,t 0, |
|
x P Tk x |
m 0 m! |
|
||
|
|
0,t 0. |
|
|
|
|
|
||
Пользуясь методами теории вероятностей можно доказать, что распределение Эрланга k-го порядка есть распределение суммы k независимых в совокупности, экспоненциально распределенных случайных величин с одним и тем же параметром . Следовательно,
ETk |
|
k |
,ETk2 |
|
|
k(k |
1) |
,DTk |
ETk2 |
ETk 2 |
k |
, k |
|
|
k |
. |
|
|
DTk |
||||||||||||||||
|
2 |
|
2 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
46
Интенсивность потока Эрланга k-го порядка совпадает с математическим ожиданием Tk .
Отметим, что обобщением пуассоновского процесса является нестационарный пуассоновский процесс или поток, для которого выполняется условие независимости приращений, 0 0, и ординарности, но вместе с тем
интенсивность t limP t h t 1 зависит от t [0; ).
h 0 h
Если для нестационарного пуассоновского потока с интенсивностью t ,
некоторый рассмотреть промежуток времени длиной 0, начинающийся с момента t0 и дискретную случайную величину X t0, — число событий,
наступающих в потоке за промежуток времени t0,t0 , то выполняется свойство, что X t0, имеет распределение Пуассона с параметром, равным
t0 |
t dt |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
t0 |
Рассмотрим теперь |
процесс |
t ,t T , |
который |
состоит из |
случайных |
|||||||
|
|||||||||||||
величин, принимающих значения из . |
|
|
|
|
|
||||||||
|
Определение. |
Винеровским |
процессом |
называется |
случайный процесс |
||||||||
t ,t T |
с непрерывным временем T 0, , удовлетворяющий условиям: |
||||||||||||
|
1). 0 0. |
|
|
|
|
|
|
|
|
|
|||
|
2). |
|
|
,t [0; ) |
– процесс |
с |
независимыми приращениями. |
|
|||||
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
3). Для любых t1 t2 |
и любого |
s приращения t |
2 |
t , |
t2 s t1 s |
одинаково |
||||||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
распределены.
4). При h 0 для моментов случайных величин h выполняются следующие свойства:
E h ah o h ,
E h2 bh o h ,
E |
|
h |
|
3 o h , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где a,b - постоянные величины, не зависящие от |
h, b 0. |
|
|
|||||||||||||||
Величины a и b являются параметрами винеровского процесса. |
||||||||||||||||||
Теорема (без доказательства). Если t ,t 0, |
- винеровский процесс, то |
|||||||||||||||||
при фиксированном |
t 0 выполняется, |
что |
t N at,bt , |
т.е t имеет |
||||||||||||||
нормальное распределение с параметрами |
|
|
at |
и |
bt, и функция распределения |
|||||||||||||
случайной величины t |
имеет вид |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
x |
|
|
|
2 |
|
|
|
|
|
F x P |
|
x |
|
|
|
|
|
exp y at |
|
dy. |
||||
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
t |
|
t |
|
|
|
|
|
|
|
|
2bt |
|
|
|
|
|
|
|
|
|
|
|
|
2 bt |
|
|
|
|
|
|
|
|
|
В отличие от пуассоновского процесса, винеровский имеет непрерывную траекторию.
47
Винеровский процесс служит математической моделью одномерного броуновского движения и часто называется процессом броуновского движения.
В случае a 0 и b 1 винеровский процесс называется стандартным винеровским процессом.
Функционалы винеровского процесса
Рассмотрим траекторию стандартного винеровского процесса t ,t 0, и
зафиксируем некоторое вещественное число a.
Пусть a - момент первого достижения уровня a траекторией случайного винеровского процесса.
Величина a является случайной величиной.
Фиксируем некоторое вещественное T 0. Рассмотрим выражение
P T a| a T P T a| a T P T a| a T 1
Утверждается, что
P T a| a T P T a| a T .
Если сформулировать нестрого, то число траекторий, входящих в левую часть, совпадает с числом траекторий, входящих в правую (т.е. происходит «отражение»). Каждой траектории мы сопоставляем отражённую траекторию.
Для винеровского процесса, обладающего непрерывной траекторией, выполняется равенство P T a| a T 0. Следовательно,
P T a| a T P T a| a T 1/ 2.
Данное выражение есть принцип отражения.
В силу непрерывности траекторий винеровского процесса верно следующее включение: T a a T . Следовательно,
48
P a| |
a |
T |
P T a, a T |
|
P T a |
|
1 |
. |
|
|
|
||||||
T |
|
P a T |
|
P a T 2 |
||||
|
|
|
|
|||||
Из этого равенства следует, что P a T 2P T |
a . |
|||||||
Рассмотрим стандартный винеровский процесс и найдём распределение его основных функционалов.
Как мы уже определили ранее, a - момент первого достижения уровня a траекторией винеровского процесса. По определению функции распределения
Fa t =P a t .
Из принципа отражения, как показано выше, следует, что
|
|
t =P a t 2P t a 2 |
|
1 |
|
|
|
|
y2 |
|
||||||||||||||||
|
Fa |
|
|
|
|
|
a e 2t dy. |
|||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||
|
2 t |
|||||||||||||||||||||||||
Так как t t 0 |
N(0,t) |
и с учётом замены y / |
|
|
x мы получаем, что это |
|||||||||||||||||||||
|
t |
|||||||||||||||||||||||||
выражение равно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
dx. |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Таким образом, |
P |
|
|
|
|
|
|
2 dx 1. |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Следовательно, |
с |
|
вероятностью |
1 |
|
|
случайная |
величина a конечна, и с |
||||||||||||||||||
вероятностью 1 траектория винеровского процесса достигает любого уровня. Теперь рассмотрим плотность a :
|
|
dF |
(t) |
|
1 |
|
|
a |
|
a2 |
|
|
|
|
|
|
|||||||
p a |
|
|
|
|
|
|
|
e 2t |
|||
|
|
|
|
|
|||||||
|
|
||||||||||
|
|
a |
|
2 |
|
|
t3/2 |
|
|
||
dt |
|
|
|||||||||
|
|
|
|
0, |
|
|
|||||
|
|
|
|
|
|
|
|
||||
,t 0. t 0
Используя эту плотность, легко показать, что E a . Таким образом, несмотря на то, что с вероятностью 1 случайная величина a конечна, ее среднее значение равно бесконечности.
Рассмотрим другой функционал винеровского процесса, равный
максимальному значению s |
по всем s |
|
|
из отрезка 0,t : |
t max s . |
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s 0,t |
|
|
|
|
|
|||
|
max |
|
|
x |
|
|
|
t |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Очевидно равенство s 0,t |
|
|
s |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
2 |
|
|
|
1 F x P |
|
|
|
|
|
|
|
|
P |
|
t |
|
|
|
|
|
t |
|
2 |
|
|
|
|
|
|||||||||||
max |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
|
|||||||||||||||||
s |
x |
x |
F |
|
|
|
2 dy |
. |
|||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||
t |
s 0,t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
|
|
|
|
|
|
dv |
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
Если сделать замену y |
|
v; |
y |
|
|
; |
|
dy |
|
, то получаем: |
|
|
|
|
|
||||||||||||||||||||
t |
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
t |
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
49
|
|
|
|
1 Ft x 2 |
|
||
|
e |
||
|
|
|
|
|
t |
|
x |
Отсюда следует, что
v2 |
|
|
2 |
|
|
|
|
|
|
v2 |
|
|
|
|
|
|
|
|
|
|
|
||||
2t dv |
|
|
e 2t dv |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
t |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
x |
|
v2 |
|
|
|
|
|
|
|
|
|
x |
|
|
v2 |
||||
|
2 |
|
|
|
|
dv 1 2 |
|
|
|
||||||||||||||||
|
|
e |
2t |
|
|
e |
2t dv. |
||||||||||||||||||
|
|
t |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
t |
|
0 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
v2 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||
Ft x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
e 2t dv при |
x 0 . |
|||||||||||||||||||
|
|
|
t |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||
Зная функцию распределения t |
, легко найти плотность |
|||||||||||
|
|
dF t |
x |
|
|
|
|
|
x2 |
|
||
|
x |
|
|
2 |
|
e |
|
|||||
p |
|
|
2t |
при |
x 0. |
|||||||
dx |
t |
|||||||||||
t |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||
Можно еще рассмотреть функционал винеровского процесса, равный минимальному значению s по всем s из отрезка 0;t :
t min s .
s 0,t
Докажите самостоятельно равенство для плотности этого функционала: p t x p t x при x<0.
Для доказательства рассмотрите новый винеровский процесс
t ,t 0, .
Тема № 5 Теория массового обслуживания
Основные понятия и классификация
Теория массового обслуживания (иначе называемая теорией очередей) изучает случайные процессы в системах массового обслуживания с приложением к рациональному построению этих систем. Она устанавливает зависимость между характером потока заявок (требований), производительностью отдельного канала (обслуживающего аппарата), числом каналов и эффективностью обслуживания.
Характерной особенностью массового обслуживания является наличие потока однородных (идентичных, однотипных) требований (событий, заявок), которые подвергаются обслуживанию. Заявкой (или требованием) называется спрос на удовлетворение какой-либо потребности. Выполнение заявки называется обслуживанием заявки.
Система массового обслуживания (СМО) – это объект, в котором выполняется последовательность (элементарных) операций или, иначе говоря, любая система для выполнения заявок, поступающих в неё в случайные моменты времени.
Множество моментов поступления в систему требований называется входным потоком требований данной системы массового обслуживания.
Будем считать, что в момент поступления требования происходит событие. Тогда множество моментов, когда происходят события, называется потоком однородных событий.
50