Глава 5 теория очередей
§ 27. Очереди к одному обслуживающему прибору
Теория очередей начала развиваться в XX веке в связи с запросами телефонии. Первые работы были сделаны Эрлангом [20 — 22], изучившим распределение времени задержки вызовов в телефонных переговорах. Значительный прогресс в математической теории очередей был достигнут в тридцатые годы благодаря работам Полячека [53, 54], А. Н. Колмогорова [45], А. Я. Хинчина [42, 43] и других. В настоящее время имеется огромное количество литературы по теории очередей и ее применениям (см. Дойг [19], Саати [66], Вольд [90]). Применяется теория очередей в основном в технике (телефония, сети связи, электронные вычислительные машины), в промышленности (обслуживающие автоматы, конвейеры, склады), на транспорте (аэропорты, гавани, железнодорожные и автобусные станции, уличное движение, почта), в торговле (рынки, сбыта, банки, билетные кассы), а также в повседневной жизни (лифты, рестораны, парикмахерские).
Мы будем пользоваться терминами: система обслуживания, требования, обслуживающие приборы, время обслуживания. На этом языке можно описать любой мыслимый процесс. Например, в телефонии термины телефонная станция, вызовы, линии, время занятости соответствуют системе обслуживания, требованиям, обслуживающим приборам, времени обслуживания.
Механизм образования очереди очень прост. Требования поступают в очередь и обслуживаются одним или более приборами. После обслуживания каждое требование выходит из системы. Время, проведенное требованием в системе, состоит из времени ожидания (возможно, равного нулю) и времени обслуживания. Время прибора складывается из чередующихся между собой периодов занятости и периодов, свободных от требований.
Наиболее важные задачи в теории очередей связаны со случайными флуктуациями длины очереди (линии ожидания) и случайными флуктуациями времени ожидания (задержки). Знание стохастических законов, управляющих этими флуктуациями, дает возможность проектировать требуемые системы обслуживания (достаточно большое помещение для ожидания, достаточное число приборов и т. д.).
В этой главе мы будем рассматривать следующую математическую модель обслуживания: в интервале времени [0, ∞) требования прибывают на обслуживание в соответствии с некоторым случайным процессом. Прибывшие требования обслуживаются одним обслуживающим прибором, причем времена обслуживания являются случайными величинами. Порядок обслуживания не задается, но предполагается, что обслуживающий прибор занят, если в системе есть хотя бы одно требование.
Мы изучим задачи, связанные с флуктуациями длины очереди и флуктуациями времени ожидания.
Длина очереди в момент t обозначается через ξ(t) и определяется как число требований в системе в момент t, включая обслуживаемое, если таковое имеется. Будем обозначать через ᶓп длину очереди непосредственно перед поступлением n-го требования, а через t,n— длину очереди непосредственно после окончания обслуживания n-го требования.
Время ожидания в момент t обозначается через η(t) и определяется как время, необходимое для завершения обслуживания всех требований, имеющихся в системе к моменту t. Если, в частности, обслуживание производится в порядке поступления, то η(t) является временем ожидания требования, поступившего в систему в момент t. При этом η(t) можно интерпретировать как виртуальное время ожидания в момент t, определяемое для всех t ≥ 0. Если требование поступает в момент t, то его действительное время ожидания равно η(t — 0). Виртуальному времени ожидания можно придавать реальный физический смысл. Например, если рассматривать поступление телеграфных сообщений, то виртуальное время ожидания равно длине ленты непрочитанной части телеграфного сообщения к моменту t. Можно даже представить себе, что используется стрелка с часовым механизмом, отсчитывающая время, и в момент поступления требования мы передвигаем стрелку вперед на длину, равную времени, необходимому для его обслуживания. Так как такие часы идут только до тех пор, пока в системе есть требования, они в любой момент показывают виртуальное время ожидания. Таким образом, на этих часах прибывающий клиент может немедленно увидеть свое действительное время ожидания. Вообще η(t) можно интерпретировать как время занятости (полной загрузки) прибора в момент t. Через ηп мы будем обозначать время ожидания непосредственно перед поступлением η-го требования. Если обслуживание производится в порядке поступления, то цп есть истинное время ожидания η-го поступившего требования.
Процесс обслуживания можно охарактеризовать с двух различных точек зрения в соответствии с тем, интересуемся ли мы флуктуациями длины очереди или флуктуациями времени ожидания.
Процесс Q. Предположим, что обслуживающий прибор начинает работу в момент времени t= 0 и к этому моменту ζ0 требований уже ожидают обслуживания. Начальная длина очереди ζ0 является случайной величиной, принимающей неотрицательные целые значения. Пусть vb v2, ..., vr, ...—числа требований, вставших в очередь во время обслуживания первого, второго, . . ., r-го, ...требований, и пусть No = 0, Nr = 4i + ... +\г для г =
= 1, 2 В этом случае будет рассматриваться процесс типа
Q = {ζo; Nn r = 0, 1, 2, ...}.
Нас будут интересовать распределения следующих случайных величин:
ᶓn (n=1, 2, ...), длины очереди непосредственно перед поступлением n-го требования;
ζn (n=1, 2, ...), длины очереди непосредственно после окончания обслуживания n-го требования; ζ0 — начальная длина очереди;
pn (n = 0, 1, 2, ...), числа требований, обслуженных в и-й период занятости;
an числа нулей среди ζ0,ζ1..., ζn_1
ßn числа положительных членов среди ζ0,ζ1..., ζn_1 . Очевидно, что аn + ßn= п.
Случайные величины ап и ßn можно также интерпретировать как число обслуживании среди первых п обслуживании, которым предшествует период бездействия прибора, и число обслуживании среди первых п обслуживании, которым не предшествует период бездействия прибора. При ζ0 = i среди первых п поступающих требований число требований, для которых в момент их поступления обслуживающий прибор занят, равно ßn+1-j.
Все эти случайные величины полностью определяются заданием
ζ0 и Nr r= 0, 1 ,2, ....
В дальнейшем мы всегда будем предполагать, что ζ0 и Nr, r = 0, 1, 2, . . ., независимы. Случайные величины v1 v2, .. ., vr, ... будут либо переставляемыми случайными величинами, принимающими неотрицательные целые значения, либо, в частности, взаимно независимыми и одинаково распределенными случайными величинами, принимающими неотрицательные целые значения.
Процесс W. Пусть обслуживающий прибор начинает работу в момент времени u = 0 и в этот момент его начальное время занятости определяется неотрицательной случайной величиной η (0). Обозначим через X(u) полное (накопленное) время обслуживания всех требований, прибывших в интервале времени [0, и]. В этом случае рассматриваемый процесс будет процессом образования очереди — процессом типа
w = η(0); x(a), 0<u<∞}.
Нас будут интересовать распределения следующих случайных величин:
η(t), времени ожидания в момент t;
θr(r = 0, 1,2, ...), длины r-го периода занятости;
a(t), полного (накопленного) времени бездействия обслуживающего прибора в интервале (0, t);
ß(t), полного (накопленного) времени занятости обслуживающего прибора в интервале (0, t). Очевидно, что а (t) + ß (t) = t.
Все эти случайные величины полностью определяются заданием
η(0) и {х(u), 0≥u<∞}.
В дальнейшем мы всегда будем предполагать, что случайная величина η(0) и случайный процесс {x(u), 0≤u< ∞} независимы. Случайный процесс {x(u), 0≤u<∞} будет иметь либо переставляемые, либо стационарные независимые приращения. Почти все его выборочные функции будут неубывающими ступенчатыми функциями, обращающимися в нуль при и = 0.