Введение
Системы массового обслуживания с ограниченным временем жизни заявок встречаются во многих областях деятельности человека: информационные технологии, финансы, телемедицина, техника, маркетинг и др. [1–3]. Характерной чертой таких систем является риск утраты заявок в процессе обслуживания, что приводит к неэффективному использованию ресурсов [4]. Ограничение на срок жизни требования может определяться случайным образом, согласно некоторому вероятностному распределению [5–7]. Случайное время пребывания заявки в системе характерно для многих важных процессов: оказание экстренной помощи тяжелым больным, затяжные процессы со сроком давности в юриспруденции, спасение людей в чрезвычайных ситуациях, обработка объектов в ограниченной зоне действия обслуживающего устройства [8].
Получение эксплуатационных характеристик таких СМО необходимо не только для оценки эффективности их функционирования, но и для прогнозирования поведения систем при различных сценариях нагрузки в процессе предпроектных исследований.
Цель исследования – получение аналитического решения задачи нахождения характеристик функционирования одноканальной системы массового обслуживания с накопителем и случайным ограничением на время жизни заявок для проверки адекватности имитационной модели.
Материалы и методы исследования
Системы массового обслуживания с ограниченным временем жизни заявок требуют особого подхода к изучению и специализированного математического аппарата, так как введение временных границ нарушает свойство отсутствия последействия, характерное для классических марковских моделей. Это происходит из-за того, что вероятность перехода заявки из состояния ожидания в состояние обслуживания начинает зависеть от времени, которое она уже провела в накопителе. Таким образом, пространство состояний перестает быть дискретным, что влечет необходимость введения функций распределения дополнительных случайных величин. Это позволяет сохранить аналитическую замкнутость и избежать расширения пространства состояний системы.
Одними из первых публикаций, посвященных моделированию СМО с ограниченным временем жизни заявок, являются работы советских математиков Б. В. Гнеденко и И. Н. Коваленко [9, 10], а также венгерского ученого Л. Такача [11]. Ими проработаны концептуальные вопросы моделирования одноканальных систем с простейшим входным потоком заявок и предложен соответствующий математический аппарат. Фокус исследований сосредоточен на нахождении функции распределения фактического времени ожидания обслуживания, которая является аналитическим базисом для определения ряда производных характеристик СМО. В частности, на ее основе находятся функция распределения реального времени пребывания заявок в системе, вероятность полного обслуживания, вероятность отказа и др. Таким образом, построение функции распределения фактического времени ожидания обслуживания обеспечивает полное вероятностное описание рассматриваемой СМО.
Результаты исследования и их обсуждение
Рассмотрим одноканальную СМО с накопителем, время жизни заявок в которой ограничено и является случайной величиной, входной поток требований простейший.
Для построения модели вводятся следующие характеристики:
λ – интенсивность простейшего входного потока заявок;
H(х) – функция распределения случайной величины – времени обслуживания;
В(х) – распределение вероятностей того, что время ожидания заявки в очереди равно некоторой случайной величине х;
Gy(x) – распределение вероятностей того, что после пребывания заявки в очереди в течение отрезка времени y, заявка попадет на обслуживание, где будет находиться в течение отрезка времени, не большего х;
ξ(t) – случайная величина – время ожидания начала обслуживания заявки, поступившей в систему в момент времени t;
F (t, x) = P(ξ(t) < x) – функция распределения фактического времени ожидания обслуживания.
Пусть заявка пришла в систему в момент времени t. За промежуток времени t + Δt могут произойти следующие несовместные события, которые обеспечат выполнение неравенства ξ(t + Δt) < x:
− за указанный промежуток времени не поступит ни одна заявка;
− за указанный промежуток времени поступает одна заявка и уходит из очереди из-за истечения срока жизни;
− за указанный промежуток времени поступает одна заявка, ожидает в очереди в течение некоторого времени, начинает обслуживание, но не успевает его завершить из-за истечения срока жизни;
− за указанный промежуток времени поступает одна заявка, ожидает в очереди в течение некоторого времени, начинает обслуживание и завершает его до истечения срока жизни.
Рассмотренные несовместные события образуют полную группу, поэтому вероятность того, что неравенство ξ(t + Δt) < x верно, равно сумме вероятностей описанных событий:

(1)

При делении уравнения (1) на Δt, Δt→0 и учете того, что в стационарном состоянии, существование которого показано в работе [12], верно равенство
(2)
получается уравнение для нахождения функции распределения ожидания обслуживания:
(3)
Функция F(x) терпит разрыв в точке 0 и абсолютно непрерывна при х > 0 [10]. Если обозначить через f(x) производную справа от функции F(x), то формула (3) преобразуется к виду

