Моделирование систем / Лекции / Q-схемы

НЕПРЕРЫВНО-СТОХАСТИЧЕСКИЕ МОДЕЛИ (Q-СХЕМЫ)

Особенности непрерывно-стохастического подхода рассмотрим на примере использования в качестве типовых математических схем систем массового обслуживания, (англ. queueing system), которые будем называть Q-схемами. Системы массового обслуживания представляют собой класс математических схем, разработанных в теории массового обслуживания и различных приложениях для формализации процессов функционирования систем, которые по своей сути являются процессами обслуживания.

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

В любом элементарном акте обслу­живания можно выделить две основные составляющие: ожидание обслужива­ния заявкой и собственно обслужива­ние заявки. Это можно изобразить в виде некоторого i-ro прибора обслужи­вания (рис. 8), состоящего из на­копителя заявок , в котором может одновременно находиться за­явок, где - емкость i-ro накопителя, и канала обслуживания за­явок (или просто канала) . На каждый элемент прибора обслу­живания поступают потоки событий: в накопитель поток за­явок , на канал - поток обслуживании .

Рис. 8. Прибор обслужи­вания заявок

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

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

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

Пример потока событий приведен на рис. 9, где обозначено Т_j – интервал между событиями (случайная величина); Т_н - время наблюдения, Т_с – момент совершения события.

Рис. 9. Поток событий

Интенсивность потока можно рассчитать экспериментально по формуле:

,

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

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

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

- ординарными, когда вероятность одновременного появления 2-х и более событий равна нулю;

- стационарными, когда частота появления событий постоянная;

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

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

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

, ,

где - величина, порядок малости которой выше, чем , т.е.

.

Стационарным потоком событий называется поток, для которого вероятность появления того или иного числа событий на интервале времени зависит лишь от длины этого участка и не зависит от того, где на оси времени 0t взят этот участок.

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

.

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

.

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

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

Заявки, обслуженные каналом , и заявки, покинувшие при­бор по различным причинам необслуженными (например, из-за переполнения накопителя ), образуют выходной поток , т.е. интервалы времени между моментами выхода заявок образуют подмножество выходных переменных.

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

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

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

Собственными (внутренними) параметрами Q-схемы будут яв­ляться количество фаз , количество каналов в каждой фазе , , количество накопителей каждой фазы , , ем­кость i-ro накопителя . Следует отметить, что в теории массо­вого обслуживания в зависимости от емкости накопителя приме­няют следующую терминологию для систем массового обслужива­ния: системы с потерями (, т.е. накопитель в приборе отсутствует, а имеется только канал обслуживания ), системы с ожиданием (, т.е. накопитель имеет бесконечную емкость и очередь заявок не ограничивается) и системы смешан­ного типа (с ограниченной емкостью накопителя ). Всю совокуп­ность собственных параметров Q-схемы обозначим как подмноже­ство H.

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

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

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

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

При ряде упрощающих предположений относительно подмно­жеств входящих потоков W и потоков обслуживания U (выполне­ние условий стационарности, ординарности и ограниченного после­действия), оператора сопряжения элементов структуры R (однофаз­ное одноканальное обслуживание в разомкнутой системе), подмно­жества собственных параметров H (обслуживание с бесконечной емкостью накопителя), оператора алгоритмов обслуживания зая­вок А (бесприоритетное обслуживание без прерываний и блокиро­вок) для оценки вероятностно-временных характеристик можно использовать аналитический аппарат, разработанный в теории мас­сового обслуживания. При принятых предположениях в обозначе­ниях Д. Кендалла будет иметь место классическая система обслу­живания типа М/М/1 (одноканальная система с марковским вхо­дящим потоком заявок и марковским потоком обслуживания). Рассмотрим на примере основные аналитические соотношения для такой элементарной Q-схемы.

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

, ;

,

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

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

, .

Перенеся влево и устремив к нулю, получим систему дифферен­циальных уравнений

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

Пусть в первом уравнении . Тогда . Подставив сюда значение из второго уравнения, находим . Повторяя эти операции, по­лучаем .

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

или .

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

Математическое ожидание числа заявок, находящихся в системе (приборе):

.

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

При этом

.

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

.

Математическое ожидание числа заявок, находящихся в накопителе:

.

Среднее время ожидания заявок в накопителе

.

Пример 2. Составить Q-схему системы.

Система передачи данных обеспечивает передачу пакетов данных из пункта А в пункт С через транзитный пункт В. В пункт А пакеты поступают через мс. Здесь они буферируются в накопителе емкостью 20 пакетов и передаются по любой из двух линий АВ1 - за время 20 мс или АВ2 – за время мс. В пункте В они снова буферируются в накопителе емкостью 25 пакетов и далее передаются по линиям BC1 (за мс) и ВС2 (за 25 мс). Причем пакеты из АВ1 поступают в ВС1, а из АВ2 - в ВС2.

Пример 3. Составить Q-схему системы.

Система обработки информации содержит мультиплексный канал и три мини-ЭВМ. Сигналы от датчиков поступают на вход канала через интервалы времени мкс. В канале они буферируются и предварительно обрабатываются в течение мкс. Затем они поступают на обработку в ту мини-ЭВМ, где имеется наименьшая по длине входная очередь. Емкости входных накопителей во всех мини-ЭВМ рассчитаны на хранение величин 10 сигналов. Время обработки сигнала в любой мини-ЭВМ равно 33 мкс.