книги / Прикладная теория систем массового обслуживания
..pdfСреднее число занятых каналов
к =Iар tR(h - 1,р,) + жхР(й,р/) х Ц ^ ----
1-Х
Вероятность того, что отдельный канал будет занят
к
ПЗ.К= - •
п
Вероятность занятости всех каналов системы
\—уп h |
, |
a-P(h,Pi)X г |
|
1-Х |
|
(п - h)(n + h +1) |
|
2 |
|
n - h + l, |
|
1_ уП-h+l
ОP(h,Pl)xl— ±----- ,
i-x
1 |
П |
Х*1, |
A < -, |
||
|
е |
|
f |
п |
Х = 1, |
A < -, |
||
|
е |
|
:s- и |
|
Х= 1, |
|
|
е Х*1-
Аналогичным образом могут быть определены и другие характери стики системы.
3.4.4.Системы массового обслуживания с отказами
инеоднородными потоками
Постановка задачи. На вход n-канальной СМО поступает неодно родный простейший поток с суммарной интенсивностью Х,ь причем
|
N |
ta~ |
, |
|
/=1 |
где Х{- интенсивность заявок в /-м источнике.
Так как поток заявок рассматривается как суперпозиция требований от различных источников, то объединенный поток с достаточной для прак тики точностью [23] можно считать пуассоновским для N = 5...20 и X,- = tai (ielfl). Интенсивность обслуживания одного прибора распределена по экспоненциальному закону и равна р = 1/г. Обслуживающие приборы для обслуживания заявки соединяются последовательно, что равносильно уве личению времени обслуживания во столько раз, сколько приборов объеди няется для обслуживания:
'обе = fo, Цобс= 1 /kt=\i/k,
где /обе - время обслуживания заявки; к - число обслуживающих приборов; Робе - интенсивность обслуживания заявки.
В рамках принятых в главе 2 допущений состояние СМО представим в виде вектора *, = >-’* т > - Л тах ), где кт- число заявок в
системе, каждая из которых обслуживается т приборами; L = qmfiX- #тт+1 - число входных потоков.
Тогда количество занятых и свободных приборов (n2m(Xj ), лСн(*, )) в состоянии Xj определяется следующим образом:
9max
I>*m >
т~Ят\
псъ(**) п ^зан(**)*
Из состояния х; система может перейти в любое другое состояние
Xj . Так как в системе действует L входных потоков, то из каждого состоя
ния потенциально возможно L прямых переходов. Однако из-за ограни ченности ресурсов системы не все эти переходы осуществимы. Пусть СМО находится в состоянии х, и приходит заявка, требующая т приборов. Если
т < исв(х, ), то заявка принимается на обслуживание и система переходит в
состояние Xj = (kqmin ,fc?min+1 ,- ,(к т + ) с интенсивностью К-
Если же заявка требует приборов больше, чем имеется свободных, то она получит отказ в обслуживании, а СМО останется в состоянии х ,. Если в состоянии Xj находятся заявки, требующие т приборов, то каждая из них
обслуживается с интенсивностью р/т , а общая интенсивность обслужи вания таких заявок (рт) определяется как р т = кт\х / т. При завершении обслуживания одной из заявок система перейдет в состояние, в котором соответствующая координата имеет значение, на единицу меньшее, чем в
состоянии X,, xv =[k4m.n,kqm.n+l,...,(km-ï),...,kqmM), т.е. произойдет об-
ратный переход. На рис. 3.9 представлен пример векторной модели СМО для п = 3, L = 3, <7mjn = 1, Ятах= 3, Р(т) = 1/3, XIt = X, интенсивность обслу живания прибора - р.
Рис. 3.9. Пример графа векторной модели СМО с отказами в обслуживании
Итак, каждое состояние х,• характеризуется числом обслуживаемых заявок определенного типа. Например, в состоянии х5(1,1,0) обслуживает ся одна заявка одним прибором и одна заявка двумя приборами. В этом со стоянии все приборы заняты, следовательно, возможны лишь обратные пе реходы (приход любой заявки в этом состоянии приводит к отказу в об служивании). Если раньше закончилось обслуживание заявки первого ти па, то система перейдет в состояние хА(0,1,0) с интенсивностью ц, если же раньше закончилось обслуживание заявки второго типа, то система перей дет в состояние х{(0,1,0) с интенсивностью р/2.
По графу состояний с нанесенными интенсивностями переходов со ставляется система линейных алгебраических уравнений. Из решения этих уравнений находятся вероятности Р(х{), по которым определяется харак теристика СМО.
Рассмотрим нахождение Pç^ (вероятность отказа в обслуживании).
где S - число состояний графа векторной модели СМО; Р(х, ) - вероят ность нахождения системы в состоянии х ,.
Число состояний согласно [11] определяется следующим образом:
S |
(3.22) |
где
Яmax
So - I
т~Яmin
|
|
|
|
|
|
N |
|
|
*7max |
W,-l |
|
ГПц-X"1 |
mN“1 |
n - 'Zmj |
|
|
|
7=1 |
|||||
S„=ml=9Il.ir+lm2I=9mn,+l |
=I |
EA |
|||||
(ОТд, + 1) |
|||||||
|
|
|
|
<7min "*"1 |
1-*<7min + 1 |
||
AL |
^ |
ffmin |
|
|
|
|
|
. ^7min 1 „ |
|
|
|
|
|||
|
|
|
|
|
|||
Определим число состояний векторной модели СМО по (3.22) для |
|||||||
примера, представленного на рис. 3.9. |
|
|
|||||
|
|
■*о=Х]3//и[ = 3 + 1 + 1 = 5 |
|
||||
|
|
т=1 |
|
|
|
|
|
|
|
Нщ* =](3-1)/(1 + !)[=!• |
|
||||
|
sn=т\~Ят\пЕ |
+1 |
3 - |
!(щ + !) = 1. |
|||
|
/=1 |
|
|
Следовательно, 5 = 1 + 5 + 1 = 7.
Для реализации реальных требований к обслуживающим приборам необходимо достаточно большое число п (40, ..., 50), а запросы на число обслуживающих приборов заявки на практике лежат в пределах 8-16. При таком соотношении приборов и запросов предложенный путь нахождения вероятностей становится чрезвычайно громоздким, т.к. векторная модель СМО имеет большое число состояний 5(50) = 1790, 5(60) = 4676, 5(70) = = 11075, а размер матрицы коэффициентов системы алгебраических урав нений пропорционален квадрату 5 [24], что требует большого объема па мяти ЭВМ и значительных затрат машинного времени. Стремление сни зить объем вычислений стимулировало поиск рекуррентных возможностей расчета P(xt ) на основе мультипликативных форм представления вероят
ностей состояний. В работе [11] представлен подход к расчету P(xt ):
|
9max |
<7max |
|
|
т~Ят\пП (РШт~ЯттE*mтml" |
||||
т ) = 5 |
9max |
*7max |
(3.23) |
|
тт\ |
||||
/Е- |
|
р к" |
||
П №) )Е |
|
|||
0m=<7mjn |
m-qm\ |
|
Использование предложенного в работе [11] критерия эквивалентно сти глобального и детального балансов цепей Маркова позволяет снижать размерность задачи и выполнять вычисления на ЭВМ средней мощности, используя рекуррентность вычислений. Кроме того, имеется возможность:
-произвести расчет для любых значений л;
-ускорить расчет и снизить затраты машинного времени. Аналогичным образом могут быть определены и другие характери
стики системы.
3.4.5. Примеры систем массового обслуживания с отказами
Пример 3.1. Рассматривается работа районной автоматической те лефонной станции (АТС), которая обеспечивает не более 120 переговоров одновременно. Средняя длительность переговоров 1/р = 1 мин. Вызовы на станцию поступают через 0,5 с, т.е. \Гк = 0,5 с. Требуется найти характери
стики работы АТС: к , J^j6c9 пЗ К, t з.к , / п.к •
Решение
АТС представляет собой 120-канальную СМО с отказами. Парамет ры системы следующие: п = 120; X= 2 сооб/с; р = 1/60 сооб/с; р = А/р = 120.
Среднее число занятых каналов
|
|
|
|
<DJ |
|
119 + 0,5-120] |
J |
Л(я-1,р) _ 12Qfi(l 19,120) |
120 |
|
У ш |
||||
* = рЛ>бс=Р Л(п,р) |
|
|
~ |
г |
120 + 0,5-120"] |
112. |
|
|
Л(120,120) |
Ф* |
J |
||||
|
|
|
|
ч |
■JÏ2Ô |
||
Вероятность обслуживания |
|
|
|
|
|||
|
|
112 «0,931. |
|
|
|
||
Робе |
к /р |
|
|
|
|||
Вероятность занятости канала |
|
120 |
|
|
|
|
|
|
|
|
|
|
|
||
7I3K=Â/w= — «0,931. |
|
|
|
||||
|
3 |
|
120 |
|
|
|
|
Среднее время занятости канала |
|
|
|
|
|
||
Г3.к= 1/р = 1 мин = 60 с. |
|
|
|
||||
Среднее время простоя канала |
|
|
|
|
|
|
|
ti |
—7 ^ |
^п.к |
~ 4,5 с. |
|
|
|
|
“ |
«з.к |
|
|
|
|
Ппи
Пример 3.2. Обслуживание заявок производится СМО с отказами, параметрами системы л, X, р. Обслуживание каждой заявки приносит среднюю прибыль Ci. С целью увеличения доходов обслуживания предла гается провести реорганизацию, состоящую в том, что система будет до пускать взаимопомощь между каналами. На преобразование СМО с отка зами в СМО с отказами и взаимопомощью требуется израсходовать стой-
мость Сг. Переоборудование занимает время тп. Определить, по истечении какого срока времени t после начала преобразования вновь организованная СМО с отказами и взаимопомощью начнет приносить прибыль.
Решение
Задачу будем решать при условии, что время установления стацио нарного режима в системах мало по сравнению со временем переоборудо вания тп и временем t (t > т„). Начало отсчета для времени тп и t одно. В этом случае можно записать следующее уравнение:
С1Х0(,)г= -С2+С,Ао(2)( г- тп)>
- абсолютная пропускная способность СМО с отка-
R{n,p)
1 Y ^
зами; Хо(2)= X-----^-г - абсолютная пропускная способность СМО с отка- 1 -х "+
зами и взаимопомощью; р = ХУр; х = ХУир.
Решая это уравнение, получим /, по истечении которого реорганиза ция начнет приносить прибыль:
/ = + С т Д 1-хя с,х i-x" |
R{n-IPŸ |
1-х ” +1 |
Д(и>р) , |
Пример 3.3. Определить, насколько увеличится вероятность обслу живания для СМО с отказами, имеющей следующие параметры: п = 10, X= = 8 сооб/мин; р = 0,8 сооб/мин, если обеспечить взаимопомощь группы из двух каналов (/ = 2).
Решение
Для СМО с отказами без взаимопомощи вероятность обслуживания
л>бс(,)= Л(я-1,Р)
R(n,p)
Для примера
р <» = |
*(9,64) |
0,795. |
•г обе |
Л(10,64) |
|
|
|
Для системы с частичной взаимопомощью имеем
N=10; х = Х/иц= 1; р/= X//р = 0,5; /г = ]«//[ = 5.
Для этой системы в случае, когда %= 1, вероятность обслуживания
Таким образом,
Робс(2)- Ллс(1) = 0,883 - 0,795 = 0,088.
Пример 3,4. Рассматривается система ПВО с нарушенным управле нием. Нарушение управления состоит в том, что каждую влетающую в зо ну обстрела цель обстреливают все свободные к этому времени каналы. Обстрел цели каждым каналом длится случайное время, распределенное по показательному закону с параметром ц. За это время каждый канал по ражает цель с вероятностью р независимо от других каналов, принимаю щих участие в обстреле. Определить характеристики работы системы, если параметры системы ПВО следующие: п = 4; X = 4 1/мин; ц = 1 1/мин; р = = 0,5.
Решение
у
В рассматриваемом примере величина а = —= 4 - целое число. Ве-
роятности различных состояний системы следующие:
а+п
Вероятности поражения цели
Вероятность полной загрузки системы
а
я,п.з а +п = 0,5.
Среднее время полной загрузки
tnз = — = 0,25 мин. n\i
Среднее время неполной загрузки
Гн з = ~ = 0,25 мин.
к
Среднее число занятых каналов
k = Z k p k =3,2.
k=0
Вероятность того, что канал занят,
л з.к = — = 0,8.
п
Среднее время занятости канала
t3к = —= 1 мин.
Среднее время простоя канала
^п.к |
1 |
1 - 7 1 , - = 0,25 мин. |
|
V |
*з.к |
Среднее время пребывания заявки в системе
1 п /^ а - [ , n-k- 1 |
1 |
X (-1 Г С П 1- , — ^ - I 0,429 мин. |
|
Ц*=1 Са+„ <п=о |
(m +1) |
Среднее число обстреливаемых целей
I =tX =1,75.
4.СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ
СОЖИДАНИЕМ
До сих пор рассматривались системы массового обслуживания с отказами; характерной особенностью таких систем было то, что любая по ступившая заявка либо немедленно принималась к обслуживанию, либо немедленно получала отказ и покидала систему.
С ростом X для достижения заданных величин вероятностно временных характеристик необходимы значительные аппаратурные затра ты и, может быть, выгоднее перейти к использованию в составе систем устройств памяти. Отсюда появляется целесообразность рассмотрения СМО с очередью.
Вэтой главе будут рассмотрены системы массового обслуживания
сожиданием, в которых заявка, заставшая все каналы занятыми, не полу чает немедленного отказа, а должна стать в очередь и ожидать освобожде ния канала, который может ее обслужить.
Системы ожидания бывают чистого или смешанного типа. В чистой СМО с ожиданием число мест в очереди и время ожидания в ней ничем не ограничены: каждая заявка рано или поздно будет обслужена. Для такой
системы понятие «отказ» не имеет смысла. Это системы массового обслу живания с бесконечной очередью.
Всистеме с ожиданием смешанного типа возможны как отказы, так
иожидание заявки в очереди. Отказы (отсутствие обслуживания) могут быть связаны или с ограниченным числом мест в очереди, или с ограни ченным временем ожидания, которым располагает заявка. Это системы массового обслуживания с ограниченной очередью.
При рассмотрении СМО с ожиданием необходимо учитывать сис тему правил, регламентирующих порядок образования и обслуживания очереди (так называемую дисциплину очереди). Необходимо указать, яв ляется ли очередь общей или образуется к каждому каналу отдельно; каков порядок вызовов заявок из очереди и т.д. Из всего многообразия известных дисциплин обслуживания [17, 25] выберем дисциплину со следующими особенностями: очередь является общей, т.е. все приходящие заявки ста новятся в одну очередь к обслуживающим приборам; порядок вызовов зая вок из очереди FIFO (first in - first out) - «кто раньше встал в очередь, тот и раньше обслуживаете*».
Поведение заявок в очереди также входит в понятие «дисциплина очереди». Заявки в очереди могут «терпеливо» ждать начала обслужива ния, а могут и уходить из системы, не дождавшись своей очереди. В этой главе будут рассматриваться только «терпеливые» заявки.
4.1.Классическая система массового обслуживания с ожиданием
Постановка задачи. В качестве обслуживающей системы рассмот рим неблокирующую коммутационную систему, которая обслуживает полнодоступный пучок емкостью я ( 1 < я < оо) линий, на которые поступа ет простейший поток с параметром X. Каждая поступившая заявка для об служивания занимает любую свободную линию пучка и обслуживается с интенсивностью р. Если все я линий пучка заняты, то она становится в очередь и «терпеливо» ждет своего обслуживания. Дисциплина очереди - FIFO, максимальное число мест в очереди - т (СМО с конечной очере дью). Если заявка застает все т мест в очереди занятыми, то она получает отказ и исключается из обслуживания. Величины я, X, р, т будем называть параметрами СМО с ожиданием.
Макросостояния рассматриваемой системы будем связывать с чис лом заявок, находящихся в системе (обслуживаемых и ожидающих в оче реди):
хк - в системе имеется к заявок (к= 0, 1,2, ..., я), они обслуживают ся к каналами, очереди нет;
хп+ г - в системе имеется я + г заявок (г = 1,2, ..., /я), я из них об служиваются в я каналах и г заявок находятся в очереди.
Таким образом, система имеет я + т + 1 макросостояний. Граф мак росостояний рассматриваемой системы показан на рис. 4.1.
М- |
2ц (л-1)ц |
пц |
n\L |
n\i |
иц |
Рис. 4.1. Граф макросостояний классической системы массового обслуживания с ожиданием (СМО с конечной очередью)
Этому графу состояний соответствует система дифференциальных уравнений для вероятностей состояний, которую обычно интегрируют для начальных условий
/> о ( 0 ) = 1 , Л ( 0 ) = 0 ( 1 * 0 ) ,
т.е. в момент времени t = 0 система свободна от заявок. Для стационарного режима работы СМО с ожиданием, когда X = const, р = const,
используя результаты главы 3, имеем
а* Рк =— Р0 (к =0,1,2,...,л); Рп+Г = хгР„ (г = 0,1,2,..,,т) ,
к\