(4)
где F(0) находится из условия нормировки
(5)
здесь Т – максимальное время жизни требования.
Наиболее значимыми эксплуатационными характеристиками функционирования любых СМО являются доли обслуженных и необслуженных заявок. В случае, когда ограничено время жизни заявок, появляется еще и множество заявок, которые начали обслуживание, но не завершили его.
Для того, чтобы получить долю заявок, совсем не попавших на обслуживание, нужно применить формулу
(6)
Для того, чтобы найти долю заявок, не получивших полное обслуживание, нужно применить формулу
Pнедо
(7)
Уравнение (4) универсально в том смысле, что обеспечивает возможность моделирования разных способов наложения ограничения на время жизни заявок: отсутствие ограничения как такового, случайное ограничение на время ожидания обслуживания в очереди, случайное ограничение на время пребывания заявок в системе.
Важно отметить, что функция Gy(x) предопределяется только в системах со сложной дисциплиной обслуживания. В случае применения стандартной дисциплины FIFO процесс функционирования системы обработки заявок целиком и полностью обусловлен законами распределения входного потока, потока обслуживания и времени жизни заявок.
В рассматриваемом авторами случае ограничение времени жизни заявок является общим, то есть накладывается на время пребывания заявок в системе. При этом истечение времени жизни требования может состояться как в ходе ожидания, так и в процессе обслуживания. Для такого характера ограничения функция B(x) определяет распределение времени пребывания заявок в системе, уравнение (4) преобразуется к виду
(8)
а функция Gy(x) находится как
(9)
Уравнение (8) относится к линейным уравнениям Вольтерра второго рода. Основными методами решения таких уравнений являются метод преобразований Лапласа или метод последовательных приближений. В отдельных случаях оно может быть сведено к дифференциальному уравнению дифференцированием по верхнему пределу.
Оценим трудоемкость аналитического решения на конкретном примере.
Пусть средняя интенсивность простейшего входного потока заявок λ = 6, средняя интенсивность обслуживания μ = 10, ограничение времени жизни заявки распределено равномерно на интервале [0, 2]. Тогда уравнение (8) принимает вид
(10)
Решим уравнение, учитывая нормировку (3) при условии х <= 2.
Введем обозначение
(11)
Тогда
(12)
Подставив выражение (12) в уравнение (10), получим

или
(13)
Это линейное дифференциальное уравнение первого порядка. Для решения введем интегрирующий множитель:

Решением уравнения (11) является функция
(14)
Тогда
(15)
Для нахождения константы F(0) воспользуемся условием нормировки (5):

откуда

График функции f(x) представлен на рисунке.

Плотность распределения фактического времени ожидания обслуживания Примечание: составлен авторами по результатам данного исследования
F(0) ≈ 0,456 означает, что примерно 45,6 % поступивших заявок попали на обслуживание без ожидания в очереди.
Получим долю заявок, совсем не попавших на обслуживание по формуле (6):
(16)
Вычислим долю заявок, не получивших полного обслуживания, по формуле (7), используя соотношение (9):
Pнедо
(17)
Тогда доля заявок, получивших полное обслуживание, находится так:
Рполн = 1 – РОТК – Рнедо = 0,8927.
Таким образом, при данных входных параметрах почти 90 % пришедших в систему заявок будет обслужено полностью.
Приведенный пример показывает, что даже в простом варианте – с экспоненциальным обслуживанием и равномерно распределенным временем жизни заявок – модель сложна для практического применения. А для вариантов с другими функциями распределения потоков или групповым входным потоком почти непригодна. В этих случаях наиболее целесообразным решением является применение имитационного моделирования [13, 14].
На основе представленной математической модели была проведена верификация имитационной модели, построенной авторами для оценки характеристик СМО с ограниченным временем жизни заявок [15, 16]. Результаты работы имитационного алгоритма с описанными входными параметрами таковы:
Рполн = 0,9, РОТК = 0,0454, Рнедо = 0,052.
Полученные данные свидетельствуют о высокой точности имитационной модели.
Выводы
1. Проведен анализ одноканальной системы массового обслуживания заявок со случайным ограничением на общее время пребывания и стандартной дисциплиной обслуживания заявок в порядке их поступления.
2. На основе концептуальной модели Б. В. Гнеденко, И. Н. Коваленко для равномерно распределенного времени жизни заявок сделана постановка задачи и получено аналитическое решение.
3. Путем сравнения с аналитическим решением проведена верификация построенной авторами имитационной модели для изучения поведения СМО с ограниченным временем жизни заявок.
4. Доказана практически (до сотых долей) точность совпадения результатов имитационного и аналитического моделирования, которая позволяет утверждать, что разработанная имитационная модель может быть использована для других функций распределения потоков заявок без поиска аналитического решения.
5. Показано, что в практических целях только имитационное моделирование позволяет разрешать задачи нахождения эксплуатационных характеристик функционирования СМО с ограниченным временем жизни заявок для любых законов распределения случайных величин и любого характера входного потока требований.
6. Обосновано, что разработанная имитационная модель является уникальной и представляет ценность для решения прикладных задач проектирования и эксплуатации СМО: позволяет модифицировать структуру системы; дает возможность анализировать поведение всей системы и отдельных ее частей на этапе проектирования; реализовать поддержку эффективных управленческих решений с привлечением искусственного интеллекта.